Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 287,925 A→C 5.2% Q379H (CAA→CAC yagH → CP4‑6 prophage; putative xylosidase/arabinosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913287,9250AC 5.2% 950.2 / 5.5 336Q379H (CAA→CACyagHCP4‑6 prophage; putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base A (157/152);  new base C (2/25);  total (159/177)
Fisher's exact test for biased strand distribution p-value = 8.01e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 337 total aligned reads displayed.

CTGGCGGGACGATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGG  >  NC_000913/287778‑288074
                                                                                                                                                   |                                                                                                                                                     
   gCGGGACGATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGAc                                                                                                                                                  >  2:306375/1‑150 (MQ=255)
   gCGGGACGATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGAc                                                                                                                                                  >  2:1630884/1‑150 (MQ=255)
     gggACGATTGCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCAGTCCGTTCGCCGCCACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCATCGACTCGCTCAATTCGACCTTCCCCCAATCGACCg                                                                                                                                                <  2:2947246/150‑1 (MQ=255)
        aCGATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTgg                                                                                                                                             >  2:1487405/1‑150 (MQ=255)
           aTTTCGACGCCTTTTCGCTTGCCCCGGAGCTGAAGACCCTGCGCATTCCGTTCGACGACACCCCCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTAGGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGgcgc                                                                                                                                          <  2:1026644/150‑1 (MQ=255)
              tCGCCGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTCCGGCTCTAGGGCAACGACTCGCTCAATGCGACCTTCCCCCCCTCGACCGTGGCGcgcc                                                                                                                                       <  2:2894640/150‑1 (MQ=255)
                gACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGcgccgc                                                                                                                                     <  2:3511007/150‑1 (MQ=255)
                gACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGcgccgc                                                                                                                                     <  2:3022648/150‑1 (MQ=255)
                gACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGcgccgc                                                                                                                                     <  2:2716258/150‑1 (MQ=255)
                   gCCTGTTCGCTTTACCCGTAGCTGCAGACCCTGCGAATTCGTTACGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTCCGGCTCTTTGGCAACGACTCGCTAAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTgg                                                                                                                                  <  1:227051/150‑1 (MQ=255)
                    ccAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGgc                                                                                                                                 >  1:3523391/1‑150 (MQ=255)
                     cAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGgca                                                                                                                                <  2:3397571/150‑1 (MQ=255)
                        ttCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGgcagca                                                                                                                             >  2:2012717/1‑150 (MQ=255)
                        ttCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGgcagca                                                                                                                             >  1:240113/1‑150 (MQ=255)
                        ttCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGgcagca                                                                                                                             >  1:2560585/1‑150 (MQ=255)
                        ttCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGgcagca                                                                                                                             >  2:1572566/1‑150 (MQ=255)
                        ttCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGgcagca                                                                                                                             >  1:1043192/1‑150 (MQ=255)
                        ttCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGgcagca                                                                                                                             >  2:665167/1‑150 (MQ=255)
                          ccagttgcccGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACt                                                                                                                           >  2:723225/8‑150 (MQ=255)
                            cTTGACCGGGTGCTGCAGACCCTGCGCTTTCCGTTCGCCGACCCCCTCGGCTGGCTCACCGCGCGCCGGGGCTTCTTACGGCTCTATGGAAACGACTCGCTCAATTCGACCTTCACCCACTCGACCGTGGCGCGCCGCTGGCAGCACTTc                                                                                                                         <  1:2088737/150‑1 (MQ=255)
                             ttGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCg                                                                                                                        >  1:3495856/1‑150 (MQ=255)
                             ttGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCg                                                                                                                        >  1:985021/1‑150 (MQ=255)
                                   cGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTcc                                                                                                                  <  2:425580/150‑1 (MQ=255)
                                           cAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGaa                                                                                                          >  2:1773440/1‑150 (MQ=255)
                                             gACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAAc                                                                                                        >  2:2632273/1‑150 (MQ=255)
                                               cccTGCGCATTCCGTTCGACGACACCCTCGGCTCGCCCACCGCGCGCCCGGGCTTCTTCGGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCACTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAAcgc                                                                                                      <  1:3038155/150‑1 (MQ=255)
                                                ccTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAAcgcg                                                                                                     >  2:1912544/1‑150 (MQ=255)
                                                  tGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCACTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGa                                                                                                   <  2:248843/150‑1 (MQ=255)
                                                  tGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGa                                                                                                   >  2:1509773/1‑150 (MQ=255)
                                                  tGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGa                                                                                                   >  2:318695/1‑150 (MQ=255)
                                                      cATTCCGTTCGACGACACCCTCGGCTCGCTCCCCGCGCGCCCGGGCTTCTTACGGCTCTAGGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGa                                                                                               <  1:2795316/150‑1 (MQ=255)
                                                      cATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGa                                                                                               <  1:3601261/150‑1 (MQ=255)
                                                      cATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTATATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGa                                                                                               <  1:321813/150‑1 (MQ=255)
                                                      cATTCCGTTCGACGACACCCTCGGCTCGCCCACCGCGCGCCCGGGCTTCTTACGGCTCTAGGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGa                                                                                               <  1:3309515/150‑1 (MQ=255)
                                                       aTTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTACTTACGGATCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAg                                                                                              <  1:1909602/150‑1 (MQ=255)
                                                         tCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGtt                                                                                            <  2:3474021/150‑1 (MQ=255)
                                                            gTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTtctc                                                                                         <  1:595962/150‑1 (MQ=255)
                                                                gacgacACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCACTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCg                                                                                     <  1:1872748/150‑1 (MQ=255)
                                                                gacgacACCCTCGGCTCGCTCACCGCGCGCCCGGGATTCTTAGGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCg                                                                                     <  2:1897331/150‑1 (MQ=255)
                                                                 acgacACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCgg                                                                                    >  2:2874646/1‑150 (MQ=255)
                                                                  cgacACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGt                                                                                   >  2:1877417/1‑150 (MQ=255)
                                                                   gacACCCGCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTg                                                                                  <  2:805801/150‑1 (MQ=255)
                                                                      acCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCAc                                                                               <  1:2953507/150‑1 (MQ=255)
                                                                        ccTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACtt                                                                             >  1:2161158/1‑150 (MQ=255)
                                                                        ccTCGGCTCGCTCACCGCGCGCCCGGGCTTCTCACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACtt                                                                             >  1:3098446/1‑150 (MQ=255)
                                                                         cTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGCCCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACGTc                                                                            >  2:850322/1‑150 (MQ=255)
                                                                          tCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTcc                                                                           >  2:3746222/1‑150 (MQ=255)
                                                                          tCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTcc                                                                           >  2:2818283/1‑150 (MQ=255)
                                                                           cGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCca                                                                          >  1:715066/1‑150 (MQ=255)
                                                                               tcgctcACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCcagcag                                                                      >  2:449575/1‑150 (MQ=255)
                                                                                 gctcACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCagag                                                                    >  2:2431653/1‑150 (MQ=255)
                                                                                   tcACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                  >  2:3102082/1‑150 (MQ=255)
                                                                                       cgcgcgCCGGGGCTTCTTACGGCTCTAGGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCggg                                                              <  1:1417972/150‑1 (MQ=255)
                                                                                        gcgcgcCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCgggg                                                             <  1:2631225/150‑1 (MQ=255)
                                                                                        gcgcgcCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCgggg                                                             <  2:2568895/150‑1 (MQ=255)
                                                                                          gcgcCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCt                                                           <  1:3885918/150‑1 (MQ=255)
                                                                                          gcgcCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCt                                                           <  1:815095/150‑1 (MQ=255)
                                                                                               cGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                      <  2:3844776/150‑1 (MQ=255)
                                                                                                   cttcttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGcta                                                  <  1:2449076/150‑1 (MQ=255)
                                                                                                    ttcttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctac                                                 <  2:930444/150‑1 (MQ=255)
                                                                                                      cttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctacta                                               >  2:3927611/1‑150 (MQ=255)
                                                                                                      cttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctacta                                               <  2:3583337/150‑1 (MQ=255)
                                                                                                         aCGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTacaa                                            >  2:2214010/1‑150 (MQ=255)
                                                                                                         aCGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTacaa                                            <  1:1859815/150‑1 (MQ=255)
                                                                                                         aCGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTacaa                                            <  1:619392/150‑1 (MQ=255)
                                                                                                          cGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTacaac                                           <  1:2892232/150‑1 (MQ=255)
                                                                                                           ggCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTacaaca                                          <  1:18061/150‑1 (MQ=255)
                                                                                                           ggCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTacaaca                                          <  1:3562174/150‑1 (MQ=255)
                                                                                                                  tGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCaaaaa                                   >  1:2058269/1‑150 (MQ=255)
                                                                                                                   ggCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAAc                                  <  2:1209078/150‑1 (MQ=255)
                                                                                                                    gCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACt                                 >  2:3262278/1‑150 (MQ=255)
                                                                                                                    gCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACt                                 >  1:1725745/1‑150 (MQ=255)
                                                                                                                     cAACTACTCGCTCATTTCGACCTTCACCCACTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTg                                <  2:1678271/150‑1 (MQ=255)
                                                                                                                     cAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTg                                <  1:3771550/150‑1 (MQ=255)
                                                                                                                       aCGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGa                              <  1:1736329/150‑1 (MQ=255)
                                                                                                                       aCGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGa                              <  1:1404177/150‑1 (MQ=255)
                                                                                                                       aCGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGa                              >  1:1004804/1‑150 (MQ=255)
                                                                                                                          aCTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCt                           >  2:1061994/1‑150 (MQ=255)
                                                                                                                            tCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTAc                         >  1:3679140/1‑150 (MQ=255)
                                                                                                                              gCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGTTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTg                       <  1:2491389/150‑1 (MQ=255)
                                                                                                                              gCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTg                       <  2:2365525/150‑1 (MQ=255)
                                                                                                                                tCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCCAAAACTGGAGCTACTGCt                     >  1:3187508/1‑150 (MQ=255)
                                                                                                                                tCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCt                     >  1:168686/1‑150 (MQ=255)
                                                                                                                                tCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCt                     >  2:1352395/1‑150 (MQ=255)
                                                                                                                                tCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCt                     >  1:1213015/1‑150 (MQ=255)
                                                                                                                                  aaTTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCttt                   >  2:1566311/1‑150 (MQ=255)
                                                                                                                                   aTTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTtg                  <  1:770137/150‑1 (MQ=255)
                                                                                                                                     tCGACCTTCACCCACTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTtgtg                <  1:1530361/150‑1 (MQ=255)
                                                                                                                                     tCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTtgtg                <  1:3037475/150‑1 (MQ=255)
                                                                                                                                        aCCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGAc             >  2:2250566/1‑150 (MQ=255)
                                                                                                                                          cTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTa           >  1:1998115/1‑150 (MQ=255)
                                                                                                                                          cTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTa           <  2:582275/150‑1 (MQ=255)
                                                                                                                                              aCCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgag       <  1:3235959/150‑1 (MQ=255)
                                                                                                                                              aCCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgag       <  1:2158038/150‑1 (MQ=255)
                                                                                                                                               cccAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgagg      >  2:3765041/1‑150 (MQ=255)
                                                                                                                                               cccAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgagg      >  1:3117363/1‑150 (MQ=255)
                                                                                                                                                ccACTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgagga     <  2:1410210/150‑1 (MQ=255)
                                                                                                                                                 cAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgaggag    <  2:1137823/150‑1 (MQ=255)
                                                                                                                                                 cAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgaggag    <  1:499222/150‑1 (MQ=255)
                                                                                                                                                 cAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgaggag    <  1:150466/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGGCGGGACGATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGG  >  NC_000913/287778‑288074

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