Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP034943 2,978,137 A→T 12.9% noncoding (2664/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,978,1370AT12.9% 1.7 / 5.7 36noncoding (2664/2928 nt)EO946_RS1543023S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (5/25);  new base T (2/2);  total (9/27)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GGGGGCTTCATGCTTAGATGCTTTCAGCACTTATCCCGTCCGCACATAGCTACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCG  >  NZ_CP034943/2978001‑2978271
                                                                                                                                        |                                                                                                                                      
ggtggCTTCATGCTTAGATGCTTTCAGCACTTATCCCCTCCGCACATAGCTACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAAttt                                                                                                                                     >  2:4198451/4‑140 (MQ=1)
               aaatccTTTCAGCACAAATCCCGTCCGCACAAAGCTACACAGAGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATACAGGTCCTCACGTACTAAGTACACCTCCTCTCAAATATCCGGCGCCCGCGACg                                                                                                                      >  2:120340/6‑140 (MQ=1)
                  tggtttCAGAACTTATCCCGTCCGCACATAGCAACCCAGCGTGGCCGTTGGAAGAACACGGGGTACACCCGGGGTGCTTCAATCCGGGGCCTATCTAACTACGCACAGCTCCTTAAAATTTTCCTGCGCCGGCGACgttt                                                                                                                   <  2:275598/137‑4 (MQ=0)
                  tgcgtgcAGCACTTCGCCCGTCCGCACATAGCGACCCCGCGAGGCCCGGGGCAGAACAACTGGGACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGAt                                                                                                                   <  1:3843185/134‑1 (MQ=1)
                            actgagccccTCCGCACCGAGCTCCCCAGCGAGGCCCCGGGCAGAACAACGCCTACACCAGCGGTGCGTCCCTCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGaa                                                                                                         <  2:3453295/134‑1 (MQ=1)
                            acgtATCCCGTCCGCACATCGCTACCCAGCGATGCCCTTCGCAGAACAACTGGTACACCAGCGGTGCGTCCAGCCCGGTCCTCTCCTACTAAGCACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGaa                                                                                                         <  2:4477630/137‑1 (MQ=1)
                            acgtATCCCGTCCGCACATAGCTACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGGACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGaa                                                                                                         <  2:1527800/137‑1 (MQ=1)
                             ctgaTCCCGGCCGCACCTAGCGACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAAc                                                                                                        <  1:2201495/137‑1 (MQ=0)
                                            cagaGCTACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTc                                                                                         <  1:3031604/137‑1 (MQ=1)
                                                gcgaCCCAGCGATGCCCTTGGCAGAACAACTGGTACACCCGCGCTGCGTCCATCCCGGTCCGCTCGTACTAAGGACAGCTCCTCTCAAAATTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGaa                                                                                     <  1:168528/137‑1 (MQ=1)
                                                gcgaCCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGaa                                                                                     <  1:4430876/137‑1 (MQ=1)
                                                         cgctGCCCGTGGCAGAACACCTGGTACACCAGCGGTGCGTCCAGCCCGGGCGTCTCGTGCTAAGGAGAGCGGCGCGCAACTTGCCTGCGCGCGCGACGGATCGGGCCCGCACTGTCTCACGACGTGCGGACCCCAGCTcg                                                                            >  2:4186972/4‑140 (MQ=0)
                                                             gcgcTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCGCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTGTGAACCCAGCTCGCGTa                                                                        <  2:3952262/137‑1 (MQ=0)
                                                             gcccggggCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTa                                                                        <  2:4388865/134‑1 (MQ=1)
                                                              ccgttGGCAGAACAACTGGTACAGCAGCGGTGCGTCCATCCCGGTCGTCTCGTACTAAGGACCGCTCCTCGCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCAAGACGTTCGGAACCCAGCTCGCGTAc                                                                       >  1:2968287/4‑140 (MQ=1)
                                                                        aaaaacTTGTACACCAGCGGTTCGTCCATAAAGGAACTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTACTGAACCCAGCTCGCGTACCGCTTTAATg                                                             >  1:640263/4‑140 (MQ=1)
                                                                              tgctACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGCCGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                                       <  2:4661319/137‑1 (MQ=1)
                                                                              tgctACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                                       <  1:5649845/137‑1 (MQ=1)
                                                                              tgctACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                                       <  2:4640963/137‑1 (MQ=1)
                                                                              tgctACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                                       <  1:467677/137‑1 (MQ=1)
                                                                              tgctACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGGACCGCTTTAATGGGCGaa                                                       <  2:4020454/137‑1 (MQ=1)
                                                                              tgcgaccccAGCGCGGCGTCCATCCCCGTCCTCACGTACCAAGGACAGCTCCTCTCAAATTTCCTGCCCCCGCGACCGATAGGGACCGAACTGTCACACGACGTTCTGAACCCAGCTCGCGAACCGCTTTAATGGGCGaa                                                       <  1:3042538/133‑1 (MQ=1)
                                                                              cgggacacCCGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATGTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                                       <  2:3300131/136‑1 (MQ=0)
                                                                              acggaccccAGCGGTGCGTCCATCCCGCTCCGCTCGTACTAAGCACACCTCCTCGCAAAGTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGaa                                                       <  2:3358332/133‑1 (MQ=1)
                                                                                     cctgCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCa                                                <  2:3512047/137‑1 (MQ=1)
                                                                                     ccggCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCa                                                <  2:3751760/137‑1 (MQ=1)
                                                                                     cccgCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCa                                                <  2:737406/137‑1 (MQ=1)
                                                                                            ggcctccctcccGCTCCTCTCGGACTCAGGACAGCTGCTCTCAAATTTCCTGAGCCCGCGACGGATAGGGACCGAACTGGCTCACGACGTTCTGACCCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTg                                         <  2:4352478/136‑1 (MQ=1)
                                                                                                tccactccGGTCCTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGaa                                                       >  2:3903664/7‑122 (MQ=0)
                                                                                                tccactccGGTCCTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGaa                                                       <  1:3903664/116‑1 (MQ=0)
                                                                                                tccactccGGTCCTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGAc                                     >  1:1733360/7‑140 (MQ=0)
                                                                                                gccctCCCGGTCCTCTCGTACTCAGGACAGCTCCTCTCAAATGTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGAc                                     <  1:541993/136‑1 (MQ=0)
                                                                                                           ccgctcGTACTAAGCACAGCTCCTCTCAAAGTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGCGACCGACTACAGcc                          <  2:4627007/137‑1 (MQ=0)
                                                                                                                           catctcctcTCTAATTTCCTGCGCCCGCGACGGATCGGGCCCGAACTGTCTCACGCCGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACGACAGCCCCAGGATGCGATGAGc          <  1:5029822/137‑1 (MQ=1)
                                                                                                                              ctgctctCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCACAGGATGCGATGAGCCGa       >  1:639076/4‑140 (MQ=1)
                                                                                                                              ctgctcgCAAGTTTCCTGCGCCCGCGACGGGTAGGGACCGAACTGTCTCACGACGGTCTGCACCCAGCTCGCGTACCGCTTTAATGGGCGACCAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGa       >  2:1447534/4‑140 (MQ=1)
                                                                                                                                   ctaaaaTTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTATCACGACGTTCTGAACCCATCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGATCCGACATCg  >  2:3707978/4‑140 (MQ=0)
                                                                                                                                        |                                                                                                                                      
GGGGGCTTCATGCTTAGATGCTTTCAGCACTTATCCCGTCCGCACATAGCTACCCAGCGATGCCCTTGGCAGAACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCG  >  NZ_CP034943/2978001‑2978271

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