Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,478,271 G→A 5.8% A16T (GCA→ACA)  ynbB → putative CDP‑diglyceride synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,478,2710GA 5.8% 274.7 / 3.6 87A16T (GCA→ACA) ynbBputative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base G (41/40);  new base A (2/3);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

GGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCG  >  NC_000913/1478123‑1478420
                                                                                                                                                    |                                                                                                                                                     
ggATGCAGTGGAATAATCTATTATGGAGCATTGCTGCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGc                                                                                                                                                      <  2:1313170/150‑1 (MQ=255)
     cAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGCTAAAATTAACAACGc                                                                                                                                                 <  2:990727/150‑1 (MQ=255)
      aGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCt                                                                                                                                                <  1:1516241/150‑1 (MQ=255)
      aGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGCTAAAATTAACAACGcg                                                                                                                                                <  1:1321055/150‑2 (MQ=255)
       gTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCTACACTTTTCGCGTTGTTAATTTTAGCTACGCTg                                                                                                                                               >  1:395092/1‑150 (MQ=255)
        tGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGa                                                                                                                                              <  2:705666/150‑1 (MQ=255)
         ggAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGCTAAAATTAACAACGcgaaa                                                                                                                                             <  2:1188040/150‑5 (MQ=255)
         ggAATAATCTTTTATGGAGAATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGAt                                                                                                                                             >  1:822508/1‑150 (MQ=255)
           aataatCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAAACTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATaa                                                                                                                                           >  2:1644671/1‑150 (MQ=255)
                    ttATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTc                                                                                                                                  >  1:1599465/1‑150 (MQ=255)
                     tATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCt                                                                                                                                 <  1:1532756/150‑1 (MQ=255)
                     tATGGAGCATTCCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAAAGCTGGAAAAATCTCTGGCAACACTTTTCGCGATGTTAATTTTAGCAACGCTGATAAAACGTTTTCt                                                                                                                                 >  1:168863/1‑150 (MQ=255)
                      aTGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTg                                                                                                                                <  2:396353/150‑1 (MQ=255)
                      aTGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTg                                                                                                                                <  2:547221/150‑1 (MQ=255)
                        ggAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCACCGCTGATAAACCGTTTTCTGCt                                                                                                                              >  2:972593/1‑150 (MQ=255)
                        ggAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCt                                                                                                                              <  1:499902/150‑1 (MQ=255)
                        ggAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCt                                                                                                                              >  1:423632/1‑150 (MQ=255)
                         gAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCtt                                                                                                                             >  1:1526726/1‑149 (MQ=255)
                         gAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTa                                                                                                                             >  2:1023215/1‑150 (MQ=255)
                         gAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTa                                                                                                                             >  1:162527/1‑150 (MQ=255)
                         gAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTa                                                                                                                             <  1:806173/150‑1 (MQ=255)
                              ttGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCTTTTTCTGCTATGGCg                                                                                                                        <  1:912727/150‑1 (MQ=255)
                                gCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGt                                                                                                                      <  1:1361814/150‑1 (MQ=255)
                                 cTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGtt                                                                                                                     >  1:1271618/1‑150 (MQ=255)
                                     aaTTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCg                                                                                                                 >  1:1101937/1‑150 (MQ=255)
                                       ttCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGga                                                                                                               >  2:710568/1‑150 (MQ=255)
                                           gCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGgagaga                                                                                                           <  2:473680/150‑1 (MQ=255)
                                              tctcTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAACCCGTTTTCTGCTATGGCGGTTACCGGAGAGataa                                                                                                        >  1:865744/1‑150 (MQ=255)
                                              tctcTGGACTGCGATTAATCGTAGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACAGTTTTCTGCTAATGCGGTTACCGGAGAGaaaa                                                                                                        >  2:1077543/1‑150 (MQ=255)
                                               ctctGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGGTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAg                                                                                                       <  1:1317082/150‑1 (MQ=255)
                                                tctGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAgg                                                                                                      <  2:407668/150‑1 (MQ=255)
                                                tctGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAgg                                                                                                      <  1:217726/150‑1 (MQ=255)
                                                 ctGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAggg                                                                                                     >  2:406151/1‑150 (MQ=255)
                                                   ggACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAggggg                                                                                                   >  1:757596/1‑150 (MQ=255)
                                                     aCTGCGATTAATCGGTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGTGGTTACCGGAGAGAAAAGGGGGTg                                                                                                 <  2:761341/150‑1 (MQ=255)
                                                      cTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGa                                                                                                >  1:1448299/1‑150 (MQ=255)
                                                        gCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGTTGAgg                                                                                              >  2:1195575/1‑150 (MQ=255)
                                                        gCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAgg                                                                                              <  2:309056/150‑1 (MQ=255)
                                                         cGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTTGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGTTTACCGGAGAGAAAAGGGGGTGAGGt                                                                                             <  1:1201493/150‑1 (MQ=255)
                                                         cGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGt                                                                                             >  1:975713/1‑150 (MQ=255)
                                                             tAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACa                                                                                         <  2:980991/150‑1 (MQ=255)
                                                              aaTCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGGGTTTCCCGAGAGAAAAAGGGGTGCGGtttaatt                                                                                       >  2:162423/1‑145 (MQ=255)
                                                                tCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTa                                                                                      <  2:975713/150‑1 (MQ=255)
                                                                  gTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGGTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACg                                                                                    <  2:864718/150‑1 (MQ=255)
                                                                   ttGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCCCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGt                                                                                   <  1:138613/150‑1 (MQ=255)
                                                                    tGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTa                                                                                  >  2:1482333/1‑150 (MQ=255)
                                                                        ggAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCCCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCAGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCg                                                                              <  2:259569/150‑1 (MQ=255)
                                                                             gtgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCt                                                                         <  1:738382/150‑1 (MQ=255)
                                                                              tgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtg                                                                        >  1:466383/1‑150 (MQ=255)
                                                                              tgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGCGGTGAGGTTACATTACGTATTCGTACCtg                                                                        >  2:577619/1‑150 (MQ=255)
                                                                                cttcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtggt                                                                      >  1:1308687/1‑150 (MQ=255)
                                                                                  tcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCTCGTTGTTAATTTTAGCAACGATGATAAACATTTTTCTGCTTTGGCGGTTAAACGATAGAAAATGGGGTGAGGTTACATTACGTATTCGTACCtggtgg                                                                    >  1:1184901/1‑150 (MQ=255)
                                                                                  tcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATTGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtggtgg                                                                    >  1:85172/1‑150 (MQ=255)
                                                                                  tcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAGCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtggtgg                                                                    >  1:291090/1‑150 (MQ=255)
                                                                                  tcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtggtgg                                                                    >  2:1339731/1‑150 (MQ=255)
                                                                                     tATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGc                                                                 >  1:1208812/1‑150 (MQ=255)
                                                                                         aGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCg                                                             >  2:1098894/1‑150 (MQ=255)
                                                                                         aGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGGTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCg                                                             >  1:858520/1‑150 (MQ=255)
                                                                                             cTGAAATATAATGCTGGAAAACTCTCTGGCAACACTTTTCGCGTTGTTAATTTAAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGTGCATCGTCAt                                                         <  1:1591296/150‑1 (MQ=255)
                                                                                             cTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCAt                                                         >  1:849073/1‑150 (MQ=255)
                                                                                                aaaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTg                                                      <  2:1371573/150‑1 (MQ=255)
                                                                                                  atatAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGtt                                                    <  2:573455/150‑1 (MQ=255)
                                                                                                     taATGCTGGAAAAATCTCTTGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGttttt                                                 <  2:822508/150‑1 (MQ=255)
                                                                                                      aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTc                                                <  1:183302/150‑1 (MQ=255)
                                                                                                         gCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAAt                                             <  1:1023215/150‑1 (MQ=255)
                                                                                                          cTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATg                                            >  1:555237/1‑150 (MQ=255)
                                                                                                          cTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATg                                            >  2:1412888/1‑150 (MQ=255)
                                                                                                          cTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATg                                            >  1:694458/1‑150 (MQ=255)
                                                                                                           tGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGGTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATgg                                           <  1:1387424/150‑1 (MQ=255)
                                                                                                            ggAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGt                                          >  2:1341965/1‑150 (MQ=255)
                                                                                                             gAAAAATCTCTGGCCACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGGGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTg                                         <  2:1101937/150‑1 (MQ=255)
                                                                                                             gAAAAATCTCTGGCAACACTTTCGCGTTGATAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAACAGTGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGa                                        <  2:727735/150‑1 (MQ=255)
                                                                                                                  aTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTc                                    >  1:661680/1‑150 (MQ=255)
                                                                                                                   tctcTGGCTAGACTTTACGCGTTGTTAATTTTATCAACGCTGATAAAGCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTAATCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCa                                   <  1:559357/150‑1 (MQ=255)
                                                                                                                     actGGCAACACTTTTCGCGTTGTTAATTTTAGCCACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTGACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAgg                                 <  2:1198443/149‑1 (MQ=255)
                                                                                                                        ggCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTcc                              <  2:1479501/150‑1 (MQ=255)
                                                                                                                         gCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCa                             <  1:406151/150‑1 (MQ=255)
                                                                                                                             cacTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCGACGCt                         >  1:200414/1‑150 (MQ=255)
                                                                                                                                 tttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAt                     <  2:613065/150‑1 (MQ=255)
                                                                                                                                 atcagcgTTGCTAAAATTAACAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAt                     >  1:706546/5‑150 (MQ=255)
                                                                                                                                   cagcgTTGCTAAAATTAACAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtga                   >  2:1376506/3‑150 (MQ=255)
                                                                                                                                        ttgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttg              <  2:441333/150‑1 (MQ=255)
                                                                                                                                         tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttga             >  1:937957/1‑150 (MQ=255)
                                                                                                                                         tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttga             <  2:1339781/150‑1 (MQ=255)
                                                                                                                                             aaTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATgacgttgacgtt         >  2:1401740/1‑150 (MQ=255)
                                                                                                                                              aTTTTAGCAACGCTGATAAACCGTATTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGttt        <  2:1311004/150‑1 (MQ=255)
                                                                                                                                                  tAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTg    <  2:1526726/150‑1 (MQ=255)
                                                                                                                                                    gCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTgcg  <  1:710568/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCG  >  NC_000913/1478123‑1478420

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: