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‑‑GTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAA  >  NC_000913/3216234‑3216722
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
aCGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCggggg‑‑ggtg                                                                                                                                                                                                                                                    >  1:467805/1‑248 (MQ=255)
  gcgACGGGGGGCGTGGTCCTGCCGGCCAGACGATGCCCTCTGCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACTATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGGCAGTTCACGGAAGCGCAGAGCATTTCGAACGCGCTGCGGGTGGCggggg‑‑ggtgtt                                                                                                                                                                                                                                                  <  2:191421/251‑6 (MQ=255)
    gACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTAt                                                                                                                                                                                                                                              <  1:10110/251‑1 (MQ=255)
      gggTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTc                                                                                                                                                                                                                                            <  1:367253/250‑1 (MQ=255)
       ggTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACAGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTc                                                                                                                                                                                                                                            >  1:44287/1‑250 (MQ=255)
        gTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCAt                                                                                                                                                                                                                                          <  2:391943/251‑1 (MQ=255)
          gggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTGCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTg                                                                                                                                                                                                                                        <  2:156450/251‑1 (MQ=255)
            gCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCATTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGtt                                                                                                                                                                                                                                      <  2:341303/251‑1 (MQ=255)
                       ccgcgcAGACGACGCCCTAATCCGTTACGGATGGCGTCACAAAGTTAGCGTAAGGGTGAGGAAAGAAGAGTTTGCTGTGATCGACAGAGCCACGCGATGTAGAGCCGGCCAGCGCCGCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAACGCAGATCATTGCGAACGGGCTGCGTGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTct                                                                                                                                                                                                                            <  2:259221/245‑2 (MQ=255)
                       ccGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAACGTGAGCGTCCGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTGACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGt                                                                                                                                                                                                                           >  2:60343/1‑251 (MQ=255)
                         ggCCAGACGATGCCCTCTTCCGTGACCGTTGGTGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGa                                                                                                                                                                                                                         <  1:13981/251‑1 (MQ=255)
                            cAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATgg                                                                                                                                                                                                                     >  1:395979/1‑251 (MQ=255)
                             aGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGt                                                                                                                                                                                                                    >  1:195809/1‑251 (MQ=255)
                               aCGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACCGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGCTGGTa                                                                                                                                                                                                                   >  2:220332/1‑250 (MQ=255)
                               aCGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTaa                                                                                                                                                                                                                  <  1:21426/251‑1 (MQ=255)
                                 gattcccTCTGCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGAAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGTCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATc                                                                                                                                                                                                                <  1:80292/247‑1 (MQ=255)
                                             cGTTACCGTTGGCGGCACAAAGGGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCGCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGGGATTGTAATCAGATATATCtaa                                                                                                                                                                                                    <  2:406586/251‑1 (MQ=255)
                                               ttACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGACGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGAGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCtaaat                                                                                                                                                                                                  >  2:213231/1‑251 (MQ=255)
                                                     ttGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGATTGCTGTGATCGTCAAAGCCACGCGATGTAACACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACt                                                                                                                                                                                            >  1:64612/1‑251 (MQ=255)
                                                           gcACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCaaa                                                                                                                                                                                      <  1:413103/251‑1 (MQ=255)
                                                               aaaGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCAtt                                                                                                                                                                                                                                         >  2:366497/1‑197 (MQ=255)
                                                               aaaGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCAtt                                                                                                                                                                                                                                         <  1:366497/197‑1 (MQ=255)
                                                                      gCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCTCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTTATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGAtt                                                                                                                                                                           <  2:411986/251‑1 (MQ=255)
                                                                      gCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGGAACACAGTCAGTTCACGGAGGCGCAGATCATTGCGAACGCGCGGCTGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCGTTCGTGATGGTAATCAGATAGATCTAAATAGAACTCGCAAATGATAATGAtt                                                                                                                                                                           <  2:480551/251‑1 (MQ=255)
                                                                           aGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATTGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTaa                                                                                                                                                                      <  1:35978/251‑1 (MQ=255)
                                                                             gTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATc                                                                                                                                                                    <  1:146489/251‑1 (MQ=255)
                                                                                 aGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATg                                                                                                                                                                 >  2:78878/1‑250 (MQ=255)
                                                                                           agTTTGCTGTGATCGTCAAAGCCCCGCGATGTAACACCGTCCAGCGCCTCGCAGCCGAGGACAATGCGCTGAAAGCCGGCGCTGATGCGCTTAACGCGTAACACAGTCAGTTCACGGAAGCGCCGATCATTGCGAACGCGCTGCGGGGGGCGGGGG‑‑GGGGTTCTTCATTGTTATCTCCTTCGTGCTGGTAATCAGCTATATCTCAATAAAACTCGCAAATGATAATGTTTGTTAATCATGCTAAATGCa                                                                                                                                                        >  2:312937/1‑249 (MQ=255)
                                                                                                gCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGAtt                                                                                                                                                 <  2:361959/251‑1 (MQ=255)
                                                                                                   gtgATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTAGTCATTTTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAAGTATAATGATTGTTAATCATGATAAATGCAAGCGAttt                                                                                                                                                <  2:478953/249‑1 (MQ=255)
                                                                                                                ccACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGt                                                                                                                                 <  1:281521/251‑1 (MQ=255)
                                                                                                                 cACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGCACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTACTCATGATAAATGCAAGCGATTTGTAGAACTGATATGTc                                                                                                                                >  1:469963/1‑251 (MQ=255)
                                                                                                                           aaaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTg                                                                                                                       <  1:175566/250‑1 (MQ=255)
                                                                                                                           aaaaCCGCCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAGCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGGGATGGTAATCAGATATATCTAAAGAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTg                                                                                                                       <  1:312937/250‑1 (MQ=255)
                                                                                                                            aaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGGAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGAGTTGTAGAACTGATATGTCTATAGTCTGAt                                                                                                                     <  2:389975/251‑1 (MQ=255)
                                                                                                                            aaaCCGGCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGAt                                                                                                                     <  2:152093/251‑1 (MQ=255)
                                                                                                                               ccGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGTGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAg                                                                                                                  <  1:202902/251‑1 (MQ=255)
                                                                                                                                   ccAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCGGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACg                                                                                                               <  1:479434/251‑1 (MQ=255)
                                                                                                                                    cAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGa                                                                                                              >  2:321629/1‑250 (MQ=255)
                                                                                                                                     aGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAAc                                                                                                            <  2:467805/251‑1 (MQ=255)
                                                                                                                                      gcgcCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAATCGATTTGGAGAACTGATATGTCTATAGGCTGATAAGACGAAcc                                                                                                           <  2:64612/251‑1 (MQ=255)
                                                                                                                                       cgcCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCg                                                                                                          <  1:392715/251‑1 (MQ=255)
                                                                                                                                             gccgccGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCtctt                                                                                                    >  2:283929/1‑251 (MQ=255)
                                                                                                                                                    aGGACAATGCGCTGAACACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGc                                                                                             >  1:407034/1‑251 (MQ=255)
                                                                                                                                                                 gAAACCCGGCTCTGATGCGCTAAACGGGTAACACAGTCAGTTCACGGAGGCGCAAATCATTGCGACCGCGCTGGGGGGGGGGGGGG‑‑GGTGTTATTCATTGTTACCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACGCGAAAATGATAATGTGTGTTAATCAGAATAAAGGCAAGCGATTTGTAGAACTGATATTTTTATGTTCTAATAAGACGAGCCGCCTCTTCTCAGGCATCATTACTCa                                                                                  <  2:125673/249‑1 (MQ=255)
                                                                                                                                                                           gcTGATGCGCTCAACGCGAAACACAGTCAGTTCACGGAAGCGCAGATAATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCACGCGATTTGTAGAACTGAAATGTCCATAGTCTGATAAGACGAACCGCCTCTTATCAGGCATCATTACTCAACCCCggcagcg                                                                      >  2:326595/1‑246 (MQ=255)
                                                                                                                                                                                                   agtcattTCGCGGAAGCGCAGATCATGGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTTTTCATTGTTATCGCTTTCTGAATGGTAATCAGAAATATCTAAAAAAACCTCGCAAATGATAATGATTGTAAAGCATGATAAATGCAACCGTTTTGTAGAACTGAAATGTTTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                              <  1:371084/251‑1 (MQ=255)
                                                                                                                                                                                                       ggtTCACGGAAGCGCAGTTTATTGCGATCGCGCTGGGGGGGGGGGGGG‑‑GGTGTTATTCTTTGTTATCGCCTTCGTTATGGTAATCAGATATATCTAAAAAAAACGCGCAAATGAAAATGATTGTTAATCAGGATAAATGCAAGCGATTGGTAGAACTGATATGTGTATAGTCGGATAAGACGAACCGCCTCTTCTCAGGCATCAGTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtgt                                          <  2:407034/250‑1 (MQ=255)
                                                                                                                                                                                                       agtTCACGGAAGCGCAGATCATGGCGAACGCGCTGCGGGGGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGGGATGGTAATCAGATATATCTAAAAAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGAAAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtgt                                          <  1:213231/251‑1 (MQ=255)
                                                                                                                                                                                                            ccgaaaGCTCAGATCATTGCGACCCCGCTGCGGGTGGCGGGGG‑‑GTTGTTATTCATGGTTAGCGCCTTCGTAATGGTAATCAAATATATCTAAATAAACCCCGCAAATGAAAATGATTGTTAATCATGATAAATCCAAGCGGTTGGTAGACCTAATATTTCTATAGTCGGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGAt                                     <  1:419865/247‑1 (MQ=255)
                                                                                                                                                                                                            aCGGAAGCTCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGTGAt                                     >  2:380169/1‑251 (MQ=255)
                                                                                                                                                                                                                    gcAGATGATTGCGAACGCGCGGCGGGTGGCGGGGG‑‑GGTTTTATTCATTGTTATCGCCTTCGTGAGGTTAATAAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCACGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCAc                             <  2:86126/251‑1 (MQ=255)
                                                                                                                                                                                                                       gATCATTGCGAACGCGCTGCGGGTGGGGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCAAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAGGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGcatcat                                                                                        <  1:74954/189‑1 (MQ=255)
                                                                                                                                                                                                                       gATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTCAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGcatcat                                                                                        >  2:74954/1‑189 (MQ=255)
                                                                                                                                                                                                                        aTCATTGCGACCGCGCTGCGGGTGGCGGGGG‑‑GGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGttt                         <  1:92089/251‑1 (MQ=255)
                                                                                                                                                                                                                        aTCATTGCGACCGCGCGGCGGGTGGCGGGGG‑‑GGTGTTATTCATTTTTATCGCCTTCGTGATGGTAACCAGATATATCTAAAAAAAACTCGCAAATGATAATGATTGTTAATCATGAAAAATGCAAGCGATTTGTAAAACTAATTTGTCTATAGTCGGATAAGACAAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGtt                          <  1:254944/250‑1 (MQ=255)
                                                                                                                                                                                                                                      gcTGCGGGTGGCGGGGG‑‑GGTGTTCTTTATTGTTATCGCCTTCGTGATGGTAATCAGATATATATAAATAAAACTCGCAACTGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGCTAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGAGAACGCCTTATCCGGCCTACGGGTGAGGTGAGTCCACGTTTAAAAATCAttttt            >  2:306438/1‑250 (MQ=255)
                                                                                                                                                                                                                                              ggcgggggg‑‑gTGTTATTCATTGTTATCCCCTTCGTGATGTTAATCAGATATATCTAAAAAAACCCCGAAAATGATAATGTTTTTTAATCATTATAAATGCAAGCGATTTGTAGACCTGATATGTATATAGTCTGATAAGACGAGCCGCCTCTCCTCAGGCATCATTACTCAACGCCTGATTCGGCGTGAACGCCTTATCCTGCCGACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCaa  <  1:306438/247‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
ACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGG‑‑GTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAA  >  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: