Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 285,887 C→A 6.4% L165M (CTG→ATG)  yagG → CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913285,8870CA 6.4% 327.7 / 8.0 95L165M (CTG→ATG) yagGCP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base C (38/50);  new base A (4/2);  total (42/53)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGG  >  NC_000913/285739‑286032
                                                                                                                                                    |                                                                                                                                                 
tctcCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCgcgc                                                                                                                                                   <  2:402404/149‑1 (MQ=255)
 ctcCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCt                                                                                                                                                  <  4:8103/149‑1 (MQ=255)
  tcCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTg                                                                                                                                                 >  1:765250/1‑149 (MQ=255)
   ccTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGc                                                                                                                                                <  1:354774/149‑1 (MQ=255)
       aCCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCAGCGCGATGCCGCt                                                                                                                                            <  2:690554/149‑1 (MQ=255)
         ccTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTgg                                                                                                                                          >  3:697747/1‑149 (MQ=255)
         ccTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCAGCGCGATGCCGCTga                                                                                                                                          <  4:374425/149‑2 (MQ=255)
            ggTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGa                                                                                                                                       <  2:259998/149‑1 (MQ=255)
            ggTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGa                                                                                                                                       <  3:286619/149‑1 (MQ=255)
              tCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGc                                                                                                                                     >  4:251666/1‑149 (MQ=255)
               cTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGca                                                                                                                                    <  4:598891/149‑1 (MQ=255)
                tACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGcat                                                                                                                                   >  3:669696/1‑149 (MQ=255)
                tACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGcat                                                                                                                                   >  2:366065/1‑149 (MQ=255)
                   acCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGcatcat                                                                                                                                >  3:690570/1‑149 (MQ=255)
                   acCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCGCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGcatcat                                                                                                                                >  3:498540/1‑149 (MQ=255)
                    ccTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCatcatc                                                                                                                               >  4:397067/1‑149 (MQ=255)
                     cTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCg                                                                                                                              >  1:752659/1‑149 (MQ=255)
                     cTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCg                                                                                                                              <  4:84680/149‑1 (MQ=255)
                      ttCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCgg                                                                                                                             >  2:728703/1‑149 (MQ=255)
                         gTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCaa                                                                                                                          >  1:236690/1‑149 (MQ=255)
                          ttAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCaaa                                                                                                                         >  2:454863/1‑149 (MQ=255)
                          ttAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCaaa                                                                                                                         >  4:281074/1‑149 (MQ=255)
                          ttAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCaaa                                                                                                                         >  3:599152/1‑149 (MQ=255)
                          ttAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCaaa                                                                                                                         <  1:366065/149‑1 (MQ=255)
                               gTGCCGTACTGCGCCATGCCGGGCGTCATCGCCGCCGACCCGAAAGAGCGTCACGCACTGCAGTCCTGGCGCTTCTTCCTGGCGGAGGCGGGCTCGATCGCTATCAGCGGAATCGCGTTGCCGCTGGTGAGCATCATCGGCAAAGGGGa                                                                                                                    <  2:306851/149‑1 (MQ=255)
                                tGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGAc                                                                                                                   <  1:71864/149‑1 (MQ=255)
                                tGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGAc                                                                                                                   <  4:38693/149‑1 (MQ=255)
                                  ccGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGa                                                                                                                 <  1:541585/149‑1 (MQ=255)
                                   cGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTAGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAg                                                                                                                <  1:350460/149‑1 (MQ=255)
                                    gTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGc                                                                                                               <  1:752566/149‑1 (MQ=255)
                                        tGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGt                                                                                                           <  2:737250/149‑1 (MQ=255)
                                        tGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGt                                                                                                           <  1:64178/149‑1 (MQ=255)
                                         gcgcCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTg                                                                                                          <  3:329140/149‑1 (MQ=255)
                                          cgcCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTgg                                                                                                         <  2:393227/149‑1 (MQ=255)
                                          cgcCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTgg                                                                                                         >  1:593783/1‑149 (MQ=255)
                                              aTGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTa                                                                                                     >  1:22142/1‑149 (MQ=255)
                                                 ccGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGc                                                                                                       <  3:732728/144‑1 (MQ=255)
                                                 ccGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGc                                                                                                       >  4:732728/1‑144 (MQ=255)
                                                  cGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTc                                                                                                 >  1:525384/1‑149 (MQ=255)
                                                     gCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGgc                                                                                              <  1:191887/149‑1 (MQ=255)
                                                        tcatcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGcc                                                                                           >  2:706455/1‑149 (MQ=255)
                                                        tcatcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGcc                                                                                           <  2:593783/149‑1 (MQ=255)
                                                            caCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAtgt                                                                                       >  1:410696/1‑149 (MQ=255)
                                                                 ccgACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATgtgcgtg                                                                                  <  1:620823/149‑1 (MQ=255)
                                                                     cccGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGCGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTgg                                                                              <  4:439322/149‑1 (MQ=255)
                                                                      ccGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAg                                                                                                             >  4:95088/1‑117 (MQ=255)
                                                                      ccGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAg                                                                                                             <  3:95088/117‑1 (MQ=255)
                                                                      ccGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTggg                                                                             >  4:603067/1‑149 (MQ=255)
                                                                       cGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTgggg                                                                            <  1:334136/149‑1 (MQ=255)
                                                                        gAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCAACATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGc                                                                           >  3:279668/1‑149 (MQ=255)
                                                                           agagCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGa                                                                        >  4:116984/1‑149 (MQ=255)
                                                                                gTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAgcggc                                                                   <  4:57016/149‑1 (MQ=255)
                                                                                    cGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCgtgg                                                               <  4:599152/149‑1 (MQ=255)
                                                                                     gCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCgtggt                                                              <  2:435093/149‑1 (MQ=255)
                                                                                     gCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCAGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCACGTGCGTGCTGGGGCTGAGCGGCgtggt                                                              <  3:397067/149‑1 (MQ=255)
                                                                                        cTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTgg                                                                              <  4:29305/130‑1 (MQ=255)
                                                                                        cTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTgg                                                                              >  3:29305/1‑130 (MQ=255)
                                                                                         tGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgctg                                                          >  4:639563/1‑149 (MQ=255)
                                                                                          gCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgctgc                                                         >  4:213043/1‑149 (MQ=255)
                                                                                          gCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgctgc                                                         >  2:109732/1‑149 (MQ=255)
                                                                                             gTCCTGGCGCGTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGctct                                                      <  1:208566/149‑1 (MQ=255)
                                                                                               ccTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGAGGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTAc                                                    <  2:609104/149‑1 (MQ=255)
                                                                                                 tGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGt                                                  >  3:467961/1‑149 (MQ=255)
                                                                                                  ggCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTc                                                 <  2:22142/149‑1 (MQ=255)
                                                                                                   gcgcTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCt                                                <  1:725523/149‑1 (MQ=255)
                                                                                                     gcTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGc                                              <  4:748962/149‑1 (MQ=255)
                                                                                                     gcTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGc                                              >  1:185943/1‑149 (MQ=255)
                                                                                                        tcttcCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGcttc                                           <  2:236690/149‑1 (MQ=255)
                                                                                                             cTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCAc                                      <  4:669696/149‑1 (MQ=255)
                                                                                                             cTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCAc                                      <  1:236089/149‑1 (MQ=255)
                                                                                                             cTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCAc                                      <  4:446169/149‑1 (MQ=255)
                                                                                                               ggcggcggcGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGa                                    <  2:765250/149‑1 (MQ=255)
                                                                                                                          gctcgctcGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAgcgc                         <  2:419707/149‑1 (MQ=255)
                                                                                                                            tcgctcgcTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTa                       >  3:419120/1‑149 (MQ=255)
                                                                                                                            tcgctcgcTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTa                       >  1:174073/1‑149 (MQ=255)
                                                                                                                                    tATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTg               >  2:568658/1‑149 (MQ=255)
                                                                                                                                    tATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTg               <  3:639563/149‑1 (MQ=255)
                                                                                                                                     aTCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGa              <  1:706455/149‑1 (MQ=255)
                                                                                                                                      tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGt                                                                                                           <  4:257239/55‑1 (MQ=255)
                                                                                                                                      tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGt                                                                                                           >  3:257239/1‑55 (MQ=255)
                                                                                                                                      tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGc                                              <  2:535844/116‑1 (MQ=255)
                                                                                                                                      tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGc                                              >  1:535844/1‑116 (MQ=255)
                                                                                                                                      tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTagag             <  2:227929/149‑1 (MQ=255)
                                                                                                                                      tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAg             <  4:145606/149‑1 (MQ=255)
                                                                                                                                      tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAg             <  1:568658/149‑1 (MQ=255)
                                                                                                                                      tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAg             <  2:301302/149‑1 (MQ=255)
                                                                                                                                      tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAg             <  3:292916/149‑1 (MQ=255)
                                                                                                                                      ccAGCGGCAGCGCGATGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAg             >  2:350645/2‑149 (MQ=255)
                                                                                                                                       cAGCGGCAGCGCGATGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAgg            >  4:302963/1‑149 (MQ=255)
                                                                                                                                        aGCGGCAGCGCGATGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGt           >  2:97014/1‑149 (MQ=255)
                                                                                                                                         gCGGCAGCGCGATGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTg          >  2:337308/1‑149 (MQ=255)
                                                                                                                                            gCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAg       >  1:262985/1‑149 (MQ=255)
                                                                                                                                            gCATCGCGCTGCCGCTGGTGAGCAACATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAg       <  4:279668/149‑1 (MQ=255)
                                                                                                                                             cATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGc      >  4:42338/1‑149 (MQ=255)
                                                                                                                                               tCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCttt                >  1:385302/1‑137 (MQ=255)
                                                                                                                                               tCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCttt                <  2:385302/137‑1 (MQ=255)
                                                                                                                                                 gcgcTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCggg  >  2:297516/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGG  >  NC_000913/285739‑286032

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