Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,216,481 (G)6→7 coding (11/765 nt) nfeF ← NADPH‑dependent ferric‑chelate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,216,4811.G92.7% 135.2 / 7.3 55coding (11/765 nt)nfeFNADPH‑dependent ferric‑chelate reductase
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (17/34);  total (18/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

ACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGG‑‑G‑T‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAA  >  NC_000913/3216234‑3216722
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
aCGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCggggg‑‑ggtg                                                                                                                                                                                                                                                    >  1:467805/1‑248 (MQ=255)
  gcgACGGGGGGCGTGGTCCTGCCGGCCAGACGATGCCCTCTGCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACTATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGGCAGTTCACGGAAGCGCAGAGCATTTCGAACGCGCTGCGGGTGGCggggg‑‑ggtgtt                                                                                                                                                                                                                                                  <  2:191421/251‑6 (MQ=255)
    gACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTAt                                                                                                                                                                                                                                              <  1:10110/251‑1 (MQ=255)
      gggTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTc                                                                                                                                                                                                                                            <  1:367253/250‑1 (MQ=255)
       ggTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACAGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTc                                                                                                                                                                                                                                            >  1:44287/1‑250 (MQ=255)
        gTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCAt                                                                                                                                                                                                                                          <  2:391943/251‑1 (MQ=255)
          gggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTGCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTg                                                                                                                                                                                                                                        <  2:156450/251‑1 (MQ=255)
            gCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCATTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGtt                                                                                                                                                                                                                                      <  2:341303/251‑1 (MQ=255)
                       ccgcgcAGACGACGCCCTAATCCGTTACGGATGGCGTCACAAAGTTAGCGTAAGGGTGAGGAAAGAAGAGTTTGCTGTGATCGACAGAGCCACGCGATGTAGAGCCGGCCAGCGCCGCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAACGCAGATCATTGCGAACGGGCTGCGTGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTct                                                                                                                                                                                                                            <  2:259221/245‑2 (MQ=255)
                       ccGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAACGTGAGCGTCCGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTGACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGt                                                                                                                                                                                                                           >  2:60343/1‑251 (MQ=255)
                         ggCCAGACGATGCCCTCTTCCGTGACCGTTGGTGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGa                                                                                                                                                                                                                         <  1:13981/251‑1 (MQ=255)
                            cAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑Tgg                                                                                                                                                                                                                     >  1:395979/1‑251 (MQ=255)
                             aGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGt                                                                                                                                                                                                                    >  1:195809/1‑251 (MQ=255)
                               aCGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACCGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGC‑TGGTa                                                                                                                                                                                                                   >  2:220332/1‑250 (MQ=255)
                               aCGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTaa                                                                                                                                                                                                                  <  1:21426/251‑1 (MQ=255)
                                 gattcccTCTGCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGAAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGTCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATc                                                                                                                                                                                                                <  1:80292/247‑1 (MQ=255)
                                             cGTTACCGTTGGCGGCACAAAGGGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCGCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGGGA‑TTGTAATCAGATATATCtaa                                                                                                                                                                                                    <  2:406586/251‑1 (MQ=255)
                                               ttACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGACGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGA‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCtaaat                                                                                                                                                                                                  >  2:213231/1‑251 (MQ=255)
                                                     ttGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGATTGCTGTGATCGTCAAAGCCACGCGATGTAACACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACt                                                                                                                                                                                            >  1:64612/1‑251 (MQ=255)
                                                           gcACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCaaa                                                                                                                                                                                      <  1:413103/251‑1 (MQ=255)
                                                               aaaGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCAtt                                                                                                                                                                                                                                         >  2:366497/1‑197 (MQ=255)
                                                               aaaGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCAtt                                                                                                                                                                                                                                         <  1:366497/197‑1 (MQ=255)
                                                                      gCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCTCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTTA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGAtt                                                                                                                                                                           <  2:411986/251‑1 (MQ=255)
                                                                      gCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGGAACACAGTCAGTTCACGGAGGCGCAGATCATTGCGAACGCGCGGCTGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCGTTCGTGA‑TGGTAATCAGATAGATCTAAATAGAACTCGCAAATGATAATGAtt                                                                                                                                                                           <  2:480551/251‑1 (MQ=255)
                                                                           aGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TTGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTaa                                                                                                                                                                      <  1:35978/251‑1 (MQ=255)
                                                                             gTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATc                                                                                                                                                                    <  1:146489/251‑1 (MQ=255)
                                                                                 aGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATg                                                                                                                                                                 >  2:78878/1‑250 (MQ=255)
                                                                                           agTTTGCTGTGATCGTCAAAGCCCCGCGATGTAACACCGTCCAGCGCCTCGCAGCCGAGGACAATGCGCTGAAAGCCGGCGCTGATGCGCTTAACGCGTAACACAGTCAGTTCACGGAAGCGCCGATCATTGCGAACGCGCTGCGGGGGGCGGGGG‑‑G‑G‑GGTTCTTCATTGTTATCTCCTTCGTGC‑TGGTAATCAGCTATATCTCAATAAAACTCGCAAATGATAATGTTTGTTAATCATGCTAAATGCa                                                                                                                                                        >  2:312937/1‑249 (MQ=255)
                                                                                                gCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGAtt                                                                                                                                                 <  2:361959/251‑1 (MQ=255)
                                                                                                   gtgATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTAGTCATTTTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAAGTATAATGATTGTTAATCATGATAAATGCAAGCGAttt                                                                                                                                                <  2:478953/249‑1 (MQ=255)
                                                                                                                ccACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGt                                                                                                                                 <  1:281521/251‑1 (MQ=255)
                                                                                                                 cACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGCACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTACTCATGATAAATGCAAGCGATTTGTAGAACTGATATGTc                                                                                                                                >  1:469963/1‑251 (MQ=255)
                                                                                                                           aaaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTg                                                                                                                       <  1:175566/250‑1 (MQ=255)
                                                                                                                           aaaaCCGCCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAGCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGGGA‑TGGTAATCAGATATATCTAAAGAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTg                                                                                                                       <  1:312937/250‑1 (MQ=255)
                                                                                                                            aaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGGAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGAGTTGTAGAACTGATATGTCTATAGTCTGAt                                                                                                                     <  2:389975/251‑1 (MQ=255)
                                                                                                                            aaaCCGGCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGAt                                                                                                                     <  2:152093/251‑1 (MQ=255)
                                                                                                                               ccGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGT‑G‑T‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAg                                                                                                                  <  1:202902/251‑1 (MQ=255)
                                                                                                                                   ccAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCGGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTG‑T‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACg                                                                                                               <  1:479434/251‑1 (MQ=255)
                                                                                                                                    cAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGa                                                                                                              >  2:321629/1‑250 (MQ=255)
                                                                                                                                     aGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAAc                                                                                                            <  2:467805/251‑1 (MQ=255)
                                                                                                                                      gcgcCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAATCGATTTGGAGAACTGATATGTCTATAGGCTGATAAGACGAAcc                                                                                                           <  2:64612/251‑1 (MQ=255)
                                                                                                                                       cgcCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCg                                                                                                          <  1:392715/251‑1 (MQ=255)
                                                                                                                                             gccgccGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCtctt                                                                                                    >  2:283929/1‑251 (MQ=255)
                                                                                                                                                    aGGACAATGCGCTGAACACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGc                                                                                             >  1:407034/1‑251 (MQ=255)
                                                                                                                                                                 gAAACCCGGCTCTGATGCGCTAAACGGGTAACACAGTCAGTTCACGGAGGCGCAAATCATTGCGACCGCGCTGGGGGGGGGGGGGG‑‑GGT‑G‑TTATTCATTGTTACCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACGCGAAAATGATAATGTGTGTTAATCAGAATAAAGGCAAGCGATTTGTAGAACTGATATTTTTATGTTCTAATAAGACGAGCCGCCTCTTCTCAGGCATCATTACTCa                                                                                  <  2:125673/249‑1 (MQ=255)
                                                                                                                                                                           gcTGATGCGCTCAACGCGAAACACAGTCAGTTCACGGAAGCGCAGATAATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCACGCGATTTGTAGAACTGAAATGTCCATAGTCTGATAAGACGAACCGCCTCTTATCAGGCATCATTACTCAACCCCggcagcg                                                                      >  2:326595/1‑246 (MQ=255)
                                                                                                                                                                                                   agtcattTCGCGGAAGCGCAGATCATGGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTTTTCATTGTTATCGCTTTC‑TGAATGGTAATCAGAAATATCTAAAAAAACCTCGCAAATGATAATGATTGTAAAGCATGATAAATGCAACCGTTTTGTAGAACTGAAATGTTTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                              <  1:371084/251‑1 (MQ=255)
                                                                                                                                                                                                       ggtTCACGGAAGCGCAGTTTATTGCGATCGCGCTGGGGGGGGGGGGGG‑‑GGT‑G‑TTATTCTTTGTTATCGCCTTCGTTA‑TGGTAATCAGATATATCTAAAAAAAACGCGCAAATGAAAATGATTGTTAATCAGGATAAATGCAAGCGATTGGTAGAACTGATATGTGTATAGTCGGATAAGACGAACCGCCTCTTCTCAGGCATCAGTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtgt                                          <  2:407034/250‑1 (MQ=255)
                                                                                                                                                                                                       agtTCACGGAAGCGCAGATCATGGCGAACGCGCTGCGGGGGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGGGA‑TGGTAATCAGATATATCTAAAAAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGAAAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtgt                                          <  1:213231/251‑1 (MQ=255)
                                                                                                                                                                                                            ccgaaaGCTCAGATCATTGCGACCCCGCTGCGGGTGGCGGGGG‑‑G‑TTG‑TTATTCATGGTTAGCGCCTTCGTAA‑TGGTAATCAAATATATCTAAATAAACCCCGCAAATGAAAATGATTGTTAATCATGATAAATCCAAGCGGTTGGTAGACCTAATATTTCTATAGTCGGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGAt                                     <  1:419865/247‑1 (MQ=255)
                                                                                                                                                                                                            aCGGAAGCTCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGTGAt                                     >  2:380169/1‑251 (MQ=255)
                                                                                                                                                                                                                    gcAGATGATTGCGAACGCGCGGCGGGTGGCGGGGG‑‑GGT‑T‑TTATTCATTGTTATCGCCTTCGTGA‑GGTTAATAAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCACGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCAc                             <  2:86126/251‑1 (MQ=255)
                                                                                                                                                                                                                       gATCATTGCGAACGCGCTGCGGGTGGGGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCAAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAGGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGcatcat                                                                                        <  1:74954/189‑1 (MQ=255)
                                                                                                                                                                                                                       gATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTCAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGcatcat                                                                                        >  2:74954/1‑189 (MQ=255)
                                                                                                                                                                                                                        aTCATTGCGACCGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGttt                         <  1:92089/251‑1 (MQ=255)
                                                                                                                                                                                                                        aTCATTGCGACCGCGCGGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTTTTATCGCCTTCGTGA‑TGGTAACCAGATATATCTAAAAAAAACTCGCAAATGATAATGATTGTTAATCATGAAAAATGCAAGCGATTTGTAAAACTAATTTGTCTATAGTCGGATAAGACAAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGtt                          <  1:254944/250‑1 (MQ=255)
                                                                                                                                                                                                                                      gcTGCGGGTGGCGGGGG‑‑GGT‑G‑TTCTTTATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATATAAATAAAACTCGCAACTGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGCTAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGAGAACGCCTTATCCGGCCTACGGGTGAGGTGAGTCCACGTTTAAAAATCAttttt            >  2:306438/1‑250 (MQ=255)
                                                                                                                                                                                                                                              ggcgggggg‑‑g‑T‑G‑TTATTCATTGTTATCCCCTTCGTGA‑TGTTAATCAGATATATCTAAAAAAACCCCGAAAATGATAATGTTTTTTAATCATTATAAATGCAAGCGATTTGTAGACCTGATATGTATATAGTCTGATAAGACGAGCCGCCTCTCCTCAGGCATCATTACTCAACGCCTGATTCGGCGTGAACGCCTTATCCTGCCGACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCaa  <  1:306438/247‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
ACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGG‑‑G‑T‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAA  >  NC_000913/3216234‑3216722

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