Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,478,266 T→A 5.7% I14N (ATT→AAT) ‡ ynbB → putative CDP‑diglyceride synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,478,2660TA 5.7% 262.0 / 4.4 87I14N (ATT→AAT) ‡ynbBputative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base T (41/41);  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 = 5.39e-01

ACAGTGGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTT  >  NC_000913/1478118‑1478414
                                                                                                                                                    |                                                                                                                                                    
acaGTGGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAAtt                                                                                                                                                     >  2:217726/1‑150 (MQ=255)
   gTGGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGATTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTa                                                                                                                                                  <  1:1161750/150‑1 (MQ=255)
     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:423632/1‑150 (MQ=255)
                             ggAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCt                                                                                                                        <  1:499902/150‑1 (MQ=255)
                              gAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCtt                                                                                                                       >  1:1526726/1‑149 (MQ=255)
                              gAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTa                                                                                                                       >  1:162527/1‑150 (MQ=255)
                              gAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTa                                                                                                                       >  2:1023215/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                                                                                                <  1:217726/150‑1 (MQ=255)
                                                     tctGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAgg                                                                                                <  2:407668/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                                      >  2:1412888/1‑150 (MQ=255)
                                                                                                               cTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATg                                      >  1:555237/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)
                                                                                                                                                    |                                                                                                                                                    
ACAGTGGATGCAGTGGAATAATCTTTTATGGAGCATTGCTTCAATTCTGCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTT  >  NC_000913/1478118‑1478414

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: