Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,157,644 A→C intergenic (+194/‑6) dkgA → / → yqhG 2,5‑diketo‑D‑gluconate reductase A/DUF3828 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,157,6440AC100.0% 294.1 / NA 76intergenic (+194/‑6)dkgA/yqhG2,5‑diketo‑D‑gluconate reductase A/DUF3828 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (35/41);  total (35/41)

GTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGG  >  NC_000913/3157416‑3157877
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         
gTCTCGGTCCCGATCCTGGCCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGaa                                                                                                                                                                                                                                 <  1:392773/239‑1 (MQ=255)
     ggTCCCGATCCTGACCAGTGCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTCTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAaataa                                                                                                                                                                                                                            <  1:811549/239‑1 (MQ=255)
          cGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACtt                                                                                                                                                                                                                       <  1:1193707/239‑1 (MQ=255)
          cGATACTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACtt                                                                                                                                                                                                                       <  2:250911/239‑1 (MQ=255)
                   ccAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTa                                                                                                                                                                                                              >  1:309447/1‑239 (MQ=255)
                       ttCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACtt                                                                                                                                                                                                                       >  2:1088102/1‑226 (MQ=255)
                       ttCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACtt                                                                                                                                                                                                                       <  1:1088102/226‑1 (MQ=255)
                       ttCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTagca                                                                                                                                                                                                           >  1:995005/1‑238 (MQ=255)
                       ttCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTagca                                                                                                                                                                                                           >  2:678898/1‑238 (MQ=255)
                                       cAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGtt                                                                                                                                                                                              <  2:615953/235‑1 (MQ=255)
                                            tctcCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATAGAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTac                                                                                                                                                                                     >  2:662485/1‑239 (MQ=255)
                                                 cGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGca                                                                                                                                                                                <  1:825868/239‑1 (MQ=255)
                                                 cGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGca                                                                                                                                                                                <  1:34424/239‑1 (MQ=255)
                                                 cGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGca                                                                                                                                                                                <  1:635998/239‑1 (MQ=255)
                                                      gCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTCCACGCACAGc                                                                                                                                                                            >  2:1194051/1‑238 (MQ=255)
                                                        ggTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGccc                                                                                                                                                                          <  2:743316/238‑1 (MQ=255)
                                                                   gCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGAcc                                                                                                                                                              <  2:520718/239‑1 (MQ=255)
                                                                             ttttCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGaacaaa                                                                                                                                                    <  2:32384/239‑1 (MQ=255)
                                                                                     aaTCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTc                                                                                                                                               <  1:1194051/236‑1 (MQ=255)
                                                                                                 cccTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAg                                                                                                                                <  1:678898/239‑1 (MQ=255)
                                                                                                 cccTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAg                                                                                                                                <  1:940308/239‑1 (MQ=255)
                                                                                                 cccTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAg                                                                                                                                <  2:995005/239‑1 (MQ=255)
                                                                                                   cTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGaa                                                                                                                              <  1:723971/239‑1 (MQ=255)
                                                                                                   cTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGa                                                                                                                               >  1:150783/1‑238 (MQ=255)
                                                                                                              cccgTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGaacaaa                                                                                                                                                    >  2:555036/1‑206 (MQ=255)
                                                                                                              cccgTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGaacaaa                                                                                                                                                    <  1:555036/206‑1 (MQ=255)
                                                                                                                    gttACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCa                                                                                                             <  1:1077381/239‑1 (MQ=255)
                                                                                                                    gttACTCTTCCTTGATCAGGAATGCCAAATATAATGACATCATCATGCAGATCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACGTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGAACTCGCACAGCCCCACAAACAAACCGTAAAACAAACAGTCCGGCATATTTATCAGAACAATAAAGCAGATGCCa                                                                                                             <  2:233484/239‑1 (MQ=255)
                                                                                                                      tACACTTCCGTGTTCAGGAATGCCAAATAGAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTGTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACt                                                                                                           <  1:353788/239‑1 (MQ=255)
                                                                                                                             ccTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGACCTATAAATCAGATGca                                                                                                              >  2:813306/1‑228 (MQ=255)
                                                                                                                                            ccAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGca                                                                                                                                                                                <  2:511811/148‑1 (MQ=255)
                                                                                                                                            ccAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGca                                                                                                                                                                                >  1:511811/1‑148 (MQ=255)
                                                                                                                                            ccAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCg                                                                                     >  1:336867/1‑239 (MQ=255)
                                                                                                                                            ccAAATAGAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCg                                                                                     >  2:981023/1‑239 (MQ=255)
                                                                                                                                                atatAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATCCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGgaga                                                                                 <  1:499070/239‑1 (MQ=255)
                                                                                                                                                 tataAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGc                                                                                                               >  2:931514/1‑208 (MQ=255)
                                                                                                                                                 tataAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGc                                                                                                               <  1:931514/208‑1 (MQ=255)
                                                                                                                                                      ggACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGc                                                                           >  2:962320/1‑239 (MQ=255)
                                                                                                                                                           tcatcatGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAg                                                                                                                                                                                                                                      >  2:1078142/1‑79 (MQ=255)
                                                                                                                                                           tcatcatGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAg                                                                                                                                                                                                                                      <  1:1078142/79‑1 (MQ=255)
                                                                                                                                                           tcatcatGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATa                                                                       >  2:989774/1‑238 (MQ=255)
                                                                                                                                                            catcatGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACtat                                                                                                                          >  1:213325/1‑186 (MQ=255)
                                                                                                                                                            catcatGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACtat                                                                                                                          <  2:213325/186‑1 (MQ=255)
                                                                                                                                                                 tGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGAt                                                                                                                                                                                                                                    <  1:175494/75‑1 (MQ=255)
                                                                                                                                                                 tGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGAt                                                                                                                                                                                                                                    >  2:175494/1‑75 (MQ=255)
                                                                                                                                                                 tGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACt                                                                    >  1:80393/1‑235 (MQ=255)
                                                                                                                                                                      agCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTcaca                                                                                                                                                                  <  2:151162/136‑1 (MQ=255)
                                                                                                                                                                      agCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTcaca                                                                                                                                                                  >  1:151162/1‑136 (MQ=255)
                                                                                                                                                                      agCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGt                                                           >  2:1071656/1‑239 (MQ=255)
                                                                                                                                                                        ccGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTAt                                                         >  1:833579/1‑239 (MQ=255)
                                                                                                                                                                        ccGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTAt                                                         >  1:157861/1‑239 (MQ=255)
                                                                                                                                                                                          acaacTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAaataa                                                                                                                                                                                                                            <  2:1141097/58‑1 (MQ=255)
                                                                                                                                                                                          acaacTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAaataa                                                                                                                                                                                                                            >  1:1141097/1‑58 (MQ=255)
                                                                                                                                                                                            aacTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGcc                                                                                                                                                                                                        <  2:1155562/76‑1 (MQ=255)
                                                                                                                                                                                            aacTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGcc                                                                                                                                                                                                        >  1:1155562/1‑76 (MQ=255)
                                                                                                                                                                                              cTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCg                                                                                                                                             <  2:688722/133‑1 (MQ=255)
                                                                                                                                                                                              cTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCg                                                                                                                                             >  1:688722/1‑133 (MQ=255)
                                                                                                                                                                                                   cTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAg                                                                                                                                                                                                <  1:764420/77‑1 (MQ=255)
                                                                                                                                                                                                   cTATATCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAg                                                                                                                                                                                                >  2:764420/1‑77 (MQ=255)
                                                                                                                                                                                                        tCCGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCtt                          >  1:635019/1‑238 (MQ=255)
                                                                                                                                                                                                         ccGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACtt                                                                                                                                                                                                                       >  2:393418/1‑48 (MQ=255)
                                                                                                                                                                                                         ccGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACtt                                                                                                                                                                                                                       <  1:393418/48‑1 (MQ=255)
                                                                                                                                                                                                          cGGGATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACg                       <  2:309447/239‑1 (MQ=255)
                                                                                                                                                                                                            ggATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGc                                                                                                                                                                                                   >  2:512136/1‑65 (MQ=255)
                                                                                                                                                                                                            ggATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGc                                                                                                                                                                                                   <  1:512136/65‑1 (MQ=255)
                                                                                                                                                                                                             gATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCata                                                                                                                                        <  2:848464/123‑1 (MQ=255)
                                                                                                                                                                                                             gATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCata                                                                                                                                        >  1:848464/1‑123 (MQ=255)
                                                                                                                                                                                                             gATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAg                                                                                                                                >  1:106131/1‑131 (MQ=255)
                                                                                                                                                                                                             gATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAg                                                                                                                                <  2:106131/131‑1 (MQ=255)
                                                                                                                                                                                                             gATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTaa                        >  2:1123216/1‑234 (MQ=255)
                                                                                                                                                                                                             gATAAACGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTAc                        >  1:74235/1‑235 (MQ=255)
                                                                                                                                                                                                                  aCGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCata                                                                                                                                        >  2:646659/1‑118 (MQ=255)
                                                                                                                                                                                                                  aCGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCata                                                                                                                                        <  1:646659/118‑1 (MQ=255)
                                                                                                                                                                                                                   cGCAACGGAGAGGTGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGc              <  2:80393/239‑1 (MQ=255)
                                                                                                                                                                                                                                tGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggctgg  <  2:157861/238‑1 (MQ=255)
                                                                                                                                                                                                                                tGAACAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggctgg  <  1:981023/238‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         
GTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGG  >  NC_000913/3157416‑3157877

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