Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,731,291 G→A intergenic (‑134/+309) rrsG ← / ← clpB 16S ribosomal RNA/ClpB chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,731,2910GA92.3% 31.1 / ‑3.4 13intergenic (‑134/+309)rrsG/clpB16S ribosomal RNA/ClpB chaperone
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (7/5);  total (7/6)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCCGGCGGCGTGT  >  NC_000913/2731238‑2731427
                                                     |                                                                                                                                        
gCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTATTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                                                    >  2:958868/1‑140 (MQ=18)
  gACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTATTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCt                                                  <  1:1010987/140‑1 (MQ=25)
      ttAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTATTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtc                                              >  2:893352/1‑140 (MQ=38)
       tAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTATTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtct                                             <  1:91919/140‑1 (MQ=255)
        aaGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTTGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                                            <  1:887693/140‑1 (MQ=255)
          gAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTATTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTg                                          <  2:872769/140‑1 (MQ=255)
                  aTCTTCGAGTGCCCACACAGATTGTCTGATAAATTATTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAtt                                  <  1:458784/140‑1 (MQ=255)
                       cGAGTGCCCACACAGATTGTCTGATAAATTATTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg                             >  2:261356/1‑140 (MQ=255)
                               cacacaGATTGTCTGATAAATTATTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAccc                     >  2:945107/1‑140 (MQ=255)
                                acacaGATTGTCTGATAAATTATTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAcccc                    <  1:478109/140‑1 (MQ=255)
                                           cTGATAAATTATTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCcggc         >  2:22343/1‑140 (MQ=255)
                                            tGATAAATTATTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCcggcg        >  1:93301/1‑140 (MQ=255)
                                            tGATAAATTATTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCcggcg        >  1:990570/1‑140 (MQ=255)
                                                  attattAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCCGGCGGCgtgt  >  1:377336/1‑140 (MQ=255)
                                                     |                                                                                                                                        
GCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCCGGCGGCGTGT  >  NC_000913/2731238‑2731427

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: