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

CTCATACCGGCATTCTCACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCCACATCCTTTTCCACTTAACGTATACTTTGGGACCTTAGCT  >  NZ_CP034943/2979495‑2979759
                                                                                                                                     |                                                                                                                                   
ctgaTACCGGCATTCTCACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCATCTTCGGTGATACTTTTAGCCCCGGTACATTATCGgcg                                                                                                                               <  1:3846256/137‑1 (MQ=0)
  cagaCCGGCATTCTCACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCGTCACAGCCCTTAGAACGCTCTCCTACCACGGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCa                                                                                                                             <  2:3670909/137‑1 (MQ=0)
    tagcGGCATTCTCACTTCTAAGCGCGCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCaga                                                                                                                           <  2:2152653/137‑1 (MQ=1)
          cagtctcACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACt                                                                                                                     <  2:2330039/137‑1 (MQ=1)
           atactcACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTATCGGCGCAGAGTCACTc                                                                                                                    <  2:4480060/137‑1 (MQ=1)
                   ttcgacgcgcgcCACCAGTCCGTCCGGTCGGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTg                                                                                                            <  1:891943/134‑1 (MQ=1)
                   ttcgacgcgcTCCACGAGTCCTTCCGGTCTGGCGTCACAGCCCTTAGAACGCTCTCCTACCACTGGTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTGAGCCCCGGTACCCTTACGGCGCAGAGTCACTCGACCAGTg                                                                                                            <  1:2594371/134‑1 (MQ=1)
                   gtcgaaGCGCCCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTCCCACTGTTCGAAGAACAGTCCGCCGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTg                                                                                                            <  2:3668919/136‑1 (MQ=0)
                     ctcaGCGCTCCACCAGTCCTTCCGGTCTGTCTTCCCAGCCCTGAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGTGCAGAGTCACTCGACCAGTGAg                                                                                                          <  2:2352050/137‑1 (MQ=1)
                     ctcaGCGCTCCACCAGTCCTTCCGGTCTGGCTTGACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTATCGGCACAGAGTCACTCGACCAGTGAg                                                                                                          <  1:2276351/137‑1 (MQ=0)
                     ctcaGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAg                                                                                                          <  2:4003368/137‑1 (MQ=1)
                     ctcaGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCGACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACCTTTTCGGCGCAGAGTCACTCGACCAGTGAg                                                                                                          <  1:4525396/137‑1 (MQ=0)
                     ctcaGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGACCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAg                                                                                                          <  1:3327366/137‑1 (MQ=1)
                     ctcaGCGCTCCACCACTCCTTCCGGTCTGGCGTCACAGCCCTTAGAACGCTCTCCGCCCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAg                                                                                                          <  1:4099222/137‑1 (MQ=1)
                             tctaccaGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACg                                                                                                  >  2:4716554/4‑140 (MQ=1)
                             gccccccgTCCTTCCGGTCGGGCTTCACAGCCCTAAGAACGCTCTCCGACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACg                                                                                                  <  2:3038688/133‑1 (MQ=1)
                              cccccaGTCCGTCCGGTCGGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGc                                                                                                 <  1:2776633/137‑1 (MQ=0)
                                ccccgtccgtccGGTCGGGCGTCACAGCCCTTAGAACGCTCTCCTACCACTGGTCGAAGACCAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCAc                                                                                               <  2:1532439/136‑1 (MQ=1)
                                acgaGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTGACTCGACCAGTGAGCTATTACGCAc                                                                                               <  1:3679789/137‑1 (MQ=1)
                                acaaGTCATTCCGGTCTGTCTTCACATCCCTTATAACGCTCTCCTACCAGTTTTATAAGAACAGTCCGCAGCTTCGGTTATACGTTTATCCCCGGTACATTTTCGTAGCAGAGTCACTCGACCAGTGAGCTATTACGCAc                                                                                               <  1:237307/137‑1 (MQ=1)
                                         tcaggTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTATCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAAt                                                                                      >  1:1938594/4‑140 (MQ=1)
                                              tcgggCTTCACAGCCCTGAGAACGCTCTCCTACCCCGGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAAtggtgg                                                                                 <  2:2996846/137‑1 (MQ=0)
                                               ctcgCTACACAGCCCTGAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCCGCTTCGGTGATCCGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGgc                                                                                <  2:292532/137‑1 (MQ=1)
                                                       accgCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATCCGTTTAGCCCCGGTACATTATCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTa                                                                        <  2:2905482/137‑1 (MQ=0)
                                                           ccattAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGcc                                                                    <  2:1803397/137‑1 (MQ=1)
                                                                       tcaccTACCACTGTTCGATGAACAGTCCGCATCATAGGTGATACGGTTATCCCCGGAACATTTTCGGCGCAGAGTCACACGACCCGAGAGCTATTACGCACACTTTAAATGGTGCCTGCTTCTAAGCCAAAATCGTGGta                                                        >  2:2208451/4‑139 (MQ=1)
                                                                         tcgtACCCCTGTTCGCAGAACAGTCCGCAGCTTCTGTGATACGTTTACCCCCGGTACATTTTCGGCGCAGTGTCACTCGACCAGTGATCTATTACGCAATCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGt                                                      <  2:2639781/137‑1 (MQ=1)
                                                                             cccccTGTTCGATCAACAGTCCGCAGCTTGGGTGATACGTTTAGCCCCGGTACATTTTCGCCGAAGAGTCACTCGACCAGTGAGCTATTACTCAATCTTTTAATGTTGGCTTCTTCTAAGCCAACATACTGGTTGTCTaa                                                  <  1:3888842/136‑1 (MQ=1)
                                                                                acggTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCa                                               <  2:5105807/137‑1 (MQ=1)
                                                                                acggTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTGAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCa                                               <  2:2166888/137‑1 (MQ=1)
                                                                                accgTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCa                                               <  2:184310/137‑1 (MQ=1)
                                                                                 cttttCTAAGAACAGTCCGCAGCTTGGGTGATACGTTTAGCCCCGGTACTTTTTCGGCGCAGTGTCACTCGACCAGTGAGCTATTACGCACTCTTTAATTGGTGGCTGCTTCTAAGCCAACTTCCTGGTTTTCTAAGCaa                                              <  2:537696/137‑1 (MQ=1)
                                                                                   gtacgcagacCAGTCCGCAGCTTCGGTGATACGTTTATCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACGCTTTAAATGGTGGCTCCTTCTAAGCCACCATCCTGGTTGTAGAAGCAACt                                            >  2:1044920/7‑140 (MQ=1)
                                                                                   gtacGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCCCTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAATCAAcg                                            >  2:2089108/4‑139 (MQ=0)
                                                                                                   gctgCTTCTGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTAAGCACTCTTTAATTTGTGTCTTCTTCTAAGCCAACATCATGGTTGTCTAAGCAACTCCACATCCTTTTCCAc                            <  2:4075046/137‑1 (MQ=1)
                                                                                                   gccgCTGCGGTGATCCGAGTCGCCCCGTTACAGTGACGGCCCAGAGTCACGCGACCAGTGAGCACTTACGCACTCTTTAAATGGTGGCTTCTTCTAACCCAACATCCAGGTTGTCTAAGCACCTCCACATCCTATTCCAc                            <  2:3020632/137‑1 (MQ=0)
                                                                                                       ctacGGTGATACGTTTAGACCCTTTAAATTATCCCCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCCACATCCTTTTCCACTTaa                        >  2:4118912/4‑140 (MQ=1)
                                                                                                           ggggATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCCACATCCTTTTCCACTTAACGta                    <  2:2667635/137‑1 (MQ=1)
                                                                                                           ggagATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGATCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCCACATCCTTTTCCACTTAACGta                    <  1:3965383/137‑1 (MQ=1)
                                                                                                                             gggaCATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCCACATCCTTTTCCACTTAACGTATACTTTGGGACCTTAGCt  >  2:4234696/4‑140 (MQ=1)
                                                                                                                                     |                                                                                                                                   
CTCATACCGGCATTCTCACTTCTAAGCGCTCCACCAGTCCTTCCGGTCTGGCTTCACAGCCCTTAGAACGCTCTCCTACCACTGTTCGAAGAACAGTCCGCAGCTTCGGTGATACGTTTAGCCCCGGTACATTTTCGGCGCAGAGTCACTCGACCAGTGAGCTATTACGCACTCTTTAAATGGTGGCTGCTTCTAAGCCAACATCCTGGTTGTCTAAGCAACTCCACATCCTTTTCCACTTAACGTATACTTTGGGACCTTAGCT  >  NZ_CP034943/2979495‑2979759

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