Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP034943 2,980,160 G→A 15.6% noncoding (641/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,980,1600GA15.6% ‑1.6 / 5.0 26noncoding (641/2928 nt)EO946_RS1543023S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (3/18);  new base A (2/2);  total (5/21)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TGTTACCTGAACTTCACCCTGGACATGGGTAGATCACCTGGTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCAC  >  NZ_CP034943/2980036‑2980291
                                                                                                                            |                                                                                                                                   
tgatACCTGAACTTCAACCTTTAAATAGGTAGAAAACCTGGTTTCGGGTCTACGACCACGAACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACAGGATCGTAACTCGcc                                                                                                                      >  2:2859346/4‑140 (MQ=1)
 gtgaCCTGAACTTCACCCTGGACATGGGTAGATCACCTGGTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCg                                                                                                                     <  1:2374147/137‑1 (MQ=1)
        gaacggcACCCTGGACATGGGTAGATCACCTGGCTGCGGGTCTACGACCACGTACTCATGCGCCCTAGTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGAACGTAACTCGCCGGTTCAtt                                                                                                              <  2:2861042/134‑1 (MQ=1)
                          ggttAGATCACCTGGTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCAt                                                                                            <  2:1372011/137‑1 (MQ=1)
                                    ccgggTTTCGGGTCTACGCCCACGGACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCGTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTaa                                                                                  <  2:2266517/137‑1 (MQ=0)
                                    ccgggTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTaa                                                                                  <  2:912968/137‑1 (MQ=1)
                                     cttgTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACAGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAAc                                                                                 <  2:2295684/137‑1 (MQ=1)
                                              ggacTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGa                                                                        <  2:2771680/137‑1 (MQ=1)
                                                             acacATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGcaca                                                         <  2:1181832/137‑1 (MQ=1)
                                                                       cctagacAGACGCGCTGTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAgg                                               <  1:2034220/134‑1 (MQ=0)
                                                                       acttttCAGACTCGCGTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAgg                                               <  1:2842532/136‑1 (MQ=0)
                                                                         tagtccgACTCGCAGTCGCTGCGGCTCCGCATCTTCTGCGTCACCCTGCACGGGATCGGACCGCGCCGGTTCATTCTACAAAACGCACGCCATCACCCGTTAACGGGCTCGGACTACTGCTAGGCCCCCGGTTTCAGCAt                                             <  2:3829857/134‑1 (MQ=1)
                                                                         tagtCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGAt                                             <  1:3562347/137‑1 (MQ=1)
                                                                         tagtCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGAt                                             <  1:2251484/137‑1 (MQ=1)
                                                                                    gcgttCGCTTCGGCTCCGCATCTTGTGCTTATCCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGGACGCGATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTc                                  <  1:2934704/137‑1 (MQ=1)
                                                                                     ctgtCGCTGCGGCTCCGCATCTTCTGGTTAAGCTTGCACAGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTCCTTGGAGGCACACGGTTTCAGGAGCTCTTTCACTcc                                 >  1:747518/4‑140 (MQ=1)
                                                                                     ctctCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTcc                                 >  1:149967/4‑140 (MQ=1)
                                                                                     ctatCGCTGCGGCTCCGCATCTTCTGCTTAACCGTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACCCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCCCACGGGTTCAGGATCTCTTTCACTcc                                 >  1:821439/4‑140 (MQ=1)
                                                                                                           tcggCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCAGGCCATCACCCGTTAACGGGCTCTGACTCCTTGTAGGCACACGGTTTCCGGATCTCTGTCCCTCCCCTTCCGGGGTGCTTTTCACCt           >  2:1257037/4‑140 (MQ=1)
                                                                                                               ctaaaCCGTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTcc       <  1:2739768/137‑1 (MQ=1)
                                                                                                                 tacccTTGCCCGCGATCGTCACCCGCCGGTTCATTCTACACAAGGCACGCCATCACCCGATAACGGGCTCTGACTACGTCTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTcccg     <  2:1293745/137‑2 (MQ=1)
                                                                                                                   acgttGCACTGGATCTTAACTCTCCGGTTAATTCTTCATAAGGCACGCCATCACACTTTAAAGGGCTCTGAATACTTGTAGGCACACGTTTTGAGGATCTCTTTCACTCCCCTTCAGTGTTGCTTTTCACCTTTCCCTCa   <  1:1349519/137‑1 (MQ=1)
                                                                                                                   acctcccACGGGCTCGTAACCCGCCCGTTCCTGCTACAAACCGCACGCCAACACCCCTTACCGGGCTCTGACTACTGGTAGGCACACGGTTTAAGCATCTCTGTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCGCa   <  2:1503769/134‑1 (MQ=1)
                                                                                                                    ccgtGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCAc  <  2:2432781/137‑1 (MQ=1)
                                                                                                                    ccgtGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCAc  <  1:2981784/137‑1 (MQ=1)
                                                                                                                    ccgtGCACAGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCAc  <  1:1721989/137‑1 (MQ=1)
                                                                                                                            |                                                                                                                                   
TGTTACCTGAACTTCACCCTGGACATGGGTAGATCACCTGGTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCAC  >  NZ_CP034943/2980036‑2980291

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: ‑