Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,817,197 G→A R103C (CGC→TGC)  yqaB ← fructose‑1‑phosphate phosphatase YqaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,817,1970GA98.5% 215.2 / NA 68R103C (CGC→TGC) yqaBfructose‑1‑phosphate phosphatase YqaB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/34);  minor base C (0/1);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGTAGG  >  NC_000913/2816952‑2817443
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
aTCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCaacgac                                                                                                                                                                                                                                                   >  1:264234/1‑251 (MQ=255)
    aCGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGATTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAAACCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTATCACTCCCCGTTCCTACCGCCATGGGGCAACTACCATg                                                                                                                                                                                                                                               >  2:434410/1‑251 (MQ=255)
         gTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAAc                                                                                                                                                                                                                                          <  1:168313/251‑1 (MQ=255)
           ccATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCCTCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCTTACAGCCATTGGGCAACGACCACGCCAACtt                                                                                                                                                                                                                                        >  2:452666/1‑251 (MQ=255)
            cATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACCGCCATTGGGCAACGACCATGCCAACttt                                                                                                                                                                                                                                       >  2:584027/1‑251 (MQ=255)
            cATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGAACATGCCAACttt                                                                                                                                                                                                                                       >  2:435947/1‑251 (MQ=255)
               gCCTGCTGCACGGGCCGCCTGAATACCGAAATCGTCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCAGACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTca                                                                                                                                                                                                                                    <  2:442290/251‑1 (MQ=255)
                         cGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCCCAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGCTCGGCAGCGACGCCGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCCCATcaac                                                                                                                                                                                                                          >  2:414887/1‑251 (MQ=255)
                                     aTACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGTCGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCgg                                                                                                                                                                                                              >  2:319530/1‑251 (MQ=255)
                                       aCCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGtt                                                                                                                                                                                                            <  1:487036/251‑1 (MQ=255)
                                         cGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCg                                                                                                                                                                                                          <  1:474994/251‑1 (MQ=255)
                                                ggCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTAt                                                                                                                                                                                                   >  1:280747/1‑251 (MQ=255)
                                                    tCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGAACGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGAGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTACACCACATCAACAAGAGGAAGCGGTTCGACGCTATcc                                                                                                                                                                                                 >  1:542445/1‑249 (MQ=255)
                                                      ttCAAAGGCCACACACTGCATCGGTTGCACGCCCATACGCTGCGCGCACCAAAAAAATGTGTCTGGCGCGGGTTTATGGCGTTTGACGAGATCGGCAGCGACGACGGCGTCAAAATAATGGAGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACCTCACCAAGAGGAAGCGGTTAGACGCTATCcatc                                                                                                                                                                                              >  2:32945/1‑248 (MQ=255)
                                                       tCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCcagc                                                                                                                                                                                              <  2:513936/249‑1 (MQ=255)
                                                       tCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCca                                                                                                                                                                                                <  1:231390/247‑1 (MQ=255)
                                                       tCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGGGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCCCCACATCAACAAGAGGAAGCGGTTCGACGCTATCcagc                                                                                                                                                                                              >  2:46604/1‑249 (MQ=255)
                                                                      tGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGc                                                                                                                                                                             >  1:572124/1‑251 (MQ=255)
                                                                            ggTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGt                                                                                                                                                                       <  2:472201/251‑1 (MQ=255)
                                                                             gTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGtt                                                                                                                                                                      >  1:78560/1‑251 (MQ=255)
                                                                                      ccATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCAcgc                                                                                                                                                             <  2:148080/251‑1 (MQ=255)
                                                                                        aTACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACgcgc                                                                                                                                                           >  2:22736/1‑251 (MQ=255)
                                                                                          aCGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCt                                                                                                                                                          <  2:275277/250‑1 (MQ=255)
                                                                                          aCGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACGCCCCGTTCCTACAGCCAGTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTGTTTCACGCGCTa                                                                                                                                                         <  2:47602/251‑1 (MQ=255)
                                                                                            gCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAAc                                                                                                                                                       <  1:434624/251‑1 (MQ=255)
                                                                                                   gcACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACGGCTTCTGTTTTTTCACGCGCGAACGCATGCg                                                                                                                                                <  1:445981/251‑1 (MQ=255)
                                                                                                       aacatcaaTGTGTCTGGCGTGGGTTTATGGTGTTTGACGTGAGCGGCAGCGACGACGGCGTCAAAATAATGGCGTACCCCCAGGTGCGCCAGCAATGCCGCAGCGATGGCTCTTTCACTCCCCGTTCCTACAGCCATTGGGCCACGACCATGCCACCTTGTCGCCTCTTCAACAAGAGTAAGCGGTTCGCCGCTATCCAGCAGCATGCTTCTTACTGCATCTGTTTTTTCACGCGCTAACGCAGGCGGGTc                                                                                                                                            <  2:58349/245‑1 (MQ=255)
                                                                                                       aacaaAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTc                                                                                                                                            <  2:160325/251‑1 (MQ=255)
                                                                                                                    cTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTg                                                                                                                                >  1:324519/1‑250 (MQ=255)
                                                                                                                    cTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTg                                                                                                                                >  2:169318/1‑250 (MQ=255)
                                                                                                                    cTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTg                                                                                                                                >  1:104880/1‑250 (MQ=255)
                                                                                                                     tGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTGCTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGAt                                                                                                                              >  1:584857/1‑251 (MQ=255)
                                                                                                                          cgGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGc                                                                                                                         <  2:99491/251‑1 (MQ=255)
                                                                                                                          cgGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGCGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGCTCGGACTGATTCAGc                                                                                                                         >  1:26092/1‑251 (MQ=255)
                                                                                                                             gTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTc                                                                                                                       >  1:330361/1‑250 (MQ=255)
                                                                                                                                tATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCaat                                                                                                                    >  1:147315/1‑250 (MQ=255)
                                                                                                                                 atgttgtTTGACGTGATCGGCAGCGACGACTGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATTGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACGTTTCACCACATCAACAATAGTAAGCGGTTCGTCGCTATCCAGCTGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCaata                                                                                                                   <  1:281835/246‑1 (MQ=255)
                                                                                                                                 aTGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCaata                                                                                                                   >  2:162831/1‑250 (MQ=255)
                                                                                                                                        ttGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCt                                                                                                            <  1:421827/250‑1 (MQ=255)
                                                                                                                                        ttGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCt                                                                                                            <  1:394624/250‑1 (MQ=255)
                                                                                                                                         tGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTg                                                                                                           <  1:215225/250‑1 (MQ=255)
                                                                                                                                             gTGATCGGCAGCGACGAGGGAGACAAAATATTGGCGTAATCCCAGAAGCGCCAGCAATGCCTAAGCGATTGAGCTATCACTCCACGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAAAAAGAGGAAGCGGTTCGACGCTATCAAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTACCGCATGCGGGTCGAGATCGGCCTGCTTCAGCTCAATCATTGCCTGAGCa                                                                                                       >  2:523759/1‑250 (MQ=255)
                                                                                                                                                aTCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATAc                                                                                                   <  2:507132/251‑1 (MQ=255)
                                                                                                                                                     cAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATCCGcc                                                                                                >  2:516650/1‑249 (MQ=255)
                                                                                                                                                     cAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGAGCCAGCAATGCCGCAGCGATGGCGCTTTCACTCCTCGTTCCCACAGCCATTGGGCAACGACCATGCCAACGTTTCACCACATCAACAAGAGGAAGCTGGTCGACGCTATCCAGCAGCATACTTCTTGCTGCTTCTGTTTTTGCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAg                                                                                              <  1:44293/251‑1 (MQ=255)
                                                                                                                                                         gacgacGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTgg                                                                                          <  1:414887/251‑1 (MQ=255)
                                                                                                                                                          acgactGCGTCAAAACAATGGCGTAAACCCAGGTGCGCCATCAATGCCGCAGACATGGCGCTGCCACACCACGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTGCTTACTGCTTCTGTTTGTTCACGCGCTAACGCATGCTGGTCGAGATCGTCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTggg                                                                                         <  2:307505/251‑1 (MQ=255)
                                                                                                                                                           cgacgGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAAGAATTGCCTGAGCAATACGCCAGGTGGGc                                                                                        <  2:511096/251‑1 (MQ=255)
                                                                                                                                                                gCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATcc                                                                                   <  1:135207/251‑1 (MQ=255)
                                                                                                                                                                   tCAAAATAATGGCGTCATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGAGCTTTCACTCCCCGTTCCTACAGCCGTTGGGCACCGACCATGCCACCTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGATTATGTTTTTTCCCGCGCTAACGCATGCGGGTCGAGATCGGCCTGCTTCAGCTCAATATTTGCCTGAGCAATACGCCAGGTGGGCGATCCAtt                                                                                >  2:579258/1‑251 (MQ=255)
                                                                                                                                                                       aataatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAgc                                                                            >  2:200688/1‑251 (MQ=255)
                                                                                                                                                                            gggCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAat                                                                       <  2:104880/250‑1 (MQ=255)
                                                                                                                                                                                  aaTCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGc                                                                 >  2:435440/1‑251 (MQ=255)
                                                                                                                                                                                       cAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGa                                                             >  2:520008/1‑250 (MQ=255)
                                                                                                                                                                                                            aGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATACAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGt                                       >  1:63670/1‑251 (MQ=255)
                                                                                                                                                                                                             gCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTg                                      <  1:304249/251‑1 (MQ=255)
                                                                                                                                                                                                             gCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTg                                      <  1:584027/251‑1 (MQ=255)
                                                                                                                                                                                                               gATGGCGCTTTCACTACCCGTTCCTCCAGCCATTGGGCAACGACCCTGCCAACTTTTCACCAAATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCCTACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGCTTCAGCTCAATAATTGCCTGAGCAATACGCAAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGcc                                    >  2:385016/1‑251 (MQ=255)
                                                                                                                                                                                                                    cgcTTTCACCCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAAt                               <  2:401932/251‑1 (MQ=255)
                                                                                                                                                                                                                     gcTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATa                              <  1:319530/251‑1 (MQ=255)
                                                                                                                                                                                                                       tttCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACg                                                                                                                                                                                                       <  2:475784/80‑1 (MQ=255)
                                                                                                                                                                                                                       tttCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACg                                                                                                                                                                                                       >  1:475784/1‑80 (MQ=255)
                                                                                                                                                                                                                        ttCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATAc                             <  1:210194/249‑1 (MQ=255)
                                                                                                                                                                                                                          cACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACTCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGTCCGTAGTGCCCTAATACTTc                          <  2:134042/250‑1 (MQ=255)
                                                                                                                                                                                                                                        acatccATTGGGCACCGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGAACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACg           >  2:84395/5‑251 (MQ=255)
                                                                                                                                                                                                                                          aGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGCCGCTATCCCGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCATCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCCTCGCCTGAATATCGTACTGAGGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTaaggg         >  2:424864/1‑247 (MQ=255)
                                                                                                                                                                                                                                          aGCCATTGGGCAAAGACCATGCCAACTTGTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGt         <  2:324519/251‑1 (MQ=255)
                                                                                                                                                                                                                                                tGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGCGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTCGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGt     >  2:48792/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 gggCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGTAgg  <  1:46604/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
ATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGTAGG  >  NC_000913/2816952‑2817443

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