Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,287,040 G→A intergenic (+56/+26) ychS → / ← rttR putative uncharacterized protein YchS/small RNA RttR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,0400GA95.2% 56.1 / NA 22intergenic (+56/+26)ychS/rttRputative uncharacterized protein YchS/small RNA RttR
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (8/12);  minor base C (1/0);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTC  >  NC_000913/1286818‑1287277
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
cATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGTTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGAtg                                                                                                                                                                                                                    <  1:73145/250‑1 (MQ=255)
              tGGGCAGAATATTTGATTGCGGATTAGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCCGCAACGCTCTCTCGCTGACGCTCGAGCCGAAGCTTCGACGAAGCTTCTCATCCTTCCCCCCATGGTCAGAATCTTTACTTGCGGATTCGCTTGAGGGTTCGGGGACTTTTGGCAGGGATGGTGGAGGGAGAAgg                                                                                                                                                                                                      >  1:328090/1‑250 (MQ=14)
                           tGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGTAAGGATTACTCAGCGCt                                                                                                                                                                                         <  1:269469/250‑1 (MQ=255)
                                        gTTTGAGAATTCCGGGCCTTTTGAAAGTGATGGTCGTCGGGGAAGGATTCCTCAGCGCTGCGCGCTTCGACCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGTGCCGAACCTTGGTCGCAGCTTCTCATCCTTCCCCGCATGAGCGGAGTATTTGATTGCGGGGTCGCTTAGGAGTTCAGGGCCTTTTGCCAGTGATGGTGGTGAGGGACGGCTTACTCAGCGCGGCGAGATTCGcact                                                                                                                                                                           >  2:313495/1‑248 (MQ=9)
                                                     ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGTATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGcc                                                                                                                                                              >  1:83364/1‑251 (MQ=255)
                                                           tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCCCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGCCTTTTGAAAGTGATGGTGGAGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGc                                                                                                                                                        >  1:280643/1‑251 (MQ=255)
                                                             tGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCa                                                                                                                                                       <  1:19166/250‑1 (MQ=255)
                                                                                                                   gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTCAGAGTTCAGGGACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGCGCTTCGCCCTTCGGGTAGTGGCCTGCGGCGACGCTCTCTCGCTGCCGCTCGCGTCGAACCTTGGTCGAAGCtaccactccttccc                                                                                                >  2:144861/1‑238 (MQ=17)
                                                                                                                       cGTTGCCTGCGGCAACGCTCGCTCGCTGACGCGCGAGCCGAACCTTAGTCGAAGCTTCGCATCCTTCCCCGCAGGGGCAGAATAGTTGATTGCGGATTCGCGTAAGAGTTCAGGGACTTTAGAAAGTGAGGGTGGTGGGGGAAGGATTACGCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCGCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                            <  1:313495/251‑1 (MQ=17)
                                                                                                                            ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAACTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGTATGGGc                                                                                       <  1:77596/251‑1 (MQ=255)
                                                                                                                                    aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGTCAGAATAtt                                                                               <  1:102177/251‑1 (MQ=255)
                                                                                                                                              cactgtctcccggacggAACCTTAGTCGAATCTTCTCATCCTTCCGGGGATGGGCAGAATATTTGATTGCGGAGTCGCTGAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTGCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCATAATATTTAATTGCGGa                                                                     <  2:280643/235‑1 (MQ=17)
                                                                                                                                                    cgcgcGAGCCGATCCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGAGTGCGGCGTCACTTAAGAGTTCTGGGTCTTTTGAGAGTGTTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGGCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                               <  2:287280/247‑1 (MQ=11)
                                                                                                                                                                         gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGTTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGCGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTg                                           >  1:129042/1‑250 (MQ=255)
                                                                                                                                                                                        tCCCCGCGTGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                             <  1:18811/249‑1 (MQ=255)
                                                                                                                                                                                           ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAt                                                                    >  1:143136/1‑207 (MQ=255)
                                                                                                                                                                                           ccGCATGGGCAGAATATTTGATTGCGGATGCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGGCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAGTGCGGAt                                                                    <  2:143136/207‑1 (MQ=255)
                                                                                                                                                                                             gCATGGGCATAATATTTGATTGCGGATTCGCGTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGCCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGt                                                                <  2:188743/209‑1 (MQ=255)
                                                                                                                                                                                             gCATGGGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGt                                                                >  1:188743/1‑209 (MQ=255)
                                                                                                                                                                                                     aGAATATTTGATTGCGGATGCGCTTAAGAGTTCAGGGACGTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGTGGCCTGCGGCAACGCTCTCTCGCCGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa              <  2:129042/251‑1 (MQ=255)
                                                                                                                                                                                                           tttGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGTGCTGCGCGCTTCGCGCTTCGGGTCGTTGCTTGCGGCAACGCTCTCTCGCGGGCGCCCGAATCGAAACTTGGTCGAAGCTTCTCATCATTACCCGCATGGGCGGGATATTTAAATTCGGATTTGTTGGGAAGTTCCGGGGCTTTTGAAAGTGATGGAGGGGGGGGAAGGATTCGatcgttc        >  1:152501/1‑245 (MQ=17)
                                                                                                                                                                                                                 tGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTc  >  1:117494/1‑251 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
CATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTC  >  NC_000913/1286818‑1287277

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