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.0 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

GCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTAT  >  NC_000913/4325806‑4326375
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                
gcCTGATAAGACGCGCGAGAGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGACACATAAACGCCTTACCCGGACTACGGAGGGAGCGGGAAGTAGTGGGCCTGATAAGACGCCCAAACGTCGCCTCAGGCGGAAGGCACCAATGCCGGATGCGGCGTAATAGCATTATACGACACACGGAGGGTGGGGGAACTTGGAGGCCTGATAGGACTGCCAAACGTCGCCTCCAGCAGGCCGCCGCGTGGGCGCGTGCGGCGTAAACGTCTTATTCGGCACACGGATGGttac                                                                                                                                                                                                                                                                                            >  2:209720/1‑286 (MQ=255)
        agacgggcaggggTCGGATGAGGCAGTGGGAACTGGTGCCGGAGGGGGCGTGAACGCCTTATCCGGCGTAGGGGGGGTGCGGGAAGTGGTAGGCCGGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTg                                                                                                                                                                                                                                                                                    <  2:1264281/277‑1 (MQ=11)
          aCGGGCGAGGGGCGCATCAGGCAGTCGGGACTGTGGCCGGGGGCGGCGTAAAAGCCTTATCCGGCGTACGGGGGGTGCGGGAATGGGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGt                                                                                                                                                                                                                                                                                   <  2:302872/288‑1 (MQ=11)
              gcgtgcgtCGCGTCAGGGAGTCGGTACGGTTGCGGGAGGCGGCGTAAGCGCCGTAGCGTGCCGACGGAGGGTGCGGGGATTTGTAGGCCTGATAAGAGGGGGAAGCGTGGCATCAGGCAGTCGGCACCATTGCCGGATGAGGCGTAAACGCCTTATCCGGCCTGCGGAGGGTGCGGGAATTTGGAGGCCGGATAAGGCGCGCGAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGc                                                                                                                                                                                                                                                                               <  2:789850/288‑1 (MQ=11)
                                   cGGCACTGTTGCCGGATGCGGAGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATTAGGCGTAAACGCCTTATCCGGCCTACGGAGGGAGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGAAGTCGGGCACGGTGCCGGATGCGGCGTAAACGCCTTATCCGGACAACGGATGGTGCGGGGATTTGTAGGCATGATAAGACGAGCACGCGTc                                                                                                                                                                                                                                                          >  1:556507/1‑288 (MQ=11)
                                       aacgtttcaggATGCGGCGTAAACGCCTTATCCGGCCTACGGGGGGTGCGGGAATTTGTAGGCCTGATAGGACGCGCAAGCGTCGCAGCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTGATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                                                                                                                                                                                                                                     <  2:1382534/280‑1 (MQ=255)
                                         gtttGCCCGATGCGGCGGAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAATGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                                                                                                                                                                                                                   <  1:1410873/287‑1 (MQ=18)
                                          ggtGCCGGTTGCGGCTTAAACGCCTTATCCGGGCTACGGAGGGTGCGGGAATGTGTAGGCCTGATAAGACGGGCAAACGTGGCATCAGGCAGTCGGCACCATGGCCGGATGCGGCGTAAACGCCTTATCCGGCCGACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAgg                                                                                                                                                                                                                                                  <  2:744929/287‑1 (MQ=11)
                                                  aTGCGGCGTAAACGCCTTATCCGGCCGACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATGGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTgg                                                                                                                                                                                                                                           <  1:1502856/288‑1 (MQ=37)
                                                   tGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCa                                                                                                                                                                                                                                         <  2:1434919/289‑1 (MQ=37)
                                                   tGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCa                                                                                                                                                                                                                                         <  1:418420/289‑1 (MQ=21)
                                                        cgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg                                                                                                                                                                                                                                                                                    <  2:1425115/241‑1 (MQ=18)
                                                                  ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAt                                                                                                                                                                                                                                                                          >  1:1425115/1‑241 (MQ=18)
                                                                    aTCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCCGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATAAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGCATTTGGAGGCCTGATAAGACGCGCACGCGTCGCATCGTGCAGTGGGCTCTGTTGCCGTACgcggc                                                                                                                                                                                                                        >  2:651641/1‑289 (MQ=11)
                                                                        ggCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTaa                                                                                                                                                                                                                    <  1:18265/289‑1 (MQ=38)
                                                                            gtcGGAGGGTGCGGGAGTGTGTAGACCGGATCAGGCGCGCAACGTGGGCATCAGGCGGTGGGAACCATTGCCGGGTGCGGGGGAAAGGCCCTATCCGGGTTAGGGAGGGTTCGGGAATTTGTCGGCCTGATAAGACGCGCAAGCGTCGAATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATTGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGc                                                                                                                                                                                                                <  2:1410812/287‑1 (MQ=11)
                                                                                  gtgggCGGGAATGTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                                                                           <  1:469356/284‑1 (MQ=38)
                                                                                  gggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCTGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                                                                           <  1:651641/288‑1 (MQ=38)
                                                                                   ggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGGAAACGCCTTATcc                                                                                                                                                                                                         >  1:909614/1‑289 (MQ=38)
                                                                                    gTGCGGGAGTTTGTAGTCCTGGTCAGACGCGAAAGCGTCCCATCAGGCATTCGGCGCCATTGCCGGATGCGGCTTAAACGCTTTAGCCGGCCTAGGGAGGGTGCGGGAATGTGTAGGCCTGATAAGAGGCGCAAGCGTCGCATCAGGAAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                                                                                         <  2:1179367/288‑1 (MQ=11)
                                                                                             tttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGa                                                                                                                                                                                                                                                                                         >  1:338085/1‑199 (MQ=37)
                                                                                                    gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                                                                                          <  1:726296/271‑1 (MQ=38)
                                                                                                    gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                                                                                          >  2:726296/1‑271 (MQ=38)
                                                                                                       tGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATGTGTAGGCCTGATAAGACGCGCTAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGTGGGTGCGGGa                                                                                                                                                                                     <  1:209720/289‑1 (MQ=255)
                                                                                                       ggATAAGACGCGAAAGCGTCGCGTCAGGCAGTCGGCACCATTGCGGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                                                                                                                                                                     <  2:1189879/288‑1 (MQ=255)
                                                                                                                 cgcAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                                                                                                                                                                                                                   <  2:1304509/217‑1 (MQ=35)
                                                                                                                 cgcAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                                                                                                                                                                                                                   >  1:1304509/1‑217 (MQ=35)
                                                                                                                   cAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGAGGGGGCGGGAATTTGTTgg                                                                                                                                                                            >  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)
                                                                                                                                                                                                           cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCggg                                                                                 >  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)
                                                                                                                                                                                                                                                    ttGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGACCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCATGACAAGACGCGCACGCATCGCACCAAGc                                        >  1:777848/1‑289 (MQ=255)
                                                                                                                                                                                                                                                         ggATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGGCTGATACGACCCGCACGCGTCGCACCAGGCATCTg                                   >  2:1580758/1‑289 (MQ=255)
                                                                                                                                                                                                                                                            tGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCa                                                                                                                                                                                                                                                                <  1:1390200/66‑1 (MQ=35)
                                                                                                                                                                                                                                                                        gCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGCTGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAAACGCACAAGCGTCCAATCAGGCAGTCACCACCGTTTCCGGATGCGGCGTACACGCCTTATCCGGCCTACGGAGGGTGCGGCACTTTTTTGTCCTGATTACACGCGCAATCGCCGCACCAGGCATCAGCAACCATTGcccgca                    >  2:988596/1‑285 (MQ=11)
                                                                                                                                                                                                                                                                                gggCGTAGGGAGGGTGCGGGTATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGTACCGTTGCCGGATGCGGCGTAAACGCCTTAGCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGcgcg            <  1:599905/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  gCCTAGGGAAGGGGCGGGAATTTGGAGGCCTGAGAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTAGCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGa           <  2:1010971/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  gCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCa                                                                                                                                                                                                                                         <  1:1451461/67‑1 (MQ=40)
                                                                                                                                                                                                                                                                                  gCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCa                                                                                                                                                                                                                                         >  2:1451461/1‑67 (MQ=40)
                                                                                                                                                                                                                                                                                  gCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAACAGGCGCGCAAGTGGCGCATCCGGCATCCGGCCCCCTTGCCGGCTGCGGCGCGa           >  1:235635/1‑288 (MQ=14)
                                                                                                                                                                                                                                                                                    cTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTTGCACCGTTGACTGATGCGGCGCGAgc         >  1:450995/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                       cGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTTCCGGATGCGGCGCGAGCGCCt     >  1:609570/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                           tGGTGCGGGAATTGCGAGGCTTGCAAAGACGCGCAAGGGTGGCATCAGGCAGTCGGAACTGTTGCCGAATGCGGCGTAAACGCCTTATCGGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGTTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTAt  <  2:1518088/288‑1 (MQ=14)
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                
GCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTAT  >  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: