Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP034943 2,978,224 C→A 14.8% noncoding (2577/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,2240CA14.8% ‑0.5 / 5.8 32noncoding (2577/2928 nt)EO946_RS1543023S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (4/23);  new base A (2/3);  total (6/26)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGCC  >  NZ_CP034943/2978086‑2978358
                                                                                                                                          |                                                                                                                                      
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)
           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)
                                                       ccagcgcCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCaa                                                                                <  2:1225802/137‑1 (MQ=1)
                                                       ccagcgcCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCaa                                                                                <  2:1462213/137‑1 (MQ=1)
                                                       ccagcgcCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCaa                                                                                <  2:1796060/137‑1 (MQ=1)
                                                                 cgccGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGt                                                                      <  1:1758396/137‑1 (MQ=1)
                                                                         agcgACCGCACGGTCTCACGACGTTCTGAACCCAGCGCGCGTACCGCTTTAATGCGCGAACAGCCCAACCCTTGGGACCGACGACACCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTgg                                                              <  2:4668958/137‑1 (MQ=1)
                                                                                 aattGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCAAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTgg                                                      >  1:1171938/4‑140 (MQ=1)
                                                                                  acggTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTACTGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTggg                                                     <  2:2578478/137‑1 (MQ=1)
                                                                                  acggTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTggg                                                     <  2:4139493/137‑1 (MQ=1)
                                                                                  acggTCTCACGACGTTCTGAACCCAGCGCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTggg                                                     <  2:1827593/137‑1 (MQ=1)
                                                                                  acggTCTCACGACGGTCTGAACCCAGCTCGCGTACCGCTTTACTGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTggg                                                     <  2:4415419/137‑1 (MQ=0)
                                                                                  accgTCTCACCCCGTTCAGCACCCAGCCGCCGGAACGCCTTTATGGGCGTACAGCCCAACCCTTGGGAACGCATACACCCCCACGCAGCGCTGACCCGCAATCGATGGTCCAAACCCGCCCGTCGTTGTTGACTCGTggg                                                     <  2:1780327/137‑1 (MQ=0)
                                                                                            gaagTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGc                                           <  2:866846/137‑1 (MQ=1)
                                                                                                ttatGAACCCAGCTCGCGTAACGCTGTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCACCATGATGCGATGACCATACATCGAGGTGACAAACCTCCCCGACGATGTGGACACTTGGGGGAGATAAGCCTGt                                       >  2:1953627/4‑140 (MQ=1)
                                                                                                     aagccAGTTCGCGTACCGCTTTAATGGGCGAACAGGCCAACCCTTGGGACCGACTACAGCCCCAGGTTGCGATGAGCCGACATCGAGGTGCCATAACTCCCCGTCGATGTGGACTCTTGGGGGAGATATGCCTGTTATcc                                  <  2:5390252/137‑1 (MQ=1)
                                                                                                           gcgcgcgTACCGCTTTAAGGGGCGAACAGCCCAACCCGTGGGACCGACGACCGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCgggg                            <  2:4054061/137‑1 (MQ=1)
                                                                                                              cgagTACCTCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAg                         <  2:765582/137‑1 (MQ=1)
                                                                                                                     cggtttAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGTTGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTAt                  <  1:2356097/137‑1 (MQ=1)
                                                                                                                     cgatttAATGGGCGAACAGCCCACCCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTAt                  <  2:587955/137‑1 (MQ=1)
                                                                                                                         ttcaTGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGt              <  2:1510555/137‑1 (MQ=1)
                                                                                                                          tactGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGtt             <  2:346860/137‑1 (MQ=1)
                                                                                                                          tactGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGAGGGGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGtt             <  2:4792378/137‑1 (MQ=0)
                                                                                                                                     caccccAACCCTTGGGACCGACTACCGCCCCAGGATGCGATTAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGATATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGcc  <  1:5762156/137‑1 (MQ=0)
                                                                                                                                          |                                                                                                                                      
CCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTGCGCCCGCGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGCC  >  NZ_CP034943/2978086‑2978358

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