Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,435,902 G→T 43.4% intergenic (‑72/+75) rpsM ← / ← rpmJ 30S ribosomal subunit protein S13/50S ribosomal subunit protein L36

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,435,9020GT43.4% 2.7 / 16.7 21intergenic (‑72/+75)rpsM/rpmJ30S ribosomal subunit protein S13/50S ribosomal subunit protein L36
Reads supporting (aligned to +/- strand):  ref base G (5/7);  new base T (6/3);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CGCCATAAATCGAAGTTAATGCGATTACGGCATGCTTATGATCAGGAATGTTAATGCCTGCTATACGGGCCACTATGCACTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAATGGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTAA  >  CP009273/3435758‑3436044
                                                                                                                                                 |                                                                                                                                              
cGCCATAAATCGAAGTTAATGCGATTACGGCATGCTTATGATCAGGAATGTTAATGCCTGCTATACGGGCCACTATGCACTCCTACTATTTAATTTGTACGTTTCATGCTGAAAAGCCCGTTTTCAGGATACTCAAATGGAAACGcaca                                                                                                                                            <  1:291777/149‑1 (MQ=255)
               ttAATGCGATTACGGCATGCTTATGATCAGGAATGTTAATGCCTGCTATACGAGCCACTATGCACTCCTACTATTTAATATGTACGCTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAATGGAAACGCACAAACATACAAAAGAtt                                                                                                                             >  1:186785/1‑149 (MQ=255)
                tAATGCGATTACGGCATGCTTATGATCAGGAATGTTAATGCCTGCTATACGGGCCACTATGCACTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAATGTAAACGCACAGACATACAAAAGATtg                                                                                                                            <  1:193922/149‑1 (MQ=255)
                            ggCATGCTTATGATCAGGAATGTTAATGCCTGCTATACGGGCCACTATGCACTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAATGGAAACTCACAGACATACAAAAGATTGGCTGGCTAATCt                                                                                                                >  2:266788/1‑149 (MQ=255)
                               aTGCTTATGATCAGGAATGTTAATGCCTGCTATACGGGCCACTATGCAGTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAATGGAAACTCACAGACATACAAAAGATTGGCTGGCTAATCTAGc                                                                                                             >  2:306468/1‑149 (MQ=255)
                                                                    gCCACTATGCACTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAATGGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGGAAAATAAGCGaaaa                                                                        >  2:280726/1‑149 (MQ=255)
                                                                      cACTATGCACTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAATGGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGaaaaaa                                                                      >  2:286115/1‑149 (MQ=255)
                                                                              aCTCCTACTATTTAATATGTACGTTCCATTCTGAAAAGCCCGTTTTCAGGATACTCAAATGGAAACTCGCAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGCAAAAATCAGCCtt                                                              >  1:80262/1‑149 (MQ=255)
                                                                                    actaTTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAATGGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAGGAAAAATATGCGAAAAAAGCAACCTTGGCGCt                                                        >  1:104079/1‑149 (MQ=255)
                                                                                         ttAATGTGTACATTCAATGATGAAAAGCCCGTTTTGAGGATATTCATATGGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCATGAAAAATATAAAAAAAAATCAGCCTTGGCGCTGTTTa                                                   <  2:22998/149‑1 (MQ=255)
                                                                                             tatGTACGTTCCATGCTGAAAAGCCCCTTTTCAGGATACTCAAATGGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCt                                               >  1:107664/1‑149 (MQ=255)
                                                                                                  acgtgccagtcggAAAAGCGCTTCTTCAGGATCTCCAATGGGAAACGCAAAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGTTCAACCCAACTTTGCAAGAAAAATATTCGAAAAAATCAGCCTTGGCGCGGTTTATGCTtcggc                                          <  2:112218/137‑1 (MQ=255)
                                                                                                              gAAAAGCCCGTTTTCAGGATACTCAAATTGAAACTCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCaa                              >  1:89840/1‑149 (MQ=255)
                                                                                                               aaaaGCCCGTTTTCAGGATACTCAAATGGAAACTCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCaaa                             >  2:9189/1‑149 (MQ=255)
                                                                                                                    cccGTTTTCAGGATACTCAAATGGAAACTCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATcaca                        <  1:306468/149‑1 (MQ=255)
                                                                                                                    cccGTTTTCAGGATACTCAAATGGAAACTCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATcaca                        <  1:266788/149‑1 (MQ=255)
                                                                                                                           tCAGGATACTCAAATGGAAACGCACAGGCATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAACTCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGa                 <  1:183609/149‑1 (MQ=255)
                                                                                                                           tCAGCACACCCAAATGGAAACGCACAGACATACAAAAGATTCGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGa                 <  2:117487/149‑1 (MQ=255)
                                                                                                                                aTACTCAAATGGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCa            <  2:110151/149‑1 (MQ=255)
                                                                                                                                    tCAAATGGAAACTCACAGACATAGAAAAGATTGGCTGGATAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCAc        >  1:172723/1‑149 (MQ=255)
                                                                                                                                           ggAAACTCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTaa  <  1:256070/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CGCCATAAATCGAAGTTAATGCGATTACGGCATGCTTATGATCAGGAATGTTAATGCCTGCTATACGGGCCACTATGCACTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAATGGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTAA  >  CP009273/3435758‑3436044

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