Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP034943 2,982,148 G→T 18.2% noncoding (372/1550 nt) EO946_RS15435 ← 16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0349432,982,1480GT18.2% ‑3.0 / 4.5 26noncoding (372/1550 nt)EO946_RS1543516S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (4/17);  new base T (2/3);  total (6/20)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCTTTCTGGTTAGGTACCGTCAAGGTACCGCCCTATTCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCATTACCTCACCAACTAGCTAATGCGCCGC  >  NZ_CP034943/2982010‑2982278
                                                                                                                                          |                                                                                                                                  
tgtcTTTCTGGTTAGGTACCGTCTCGGTACCGCCCTAGTCGCACGGTACTTGGTCTTCCCTACCAACAGAGCTTTCCGATCCGAAAACCTTCAGCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGa                                                                                                                                   <  2:943281/137‑1 (MQ=1)
             agttACCGTCAAGGTACCGCCCTATTCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAATATTCCCTActgctg                                                                                                                      <  2:1485927/137‑1 (MQ=1)
               gttccGTGAAGGTACCTCCCTATTAGAACGTTACTTGTTCTTCCCGAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGcc                                                                                                                    <  2:1568108/137‑1 (MQ=1)
                      cacggTACCGCCCTATTCGACCGGGACTTGGTCTTCCCTAACAACACAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTa                                                                                                             <  2:3574077/137‑1 (MQ=1)
                                  ctcttCGCACCGTACTTGGTCGTCCCGAACAACAGAGCTTTACGATCCGAAAACCGTCATCACTCACGCGGCGTTGCTCCGTCAGACTGTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGcc                                                                                                 <  2:2555424/137‑1 (MQ=1)
                                  ctcttCGAACGGTACTTGTTCTTCCCTAACACCAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGcc                                                                                                 <  2:2613346/137‑1 (MQ=1)
                                  ctcttCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTGTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTGTCGTCCATAGCGGAAGATTCCCTACTGCTGCCTCCCGGAGGAGTCTGGGcc                                                                                                 <  2:3910735/137‑1 (MQ=1)
                                     ttggAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAACACCTTCATCACTCACGCGGCGTTGCTCCATCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTCGGAGTCTGGGCCgtg                                                                                              >  1:17176/4‑140 (MQ=0)
                                        gatcGGTACTTGTTCGTCCCTAACAACAGAGCTTTACGAACCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCAGTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                           <  1:2430134/137‑1 (MQ=0)
                                        gaccGGTACTTGTTCTTCCCTCACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                           <  1:2850585/137‑1 (MQ=1)
                                        gaccGGTACTTGTTCTTCCCTACCAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                           <  2:4096253/137‑1 (MQ=1)
                                        gaccGGTACTTGTTCTTCCCTACCAACAGAGCTTGACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                           <  1:3458827/137‑1 (MQ=0)
                                        gaccGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGAACCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAATATTCCCTACTGCTGCCGCCCGTAGGAGTCTGGGCCGTGtct                                                                                           <  2:3798834/137‑1 (MQ=0)
                                        gaccGGTACTTGGTCTACCCTACCCACAGAGCCGTCCGAGCCGAAACCCTTCATCACACACGCGGCGTTGCTCCGTCAGACGCTCGTCCATGGCGGAAGAGGCCCTCCTGATGCCCCCCGTAGGAGTCTGGGCCGTGtct                                                                                           <  2:2842996/137‑1 (MQ=0)
                                            gggaCTTGTTCTTCCCTCACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCAGTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGt                                                                                       <  1:1279094/137‑1 (MQ=0)
                                              gacttcgtcttcCCGACCAACAGAGCTTGACGATCCGACAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGCCGTTCGTCCATTGCGGACGAGGGCCTACTGCTGCCTCGCGGCGGAGTCTTGGCCGTGTCTCAGTcc                                                                                     <  2:2176071/133‑1 (MQ=1)
                                                        tcgctcacaacaGAGCTTTACGATCCGAAAACCTTCATAACTCACGCGGCGTTGCTCCGTCAGTCTTTCGTCCAGTGCGGAAGATTCGCTAGTGCTGCGTCCCGGCGGTGTCTGGGCCGTGTCTCAGTCCCAGTGTGGcc                                                                           >  2:2438879/7‑140 (MQ=1)
                                                        tcacTAACAACAGAGCTTTACGATCCGAAAACATTCATCACTCAATCGGATTTGCTCCTACAGACTTTCTTCAATTGCGGAAGATAACCTACAGCTGACTACAGTAGCAGTCTGGTCCGTGTCTCAGTCCCAGTCTggaa                                                                           >  2:4046431/4‑138 (MQ=0)
                                                        tcacTAACAACAGAGCTTTACGATCCGAAAACATTCATAACTCACGCGGCGTTGCTCCGTAAGACTTTCGTCAATTGCGGAAGATTCCCTACTGCTTCATCCCGAAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGcc                                                                           >  1:2609457/4‑140 (MQ=1)
                                                            taccaacaGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGAGGATGATTCCCTACTGCTGCCTCCCGTAGGAGTATGTGCCGTGTCTCAGTCCCAGTGTGGCCGATc                                                                       <  2:905527/137‑1 (MQ=1)
                                                                 aaatagatttACCATCCGAAAAAATTAATAACACACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAATATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCct                                                                  >  2:3202208/8‑140 (MQ=1)
                                                                       ctgtACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAgg                                                            <  2:1778658/137‑1 (MQ=1)
                                                                                   aataCCGTCATCACTCCCGCGGCGTTGCTCCGTCAGTCTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCAt                                                <  2:1384926/137‑1 (MQ=0)
                                                                                               ccacaggcggcgTGGCTCCCTCCGACTGTCGTCCCGTGCGGAATATGCCCTCCGGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCAACGTTGCCTTGGt                                    <  2:2451420/134‑1 (MQ=1)
                                                                                                                        ctatCGTGCGTTGCGGAATATTCAATACTGCTCCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCATTACCTCACCAACTAGCTa           >  2:2670337/4‑140 (MQ=1)
                                                                                                                                 atggCGGAAGATTCCCTCCTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCATTACCTCACCAACTAGCTAATGcgccgc  <  1:196847/137‑1 (MQ=1)
                                                                                                                                          |                                                                                                                                  
TGGCTTTCTGGTTAGGTACCGTCAAGGTACCGCCCTATTCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCATTACCTCACCAACTAGCTAATGCGCCGC  >  NZ_CP034943/2982010‑2982278

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