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,0400GA94.1% 38.9 / NA 17intergenic (+56/+26)ychS/rttRputative uncharacterized protein YchS/small RNA RttR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/6);  minor base C (1/0);  total (11/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAG  >  NC_000913/1286827‑1287275
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
ccgtttGGGCAGAATATTTGATTGGGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                                                                                        <  1:321288/247‑1 (MQ=255)
                                                     gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTTGCCCTTAGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTACTCAGCGCTgcgcg                                                                                                                                                                                  >  2:272110/1‑220 (MQ=255)
                                                     gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcg                                                                                                                                                                                  <  1:272110/220‑1 (MQ=255)
                                                                           ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTACGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGt                                                                                                                             >  1:75398/1‑251 (MQ=255)
                                                                              ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGAATTTTGAACGTGATGGTGGTGGGGGAAGGATTACTCAGAGCGGCGCGCTTCGCCCTTCGGGTCGATGCCTGAGGCAACGCTCTCTCGCTGACGCTCCAGTCgg                                                                                                                          >  1:247028/1‑250 (MQ=14)
                                                                                     ccgcTGCGCCCGTCGTCCTTCGGTTCGGTGCCTGCGGCAACGCTCTCCCGCTGACGCTCGAGTCGATCCTTAGTCGAAGCCTGTCATCCTTGACCGC‑TGCGGCAGAACATTTGATTGCGGATTCGCGAAATAGTTCGGGTACTTTGGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgt                                                                                                                   <  2:24987/250‑2 (MQ=14)
                                                                                                               gTTGCATGCGGCAACGCACTCTCGCTGACGCTCGCGCCGAACATTAGTCGAAGCTTCTCATCCTTCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGTTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCATTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTGCGCTCGCGTCGAACCTTGGTAGAAGCTTCTCATACTTCCCCGCAt                                                                                         >  2:443754/1‑251 (MQ=17)
                                                                                                                                 ctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTGCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGGGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAt                                                                         <  1:443754/249‑1 (MQ=255)
                                                                                                                                                gAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGa                                                         >  1:83637/1‑250 (MQ=255)
                                                                                                                                                    cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                                                                                                 <  1:144944/174‑1 (MQ=255)
                                                                                                                                                    cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                                                                                                 >  2:144944/1‑174 (MQ=255)
                                                                                                                                                     gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTAAGCGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTACGCCCTTCGCCCTTGGTGTCGTTGCTTGCGGCAACGCTCTCTCGCTGGCGCTCGAGGGGAACCTTGGTCGTAGCTTCTCTTCCTTCCCCGCGTGGGGAGGATATTTACTTgtcgattggtgcggttcttcc                                                   >  2:188814/1‑231 (MQ=14)
                                                                                                                                                                gAAGCTTCTCATCCTTCCCCGCATG‑CGCAGAATCTTTGATAGCGGATTCGCTTCAGCGTCCAGGGCCTTTTGAAAGTGATGGTGGTGGGGGAAGGCTTACTCAGCGCTGCGCGCTTCGCACTTCGGGTCGtt                                                                                                                                                               >  2:232953/1‑132 (MQ=14)
                                                                                                                                                                                      aTG‑GGCAGAATATTTGATTGCGGCTTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCGGGCGCTCGCGGCGGACCTTGGTCGACGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTCAGTGCGGAATCGTTGGGTAGTTCAGGGCATTTTGTACGTGATGGTGGGGGGGGaggg                  >  2:123513/1‑248 (MQ=14)
                                                                                                                                                                                                      gATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGCTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGa    >  2:467356/1‑250 (MQ=255)
                                                                                                                                                                                                      gATTGCGGATTCGCTTAAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGATTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAAACTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTATTGGCGGATTCGTTGGGAAGTTCAGGGACTTTTGACAGTGATGGTGGTGGGGGAGGGATTCGCACCTTCGa    >  2:341508/1‑250 (MQ=17)
                                                                                                                                                                                                       attggtgATTCGCATAAGATTTCAGGTACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAg  <  1:306592/245‑1 (MQ=17)
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
CCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATG‑GGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAG  >  NC_000913/1286827‑1287275

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