Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 675,414 C→A 24.3% H133N (CAT→AAT)  ybeL → DUF1451 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913675,4140CA24.3% 45.1 / 5.1 27H133N (CAT→AAT) ybeLDUF1451 family protein
Reads supporting (aligned to +/- strand):  ref base C (16/4);  new base A (2/5);  total (18/9)
Fisher's exact test for biased strand distribution p-value = 2.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCAGGAGCTGGCAGACATCACCGATAAAACGCAGCTTGAATGGCGCGAAGTTTTCCAGGACCTCAATCATCATGGGGTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGTCTGCGAGAAATGTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGT  >  NC_000913/675276‑675561
                                                                                                                                          |                                                                                                                                                   
tgGCAGGAGCTGGCAGACATCACAGATAAAACGCAGCTTGAATGGCGCGAAGTTTTCCAGGACCTCAATCATCATGGGGTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGGCTGCGAGAAATGTCACTTCCATCTCCCGAt                                                                                                                                           <  4:2310/149‑1 (MQ=255)
   cAGGAGCTGGCAGACATCACCGATAAAACGCAGCTTGAATGGCGCGAAGTTTTCCAGGACCTCAATCATCATGGGGTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGTCTGCGAGAAATGTCACTTCCATCTCCCGATCTa                                                                                                                                        >  3:194431/1‑149 (MQ=255)
                  aTCACCGATAAAACGCAGCTTGAATGGCGCGAAGTTTTCCAGTACCTCAATCATCATGTGGTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGTATTCGAGCAATGTCACTTACATCTCCCGATCTACACACCGGAAGTGCt                                                                                                                         >  4:308173/1‑149 (MQ=255)
                     aCCGATAAAACGCAGCTTGAATGGCGCGAAGTTTTCCAGGACCTCAATCATCATGGGTTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGTCTGCGAGAAATGTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGAc                                                                                                                      >  6:348096/1‑149 (MQ=255)
                     aCCGATAAAACGCAGCTTGAATGGCGCGAACTTTTCCAGGACCTCAATCATCATGGGGTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGTCTGCGAGAAATTTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGAc                                                                                                                      >  7:187501/1‑149 (MQ=255)
                                gCAGCTTGAATGGCGCGAAGTTTTCCAGGACCTCAATCATCATGGGGTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGTCTGCGAGAAATGTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGACGCTATGCCCg                                                                                                            >  4:205470/1‑148 (MQ=255)
                                   gCTTGAATGGCGCGAAGTTTTCCAGGACCTCAATCATCATGGGGTTTATCACAGCGGAGAAGTGGTCGCGCGGGGAAATATGGTCTGCGAGAACTGTCACTTCCATCTCCACATATACACACCGAAAGTCCAGACGCAATCCCCGaaag                                                                                                        >  7:302709/1‑148 (MQ=255)
                                    gttGAATGGCGCGAAGTTTTCCAGGACCTCAATCATCATGGGGTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGTCTGAGAGAAATGTCACTTCAATCTCCCGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAAtg                                                                                                       >  5:223850/2‑149 (MQ=255)
                                       taaTGGCGCGAAGTTTTCCAGGAGCTCAATCATCATGGGGTTTATCACAGCGGAGAAGGGGTCGGGCTGGGAAATCTGGTCTGCGAGAAATGTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAAtgtg                                                                                                     <  3:308173/147‑1 (MQ=255)
                                             ctcgAAGCTATCCAGGACCTCAATCATCAAGGGCTTTATCACAGCGGAGAAGTGGTCTGGCTGGGAAATATTTTTTTTTTTAATGTTAAGTTGAATATCGAGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAATGTGGTCATGa                                                                                              <  8:187501/147‑1 (MQ=255)
                                                gAAGTTTTCCAGGACCTCAATCATCATGGGGTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGTCTGCGAGAAATGTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGACGATATGCCCTAAATGTGGTAATGACCa                                                                                           >  5:227055/1‑149 (MQ=255)
                                                      ttCCAGGACCTCAATCATCATGGGGTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGTCAGCGAGAAATGTCACTTCCATCTCCCGAACTACACACCGCAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACAACTTCCa                                                                                     >  1:163867/1‑149 (MQ=255)
                                                         cAGGACCTCAATCATCATGGGGTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGTCTGCGAGAATTGTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTcaa                                                                                     >  1:313591/1‑144 (MQ=255)
                                                          aGGACCTCAATCATCATGGGTTTTATCACAGCTGAGAAGTTGTCGGGCTTGGAAATATTGTATGTTTTAAATTTAACTTGAATATCGAGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAgaga                                                                                 <  4:194431/149‑1 (MQ=255)
                                                               cTCAATCATCATGGGGTTTATCACAGCGGAGAAGTGGTCGGGCTGCGAAATCTGGGATGCCAGAAATGTCACTTCAATATCCCGATCTACACACCGGAAGTGCTGACGATATGCCCGAAATCTGGTCATGACCAGTTCAAGAGActcac                                                                            >  1:126567/1‑145 (MQ=255)
                                                                           ggggTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGTCTGCGAGAAATGTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGACGCTATGACCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGt                                                                 >  1:175185/1‑148 (MQ=255)
                                                                             gttttATCACAGCGGAGAAGTGGTCGGGCTTGTAAATCTTTTTTTTTTTAATTTTAAGTTGAATATCGAGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTaa                                                               <  2:175185/146‑1 (MQ=255)
                                                                              gTTTATCACAGCGTAGAAGTGGTCGGGATGTGAAATCTTTTTTTTTTTAATTTTAAGTTGAATATCGAGATCTACACACCGGAAGTGCTGACGCTATGACCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATa                                                             <  2:313591/149‑1 (MQ=255)
                                                                                            agaAGTGGTCGGGCTGGGAAATCTGGTCTGCGAGAAATGTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGACGCAATGCCCGAAATGTCGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTATGATTTAGca                                                >  4:147754/1‑148 (MQ=255)
                                                                                                 tGGTCGGGCTGGGAAATCTGGTCTGCGAGAAATGTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGACTCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGGCTGATTTCGCACAGCaa                                          <  5:348096/149‑1 (MQ=255)
                                                                                                   gTCGGGCTGGGAAATCTGTTTTTTGTTAAATTTAAGTTGAATATCGAGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAAt                                         <  6:227055/148‑1 (MQ=255)
                                                                                                             gAAATCTGGTCTGCGAGAAATGTCACTTCCATCTCCCg                                                                                                                                             <  2:94256/38‑1 (MQ=255)
                                                                                                             gAAATCTGGTCTGCGAGAAATGTCACTTCCATCTCCCg                                                                                                                                             >  1:94256/1‑38 (MQ=255)
                                                                                                                        tGCGAGAAATGTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAATGTAGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATg                    >  3:194639/1‑148 (MQ=255)
                                                                                                                           gagaAATGTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAAt                >  2:282189/1‑149 (MQ=255)
                                                                                                                                aTGTCCCTTCCATCTCCCGATCTACACACCGGAAGTCCAGACGCTATGCCCGAAATGTGCTGATGACCAGTTCCAGCGACGCCCGGTTCAGCCGTAATAGACTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTcccc           >  5:46958/1‑149 (MQ=255)
                                                                                                                                   tCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGc        >  4:112265/1‑149 (MQ=255)
                                                                                                                                         ccATCTCCCGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGt  >  2:82491/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
TGGCAGGAGCTGGCAGACATCACCGATAAAACGCAGCTTGAATGGCGCGAAGTTTTCCAGGACCTCAATCATCATGGGGTTTATCACAGCGGAGAAGTGGTCGGGCTGGGAAATCTGGTCTGCGAGAAATGTCACTTCCATCTCCCGATCTACACACCGGAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGT  >  NC_000913/675276‑675561

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑