Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,147 T→C 9.8% intergenic (‑406/+252) 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,1470TC 9.8% 142.2 / 8.0 60intergenic (‑406/+252)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base T (31/23);  new base C (3/3);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTT  >  NC_000913/4325870‑4326431
                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                
acgtgtCCGGCGTTCGGAGGGTGCGGGAATGTGTAGGCCTGGTAAGAGGCGCAAGCGTCGCATCAGGCAGTCGGGACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGAt                                                                                                                                                                                                                                                                                       <  2:756154/283‑1 (MQ=255)
  ttgtctgtgcgaCGGGGGGGGCGGGAAGTTGTGTGCCGGAAAAGACGCGCAAGCTGCGGGTCAGGCTGTCGGAACGCTTGCCGGATGCGGCGTAAACGCCTCATCCGGACTACGGAGGGGGCGGGAATTTGTAGGCCTGATAAGACGGGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATgcg                                                                                                                                                                                                                                                                                    <  2:802499/278‑1 (MQ=11)
                    gTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCCTCAGGCAGTTGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                                                                                                                                     >  2:203203/1‑286 (MQ=25)
                    gTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                                                                                                                                     <  1:203203/286‑1 (MQ=25)
                       cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCCGGCAGTTGGCACT‑TGTGCCGGATGCAGCGGTAACGCCTTACCCGGc                                                                                                                                                                                                                                                                >  2:875005/1‑288 (MQ=255)
                        gggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTCAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATACGACGCGCAAGCAGCCCATCAGGGAGTTTGCAAT‑GTTGCCGGATTGTGCGGAAACGCCCTaacggcgct                                                                                                                                                                                                                                                              >  1:772742/1‑281 (MQ=11)
                             tttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGGCGCCTCAGGCAGCTTGCACT‑TTTTCCGGCTGCGGCGTCAACGCCGTATCCGGCCTCCggc                                                                                                                                                                                                                                                         >  2:430435/1‑288 (MQ=11)
                             tgtGTAGGCCGGAGAAGAGGAGGAAGCGTGGCATCAGGCGGTAGGCACCATTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATGTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                                                                                                                                                                         <  2:701122/287‑1 (MQ=11)
                                       tGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                                                                                                                                                                                                                               <  1:528477/289‑1 (MQ=38)
                                       agaaaacaCGCGCAAGCGGCGTACGAGGGAGTGGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGAATCAGGCGGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGGTGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                                                                                                                                                                                                                               <  2:727482/282‑1 (MQ=11)
                                        gATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGaa                                                                                                                                                                                                                                              >  1:74779/1‑289 (MQ=38)
                                             gACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACCCCTTATCCGGCCTACGGACGGTGCGGGAATTTTTAGGCCTGCTAAAACCCCCACGCGTCGCATCCGGACGTTCCCCCTCTTTCCCGGAT‑CAGCCTAAACGCCTTACACCGCCCACGGCGGGTGcacgcacatgt                                                                                                                                                                                                                                         >  2:510178/1‑279 (MQ=11)
                                                  gcaagcggcggcgaaGGCAGTGGGGGCAATTGCCGTATGCGGCGGAAACGCCGTATCCGGGCTACGGAGGGTGCGGGAATTGGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc                                                                                                                                                                                                                                     <  2:94496/274‑1 (MQ=11)
                                                   cAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGGGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGCAATTTTAGGcc                                                                                                                                                                                                                                    >  1:727482/1‑288 (MQ=255)
                                                                           aCCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCACGCAGTTGGCACT‑GTTGCCGGCAGCGGCGTAAACGCCTTATCCGGCCTACGGCGGGTACGGGAATTTGTTGGCCTGATAAGACGCGCCAGCGTccgac                                                                                                                                                                                                            >  2:320375/1‑284 (MQ=255)
                                                                              aTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGGATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                                                                                                                                                                         <  2:831491/288‑1 (MQ=9)
                                                                                                           cGGCCTACGGAGGGTGCGGGGATTTGTAGGCCGGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATAAGGCGtt                                                                                                                                                                           <  2:161774/289‑2 (MQ=255)
                                                                                                                                   tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCTGGAATTTGTCGGCCTGATAAGACGCGCAAGCCTCGCATCAGGCACTCGGCACCGTTGTCGGACGCGGCGGTACCGCCTTCTCCGGCCtcccgcgg                                                                                                                                                   >  2:381732/1‑282 (MQ=11)
                                                                                                                                       ggCCTGATAAGACGCGCAAGCGTGGCAGCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGAGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTg                                                                                                                                                <  2:80557/288‑1 (MQ=255)
                                                                                                                                                   cgcgCAAGCGTGTGACCAGACCGTCGGAACGGTTGCCGGAGGGGGCGTAAACGCCGTATCCGGCCTCCGGAGGGGGGGGGAGTTTGTAGGCCGGATAAGACGCGCAAGCGGCGCATCAGTCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCGTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                                                                                                                                    <  2:727208/288‑1 (MQ=11)
                                                                                                                                                                               aCCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCTGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGGAGTCGGCACCGCTGCCGGGTGCGGCGTAATCGGCTTACCCGGCCTCCGGCGGGGTCGGGGATTTGTAGGGCCTATCAGATCGCG‑Caaaggcgtct                                                                                                       >  2:297341/1‑280 (MQ=11)
                                                                                                                                                                               aCCGTGGCCGGATGGGGCGTAAACGCCTTATCCGGCCTACGGATGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTGGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATc                                                                                                     <  2:78636/289‑1 (MQ=25)
                                                                                                                                                                                ccGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCa                                                                                                    >  1:809283/1‑289 (MQ=37)
                                                                                                                                                                                 cGTTGCCGGGTGCGGCGAAAACGGCTTGTCCGGCTGACGGAGGGTGCGGGAATGTGTAGGCCTGAGAGGACGCGCAAGCGGAGCGTCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTAGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAg                                                                                                   <  2:422251/289‑1 (MQ=11)
                                                                                                                                                                                 cGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATccg                                                                                                   >  1:702403/1‑287 (MQ=25)
                                                                                                                                                                                   ttGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGACGCGGCGTAAACGCCTTATCCGGCCTAAGGAGGGTGCGGGAATTTGTAGGCCTGATAAGT‑CGCA‑CCAACGTcccgcctgt                                                                                                  >  2:743854/1‑280 (MQ=255)
                                                                                                                                                                                    tGCCGGATGCGGCGTAAACGCCTCATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCa                                                                                                >  1:227614/1‑289 (MQ=25)
                                                                                                                                                                                             cggcgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                                             <  2:261639/221‑1 (MQ=35)
                                                                                                                                                                                             cggcgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                                             >  1:261639/1‑221 (MQ=35)
                                                                                                                                                                                                cgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCACGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCGGGTGCGGGCATTTGTCGGCCTGATAAGA‑CGCGACACGC‑CCGCCTCAGGCATCTGGCACCGtt                                                                                    >  2:715962/1‑289 (MQ=14)
                                                                                                                                                                                                      cGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGa                                                                              <  2:809283/289‑1 (MQ=39)
                                                                                                                                                                                                            aTCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTCAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATgcggc                                                                        >  1:536850/1‑289 (MQ=17)
                                                                                                                                                                                                                 gCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCg                                                                                                         >  1:189373/1‑251 (MQ=35)
                                                                                                                                                                                                                 gCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGa                                                                   >  1:550154/1‑288 (MQ=40)
                                                                                                                                                                                                                 gCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTGGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCg                                                                                                         <  2:189373/251‑1 (MQ=21)
                                                                                                                                                                                                                  ccTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCCGTCGGCACCTGTGCCGGAAGCGGCGTAAACGGCTTATCCGGACTACGGGGGGGGCGGGTATTTGTTGGCCTGATAAGA‑CGGC‑CCAGAGTCGCACACGGCACTTGGACCCGTTGCCGGATgggcatgcag                                                                   >  2:310307/1‑279 (MQ=11)
                                                                                                                                                                                                                          tGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGAGCGCCTTAt                                                          <  1:221680/288‑1 (MQ=255)
                                                                                                                                                                                                                            gTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGAGCGCCTTATCCg                                                       >  1:711940/1‑289 (MQ=255)
                                                                                                                                                                                                                                 ggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                                            <  1:726227/186‑1 (MQ=35)
                                                                                                                                                                                                                                 ggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                                            >  2:726227/1‑186 (MQ=35)
                                                                                                                                                                                                                                  gAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGAGCGCCTTATCCGGCTTa                                                  >  1:879729/1‑288 (MQ=255)
                                                                                                                                                                                                                                   aaTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCCGATAAGA‑AGCG‑CCAGCGGCGCAACAGGCATCTGGCACCGTTTCCGGATCCGGCG‑CGGAGGCCTTCTCCCGCtctc                                                 >  2:2543/1‑285 (MQ=14)
                                                                                                                                                                                                                                    aTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGAGCGCCTTATCCGGCTTAcc                                                <  1:313819/288‑1 (MQ=255)
                                                                                                                                                                                                                                    aTTTGTAGGCCGGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGAGCGCCTTATCCGGCTTAcc                                                <  1:2543/288‑1 (MQ=255)
                                                                                                                                                                                                                                     tttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACt                                                                                                                                                                                                                                                                                                  <  1:634717/49‑1 (MQ=40)
                                                                                                                                                                                                                                       tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc                                                                                                                                 >  1:621519/1‑207 (MQ=35)
                                                                                                                                                                                                                                       tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc                                                                                                                                 <  2:621519/207‑1 (MQ=35)
                                                                                                                                                                                                                                       tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCCCATTAGGCATCACGCCCCGTTGTCGGATGCTGCG‑CGAGCGCGTTATCCgccctcctcat                                             >  2:774046/1‑278 (MQ=14)
                                                                                                                                                                                                                                            gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACg                                                                                                                                                                                                                                                           <  1:676312/80‑1 (MQ=39)
                                                                                                                                                                                                                                            gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                                                                                                                                                                  <  1:659506/89‑1 (MQ=38)
                                                                                                                                                                                                                                              cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGAGCGCCTTATCCGGCTTACCGATCacaatac                                     >  1:709474/1‑289 (MQ=255)
                                                                                                                                                                                                                                                 aTAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgc                                                                                                                                                                                 <  1:53252/149‑1 (MQ=37)
                                                                                                                                                                                                                                                 aTAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTCCGGCGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGAGCGCCTTATCCGGCTTACCGATCacaatacca                                   >  2:790943/1‑286 (MQ=255)
                                                                                                                                                                                                                                                 aTAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCT‑CGAGCGCCTTGTCCGGCTTCCCGATCacattacaa                                   >  2:135899/1‑288 (MQ=255)
                                                                                                                                                                                                                                                 aTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGAGCGCCTTATCCGGc                                                     >  1:27266/1‑270 (MQ=21)
                                                                                                                                                                                                                                                 aTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGAGCGCCTTATCCGGc                                                     <  2:27266/270‑1 (MQ=21)
                                                                                                                                                                                                                                                      aCGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCACCCGGCATCTGGCACCGTTGCCTGATGCGGCG‑CGAGCGCCTTATCCGGCTTACCCCTCcaaatacaacaagt                              >  2:532733/1‑283 (MQ=14)
                                                                                                                                                                                                                                                           cAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCAcca                                                                                                                                                                                            <  1:554588/128‑2 (MQ=18)
                                                                                                                                                                                                                                                           cAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCAGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCtttt                        >  1:20662/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                               tGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCAGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGGACCGTTGCCGGATGCGGCG‑CGAGCGCCGTATCCGGCTTACCGAACACAATACAATCAGTTCTTTTTCACCAACTCAGATTTCAGt    >  1:502423/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                 gCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTCACGGATGCGGCGCCGAG‑GCCTTCTCCGGCCTACGCATCACAATACTAACCGTCTTTTTTCCCAACATCAGGTTTCCGttt  >  2:818795/1‑289 (MQ=11)
                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                
CCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGA‑CGCG‑CAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCG‑CGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTT  >  NC_000913/4325870‑4326431

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