Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,637,439 C→A 5.3% intergenic (‑83/+35) der ← / ← bamB GTPase; multicopy suppressor of ftsJ/BamABCDE complex OM biogenesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,637,4390CA 5.3% 1379.0 / 42.5 375intergenic (‑83/+35)der/bamBGTPase; multicopy suppressor of ftsJ/BamABCDE complex OM biogenesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (188/167);  new base A (15/5);  total (203/172)
Fisher's exact test for biased strand distribution p-value = 6.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 381 total aligned reads displayed.

TTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAA  >  NC_000913/2637290‑2637588
                                                                                                                                                     |                                                                                                                                                     
ttAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTc                                                                                                                                                       <  2:825380/150‑1 (MQ=255)
  aGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCa                                                                                                                                                     >  1:5649979/1‑150 (MQ=255)
             tAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGtttt                                                                                                                                          >  2:2359315/1‑150 (MQ=255)
              aaCGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGATATATTCGCGTCGTTGTTAAAAACAGTAAAACGGCTCCTGGACAGGGGCCGTTTTc                                                                                                                                         <  2:72623/150‑1 (MQ=255)
                 gTGGATTTTCCTACGTTTGGGCGCCCGACACGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCAACAGACACCGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCaaa                                                                                                                                      <  1:361793/150‑1 (MQ=255)
                 gTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCaaa                                                                                                                                      >  2:5474105/1‑150 (MQ=255)
                 gTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCaaa                                                                                                                                      >  2:4453147/1‑150 (MQ=255)
                 gTGGATGGTCCTACGTTAGGGCGCCCGACAAGCGCGACCACCGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCaaa                                                                                                                                      <  2:5688017/150‑1 (MQ=255)
                      ttttCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGaa                                                                                                                                 >  1:3919284/1‑150 (MQ=255)
                          ccTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGAc                                                                                                                             <  2:5180867/150‑1 (MQ=255)
                                  aGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGa                                                                                                                     <  2:1611321/150‑1 (MQ=255)
                                    ggCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGAtt                                                                                                                   >  1:5632066/1‑150 (MQ=255)
                                    ggCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGAtt                                                                                                                   >  2:1137472/1‑150 (MQ=255)
                                         ccGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACgt                                                                                                              >  1:5793930/1‑150 (MQ=255)
                                         ccGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACgt                                                                                                              >  1:2608181/1‑150 (MQ=255)
                                          cGACAAGCGCGACCACAGGTTCCATGTTTATAGGCTCATTTTTATAAATCATCAGTCAACGCTCGATATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACgtg                                                                                                             <  2:2322809/150‑1 (MQ=255)
                                          cGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACgtg                                                                                                             >  1:4480554/1‑150 (MQ=255)
                                          cGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACgtg                                                                                                             >  2:1231100/1‑150 (MQ=255)
                                          cGACAAGCGCGACCACAGGTACATTTTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTTTTAAAAACAGTAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACgtg                                                                                                             <  1:1716354/150‑1 (MQ=255)
                                           gACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACgtgt                                                                                                            <  1:5178593/150‑1 (MQ=255)
                                           gACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACgagt                                                                                                            >  2:912345/1‑150 (MQ=255)
                                             cAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGATATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTaa                                                                                                          <  1:5544812/150‑1 (MQ=255)
                                               cgcgcGACCACTGGTACCTTGTTTAAAGCCTCATTTTTATAAATAATCACACAACCCACGCTATATTCGCGTCTTTGTTAAAAACAGCAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATa                                                                                                        <  1:4887296/149‑1 (MQ=255)
                                               atctggACCACAGGTACCATGTTTAACGCCTCTTTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAATACCAGGAAAACGGCGCCGGTCCAGGAGACGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTtata                                                                                                        <  2:3504277/145‑1 (MQ=255)
                                                 cgcgACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATaga                                                                                                      >  1:2842279/1‑150 (MQ=255)
                                                    gACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTa                                                                                                   >  1:2714686/1‑150 (MQ=255)
                                                    gACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTa                                                                                                   >  2:1557304/1‑150 (MQ=255)
                                                      ccACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTaca                                                                                                 <  2:4759721/150‑1 (MQ=255)
                                                          aGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGt                                                                                             >  2:2091461/1‑150 (MQ=255)
                                                            gTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTc                                                                                           >  1:1388952/1‑150 (MQ=255)
                                                             tACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGAGTAACGTGTACTAGAGTACACGGTTcc                                                                                          >  2:1926082/1‑150 (MQ=255)
                                                               ccATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTTTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGt                                                                                        >  2:1785367/1‑150 (MQ=255)
                                                                 aTGTTTAAAGCCTCATTTTTATAATTCATCAGCCAACGCACGCTCTATTCGCGTCGTTGTTAAACACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACCGAGACGATTAACGTGTTATAGAGTACACGGTTCCGTCt                                                                                      <  1:1161992/150‑1 (MQ=255)
                                                                  tGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCtt                                                                                     >  1:789452/1‑150 (MQ=255)
                                                                  tGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCtt                                                                                     <  1:545106/150‑1 (MQ=255)
                                                                        aaGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCt                                                                               >  2:4552907/1‑150 (MQ=255)
                                                                           ccTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGa                                                                            >  1:2039037/1‑150 (MQ=255)
                                                                             tCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATc                                                                          >  1:1752124/1‑150 (MQ=255)
                                                                               ctttttATAAATCATCAGACACCGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATcag                                                                        <  2:2006296/149‑1 (MQ=255)
                                                                               aTTTTTATAAATCATCAGACACCGCACGCTATAGTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATcag                                                                        <  1:3998370/150‑1 (MQ=255)
                                                                               aTTTTTATAAATCATCAGACAACGCACGCTATAATCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGATTACACGGTTCCGTCTTTTGCCTGGATcag                                                                        >  1:3820414/1‑150 (MQ=255)
                                                                                 ttttATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATcagca                                                                      <  2:1942465/150‑1 (MQ=255)
                                                                                 ttttATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATcagca                                                                      <  2:4903873/150‑1 (MQ=255)
                                                                                 ttttATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATcagca                                                                      <  2:4849430/150‑1 (MQ=255)
                                                                                      taAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTg                                                                 >  2:3272502/1‑150 (MQ=255)
                                                                                      taAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTg                                                                 >  2:3470617/1‑150 (MQ=255)
                                                                                        aaTCATCAGACCCCGCACGCTATATTCGCGTCGTTGTTAAACACCGGAACACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGcc                                                               <  2:1443195/150‑1 (MQ=255)
                                                                                        aaTCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTTAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGcc                                                               <  1:4783031/150‑1 (MQ=255)
                                                                                        aaTCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGcc                                                               <  1:2330891/150‑1 (MQ=255)
                                                                                        aaTCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGcc                                                               >  1:3365132/1‑150 (MQ=255)
                                                                                        aaTCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGcc                                                               >  1:3710048/1‑150 (MQ=255)
                                                                                             tcaGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAg                                                          >  2:1207916/1‑150 (MQ=255)
                                                                                             tcaGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAg                                                          >  2:4403187/1‑150 (MQ=255)
                                                                                               aGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCg                                                        >  2:2911861/1‑150 (MQ=255)
                                                                                                gACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCCGTTTGCCGTCAGCgg                                                       >  1:1011893/1‑150 (MQ=255)
                                                                                                gACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCgg                                                       >  2:529279/1‑150 (MQ=255)
                                                                                                gACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCgg                                                       >  1:5677075/1‑150 (MQ=255)
                                                                                                 aCAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGc                                                      >  1:5040604/1‑150 (MQ=255)
                                                                                                    ccccacgcTAGACTCGCGTCGCTGCCACCAACAGGAAACCGGCCCCAGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTCACGTGCAAGAGAGTCCACGGTTCCGTCTTGTGCCTGGATCAGCAGTGGGCCGTCAGCGGCAAc                                                   <  2:4087961/147‑1 (MQ=255)
                                                                                                     cgcacgcTATATTCGCGTCGTTGATAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAAcc                                                  <  1:620906/150‑1 (MQ=255)
                                                                                                      gcacgcTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCg                                                 >  1:5078137/1‑150 (MQ=255)
                                                                                                         cgcTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGtt                                              >  1:4280821/1‑150 (MQ=255)
                                                                                                         cgcTATAGTCGCGTCGTTGTTAACAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGtt                                              <  2:1641594/150‑1 (MQ=255)
                                                                                                           cTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCa                                            <  1:1463881/150‑1 (MQ=255)
                                                                                                            tataTTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTATCAAAGTGAACGACAGAGACGATAAACGTGTAATAGAGTACACGGTTCCGTCTTTTTCCTGGATCAGCAGTTTGCCGTAAGCGGCAAACGGTTCAg                                           >  1:3530443/1‑150 (MQ=255)
                                                                                                             atagTCGCGTCGTTGTTAAAAACAGGAACACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGt                                          <  2:5780635/150‑1 (MQ=255)
                                                                                                              tatTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTaca                                                                                                 <  1:241556/94‑1 (MQ=255)
                                                                                                              tatTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTaca                                                                                                 >  2:241556/1‑94 (MQ=255)
                                                                                                                ttCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTg                                       <  1:1917114/150‑1 (MQ=255)
                                                                                                                  cgcgTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGa                                     <  2:3739555/150‑1 (MQ=255)
                                                                                                                    cgtcgtTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGaaa                                   <  2:560124/150‑1 (MQ=255)
                                                                                                                       ctttgttAAAATCAGGAAAACTGACCCTGTCCAGGAGCCGTTTTCTATGTGAACTACAGAGACGATTAACGTGTAATAGTGTACACTGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCg                                <  1:2093551/148‑1 (MQ=255)
                                                                                                                        gttgttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCgg                               <  1:2954031/150‑1 (MQ=255)
                                                                                                                         ttgttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGa                              >  1:107828/1‑150 (MQ=255)
                                                                                                                          tgttAAACACAGGACCACGGCCCCGGTCCAGGAGCCGGTTTCTAAGTGAACGACTGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGaa                             <  2:2264576/150‑1 (MQ=255)
                                                                                                                            ttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACt                           >  2:2183626/1‑150 (MQ=255)
                                                                                                                              aaaaaCAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTAt                         <  2:526416/150‑1 (MQ=255)
                                                                                                                                aaaCAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCa                       <  2:3634524/150‑1 (MQ=255)
                                                                                                                                 aaCAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCaa                      <  2:5002496/150‑1 (MQ=255)
                                                                                                                                    aGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCACCCGGTTCAGTCTGGAAACCGGACCTATCAACtt                   >  2:3322421/1‑150 (MQ=255)
                                                                                                                                     ggAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACttt                  >  1:4632192/1‑150 (MQ=255)
                                                                                                                                     ggAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACttt                  >  2:3139373/1‑150 (MQ=255)
                                                                                                                                       caaaCGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGACCGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGACCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACtttta                <  1:5380849/149‑2 (MQ=255)
                                                                                                                                        aaaCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGGACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCCACTTTTTg               >  2:50606/1‑150 (MQ=255)
                                                                                                                                        aaaCGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTAAGCGGCAACCGGTTCAGTCTGGAAACCGGAACAATCAACTTTTTg               >  1:4243339/1‑150 (MQ=255)
                                                                                                                                         aaCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGc              >  2:5425460/1‑150 (MQ=255)
                                                                                                                                         aaCGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGc              <  1:3230825/150‑1 (MQ=255)
                                                                                                                                          aCGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCt             >  2:1478070/1‑150 (MQ=255)
                                                                                                                                           cGGCCCCTGTCCAGGAGCCGTTTTGAATGTGATCGACAGAGACGATTATCGTGTAATAGAGTACACGGTTGCGTCTTTTGCCTGGATCAGCAGTTTGGCGTCAGCTGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTg            <  2:3491393/150‑1 (MQ=255)
                                                                                                                                           cGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTg            <  2:2669842/150‑1 (MQ=255)
                                                                                                                                           cGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTg            <  2:5040403/150‑1 (MQ=255)
                                                                                                                                            ggCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATATATTACACGGTTACGTCTTTTTCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGCACAATCAACTTTTTGCTgg           >  1:5198392/1‑150 (MQ=255)
                                                                                                                                            ggCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCTGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTgg           <  2:4223782/150‑1 (MQ=255)
                                                                                                                                            ggCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTgg           <  2:1167272/150‑1 (MQ=255)
                                                                                                                                             gCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTggg          <  1:4074515/150‑1 (MQ=255)
                                                                                                                                              ccccTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGc         <  2:3358266/150‑1 (MQ=255)
                                                                                                                                               cccTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCa        >  2:1350941/1‑150 (MQ=255)
                                                                                                                                                 cTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCaac      <  2:389203/150‑1 (MQ=255)
                                                                                                                                                     ccAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCaacgaaa  <  2:564728/150‑1 (MQ=255)
                                                                                                                                                     ccAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCaacgaaa  <  1:735809/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAA  >  NC_000913/2637290‑2637588

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: