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

TGTTACCTGAACTTCACCCTGGACATGGGTAGATCACCTGGTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCACGGTA  >  NZ_CP034943/2980036‑2980295
                                                                                                                            |                                                                                                                                       
tggtACCTGAACTTCACCCTGGACATGGGTAGATCACCTGGTTACGGGTCTACGACCACGTACTCATGCGCCCTATTAAGACACGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAAGCTTGCACAGGATCGTAAATCGca                                                                                                                          >  1:480244/4‑139 (MQ=1)
   taacTGAACTTCACCCTGGACATGGGTAGATCACCTGGTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGt                                                                                                                       <  2:3958303/137‑1 (MQ=1)
                  cttgACATGGGTAGATCACCTGGTTTCGGGTCTCCGACCTCGTACTCATGCGCCCTATTCAGACTCTCTTTAGCTGCGGCTCCGCATCTTCATCTTATAATTTCACGTGCTAGTAACTAGCCTTTTCTTTATACaattgt                                                                                                        <  2:1319426/137‑5 (MQ=0)
                  cttgACATGGGTAGATCACCTGGTTTCGGGTCTACGACCACGTACTCTTGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCTTTTTATTATAAAAAAgg                                                                                                        <  2:1760880/137‑1 (MQ=1)
                    ggccATGGGTAGATCACCTGGTTGCGGGTCTACGACCACGTACTCATGCGCCCTCTTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACAGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCa                                                                                                      <  2:1696331/137‑1 (MQ=0)
                          ggttAGATCACCTGGTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCAt                                                                                                <  1:4970631/137‑1 (MQ=1)
                                    ccgggTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTaa                                                                                      <  2:1178545/137‑1 (MQ=1)
                                              ggacTACGACCAGGTACTTATGCGCACTTTTCAGACTCGATTTCGCTGCGGGTCCGCATCTTCTGCTTTACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCTTTAAGGGGCTCTGa                                                                            <  1:1383161/137‑1 (MQ=1)
                                              ggacTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGa                                                                            <  2:2306780/137‑1 (MQ=1)
                                              ggacTACGACCACGTACACATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGa                                                                            <  2:3567091/137‑1 (MQ=1)
                                               gtatACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACACAAGGCACGCCATCACCCGTTAACGGGCTCTGac                                                                           <  2:3667832/137‑1 (MQ=1)
                                               gtatACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGac                                                                           <  1:1286222/137‑1 (MQ=1)
                                               gtatACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGac                                                                           <  2:2083090/137‑1 (MQ=1)
                                               gtatACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGac                                                                           <  1:2597509/137‑1 (MQ=1)
                                               gtatACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGAATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGac                                                                           <  1:2018903/137‑1 (MQ=1)
                                                         acctACTCATTCGCCCTATTCAGACTCGCTTTCGCTGCGGATCCGCATCTTCTGCTTAAACTTGAACCGGATCGTAACTCGCCGGTTCAATCAACAAAAAGCACGCAATCACCCGTTAACGGGCTCTGACTACATGTAgg                                                                 >  2:4744693/4‑140 (MQ=1)
                                                             tctgaTGCGCCCTATTCAGACTCGCTTTCGCTGGGGGTCCGCTTCTTCTGCTTAACGTTGCACGGTATCTTAAGTCGCCGGTTCATTCTTCAAAAGGCACGCCATCACCCGTTAACGGGATCTGACTACTTGTAGGcaga                                                             <  2:5070462/136‑3 (MQ=1)
                                                             acacATGCGCCCTATTCAGACTCGGTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACCAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGcaca                                                             <  2:3775328/137‑1 (MQ=0)
                                                                 atccgcCCAATTAAGACTCGCTTTCGCTGCGCATCCGCAACAGCTTCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATGCTACAAAACGCACGCCATCACCCGTAAAAGGGCTCTGACTACTTTTAGGCACACGGt                                                         >  2:2778952/4‑140 (MQ=1)
                                                                       ccgaTTCCGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACAGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAgg                                                   <  1:3499457/137‑1 (MQ=0)
                                                                         tagtCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGAt                                                 <  2:1831112/137‑1 (MQ=1)
                                                                              agtcTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACAGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCtt                                            >  1:2493679/4‑140 (MQ=1)
                                                                                   cgatttCGCTGCGGCTCCGCATCTTCTTCTTAACCTTGCACAGGATCGTAACACGCCGGTTCATTATACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGAAGGCACACGGTTTCAGGATCTCTTTCaat                                       >  2:1158271/4‑138 (MQ=1)
                                                                                    gcgttCGGTGCGGCTCCGCATCTTCTGCTTAACCTTGGACGGGATCGGAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTc                                      <  1:2987002/137‑1 (MQ=0)
                                                                                    gcgttCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTc                                      <  1:4085982/137‑1 (MQ=1)
                                                                                    gcgttCGCTGCGGCTCCGCATCTTCTGCTTAACCGTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTACCGGGCTCTGACTACTTGTAGGCACACGGTGTCAGGATCTCTTTCCCTc                                      <  2:4091675/137‑1 (MQ=1)
                                                                                     ctatCGCTGCGGCTCCGCATCTTCTTATTAAAATTGAAAGGGATCGTAACGCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTcc                                     >  2:4185826/4‑140 (MQ=1)
                                                                                            tgctgatCCGCATCTTCTGCTTAACGTTGCACGGGATCGTAACTGGCCGGTTCATTCTACAAAAGGCACGCCATCACCGGTTAACGGGCTCTGACTACTTGTAGGCTCACGGTTTCAGGATCTCTTTCACTCCCCTTCCg                              <  1:2739059/134‑1 (MQ=1)
                                                                                              ccgcaccGCATCTTCTGCTTAACCGTGCACCGGATCGTAACACGCCGGTTCATTCTACAAAAGGCACCCCATCACCCGTTAACGGGCACTGACTACTTGTAGGCACACGGTTTCACGATCTATTACACTCCCCTTCCggg                            <  1:5104408/135‑1 (MQ=1)
                                                                                               ggatCCGCATCTTCTTCTTAACCTTGCAAGTGATAGTATCTCGCCGGTTCATTATAAAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCgggg                           <  1:3555152/137‑1 (MQ=1)
                                                                                                           tcggCTTAACCTGGCACGGGATCGTACCTCGCCGGGTCAGTCGACCAAAGGCACGCCAACACCCCTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCt               <  2:5065054/137‑1 (MQ=1)
                                                                                                                  aaacTTGCTCGGGATCGTATCTCGCCTGTTCATTCTACACAAGGCACGCCATCTCCAGTTAACGTGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTGTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTc        <  1:3120986/137‑1 (MQ=1)
                                                                                                                   acgttGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCa       <  2:2230658/137‑1 (MQ=1)
                                                                                                                        gcgcGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCACGGTa  <  1:404090/137‑1 (MQ=1)
                                                                                                                        gcccGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATGACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCACGGTa  <  1:3053981/137‑1 (MQ=1)
                                                                                                                            |                                                                                                                                       
TGTTACCTGAACTTCACCCTGGACATGGGTAGATCACCTGGTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCACGGTA  >  NZ_CP034943/2980036‑2980295

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