Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,088 T→G 8.4% intergenic (‑347/+311) yjdN ← / ← yjdM metalloprotein superfamily protein/zinc‑ribbon family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,326,0880TG 8.4% 138.3 / 9.1 59intergenic (‑347/+311)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base T (22/30);  new base G (3/3);  total (25/34)
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.76e-01

GCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCT‑ACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTAT  >  NC_000913/4325806‑4326375
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                
gcCTGATAAGACGCGCGAGAGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGACACATAAACGCCTTACCCGGACTACGGAGGGAGCGGGAAGTAGTGGGCCTGATAAGACGCCCAAACGT‑CGCCTCAGGCGGAAGGCACCAATGCCGGATGCGGCGTAATAGCATTATACGACACACGGAGGGTGGGGGAACTTGGAGGCCTGATAGGACTGCCAAACGTCGCCTCCAGCAGGCCGCCGCGTGGGCGCGTGCGGCGTAAACGTCTTATTCGGCACACGGATGGttac                                                                                                                                                                                                                                                                                            >  2:209720/1‑286 (MQ=255)
        agacgggcaggggTCGGATGAGGCAGTGGGAACTGGTGCCGGAGGGGGCGTGAACGCCTTATCCGGCGTAGGGGGGGTGCGGGAAGTGGTAGGCCGGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTg                                                                                                                                                                                                                                                                                    <  2:1264281/277‑1 (MQ=11)
          aCGGGCGAGGGGCGCATCAGGCAGTCGGGACTGTGGCCGGGGGCGGCGTAAAAGCCTTATCCGGCGTACGGGGGGTGCGGGAATGGGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGt                                                                                                                                                                                                                                                                                   <  2:302872/288‑1 (MQ=11)
              gcgtgcgtCGCGTCAGGGAGTCGGTACGGTTGCGGGAGGCGGCGTAAGCGCCGTAGCGTGCCGACGGAGGGTGCGGGGATTTGTAGGCCTGATAAGAGGGGGAAGCGT‑GGCATCAGGCAGTCGGCACCATTGCCGGATGAGGCGTAAACGCCTTATCCGGCCTGCGGAGGGTGCGGGAATTTGGAGGCCGGATAAGGCGCGCGAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGc                                                                                                                                                                                                                                                                               <  2:789850/288‑1 (MQ=11)
                                   cGGCACTGTTGCCGGATGCGGAGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATTAGGCGTAAACGCCTTATCCGGCCTACGGAGGGAGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGAAGTCGGGCACGGTGCCGGATGCGGCGTAAACGCCTTATCCGGACAACGGATGGTGCGGGGATTTGTAGGCATGATAAGACGAGCACGCGTc                                                                                                                                                                                                                                                          >  1:556507/1‑288 (MQ=11)
                                       aacgtttcaggATGCGGCGTAAACGCCTTATCCGGCCTACGGGGGGTGCGGGAATTTGTAGGCCTGATAGGACGCGCAAGCGT‑CGCAGCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTGATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                                                                                                                                                                                                                                     <  2:1382534/280‑1 (MQ=255)
                                         gtttGCCCGATGCGGCGGAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAATGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                                                                                                                                                                                                                   <  1:1410873/287‑1 (MQ=18)
                                          ggtGCCGGTTGCGGCTTAAACGCCTTATCCGGGCTACGGAGGGTGCGGGAATGTGTAGGCCTGATAAGACGGGCAAACGT‑GGCATCAGGCAGTCGGCACCATGGCCGGATGCGGCGTAAACGCCTTATCCGGCCGACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAgg                                                                                                                                                                                                                                                  <  2:744929/287‑1 (MQ=11)
                                                  aTGCGGCGTAAACGCCTTATCCGGCCGACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATGGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTgg                                                                                                                                                                                                                                           <  1:1502856/288‑1 (MQ=37)
                                                   tGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCa                                                                                                                                                                                                                                         <  2:1434919/289‑1 (MQ=37)
                                                   tGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCa                                                                                                                                                                                                                                         <  1:418420/289‑1 (MQ=21)
                                                        cgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg                                                                                                                                                                                                                                                                                    <  2:1425115/241‑1 (MQ=18)
                                                                  ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAt                                                                                                                                                                                                                                                                          >  1:1425115/1‑241 (MQ=18)
                                                                    aTCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCCGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATAAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGCATTTGGAGGCCTGATAAGACGCGCACGCGTCGCATCGTGCAGTGGGCTCTGTTGCCGTACgcggc                                                                                                                                                                                                                        >  2:651641/1‑289 (MQ=11)
                                                                        ggCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTaa                                                                                                                                                                                                                    <  1:18265/289‑1 (MQ=38)
                                                                            gtcGGAGGGTGCGGGAGTGTGTAGACCGGATCAGGCGCGCAA‑CGTGGGCATCAGGCGGTGGGAACCATTGCCGGGTGCGGGGGAAAGGCCCTATCCGGGTTAGGGAGGGTTCGGGAATTTGTCGGCCTGATAAGACGCGCAAGCGTCGAATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATTGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGc                                                                                                                                                                                                                <  2:1410812/287‑1 (MQ=11)
                                                                                  gtgggCGGGAATGTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                                                                           <  1:469356/284‑1 (MQ=38)
                                                                                  gggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCTGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                                                                           <  1:651641/288‑1 (MQ=38)
                                                                                   ggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGGAAACGCCTTATcc                                                                                                                                                                                                         >  1:909614/1‑289 (MQ=38)
                                                                                    gTGCGGGAGTTTGTAGTCCTGGTCAGACGCGAAAGCGT‑CCCATCAGGCATTCGGCGCCATTGCCGGATGCGGCTTAAACGCTTTAGCCGGCCTAGGGAGGGTGCGGGAATGTGTAGGCCTGATAAGAGGCGCAAGCGTCGCATCAGGAAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                                                                                         <  2:1179367/288‑1 (MQ=11)
                                                                                             tttGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGa                                                                                                                                                                                                                                                                                         >  1:338085/1‑199 (MQ=37)
                                                                                                    gCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                                                                                          <  1:726296/271‑1 (MQ=38)
                                                                                                    gCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                                                                                          >  2:726296/1‑271 (MQ=38)
                                                                                                       tGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATGTGTAGGCCTGATAAGACGCGCTAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGTGGGTGCGGGa                                                                                                                                                                                     <  1:209720/289‑1 (MQ=255)
                                                                                                       ggATAAGACGCGAAAGCGT‑CGCGTCAGGCAGTCGGCACCATTGCGGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                                                                                                                                                                     <  2:1189879/288‑1 (MQ=255)
                                                                                                                 cgcAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                                                                                                                                                                                                                   <  2:1304509/217‑1 (MQ=35)
                                                                                                                 cgcAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                                                                                                                                                                                                                   >  1:1304509/1‑217 (MQ=35)
                                                                                                                   cAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGAGGGGGCGGGAATTTGTTgg                                                                                                                                                                            >  1:82353/1‑286 (MQ=255)
                                                                                                                             cATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                                                                                                                                                                                                                                                                         >  2:11500/1‑168 (MQ=18)
                                                                                                                                    cAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCa                                                                                                                                                            <  1:977525/286‑1 (MQ=38)
                                                                                                                                                ttGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGc                                                                                                                                             <  1:991542/289‑1 (MQ=35)
                                                                                                                                                     ggAGGCGGGGTAGACGCCTGATCCGGCCTGAGGGGGTTGCGGGGAGGTGTCGGCCGGAGAAGACGGGCAAGCGTCGCATCAGGCAGTCGGCACCGTGGCCGGATGGGGCGTAAACGCCTTATCCGGCCTACGGATGGTACGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTc                                                                                                                                         <  2:574598/288‑1 (MQ=11)
                                                                                                                                                        tGCGGCGCAAACGCCTTAGCCGGCGTACGGAGGCGACGAGGATGTGGAGGCCGGAGAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGGTGCCGGAGGCGGCGTAAACGCCTTATCCGGCCTGCGGATGGTGCGGGCAGTTGTAGGCCTGATAAGACGCGGAAGTGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCGACGGAGGGTGCGGGAATTGGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCa                                                                                                                                     <  2:319889/289‑1 (MQ=11)
                                                                                                                                                                     gtttggcgggCTTACGGAGGGAGCGGGAATGGGTAGGCCTGATAAGACGGGCAAGCGTCGGATCAGGCAGTCGGCACCGTTGCCGGGTGCGGCGTAACCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGAt                                                                                                                         <  2:1188273/280‑1 (MQ=11)
                                                                                                                                                                       ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTCGGCCTGATAAGACGCGCAAGTGTCGCACCCGGCCTCTGGCACCGTTTCCGGATgcg                                                                                                                      >  2:1131423/1‑289 (MQ=11)
                                                                                                                                                                                           gCGGGGAGTTGTGGGCCTGCTAAGACGCGCAAGGGTCGCATCAGGCAGTGGGCACCGTTGCCGGATGCGGCGGAAACGCCGTGTCCGGGCTACGGATGGTGCGGGAATTTGGAGGCCGGATAAGACGCGCAAGCGTCGCATCAGGGAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGGGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                   <  2:95772/288‑1 (MQ=11)
                                                                                                                                                                                            cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATg                                                                                                                                                                                                                            <  1:1096774/166‑1 (MQ=38)
                                                                                                                                                                                             gggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCCGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGAGGGAATTTGTGGGCCTGATAAGACGCGCAAGCGGAGCATCCGGCGTTTGGCCCCGTTTCCGGATGCCGCGCCAGCGCCTTATCCGGct                                                                                                 >  2:408873/1‑287 (MQ=11)
                                                                                                                                                                                                    tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAg                                                                                                                                                                                                                                                                                 <  1:886033/105‑1 (MQ=32)
                                                                                                                                                                                                    tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCa                                                                                                                                                                                                                                         <  1:875079/145‑1 (MQ=37)
                                                                                                                                                                                                      tAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAAc                                                                                                                                                                                                                  >  2:1096774/1‑166 (MQ=38)
                                                                                                                                                                                                         gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGa                                                                                                                                                                                                                                                                                         <  2:322146/92‑1 (MQ=33)
                                                                                                                                                                                                           cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCT‑ACGGAGGGTGCggg                                                                                 >  1:1425517/1‑289 (MQ=38)
                                                                                                                                                                                                              attaGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTAACCGTCCTACGGATGGTGCGGGAAGTTGTAGGACTGATagt                                                                                                                                                                                                                                                                       >  2:886033/4‑103 (MQ=11)
                                                                                                                                                                                                              aTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCgg                                                                                                                                                                                                                               >  2:875079/1‑145 (MQ=39)
                                                                                                                                                                                                                    cgcgCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAttt                                                                                                                                                                                                                                                                                     <  2:974859/85‑1 (MQ=33)
                                                                                                                                                                                                                    cgcgCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAttt                                                                                                                                                                                                                                                                                     >  1:974859/1‑85 (MQ=33)
                                                                                                                                                                                                                        cAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCTTGATAAGACTGGCAAGCCTCGCATCAGGCAGTCGGCAACGTTGCCGGATGCGGCGTACACGCACTACCCGGCCTAACGGTGGGTCCggacattgggacgtg                                                                     >  2:30073/1‑276 (MQ=11)
                                                                                                                                                                                                                                                    ttGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGACCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCT‑ACGGATGGTGCGGGAATTTGTAGGCATGACAAGACGCGCACGCATCGCACCAAGc                                        >  1:777848/1‑289 (MQ=255)
                                                                                                                                                                                                                                                         ggATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGGAAACGCCTTATCCGGCCT‑ACGGAGGGTGCGGGAATTTGTAGGGCTGATACGACCCGCACGCGTCGCACCAGGCATCTg                                   >  2:1580758/1‑289 (MQ=255)
                                                                                                                                                                                                                                                            tGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCa                                                                                                                                                                                                                                                                <  1:1390200/66‑1 (MQ=35)
                                                                                                                                                                                                                                                                        gCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGCTGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAAACGCACAAGCGTCCAATCAGGCAGTCACCACCGTTTCCGGATGCGGCGTACACGCCTTATCCGGCCT‑ACGGAGGGTGCGGCACTTTTTTGTCCTGATTACACGCGCAATCGCCGCACCAGGCATCAGCAACCATTGcccgca                    >  2:988596/1‑285 (MQ=11)
                                                                                                                                                                                                                                                                                gggCGTAGGGAGGGTGCGGGTATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGTACCGTTGCCGGATGCGGCGTAAACGCCTTAGCCGGCCT‑ACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGcgcg            <  1:599905/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  gCCTAGGGAAGGGGCGGGAATTTGGAGGCCTGAGAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTAGCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCT‑ACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGa           <  2:1010971/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  gCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCa                                                                                                                                                                                                                                         <  1:1451461/67‑1 (MQ=40)
                                                                                                                                                                                                                                                                                  gCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCa                                                                                                                                                                                                                                         >  2:1451461/1‑67 (MQ=40)
                                                                                                                                                                                                                                                                                  gCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCT‑ACGGAGGGTGCGGGAATTTGTAGGCCTGAACAGGCGCGCAAGTGGCGCATCCGGCATCCGGCCCCCTTGCCGGCTGCGGCGCGa           >  1:235635/1‑288 (MQ=14)
                                                                                                                                                                                                                                                                                    cTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCT‑ACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTTGCACCGTTGACTGATGCGGCGCGAgc         >  1:450995/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                       cGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCT‑ACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTTCCGGATGCGGCGCGAGCGCCt     >  1:609570/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                           tGGTGCGGGAATTGCGAGGCTTGCAAAGACGCGCAAGGGTGGCATCAGGCAGTCGGAACTGTTGCCGAATGCGGCGTAAACGCCTTATCGGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGTTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCT‑ACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTAt  <  2:1518088/288‑1 (MQ=14)
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                
GCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGT‑CGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCT‑ACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTAT  >  NC_000913/4325806‑4326375

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