Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP034943 2,978,569 C→T 14.9% noncoding (2232/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,5690CT14.9% 0.1 / 7.3 34noncoding (2232/2928 nt)EO946_RS1543023S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (4/24);  new base T (3/2);  total (7/27)
Fisher's exact test for biased strand distribution p-value = 5.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCCCTTGTGCCTTTACACTCTGCGAATGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCAACATTCAATATCAGGCTG  >  NZ_CP034943/2978430‑2978691
                                                                                                                                           |                                                                                                                          
gcgcccTTGTGCCGTTACACTCGGCGAAGGATTGCCAACCAGTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCACTCCAACGGCCCACCTGACACTGTCTCCCCGCCCGATAAGggcggc                                                                                                                            <  2:2924609/137‑1 (MQ=1)
  tcgattttgCCTGTACACTGGGCGAATGTTTTCCAACCATTCTGAGGGAAGCTTTGGGCGCCTCCGTTAGCTTTTTGGAGGCGAGCGCCCCAGTAAAACTGCCCACCTGACACTGTGTCCCCGCCGGTTAAGGGcggcgg                                                                                                                          <  2:171243/133‑1 (MQ=1)
       tgggCCTTTACACTCTGCGAATGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGCGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATACGGGCGGCGGGTTAg                                                                                                                     <  1:1516517/137‑1 (MQ=0)
       tgggCCTTTACACTCTGCGAATGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAg                                                                                                                     <  2:4789185/137‑1 (MQ=1)
       tgggCCTTTACACTCTGCGAATGAGTTCCAACCATTCTGACGGAACCTTGGGGCGCCGCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACGGTCTCCCCGCCCGATAAGGGCGGCGGGTTAg                                                                                                                     <  2:330593/137‑1 (MQ=1)
       tgggCCCTTACACTCGGCGAATGCTTTCCAACCATTCTGAGGGAACCTTTGGGCGCCGCCGGTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAg                                                                                                                     <  1:2638467/137‑1 (MQ=1)
              ttccacTCTGCGAATGATTTCCAACCATTCTGAGGGCACCTTTGGGCGCCTCCGTTACCGTTTAGGAGGCGCCCGCCGCAGTCAACCTGCCCCCCTGACCCTGGGTCCCCGCCCGCTAAGGGCGGTGGGTTAGCAGGTCa                                                                                                              >  2:3778744/4‑140 (MQ=1)
                  accctGCGAATGATTTCCAACCAGTCTGAGGGAACCTTTGGGCGCCTCCGTTACCATTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACCCTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATAc                                                                                                          <  1:2238518/137‑1 (MQ=0)
                     cttcGAATGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGTGGGTTAGAAGGTCAATACAGc                                                                                                       >  1:1226545/4‑140 (MQ=1)
                       gctaaTGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCa                                                                                                     <  1:2594480/137‑1 (MQ=1)
                       gctaaTGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTGAGGAGGCGCCCGCCCCAGTCACACTGCCCACCTGACTCTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCa                                                                                                     <  2:4035514/137‑1 (MQ=1)
                       gctaaTGATTTCAAACCATTCTGAGGGATCCTTTGGGCGCCTCCTTTACGTTTTAGGAGGCGACCTCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCGGATAAGGGCGGTTGGTTAGAAGGTCAATACAGCCa                                                                                                     <  1:493279/137‑1 (MQ=1)
                          aaggATGTCCTACCATTCTGAGGGAACCTTTGGGCGCCGCCGTTCCCCGGGAGGAGGCGACCGCCCCAGTCAAACGGCCCACCTGACACTGGCTCCCCGCCCGATAAGGGCGGCGGGTTACAAGGTCAATACAGCCAggg                                                                                                  <  2:4583979/137‑1 (MQ=1)
                              atatcaaaCCATTATGAGGGAACCTTTGGGCGCCTCCGTTAGCTTTTAGGAGGCGACCGCCCCCGACAAACTGCCCACCACACACTGTCTCCCCCCCCGATCAGGGCGGCGCGTTAGAAGGTCAATACAGCCAGGgtagt                                                                                              >  2:3663804/7‑140 (MQ=1)
                                        cttttGATGGTACCGTTTGGCGCCTCCGTTTCCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACCGCCAGGGTAGTATCCCACCGa                                                                                    <  2:1453836/136‑1 (MQ=1)
                                                               ctacGTTACCGGTTAGGAGGCGACCGCCCCCGTCTAACTGCCCACCTTACACTGTCTCCCCTCCCGATAAGGGCTGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGACGCTGGCGCTc                                                             >  1:4967750/4‑140 (MQ=1)
                                                               atcagTTACCTTTTAGGAGGCGACCGCGACAGTCAAACTGCCCACGTGAAACTGTCTCCCCGCCCGATAAGTGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTc                                                             <  1:4917014/136‑1 (MQ=1)
                                                                 cctttACCTTTTAGTAGGCGACCGCCCCAGTCAAACTTCCCACCTGACACTGTCTCCCCGCCCGATCAGGGCGGCGTGTTAGAAGGTCATTTCAGACATGGTTGTATCCGACCGATTCCTCCACCGAAGCTGGCGCTCCg                                                           <  2:1590113/137‑1 (MQ=1)
                                                                         ttttttgaggCGTACGCCCCTTTCCAACTGCCCTCGTGACACTGTCTCCCCGCGCGATAAGGGCGGCGGGTTAGAATGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCAGGTTTCCaa                                                   <  1:2259267/134‑1 (MQ=1)
                                                                         ttgtAGGAGGCGACCGCCCCAGTCAAACGGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCaa                                                   <  2:3272655/137‑1 (MQ=1)
                                                                         ggttcggaggCGCCCGCCCCAGGCACACGGCCCACCAGACACTGTCTCCCCGCCCCATACGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTAGCCCACCGAGGCCTCCACCGAAGCGGGCGCTCCGGTTTCCaa                                                   <  2:3850947/135‑1 (MQ=1)
                                                                             ccgacgCGACCGCCCCAGTCACACGGCCCACCCGACCCGCTCGCCCCGCCCGATCACGGCGGCGGGTTAGAAGGTCAATCCAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCt                                               <  2:3336179/135‑1 (MQ=1)
                                                                                       cggcccAGTCATACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGAGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATTACTCCACCGAAGCTGGCGCTCAGGTTTCCAAGGCTCCTACCTATc                                     <  1:1093880/137‑1 (MQ=1)
                                                                                                cacaCTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAAGACAGCCAGGGTAGTATCCCACCGAGGCCTCCACCGAAGCTGGCGCTCCGGGTTCCAAGGCTCCTACCTATCCTGTACAAg                            <  2:4891280/137‑1 (MQ=0)
                                                                                                 acacggCCCACCGGACACGGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGATGCTGGGGCTCCGGGTTCCAAGGCTCCTACCTATCCTGTACAAGc                           <  2:846146/135‑1 (MQ=1)
                                                                                                 aagcTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGc                           <  1:531731/137‑1 (MQ=1)
                                                                                                   acggCCCACCTGACCCTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTg                         <  2:4303301/137‑1 (MQ=1)
                                                                                                        cccccTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCa                    <  2:682216/137‑1 (MQ=1)
                                                                                                        cccccTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCa                    >  1:4033521/4‑140 (MQ=1)
                                                                                                         cacgttacagTGTCTCCCCGACAGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCAGCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTTTACCAGCTGTACCaa                   <  2:4205503/134‑1 (MQ=1)
                                                                                                         caacTGACACTGTCTCCCCGGCCGATACGGGCGGCGGGTTATAAGGTCAATACAGCCATGGTAGTATGCCACCGATGCCTCCTCATAAGCTGGCGATCAGTTTTCCAAGGATCCTACCTATCCTGTACAAGCTGTACCaa                   <  2:2292352/137‑1 (MQ=1)
                                                                                                               acccTGTCTCCCCGCCCGATAAGGGCGGTGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCAACAGTCa             >  1:4472623/4‑140 (MQ=0)
                                                                                                                 acggTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCt                                               <  2:1945911/101‑1 (MQ=1)
                                                                                                                 acggTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTCTCCTGTACAAGCTGTAGCAACATTCAat           >  2:1880744/4‑140 (MQ=0)
                                                                                                                        ccgcgcccgATAAGGGCGGTGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCAACATTCAATATCAGGc    <  2:4444875/137‑1 (MQ=1)
                                                                                                                          gagcacgATAAGGGCGGCGGGTTAGAAGTTCAATACACCCAGTTTATTATCCAAACGATGCCTCCACAGAAGCTGTAGATCCGGTTTGCAAGGATCATACCTATCCTGTACTAGCTGTACCAACATTAAATATAAGGCTg  <  2:3913481/135‑1 (MQ=0)
                                                                                                                                           |                                                                                                                          
GCTCCCTTGTGCCTTTACACTCTGCGAATGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCAACATTCAATATCAGGCTG  >  NZ_CP034943/2978430‑2978691

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