Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,923,978 C→G 9.6% E46Q (GAA→CAA)  yqcC ← DUF446 domain‑containing protein YqcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,923,9780CG 9.6% 458.0 / 26.6 143E46Q (GAA→CAA) yqcCDUF446 domain‑containing protein YqcC
Reads supporting (aligned to +/- strand):  ref base C (60/69);  new base G (3/11);  total (63/80)
Fisher's exact test for biased strand distribution p-value = 9.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 149 total aligned reads displayed.

GCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAAT  >  NC_000913/2923830‑2924127
                                                                                                                                                    |                                                                                                                                                     
gCCAGAATAATCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTc                                                                                                                                                       <  2:627442/149‑1 (MQ=255)
gCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTcc                                                                                                                                                      >  1:2246635/1‑150 (MQ=255)
gCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTcc                                                                                                                                                      >  2:467394/1‑150 (MQ=255)
  cAGAATAAGCGCGCGTTGCGGATGGTCCGTCTCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAAACACTGGAGTGGTTGCag                                                                                                                                                    <  2:133702/150‑2 (MQ=255)
  cAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGTTGCag                                                                                                                                                    <  2:133695/150‑2 (MQ=255)
   aGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGTTGCagt                                                                                                                                                   <  1:36906/150‑3 (MQ=255)
   aGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATg                                                                                                                                                   <  1:1898212/150‑1 (MQ=255)
    gAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGTTGCagtg                                                                                                                                                  <  1:2289328/150‑4 (MQ=255)
    gAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATgg                                                                                                                                                  >  2:2469823/1‑150 (MQ=255)
    gAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATgg                                                                                                                                                  <  1:217244/150‑1 (MQ=255)
    gAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATgg                                                                                                                                                  <  1:503048/150‑1 (MQ=255)
     aaTAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGTTGCagtgg                                                                                                                                                 <  2:2021778/150‑5 (MQ=255)
      aTAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGgtg                                                                                                                                                >  2:347561/1‑150 (MQ=255)
       tAAGCGCGCGTTGCTTATGGTCGGTCTCCAGCGCCATTTCATAATAGGGTGCAACGGAAAAAGCGCCCGGTAACGGCTTTTTTTTATCTTTTAATTGGTGGATTCGATGGATCAGGAGCCAGTTCAACCAATCCAGTGGTTCCATGgtgt                                                                                                                                               <  2:512355/150‑1 (MQ=255)
         aGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTcc                                                                                                                                             >  1:1430129/1‑150 (MQ=255)
         aGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTcc                                                                                                                                             >  1:884148/1‑150 (MQ=255)
             cgcgTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATaa                                                                                                                                         <  1:1974780/150‑1 (MQ=255)
                 ttGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGaa                                                                                                                                     >  1:303769/1‑150 (MQ=255)
                  tGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAAc                                                                                                                                    <  2:1402156/150‑1 (MQ=255)
                       aTGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTg                                                                                                                               <  2:1806873/150‑1 (MQ=255)
                        tGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTAACTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTgg                                                                                                                              >  1:2103760/1‑150 (MQ=255)
                         ggTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATTGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTggg                                                                                                                             <  1:1677150/150‑1 (MQ=255)
                         ggTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTggg                                                                                                                             <  2:489962/150‑1 (MQ=255)
                         ggTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTggg                                                                                                                             <  1:2469823/150‑1 (MQ=255)
                          gTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGt                                                                                                                            >  2:1594341/1‑150 (MQ=255)
                           tCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTa                                                                                                                           >  1:1536728/1‑150 (MQ=255)
                            ccGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTAc                                                                                                                          >  2:85535/1‑150 (MQ=255)
                              gTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTa                                                                                                                        >  1:169704/1‑150 (MQ=255)
                               tCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTAt                                                                                                                       <  1:1513449/150‑1 (MQ=255)
                                   cAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTaaa                                                                                                                   >  2:1544447/1‑150 (MQ=255)
                                   cAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTaaa                                                                                                                   >  2:173709/1‑150 (MQ=255)
                                    aGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAAt                                                                                                                  <  1:1270353/150‑1 (MQ=255)
                                          aTTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGc                                                                                                            <  2:289819/150‑1 (MQ=255)
                                           tttCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCg                                                                                                           >  1:2131748/1‑150 (MQ=255)
                                              cATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCt                                                                                                        >  2:2430614/1‑150 (MQ=255)
                                                taataGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTgg                                                                                                      >  1:2122887/1‑150 (MQ=255)
                                                   taGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGt                                                                                                   <  2:2229521/150‑1 (MQ=255)
                                                   taGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGt                                                                                                   <  1:1085948/150‑1 (MQ=255)
                                                    aGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACGAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGtt                                                                                                  >  1:1348876/1‑150 (MQ=255)
                                                      ggTGCAACGGCAAAAGCGCCCGGTAACGGCGGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCGATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCa                                                                                                <  2:1598893/150‑1 (MQ=255)
                                                         gCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCg                                                                                             <  1:1594341/150‑1 (MQ=255)
                                                           aaCGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGtt                                                                                           <  2:1097216/150‑1 (MQ=255)
                                                            aCGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGttt                                                                                          <  1:2137172/150‑1 (MQ=255)
                                                              ggCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCTGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCg                                                                                        <  2:1018621/150‑1 (MQ=255)
                                                              ggCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCg                                                                                        >  1:1614530/1‑150 (MQ=255)
                                                               gCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGc                                                                                       >  1:2276053/1‑150 (MQ=255)
                                                               gCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGc                                                                                       >  2:1092380/1‑150 (MQ=255)
                                                                 aaaaGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCa                                                                                     >  1:991782/1‑150 (MQ=255)
                                                                  aaaGCGCCCGGTAACGGCTGGTTGTTATCTAACAAATCGTGCATTCGCGGGATCACCACCCACTGCAACCACTCCAGTGGTTCCATGGTTTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAg                                                                                    <  2:717626/150‑1 (MQ=255)
                                                                    aGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATAAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTACATAAAGAACGGTTGGGTACTATTAAATTGATCCGGCTGGGGTTCATCGTTTCGCCAGTg                                                                                  >  2:2158929/1‑150 (MQ=255)
                                                                     gcgcCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGc                                                                                 <  1:1452330/150‑1 (MQ=255)
                                                                        cccGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGa                                                                              <  1:2193576/150‑1 (MQ=255)
                                                                             tAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGTCTGGGGTTCATCGTTTCGCCAGTGGTGATGTTc                                                                         <  2:2030629/150‑1 (MQ=255)
                                                                               aCGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCAc                                                                       <  2:1536728/150‑1 (MQ=255)
                                                                                cGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACg                                                                      <  1:13763/150‑1 (MQ=255)
                                                                                 ggCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGGTCATCGTTTCGCCAGTGCTGATGTTCACGc                                                                     <  2:2331284/150‑1 (MQ=255)
                                                                                   cgatttgttAGATAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAg                                                                   >  1:1892898/4‑150 (MQ=255)
                                                                                    tgtttgttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGt                                                                  >  2:678691/1‑150 (MQ=255)
                                                                                    tgtttgttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGt                                                                  >  1:1152402/1‑150 (MQ=255)
                                                                                         gttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGc                                                             <  1:1992867/150‑1 (MQ=255)
                                                                                             tCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCa                                                         >  2:1225530/1‑150 (MQ=255)
                                                                                             tCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCa                                                         <  2:2123356/150‑1 (MQ=255)
                                                                                              cTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTAGTATTAAATTGATGCGGCTGGGGTTCATCGTTTCTCCAGTGCTTATGTTCACGCAGTAACGCTTCaa                                                        <  1:2022459/150‑1 (MQ=255)
                                                                                               tAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAg                                                       >  2:2306664/1‑150 (MQ=255)
                                                                                                  cAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAgcgc                                                    >  2:1135435/1‑150 (MQ=255)
                                                                                                   aaaTCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGcc                                                   <  2:303769/150‑1 (MQ=255)
                                                                                                   aaaTCGTGCATTCGCGGGAGCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGGCCATAAAGAACGGGTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACTCTTCAAGCGcc                                                   <  1:1016228/150‑1 (MQ=255)
                                                                                                    aaTCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCt                                                  >  1:1622073/1‑150 (MQ=255)
                                                                                                      tCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAATTTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGc                                                <  1:815581/150‑1 (MQ=255)
                                                                                                       cGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCa                                               >  1:246281/1‑150 (MQ=255)
                                                                                                            aTTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTgg                                          >  1:591214/1‑150 (MQ=255)
                                                                                                             ttCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTGCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGga                                         <  2:169704/150‑1 (MQ=255)
                                                                                                               cgcgGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGgaga                                       >  1:275109/1‑150 (MQ=255)
                                                                                                                 cgGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACg                                     <  1:1983581/150‑1 (MQ=255)
                                                                                                                  gggATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGa                                    >  1:1647006/1‑150 (MQ=255)
                                                                                                                     aTCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAAcg                                 >  1:524174/1‑150 (MQ=255)
                                                                                                                        aGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCgg                              >  1:382343/1‑150 (MQ=255)
                                                                                                                        aGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCgg                              <  2:354899/150‑1 (MQ=255)
                                                                                                                          cACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTc                            >  1:664562/1‑150 (MQ=255)
                                                                                                                           aCCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCa                           >  2:1719076/1‑150 (MQ=255)
                                                                                                                              ccatggaaCCACTGGAGTGGTTGCATGGTGTCCATAATGCACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATgc                        >  2:742138/7‑149 (MQ=255)
                                                                                                                              ccatggaaCCACTGGAGTGGTTGCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGa                        >  1:1808567/7‑150 (MQ=255)
                                                                                                                                atggaaCCACTGGAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGt                      >  1:1220756/5‑150 (MQ=255)
                                                                                                                                 tGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTg                     <  2:1622073/150‑1 (MQ=255)
                                                                                                                                  gCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTgg                    <  2:2125629/150‑1 (MQ=255)
                                                                                                                                   gaaCCACTGGAGTGGTTGCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGt                   >  2:1257356/2‑150 (MQ=255)
                                                                                                                                     aCCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCa                 >  1:707467/1‑150 (MQ=255)
                                                                                                                                      ccACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCAt                >  2:810072/1‑150 (MQ=255)
                                                                                                                                       cACTCCGGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATa               <  1:582752/150‑1 (MQ=255)
                                                                                                                                       cACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATa               <  1:1711981/150‑1 (MQ=255)
                                                                                                                                        aCTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATaa              <  1:85535/150‑1 (MQ=255)
                                                                                                                                         cTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATaaa             >  1:1806021/1‑150 (MQ=255)
                                                                                                                                          tCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAg            >  2:527992/1‑150 (MQ=255)
                                                                                                                                              gTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAAc        <  2:1143454/150‑1 (MQ=255)
                                                                                                                                                ggTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCt      >  1:1624219/1‑150 (MQ=255)
                                                                                                                                                ggTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCt      >  2:1623876/1‑150 (MQ=255)
                                                                                                                                                 gTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTc     <  2:1502324/150‑1 (MQ=255)
                                                                                                                                                 gTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGATTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTc     <  1:2037226/150‑1 (MQ=255)
                                                                                                                                                    ccATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAAt  <  2:1451595/150‑1 (MQ=255)
                                                                                                                                                    ccATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAAt  <  2:1430129/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAAT  >  NC_000913/2923830‑2924127

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