Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 892,443 G→C 9.3% A159A (GCG→GCC) ‡ rimK → ribosomal protein S6 modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913892,4430GC 9.3% 146.4 / 6.4 43A159A (GCG→GCC) ‡rimKribosomal protein S6 modification protein
Reads supporting (aligned to +/- strand):  ref base G (17/22);  new base C (2/2);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTC  >  NC_000913/892297‑892537
                                                                                                                                                   |                                                                                              
ggCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAggcgg                                                                                               >  4:236124/1‑148 (MQ=255)
          tgcctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGt                                                                                     >  4:385770/1‑148 (MQ=255)
           gcctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTg                                                                                    <  1:910520/148‑1 (MQ=255)
           gcctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTg                                                                                    <  4:121011/148‑1 (MQ=255)
           gcctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTg                                                                                    <  3:837685/148‑1 (MQ=255)
               gTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTg                                                                                <  3:385770/148‑1 (MQ=255)
                tCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGa                                                                               <  3:236124/148‑1 (MQ=255)
                tCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGa                                                                               <  1:756932/148‑1 (MQ=255)
                              cATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTCAATCACGCttt                                                                         >  2:855218/1‑140 (MQ=255)
                              cATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTCAATCACGCttt                                                                         <  1:855218/140‑1 (MQ=255)
                              cATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGACCGCGGAAAGCGtca                                                                                   >  2:537192/1‑128 (MQ=255)
                              cATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGACCGCGGAAAGCGtca                                                                                   <  1:537192/130‑3 (MQ=255)
                                                aGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGa                                                                               >  2:281725/1‑116 (MQ=255)
                                                aGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGa                                                                               <  1:281725/116‑1 (MQ=255)
                                                 gCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGaa                                                                                          >  3:259593/1‑104 (MQ=255)
                                                 gCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGaa                                                                                          <  4:259593/104‑1 (MQ=255)
                                                 gCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGAc                                                                              >  2:937270/1‑116 (MQ=255)
                                                 gCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGAc                                                                              <  1:937270/116‑1 (MQ=255)
                                                     tttAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGtcaat                                                                                 >  2:300746/1‑106 (MQ=255)
                                                     tttAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGtcaat                                                                                 <  1:300746/110‑5 (MQ=255)
                                                           cGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAatat                                    >  4:193288/1‑148 (MQ=255)
                                                             aCATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGAt                                                                                  <  1:1108976/100‑1 (MQ=255)
                                                             aCATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGAt                                                                                  >  2:1108976/1‑100 (MQ=255)
                                                              cATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTg                                                                                    <  2:314044/97‑1 (MQ=255)
                                                              cATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTg                                                                                    >  1:314044/1‑97 (MQ=255)
                                                                    cGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGAt                                                                                  <  3:1014807/93‑1 (MQ=255)
                                                                    cGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGAt                                                                                  >  4:1014807/1‑93 (MQ=255)
                                                                       tggtgCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGa                                                                               >  2:628157/1‑93 (MQ=255)
                                                                       tggtgCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGa                                                                               <  1:628157/93‑1 (MQ=255)
                                                                                 caggtggtCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTc                                                                        <  3:897677/88‑1 (MQ=255)
                                                                                 cTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTc                                                                        >  4:897677/1‑90 (MQ=255)
                                                                                        tCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCg                                                                                      >  1:928069/1‑69 (MQ=255)
                                                                                        tCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCg                                                                                      <  2:928069/69‑1 (MQ=255)
                                                                                        tCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCaa                               <  3:356399/124‑1 (MQ=255)
                                                                                        tCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCaa                               >  4:356399/1‑124 (MQ=255)
                                                                                         cAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGaaa                                                                                         <  1:96126/65‑1 (MQ=255)
                                                                                         cAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGaaa                                                                                         >  2:96126/1‑65 (MQ=255)
                                                                                             ttggttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTc  <  3:193288/148‑1 (MQ=255)
                                                                                                 ttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTGTGGTGCAGGAatat                                    <  4:844384/110‑1 (MQ=255)
                                                                                                 ttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAatat                                    >  3:844384/1‑110 (MQ=255)
                                                                                                    aaGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCat                                                      <  3:816066/89‑1 (MQ=255)
                                                                                                    aaGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCat                                                      >  4:816066/1‑89 (MQ=255)
                                                                                                                             gCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGaatg                                                           <  2:92923/59‑3 (MQ=255)
                                                                                                                             gCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAAcg                                                           >  1:92923/1‑59 (MQ=255)
                                                                                                                                                   |                                                                                              
GGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTC  >  NC_000913/892297‑892537

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