Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,539,659 A→C 6.2% S759R (AGC→CGC)  yhgF → putative transcriptional accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,539,6590AC 6.2% 1872.3 / 4.7 607S759R (AGC→CGC) yhgFputative transcriptional accessory protein
Reads supporting (aligned to +/- strand):  ref base A (291/254);  new base C (2/56);  total (293/314)
Fisher's exact test for biased strand distribution p-value = 4.14e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 615 total aligned reads displayed.

ATCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGC  >  NC_000913/3539510‑3539808
                                                                                                                                                     |                                                                                                                                                     
    tCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAgcgc                                                                                                                                                   >  1:6010329/1‑150 (MQ=255)
    tCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAgcgc                                                                                                                                                   <  2:4916590/150‑1 (MQ=255)
       gCGTAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCgat                                                                                                                                                <  2:7740300/150‑1 (MQ=255)
       gCGTAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCgat                                                                                                                                                <  2:133147/150‑1 (MQ=255)
       gCGTAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCgat                                                                                                                                                >  2:3208402/1‑150 (MQ=255)
        cGTAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCgatg                                                                                                                                               >  2:7503844/1‑150 (MQ=255)
         gTAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCgatga                                                                                                                                              <  1:7809052/150‑1 (MQ=255)
          tAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCgatgat                                                                                                                                             <  1:5100834/150‑1 (MQ=255)
          tAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCgatgat                                                                                                                                             <  1:3759856/150‑1 (MQ=255)
           aaaCGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGatgatg                                                                                                                                            <  2:1435828/150‑1 (MQ=255)
            aaCGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACTCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAATCAACCGCCCGGCAGCCAAACCACGCTGTCGTGAAGCGCAGACTGCCTGTAATAGCGCGATgatgg                                                                                                                                           <  1:7226016/150‑1 (MQ=255)
            aaCGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATgatgg                                                                                                                                           <  2:1098233/150‑1 (MQ=255)
             aCGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATgatgga                                                                                                                                          <  1:8964951/150‑1 (MQ=255)
              cGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATgatggat                                                                                                                                         <  1:7588197/150‑1 (MQ=255)
              cGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATgatggat                                                                                                                                         <  1:1595703/150‑1 (MQ=255)
                tATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATgc                                                                                                                                       >  2:176725/1‑150 (MQ=255)
                  tCGCCCTGACTCTGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATgcgc                                                                                                                                     <  2:4705085/150‑1 (MQ=255)
                  tCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATgcgc                                                                                                                                     <  1:2768399/150‑1 (MQ=255)
                     cccTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTgg                                                                                                                                  >  2:7352327/1‑150 (MQ=255)
                     cccTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGAGGATGGATGCGCTgg                                                                                                                                  >  2:1901139/1‑150 (MQ=255)
                      ccTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTggc                                                                                                                                 <  2:6756884/150‑1 (MQ=255)
                           cTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTggcggcgg                                                                                                                            <  2:7896477/150‑1 (MQ=255)
                            tATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTggcggcggc                                                                                                                           <  2:3037513/150‑1 (MQ=255)
                            tATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTggcggcggc                                                                                                                           >  2:5632538/1‑150 (MQ=255)
                            tATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGGCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGCTGCGCTggcggcggc                                                                                                                           >  2:5044092/1‑150 (MQ=255)
                                cgcCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCCAACCACGCGGTCGTGCAGCGCAGCCTGCCGGTAATAGCGCGATGAGGGATGCGCTGGCGGCGGCacgg                                                                                                                       >  1:5188206/1‑147 (MQ=255)
                                 gcCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATgg                                                                                                                      <  2:313496/150‑1 (MQ=255)
                                  ccTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATggg                                                                                                                     >  2:9078998/1‑150 (MQ=255)
                                   cTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGACAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGc                                                                                                                    <  1:437804/150‑1 (MQ=255)
                                    tGGATGAGCAGCCTGGCGACACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCa                                                                                                                   <  1:7359785/150‑1 (MQ=255)
                                      gATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCaaa                                                                                                                 >  2:4499451/1‑150 (MQ=255)
                                          agcagcCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAAc                                                                                                             <  2:5697849/150‑1 (MQ=255)
                                            cagcCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAGACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGt                                                                                                           >  1:2327039/1‑150 (MQ=255)
                                            cagcCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGt                                                                                                           >  2:3480011/1‑150 (MQ=255)
                                              gcCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCACAAAACGTTa                                                                                                         >  1:8355263/1‑150 (MQ=255)
                                              gcCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTa                                                                                                         <  1:8856863/150‑1 (MQ=255)
                                               ccTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCCAAAAACGTTaa                                                                                                        >  2:8943882/1‑150 (MQ=255)
                                               ccTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTaa                                                                                                        <  1:3239941/150‑1 (MQ=255)
                                               ccTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTaa                                                                                                        >  2:8227665/1‑150 (MQ=255)
                                                 tGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAAc                                                                                                      >  2:3809110/1‑150 (MQ=255)
                                                   gCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCACAAAACGTTAAACGc                                                                                                    >  1:1146013/1‑150 (MQ=255)
                                                   gCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGc                                                                                                    >  1:40062/1‑150 (MQ=255)
                                                     gAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGTGGCAATGGGCAAAAAACGTTAAACGccc                                                                                                  <  2:8627171/150‑1 (MQ=255)
                                                     gAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGccc                                                                                                  <  2:296705/150‑1 (MQ=255)
                                                        aCCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTa                                                                                               >  1:5808351/1‑150 (MQ=255)
                                                          cAACGCTCGTCGCGGCGGCGTTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCCGCCTGCCGGTAATCGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACt                                                                                             <  2:1259011/150‑1 (MQ=255)
                                                           aaCGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTg                                                                                            >  1:1346456/1‑150 (MQ=255)
                                                              gCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGcc                                                                                         >  2:7625682/1‑150 (MQ=255)
                                                                tcgtcgCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCACAAAACGGGAAACGCCCGTACTGGCCtg                                                                                       >  1:1866910/1‑149 (MQ=255)
                                                                  gtcgCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATATCGCGATGTTGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACg                                                                                     >  1:2944879/1‑150 (MQ=255)
                                                                  gtcgCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACg                                                                                     >  1:335847/1‑150 (MQ=255)
                                                                  gtcgCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACg                                                                                     <  2:8386662/150‑1 (MQ=255)
                                                                  gtcgCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACg                                                                                     <  2:7205749/150‑1 (MQ=255)
                                                                      cggcggcggTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGGTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTc                                                                                 <  2:6010329/150‑1 (MQ=255)
                                                                         cggcggTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGaa                                                                              <  2:616131/150‑1 (MQ=255)
                                                                             ggTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTg                                                                          >  1:4992307/1‑150 (MQ=255)
                                                                             ggTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGAAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTg                                                                          <  2:2098236/150‑1 (MQ=255)
                                                                               tAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCa                                                                        <  2:4357640/150‑1 (MQ=255)
                                                                                   gAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATggg                                                                                                                     <  1:5418523/101‑1 (MQ=255)
                                                                                   gAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATggg                                                                                                                     >  2:5418523/1‑101 (MQ=255)
                                                                                   gAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGaa                                                                    >  2:2637593/1‑150 (MQ=255)
                                                                                     aCGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAAt                                                                  >  2:8868381/1‑150 (MQ=255)
                                                                                        cccGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGt                                                               >  1:6398333/1‑150 (MQ=255)
                                                                                          cGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAg                                                             >  1:8695907/1‑150 (MQ=255)
                                                                                                     cGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAgg                                                  >  2:360359/1‑150 (MQ=255)
                                                                                                     cGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAgg                                                  >  1:810745/1‑150 (MQ=255)
                                                                                                         cGGCAGCCAAACCACGGGGTTGTGAAGCGCCGCCTGCCGGTAATAGCGCTATTAGTGATGCGCTGGCGGCGGCTTTGGGCAAAAAACGTTACACGCACGTACTGTCATACGGTTCGAATTTGGACGAAATCGTAGGGCAGATAAGGCGtt                                              <  2:5504424/150‑1 (MQ=255)
                                                                                                          ggCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTc                                             >  2:1652999/1‑150 (MQ=255)
                                                                                                             aGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCAcg                                          >  1:8513971/1‑150 (MQ=255)
                                                                                                              gCCAACCCCCGCGGTCGTGACGCGCCGCCTGCGGGTAATCGCGCGTTGATGGTGGCGCGGGGGGCGGCATTGGGCAAAAAACGTTACACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCAcgc                                         <  2:4144865/150‑1 (MQ=255)
                                                                                                              gCCAAACCCCGCGGTCGTGACGCGCAGCCTGCCGGTAATTGCGCGTTGATGGATGCGCTGGCGGCGGAATTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCAcgc                                         <  1:8223135/150‑1 (MQ=255)
                                                                                                              gCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCATTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCAcgc                                         <  2:2408335/150‑1 (MQ=255)
                                                                                                              gCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCGGGCGGCGGCATTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCAcgc                                         <  1:2189106/150‑1 (MQ=255)
                                                                                                               ccAAACCACGGGGTCGTGACGCGCAGCCTGCCGGTAATCGCGCGTTGATGGAGGCGCGGGCGGCGGCATTGGGAAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCAcgcc                                        <  1:4229628/150‑1 (MQ=255)
                                                                                                                cAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCAcgccg                                       >  1:1047683/1‑150 (MQ=255)
                                                                                                                cAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCAcgccg                                       >  1:6213316/1‑150 (MQ=255)
                                                                                                                 aaaCCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCAcgccgc                                      >  1:6967444/1‑150 (MQ=255)
                                                                                                                     cACGCGGTCGTGAAGCGCAGCCTGCGGGTAATAGCGCGATGATGGATGCGGTGGGGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCt                                  <  1:468646/150‑1 (MQ=255)
                                                                                                                      aCGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTg                                 <  2:5387306/150‑1 (MQ=255)
                                                                                                                       cggggTCGTGACGCGCCGCCTGCGGGTAATCGCGCGATGATGTTGGCGCTGGCGGCGGCATTGGCAAAAAAACTTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTgg                                <  1:3680269/147‑1 (MQ=255)
                                                                                                                       cgcgGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTgg                                <  2:5199754/150‑1 (MQ=255)
                                                                                                                        ccgGGCGTGAAGCGCCGCCTGCCGGTAATCGCGCGTTGATGGAGGCGCTGGGGGCGGCAATGGGCAACACACGTTACACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAATTCGTAGGGCAGATAAGGAGTTCACGCCGCATCTgcc                               <  2:8543616/149‑3 (MQ=255)
                                                                                                                         cgGTCGTGCCGCGCCGCCTGCCGGTACTCGCGCGAGGCTGGAGGCGCGGGGGGCGGCAAGGGGCAAAACACGTTACACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCa                              <  2:3338131/150‑1 (MQ=37)
                                                                                                                          ggTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaa                             >  1:1430927/1‑150 (MQ=255)
                                                                                                                             cGTGTAGCGCAGCCTGCCGGTAATAGCGCGATGATGTATGCGCTGGCGGCGGCATTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacga                          <  2:7040121/150‑1 (MQ=255)
                                                                                                                              gTGAAGCGCAGCCTGCCGGTAATCGCGCGATGATGGTTGCGCTGGCGGCGGCAAGGGGCAAAAAACGTAAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaa                         <  1:7920200/150‑1 (MQ=255)
                                                                                                                               tGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaac                        >  1:1526849/1‑150 (MQ=255)
                                                                                                                               tGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaac                        >  2:6412317/1‑150 (MQ=255)
                                                                                                                                gAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaacg                       >  1:3070448/1‑150 (MQ=255)
                                                                                                                                   gcgcAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCt                    >  2:8188873/1‑150 (MQ=255)
                                                                                                                                    cgcAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACTCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCtt                   <  2:831233/150‑1 (MQ=255)
                                                                                                                                        gCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAACCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCt               <  2:6658308/150‑1 (MQ=255)
                                                                                                                                          cTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGa             <  2:8118111/150‑1 (MQ=255)
                                                                                                                                          cTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGAAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGa             <  1:5546839/150‑1 (MQ=255)
                                                                                                                                            gCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACa           <  2:3698183/150‑1 (MQ=255)
                                                                                                                                                 tAATCGCGCGATGATGTAGGCGGTGGCGGCGGCAATGTGCAAAAAACGTCACACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaaa      <  1:941342/150‑1 (MQ=39)
                                                                                                                                                 tAATAGCGCGATGATGGATGCGCTGGCGGCGGCATTGTGCAAAAAACTTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaaa      <  1:6786897/150‑1 (MQ=255)
                                                                                                                                                    tAGCGCGATGATGGATGCGCTGGCGGCGGCACTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAAGTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg   >  2:2774258/1‑150 (MQ=255)
                                                                                                                                                    tAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg   >  2:8770262/1‑150 (MQ=255)
                                                                                                                                                    tAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg   >  2:4277240/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTGGCGAAACCAACGCTCGTCGCGGCGGCGGTAATGAACGCCCGCAAAACAACCGCCCGGCAGCCAAACCACGCGGTCGTGAAGCGCAGCCTGCCGGTAATAGCGCGATGATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGC  >  NC_000913/3539510‑3539808

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