Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,210,839 A→T 17.2% K20M (AAG→ATG)  rpsU → 30S ribosomal subunit protein S21

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,210,8390AT17.2% 64.2 / 8.2 35K20M (AAG→ATG) rpsU30S ribosomal subunit protein S21
Reads supporting (aligned to +/- strand):  ref base A (8/20);  new base T (4/2);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGC  >  NC_000913/3210697‑3210978
                                                                                                                                              |                                                                                                                                           
aaTTGCAGGAAAACTCCGCACCAGTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCCAACATGAGACGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCaa                                                                                                                                             <  2:290213/143‑1 (MQ=255)
  ttGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGc                                                                                                                                           >  1:320802/1‑143 (MQ=255)
      agggaaaTTGCGCACCAGTTTGAAGTAAGCTGGCGGTGCTGCCAGCGGGAAACCGAGTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTc                                                                                                                                       <  2:594978/139‑1 (MQ=255)
        tAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGACTTAATCAAAGGTGAGAGGCACATGCCGGCAAGTAACGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCCCGCGTTCCt                                                                                                                                     >  1:774802/1‑143 (MQ=255)
                 gCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCATGCGAAAAAg                                                                                                                            >  1:383740/1‑143 (MQ=255)
                  cACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCATGCGTTCCTGCGAAAAAGc                                                                                                                           >  1:96922/1‑143 (MQ=255)
                  cACCATTTTGAAATAAGCGGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGc                                                                                                                           <  2:759295/143‑1 (MQ=255)
                    ccATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAg                                                                                                                         <  1:510647/143‑1 (MQ=255)
                    ccATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAACCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAg                                                                                                                         <  1:290375/143‑1 (MQ=255)
                              aTAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTTCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGc                                                                                                               >  1:672459/1‑143 (MQ=255)
                                     ggCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGtt                                                                                                        >  2:775412/1‑143 (MQ=255)
                                     ggCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGtt                                                                                                        <  1:597582/143‑1 (MQ=255)
                                         ttGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtc                                                                                                    <  2:444687/143‑1 (MQ=255)
                                          tGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGGGGCTCCGCCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcg                                                                                                   <  2:580591/143‑1 (MQ=255)
                                                cagcagCAAAGCTAATCAATCAAAGCTGAGAGCCACATGCCGGTAAGTAAACTACGTGCAAACGAGCCGTTCGACGTAGCTCTCCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGa                                                                                             <  1:56472/143‑1 (MQ=255)
                                                      cAAACAGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTa                                                                                       <  2:619444/143‑1 (MQ=255)
                                                       aaaCCGAATTCATCAAAGTGGAGGGGCACATGCCGGTAATTAAAGTACGGGAAAACGAGCCGGTCGACGGAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAGAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTAt                                                                                      <  2:308141/143‑1 (MQ=255)
                                                         cccGACTTAAACAAAGGCGACAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGa                                                                                    <  1:11039/142‑1 (MQ=255)
                                                         aCCGAATTAATCAAAGGGGAGAGGCCCATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGGAGCTCTGCGTCGCTGCAAGCGTTTCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGa                                                                                    <  1:39353/143‑1 (MQ=255)
                                                             aaTGAATCCAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTGCGTCGTCGTGAGTTGTATGagaaa                                                                                >  2:504463/1‑143 (MQ=255)
                                                              aTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAc                                                                               <  1:750444/143‑1 (MQ=255)
                                                                              agGCACCTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCATGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTaa                                                               <  2:734203/143‑1 (MQ=255)
                                                                                    aTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTATGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAgcgcgc                                                         <  2:692775/143‑1 (MQ=255)
                                                                                             aTTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCATGCGTTCCTGCGAAAATGCAGGAGTTCTGGCGTAAGGTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTACGCGCGCTAAAGCtta                                                >  2:445262/1‑142 (MQ=255)
                                                                                             aTTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTc                                                <  2:654957/143‑1 (MQ=255)
                                                                                                        tGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCATGCGTTCCTGCGAAAAAGCAGGTGGTCTGGAGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGaagc                                     >  1:39405/1‑141 (MQ=255)
                                                                                                        tGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGACAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGCACGTAAGCGCGCTAAAGCTTCTGCAGTGaagc                                     >  2:749169/1‑141 (MQ=255)
                                                                                                                     ttCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCGAAAGCTTCTGCAGTGAAACGACACGCgaagaa                        >  1:141778/1‑143 (MQ=255)
                                                                                                                      gtggcGTAGCTCTGTGTAGCTTCAAGCGGTCCTGCGAAAAAGCAGCTGTTCTGGCGGAAGTTTGTCGTCGGGAGTCCTATGAAAGACCGACTACCGAACGTAAGCGCTCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGaaa                       <  2:633866/139‑1 (MQ=255)
                                                                                                                             aactcgGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTc                <  1:725718/141‑1 (MQ=255)
                                                                                                                             aGCGCTGCGGCGCTTCTAGCGTTCCTGCGAAAAATGAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGATCGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTc                <  2:570384/143‑1 (MQ=255)
                                                                                                                               ctctGCGTCGCTTCATGCGTTCCTGCGAAAAAGCAGATGTTCTGGCGGAAGTTCGCCGGCGAGAGGTCTATGAAAAACCGACTACCGAACGTACGGGCGCTACAGCTGCTGCAGTGAAACGTGACGCGAAGAAACTggagaac              >  2:78608/1‑138 (MQ=255)
                                                                                                                                    cGTCGCTTCATGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGaaaa         <  1:492633/143‑1 (MQ=255)
                                                                                                                                    cGTCGCGTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGaaaa         <  2:740502/143‑1 (MQ=255)
                                                                                                                                           tCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTGTGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAACCGTCACGCGAAGCAACTGGCTCGCGAAAaagcacgc  >  2:603098/1‑143 (MQ=255)
                                                                                                                                           tCAAGCGTGCCTGCGAAAAAGCAGGTGGTCTGGCGGCAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAacgcacgc  <  2:732289/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGC  >  NC_000913/3210697‑3210978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: