Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP034943 2,979,248 T→C 27.6% noncoding (1553/2928 nt) EO946_RS15430 ← 23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0349432,979,2480TC27.6% ‑3.9 / 9.6 20noncoding (1553/2928 nt)EO946_RS1543023S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (2/12);  new base C (2/3);  total (4/16)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

TCTCCTCGCCTACCTGTGTCGGTTTGCGGTACGGGCACCTCTCACCTCGCTAGAGGCTTTTCTTGGCAGTGTGGAATCAGGAACTTCGCTACTACATTTCGCTCGCCATCACAGCTCAGCCTTATGGGAAACGGATTTGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCAATACCGCGCTTACCCTATCCTCCTGCGTCCCCCCATTGCTCAAATGGTGAGGAGGTGGTACAGGAATATCAACCTGTTGTCCA  >  NZ_CP034943/2979124‑2979383
                                                                                                                            |                                                                                                                                       
tcgcctcGCCTACCTGTGTCGGGTTGCGGTACGGGCACCTCTCACCTCGCTAGAGGCTTTGCTTGGCAGTGTGGAATCAGGAACGTCGCTACTAGATTTCGCGCGCCATCACAGCGCGGCCGGATGGGCAAGGGATTTGc                                                                                                                          >  2:288032/4‑140 (MQ=1)
              tgggtCGGTTTGCGGTACGGGCACCTCTCACCTCGCTAGAGGCTTTTCTTGGCAGTGTGGAATCAGGAACTTCGCTACTACATTTCGCTCGCCATCACAGCTCAGCCTTACGGGAAACGGATTTGCCTATTTCCCAGCCt                                                                                                            <  2:1695444/137‑1 (MQ=1)
                              actggCACCTCTCACCTCGCTAGAGGCTTTTCTTGGCAGTGTGGAATCAGGAACTTCGCTACTACATTTCGCTCGCCATCACAGCTCAGCCTTATGGGAAACGGATTTGCCTATTTCCCAGCCTAACTGCTTGGACGCgg                                                                                            <  2:12076/137‑1 (MQ=1)
                                       tcgcACCTCGCTAGAGGCTTTTCTTGGCAGTGTGGAATCAGGAACTTCGCTACTACAGTTCGCTCGCCATCACAGCTCAGCCTTACGGGAAACGGATTTGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCaaa                                                                                   >  2:1835194/4‑139 (MQ=0)
                                         tccccTCGCTAGAGGCTTTTCTTGGCAGTGTGGAATCAGGAACTTCGCTACTACATTTCGCTCGCCATCACAGCTCAGCCTTATGGGAAACGGATTTGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCAATAc                                                                                 <  2:2800051/137‑1 (MQ=1)
                                         tccccTCGCTAGAGGCTTTTCTTGGCAGTGTGGAATCAGGAACTTCGCTACTACATTTCGCTCGCCATCACAGCTCAGCCTTACGGGAAACGGATTTGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCAATAc                                                                                 <  2:2455510/137‑1 (MQ=1)
                                                 cgagcggCTTTTCTGGGCAGTGTGGACTCAGGAACGTCGCTCCTACATTTCGCTCGCCATCACCGCTGAGCCTTATGGGAAACGGATTGGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCAATACCGCGCTTa                                                                         <  2:3389207/135‑1 (MQ=1)
                                                                                    ttagCTACTACATTTCGCTCGCCATCACAGCTCAGCCTTACGGGAAACGGATTTGCCTCTTTCCCAGCCTAACTCCTTGGACGCGGAAATCCAATACCGCGCGTACCCTCTCCTCCTGCGTCCCCCCATTGCTCAACTgg                                      >  2:2529572/4‑140 (MQ=1)
                                                                                      cgttactacATTTCGCTCGCCATCACAGCTCAGCCTTATGGGAAACGGATTTGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCAATACCGCGCTTACCCTATCCTCCTg                                                            <  1:2282202/113‑1 (MQ=1)
                                                                                      cgttactacATTTCGCTCGCCATCACAGCTCAGCCTTATGGGAAACGGATTTGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCAATACCGCGCTTACCCTATCCTCCTg                                                            >  2:2282202/4‑116 (MQ=1)
                                                                                       gcgactacCTGTCGCTCGCCATCACAGCTCAGCCTTATGGGCAACGGAGTTGCCTAGGTCCCAGCCTAACTGCTGGGACGCGGATATCCACTACCGCGCTTACCCTATCCTCCGGCGTCCCCCCATTGCTCAAATGGTga                                   <  2:1431514/137‑1 (MQ=1)
                                                                                           cgacctttCGCTCGCCATCACAGCTCAGCCTTATGGGAAACGGATTTGCCTATTTCCCAGCCTACCTGCTTGGACGCGGATATCCAATACCGCGCTTACCCTATCCTCCTGCGTCCCCCCATTGCTCAAATGGTgaggag                               <  2:1073128/135‑1 (MQ=0)
                                                                                            tacctgtcgctcgcCATCACAGCTCAGCCTTATGGGAAACGGATTTGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCAATACCGCGCTTACCCTATCCTCCTGCGTCCCCCCATTGCTCAAATGGTGaggagg                              <  1:2404675/134‑1 (MQ=1)
                                                                                            tacctgtcgctcgcCATCACAGCTCAGCCGTATGGGAAACGGATTTGCCGAGTTCCCAGCCTAACTGCTTGGACGCGGAGATCCAATACCGCGCTTACCCTATCCTCCTGCGTCCCCCCAGTGCTCAAATGGTGaggagg                              <  2:2217476/134‑1 (MQ=1)
                                                                                            taaaTTTCGCTCGCCATCACAGCTCAGCCTTAAGGGAAACGGATTTGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCAATACCGCGCTTACCCTATCCTCCTGCGTCCCCCCATTGCTCAAATGGTGaggagg                              <  1:1538374/137‑1 (MQ=1)
                                                                                               atctcgcgcgcCATCACAGCTCAGCCGTATGGGAAACGGATTTGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCAATACCGCGCTTACCCTATCCGCCTGCGTCCCCCCATTGCTCAAATGGTGAGGAggtgg                           <  1:3119351/137‑1 (MQ=0)
                                                                                                                 gcacAGCCTTACGGGAAACGGAGTTGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCAATACCGCGCTTACCCTATCCTCCTGCGTCCCCCCATTGCTCAAATGGTGAGGAGGTGGTACAGGAATATCAACCtg         <  2:1617078/137‑1 (MQ=0)
                                                                                                                       ccgtATGGGAAACGGATTTGCCTATTTCCCAGCCTTACTGCTTTGACGCGGATATCCAATACCGCGCTTACCCTATGCTCCGGCGTCCCCCCATTGCTCAAATGGTGATGATGTGGTACAGGAATATCAACCTGTTGTcc   <  1:3120816/137‑1 (MQ=1)
                                                                                                                        ctgaTGGGAAACGGATTTGCCTATGTCCCAGCCTAACTGCTTGGACGCGGATATCCAATACCGCGCTTACCCTATCCTCCTGCGTCCCCCCATTGCTCAACTGGTGAGGAGGTGGTACAGGAATATCAACCCGTTGTCCa  <  2:3468276/137‑1 (MQ=0)
                                                                                                                        ctgaTGGGAAACGGATTGGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCAATACCGCGCTTACCCTATCCTCCTGCGTCCCCCCATTGCTCAAATGGTGAGGAGGTGGTACAGGAATATCAACCTGTTGTCCa  <  1:2965360/137‑1 (MQ=1)
                                                                                                                            |                                                                                                                                       
TCTCCTCGCCTACCTGTGTCGGTTTGCGGTACGGGCACCTCTCACCTCGCTAGAGGCTTTTCTTGGCAGTGTGGAATCAGGAACTTCGCTACTACATTTCGCTCGCCATCACAGCTCAGCCTTATGGGAAACGGATTTGCCTATTTCCCAGCCTAACTGCTTGGACGCGGATATCCAATACCGCGCTTACCCTATCCTCCTGCGTCCCCCCATTGCTCAAATGGTGAGGAGGTGGTACAGGAATATCAACCTGTTGTCCA  >  NZ_CP034943/2979124‑2979383

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