Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,234,667 Δ1 bp intergenic (‑1434/‑499) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,234,6660A.96.0% 301.2 / 11.6 79intergenic (‑1433/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base . (37/34);  minor base G (0/3);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_pae_pgi/4234523‑4234809
                                                                                                                                               |                                                                                                                                               
atataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑atcatc                                                                                                                                           <  1:117515/149‑1 (MQ=255)
  ataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAAATGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcgt                                                                                                                                         <  1:167806/149‑1 (MQ=255)
   taCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCCTTGAAGAGGCGCTTGCAGTG‑ATCAtcgtc                                                                                                                                        >  1:468761/1‑149 (MQ=255)
      gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGat                                                                                                                                     <  1:888905/149‑1 (MQ=255)
      gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGag                                                                                                                                     >  1:128030/1‑148 (MQ=255)
         gTAGAGATGCACCCCAAGGTCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGGTTAGAGGCGTAAATCCCCTGGCCACCCCGAAAATAGATGCATTGAAGCGGCGCTTGCGGTG‑TTCATCGTCGatgaa                                                                                                                                  >  2:1023616/1‑146 (MQ=255)
          tAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGt                                                                                                                                 >  2:847917/1‑149 (MQ=255)
              gaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGc                                                                                                                             <  1:1164232/149‑1 (MQ=255)
                tGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATATGGGCATCGCATCTGTAAAAATGACAAAGATCCAAGGCGTAAATCTCCTGCCCACCCCGAGAATCGACGCACTGCAGCGGCGCATGGAGTG‑ATCTTCGTCGATTTGTaagtat                                                                                                                           >  2:648870/1‑144 (MQ=255)
                   cacccAAACGCGTTTAATCCTTATCACGAGATGTGCGCACCGCCTCTGCAAAACTAGCAGAGATTAGAGGCGTAAACCCCCCGGTCACCCCGAGAATAGATGCATTGAATAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGAt                                                                                                                        <  2:407996/147‑1 (MQ=255)
                     cctaaGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATaa                                                                                                                      <  1:1053039/146‑1 (MQ=255)
                     cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAAATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGACGCTTGCAGTG‑ATCACCGTCGCTATGTAGGCCGGATaa                                                                                                                      >  1:229193/1‑149 (MQ=255)
                     accAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGT‑CATCATCGTCGCTATATAGGCCGGATaa                                                                                                                      >  2:64028/2‑149 (MQ=255)
                      ccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                     >  1:511878/1‑149 (MQ=255)
                      ccAAGGGCGTTTTATCCTTGTTACGAGATGGGGGACTCGCCCCTGTAAAACTGGAAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGTATTTAAGCGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                     <  1:225495/149‑1 (MQ=255)
                           ggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                <  1:847917/149‑1 (MQ=255)
                           ggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                <  2:1521/149‑1 (MQ=255)
                           ggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTCTAAAAATGGCATAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGT‑CGTCATGGTTGATATGTAGGCAGGATATGGCGtt                                                                                                                >  2:707786/1‑149 (MQ=255)
                               ttttGTCCGTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcg                                                                                                            <  2:792432/149‑1 (MQ=255)
                               ttttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcg                                                                                                            >  1:8020/1‑149 (MQ=255)
                                    tCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGTGATTAGAGGCGTAAATCCCCTGGCCACTCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGTCGGTTTAGGCGTTCACGGCGCa                                                                                                       >  1:866680/1‑149 (MQ=255)
                                    tCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                                                       >  1:315058/1‑149 (MQ=255)
                                      cTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                     >  1:837109/1‑149 (MQ=255)
                                             cGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                              <  1:1210022/149‑1 (MQ=255)
                                                gaTGTGGCGATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGa                                                                                           <  1:409616/149‑1 (MQ=255)
                                                 aTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTTGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCTTCGATATGTAGACCGGATTAGTCGTTCACGCCGCATCCGGCAACCGAt                                                                                          >  2:409616/1‑149 (MQ=255)
                                                   gtgGGCATCGCCTCTGTAAAACTGGCAAAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGTTGc                                                                                        >  2:1189448/1‑149 (MQ=255)
                                                     gggCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCTGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                      <  1:1189837/149‑1 (MQ=255)
                                                       gagtcctctctGTAAAACTGGCATGGGTTTGAGGCGTAAATCCCCTGGCCACCTCGAAAATAGATGCATTGAAGAGGCGTTTGCAGT‑CATCATCGTCGATATGTAGGCCGGATAAGACGTTCACGCCGCGTCCGGCAACCGATGCCTGa                                                                                    <  2:81325/142‑1 (MQ=255)
                                                         aTCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                                                                  <  2:314717/149‑1 (MQ=255)
                                                          tCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGc                                                                                 >  2:143920/1‑149 (MQ=255)
                                                             ccTCTGTAAACCTGGCAGAGATTAGAGGCGAAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                              <  1:1014929/149‑1 (MQ=255)
                                                             ccTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTGGCAGTG‑ATCATGGCCGCTATGTTGGCCGGAGCAGGCGTTCACGCCACATACGGCAACCGATGACCGAttctga                                                                              >  2:459198/1‑144 (MQ=255)
                                                              ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcgtc                                                                                                                                        <  1:315898/90‑1 (MQ=255)
                                                              ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcgtc                                                                                                                                        >  2:315898/1‑90 (MQ=255)
                                                               acagaaaaaCCGGCAAAGATAAGACGCGTAAATCCCCTGACCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGAAGTG‑ATCATCGACCATATGAAGGCCAGCAAATTCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                            <  2:866680/144‑1 (MQ=255)
                                                                   tAAAACTGGCAGGTATTAGAGGCGTAAATCTCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                        <  1:674001/149‑1 (MQ=255)
                                                                        cTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                                                                   <  2:315058/149‑1 (MQ=255)
                                                                           gCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCTTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  2:381201/1‑149 (MQ=255)
                                                                           gCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  1:468654/1‑149 (MQ=255)
                                                                             agagATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTcgcatc                                                              <  2:814669/149‑6 (MQ=38)
                                                                             agagATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCAGCAACCGATGCCTGATGCAACGCGGTCGCGCCTTATc                                                              >  2:684735/1‑149 (MQ=38)
                                                                             agagATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                              >  1:1036142/1‑149 (MQ=255)
                                                                             agagATTAGAAGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAATGGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                              <  2:740595/149‑1 (MQ=255)
                                                                              gagaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                             >  1:1001013/1‑149 (MQ=255)
                                                                                gaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg                                                           <  1:143920/149‑1 (MQ=255)
                                                                                 aTTAGAGGCGTTAATCCCTTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                          <  1:918439/149‑1 (MQ=255)
                                                                                 aTTAGAGGCGTAAACCCCCTGGGCACCCCGAGAATCGGTGCATTGAAGAGGCGCTTGCAGT‑TATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                          <  2:1036142/149‑1 (MQ=37)
                                                                                          gTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc                                                 >  1:788164/1‑149 (MQ=255)
                                                                                            aaaTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                               <  2:837109/149‑1 (MQ=255)
                                                                                            aaaGCCCCTGGCCACTCCGATAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGACAACCGAGGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                               <  1:474580/149‑1 (MQ=37)
                                                                                                  ccTGGCGACCCCGGGAATAGATGCATTTCAGAGTCTCCTGCAG‑GGGTCATCGTCGATATGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                         <  1:707786/149‑1 (MQ=21)
                                                                                                    tGGCCACCGCGAGAATAGATGCATTGAAGAGGCGTTTGCAGTG‑ACCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                       <  2:534930/149‑1 (MQ=255)
                                                                                                    tGGCCACCCCGAGAATAGATGCATTGGAGAGGCGCTGGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                       <  2:591351/149‑1 (MQ=255)
                                                                                                       ccACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGAAAAGGAGTTCACGCCGCATCCGGCAACCGATGCCTTATGCGACGCGGCCGCGTCTTACCAGGCCTACAGGTCGATGCCCattttt                                    >  2:428256/1‑145 (MQ=255)
                                                                                                          ccccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACa                                 >  2:1138773/1‑149 (MQ=255)
                                                                                                           cccGAGAATAGATGCATTGACGAGGCGCTTGCAG‑GGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                <  2:637248/149‑1 (MQ=255)
                                                                                                           cccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                >  2:1124581/1‑149 (MQ=255)
                                                                                                            ccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                               <  2:932975/85‑1 (MQ=255)
                                                                                                            ccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                               >  1:932975/1‑85 (MQ=255)
                                                                                                             cGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCg                              >  1:334404/1‑149 (MQ=255)
                                                                                                               agaATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTa                            <  1:1138773/149‑1 (MQ=255)
                                                                                                                    aGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTCTCAGGCCTACAGGTCGATGCCCATATGTCAATCGTATGCgg                       >  1:498256/1‑149 (MQ=25)
                                                                                                                      aTGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGa                                                                               <  2:356470/91‑1 (MQ=255)
                                                                                                                      aTGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGa                                                                               >  1:356470/1‑91 (MQ=255)
                                                                                                                      aTGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGGCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCTTGATGCGACGCGGTCGCGTCTTATCACGCTTAAAGGTTGATGACATTTTGTAAATCGGATTTggga                     >  2:683551/1‑147 (MQ=21)
                                                                                                                       cgCATTGAAGGGGTTCTTGCAG‑GGATCATCGGCGATATGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                    <  1:356309/148‑1 (MQ=25)
                                                                                                                        gCATTGAAGGGGCGCTTGCGGTG‑ATCATAGTCGATGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCGGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAt                   <  1:64821/149‑1 (MQ=25)
                                                                                                                         cATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTAGGTAAtt                  >  1:152131/1‑149 (MQ=255)
                                                                                                                         cATTGAAGAGGCACTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                  <  1:811693/149‑1 (MQ=255)
                                                                                                                          aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                 >  2:619311/1‑149 (MQ=255)
                                                                                                                          aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                 >  1:575542/1‑149 (MQ=255)
                                                                                                                          aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTAACGCCGCCTCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                 >  1:480897/1‑149 (MQ=255)
                                                                                                                              aaGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTCCATCGTATTCGGCAATTAATAc             >  1:355758/1‑149 (MQ=255)
                                                                                                                                 gggCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa          <  2:1001013/148‑1 (MQ=255)
                                                                                                                                 agggccTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa          <  2:137213/144‑1 (MQ=255)
                                                                                                                                  ggCGCTTGCAGTGAA‑CATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAg         <  1:1065360/149‑1 (MQ=255)
                                                                                                                                   gcgcTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGc        >  1:275090/1‑149 (MQ=255)
                                                                                                                                        tGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGGGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat   <  2:509539/149‑1 (MQ=255)
                                                                                                                                         gCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgatt  >  1:659250/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                               
ATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_pae_pgi/4234523‑4234809

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑