Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP034943 2,979,628 T→A 15.4% noncoding (1173/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,979,6280TA15.4% ‑0.9 / 9.7 39noncoding (1173/2928 nt)EO946_RS1543023S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (4/29);  new base A (2/4);  total (6/33)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

CTACTCATACCGGCATTCTCACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCCACATCCTTTTCCACTTAACGTATACTTTGGGACCTTAGCTG  >  NZ_CP034943/2979492‑2979760
                                                                                                                                        |                                                                                                                                    
ctccTCATACCGGCATTCTCACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGTACAGTCCGCAGCTTCCGTGATACGTTAAGCCCCGGTACATTTTCg                                                                                                                                   >  2:4267936/4‑140 (MQ=0)
        acaggCATTCTCACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTTGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCagag                                                                                                                           >  1:2836777/4‑140 (MQ=1)
         gagacagtctcACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGt                                                                                                                          <  2:3086346/133‑1 (MQ=0)
             cagtctcACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACt                                                                                                                      <  2:4384747/137‑1 (MQ=1)
             cagtctcACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACt                                                                                                                      >  2:3762330/4‑140 (MQ=1)
             cagtctcACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACt                                                                                                                      <  1:4918683/137‑1 (MQ=1)
             cagtctcACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTATCGGCGCAGAGTCACt                                                                                                                      <  2:5759710/137‑1 (MQ=1)
             cagtctcACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCGACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACt                                                                                                                      <  2:3284624/137‑1 (MQ=1)
                      gtcgaaGCGCTCCACCAGTCCTTCCGGTCGGGCTTCCCCGCCCTTAGAACGCTCTCCTACCACGGTTCGACGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTg                                                                                                             <  2:2137134/136‑1 (MQ=1)
                        ctcaGCGCTCCCCCAGTCCTTCCGCTCTCGCGACCCACCCCTTAGCACGCTCTCCTCCCACTGTTCGACGAACAGTCCGCAGCTTCGGTGATACCTTTAGCCCCGGTACATTTTCGGCGAAGAGTCACTCGACCAGTGAg                                                                                                           <  1:5031721/137‑1 (MQ=1)
                        ctcaGCGCTCCACCAGTCCTTCCGGTCTGGCTTCCCCGCCCTTAGAACGCTCTCCTCCCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAg                                                                                                           <  2:4251986/137‑1 (MQ=0)
                                tcaaccaGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTGCATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACg                                                                                                   <  2:4704637/137‑1 (MQ=1)
                                tcaaccaGTCCTTCCGGTCTGGCTTCACAGCCATTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACg                                                                                                   <  1:375657/137‑1 (MQ=1)
                                cccccccgTCCTGCCGGTCTGGCCTCACAGCCCGTAGAACGCTCTCCTCCCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACg                                                                                                   <  2:5370674/133‑1 (MQ=1)
                                 cccccaGTCCGTCCGGTCTGGCTTCACAGCCCTGAGAACGCTCTCCGACCCCTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGc                                                                                                  <  1:3236338/137‑1 (MQ=1)
                                  cccccgTCCTTCCGGTCTGGCTTCACCGCCCTTAGAACGCTCTCCTACCCCTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCa                                                                                                 <  2:2970076/135‑1 (MQ=0)
                                                  cttgCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGgc                                                                                 <  2:4906123/137‑1 (MQ=1)
                                                     gcgtCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTATCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGgctgc                                                                              >  1:2169617/4‑140 (MQ=1)
                                                        tcccaGCCCTGAGAACGCTCTCCTACCACTGTTCGACGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTc                                                                           <  2:3596586/137‑1 (MQ=0)
                                                          actgCCCTTAGAACTCTATCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTAAATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTa                                                                         <  1:2729439/137‑1 (MQ=0)
                                                                ctgaGAACGCTCTGCTACCACTGTTCGTAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCaa                                                                   <  1:3748977/137‑1 (MQ=0)
                                                                        gcgctcCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTTTTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTgg                                                           <  2:3796911/137‑1 (MQ=1)
                                                                            gccgaCCACTCTTCGAACCACAGTCCGCAGCTTCGGCGATACGTTTAGCCCCGGTACAGTATCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTCGTTGt                                                       <  2:4354222/136‑1 (MQ=1)
                                                                             ccgaCCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCGATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTc                                                      <  2:5069141/137‑1 (MQ=1)
                                                                                  ccctcttCGAAGAACAGTCCGCCGCTTCCGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGc                                                 <  1:3357913/135‑1 (MQ=0)
                                                                                   acggTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCa                                                <  1:640280/137‑1 (MQ=1)
                                                                                    ctcttCGAAGAACAGTCCGCAGCATCGGTGATACGTTTAGCCCCGGTACATTAACGGCGCAGAGTCACTCGACCAGTGAGATATTACGCACTCTTTAAATGGTGGCTTCTTCTAAGCCAACATCGTGGTTGTCTTAGCaa                                               >  2:2998368/4‑140 (MQ=1)
                                                                                          gacgaaCAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCcac                                         <  1:3201416/137‑1 (MQ=1)
                                                                                           aacaaCAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCcaca                                        <  2:4876404/137‑1 (MQ=1)
                                                                                           aacaaCAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCcaca                                        <  1:2568286/137‑1 (MQ=1)
                                                                                           aacaaCAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTATCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCcaca                                        <  1:4793758/137‑1 (MQ=1)
                                                                                                      gctgCTTGGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCCACATCCTTTTCCAc                             <  1:2312760/137‑1 (MQ=1)
                                                                                                      gctgCTTCGGTGATACGTTTAGCCCCGGTACATTTTGGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTACATGGTGGCTGCTTCTAAGCCAACATCCTGGTGGTCTAAGCAACTCCACATCCTTTTCCAc                             <  2:4701323/137‑1 (MQ=0)
                                                                                                      gccgcgtCGGTGATACGTTTAGCCCCGCTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCCACATCCTTTTCCAc                             <  2:3535541/134‑1 (MQ=0)
                                                                                                      gccgCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCCACATCCTTTTCCAc                             <  1:5101903/137‑1 (MQ=1)
                                                                                                      gccgCTTCGGTGATACGTTTAGCCCCGGTACATGTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTCAGCAACTCCACATCCTTTTCCAc                             <  2:1092944/137‑1 (MQ=0)
                                                                                                              ggagATACGTTTAGCCCCTGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATTGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAATCAACTCCACATCCTTTTCCACTTAACGta                     <  1:5761789/137‑1 (MQ=0)
                                                                                                                    accttaaGCCCCGGTACCTTTTAGGCGCAGAGTCACTCGAACAGTGAGCAATTACGCACTCTTAAAATGGTGGCTGCTTCTAAGCCCACATCCTGGTAGTCCCAGCACCTCCACATCCTTATCCACTTAAAGTAAACttt               >  2:154322/7‑140 (MQ=1)
                                                                                                                                 gttcATTATCGGCTCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCGTGGTTGTCTAAGCAACTCCACATCCTTTTCCACTTAACGTATACTTTGGGACCTTAGCTg  <  1:3020372/137‑1 (MQ=0)
                                                                                                                                        |                                                                                                                                    
CTACTCATACCGGCATTCTCACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCCACATCCTTTTCCACTTAACGTATACTTTGGGACCTTAGCTG  >  NZ_CP034943/2979492‑2979760

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