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

CGCCATAAATCGAAGTTAATGCGATTACGGCATGCTTATGATCAGGAATGTTAATGCCTGCTATACGGGCCACTATGCACTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAAT‑GGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTAA  >  CP009273/3435758‑3436044
                                                                                                                                                 |                                                                                                                                              
cGCCATAAATCGAAGTTAATGCGATTACGGCATGCTTATGATCAGGAATGTTAATGCCTGCTATACGGGCCACTATGCACTCCTACTATTTAATTTGTACGTTTCATGCTGAAAAGCCCGTTTTCAGGATACTCAAAT‑GGAAACGcaca                                                                                                                                            <  1:291777/149‑1 (MQ=255)
               ttAATGCGATTACGGCATGCTTATGATCAGGAATGTTAATGCCTGCTATACGAGCCACTATGCACTCCTACTATTTAATATGTACGCTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAAT‑GGAAACGCACAAACATACAAAAGAtt                                                                                                                             >  1:186785/1‑149 (MQ=255)
                tAATGCGATTACGGCATGCTTATGATCAGGAATGTTAATGCCTGCTATACGGGCCACTATGCACTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAAT‑GTAAACGCACAGACATACAAAAGATtg                                                                                                                            <  1:193922/149‑1 (MQ=255)
                            ggCATGCTTATGATCAGGAATGTTAATGCCTGCTATACGGGCCACTATGCACTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAAT‑GGAAACTCACAGACATACAAAAGATTGGCTGGCTAATCt                                                                                                                >  2:266788/1‑149 (MQ=255)
                               aTGCTTATGATCAGGAATGTTAATGCCTGCTATACGGGCCACTATGCAGTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAAT‑GGAAACTCACAGACATACAAAAGATTGGCTGGCTAATCTAGc                                                                                                             >  2:306468/1‑149 (MQ=255)
                                                                    gCCACTATGCACTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAAT‑GGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGGAAAATAAGCGaaaa                                                                        >  2:280726/1‑149 (MQ=255)
                                                                      cACTATGCACTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAAT‑GGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGaaaaaa                                                                      >  2:286115/1‑149 (MQ=255)
                                                                              aCTCCTACTATTTAATATGTACGTTCCATTCTGAAAAGCCCGTTTTCAGGATACTCAAAT‑GGAAACTCGCAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGCAAAAATCAGCCtt                                                              >  1:80262/1‑149 (MQ=255)
                                                                                    actaTTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAAT‑GGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAGGAAAAATATGCGAAAAAAGCAACCTTGGCGCt                                                        >  1:104079/1‑149 (MQ=255)
                                                                                         ttAATGTGTACATTCAATGATGAAAAGCCCGTTTTGAGGATATTCATAT‑GGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCATGAAAAATATAAAAAAAAATCAGCCTTGGCGCTGTTTa                                                   <  2:22998/149‑1 (MQ=255)
                                                                                             tatGTACGTTCCATGCTGAAAAGCCCCTTTTCAGGATACTCAAAT‑GGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCt                                               >  1:107664/1‑149 (MQ=255)
                                                                                                  acgtgccagtcggAAAAGCGCTTCTTCAGGAT‑CTCCAATGGGAAACGCAAAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGTTCAACCCAACTTTGCAAGAAAAATATTCGAAAAAATCAGCCTTGGCGCGGTTTATGCTtcggc                                          <  2:112218/137‑1 (MQ=255)
                                                                                                              gAAAAGCCCGTTTTCAGGATACTCAAAT‑TGAAACTCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCaa                              >  1:89840/1‑149 (MQ=255)
                                                                                                               aaaaGCCCGTTTTCAGGATACTCAAAT‑GGAAACTCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCaaa                             >  2:9189/1‑149 (MQ=255)
                                                                                                                    cccGTTTTCAGGATACTCAAAT‑GGAAACTCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATcaca                        <  1:306468/149‑1 (MQ=255)
                                                                                                                    cccGTTTTCAGGATACTCAAAT‑GGAAACTCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATcaca                        <  1:266788/149‑1 (MQ=255)
                                                                                                                           tCAGGATACTCAAAT‑GGAAACGCACAGGCATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAACTCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGa                 <  1:183609/149‑1 (MQ=255)
                                                                                                                           tCAGCACACCCAAAT‑GGAAACGCACAGACATACAAAAGATTCGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGa                 <  2:117487/149‑1 (MQ=255)
                                                                                                                                aTACTCAAAT‑GGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCa            <  2:110151/149‑1 (MQ=255)
                                                                                                                                    tCAAAT‑GGAAACTCACAGACATAGAAAAGATTGGCTGGATAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCAc        >  1:172723/1‑149 (MQ=255)
                                                                                                                                           ggAAACTCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTaa  <  1:256070/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CGCCATAAATCGAAGTTAATGCGATTACGGCATGCTTATGATCAGGAATGTTAATGCCTGCTATACGGGCCACTATGCACTCCTACTATTTAATATGTACGTTCCATGCTGAAAAGCCCGTTTTCAGGATACTCAAAT‑GGAAACGCACAGACATACAAAAGATTGGCTGGCTAATCTAGCCAGCTCAACCCAACTTTGCAAGAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTAA  >  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: ‑