Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,555,945 A→G V302A (GTC→GCC)  rtcA ← RNA 3'‑terminal phosphate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,555,9450AG97.7% 133.0 / NA 44V302A (GTC→GCC) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (26/17);  minor base C (0/1);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCG  >  NC_000913/3555798‑3556084
                                                                                                                                                   |                                                                                                                                           
cGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCg                                                                                                                                          >  1:740451/1‑151 (MQ=255)
cGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCg                                                                                                                                          >  1:525648/1‑151 (MQ=255)
    tAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTa                                                                                                                                        >  1:180415/1‑149 (MQ=255)
       cgcgCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTc                                                                                                                                   >  2:117357/1‑151 (MQ=255)
             aGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCcgc                                                                                                                             <  2:426346/151‑1 (MQ=255)
                    gCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGcgcccgcc                                                                                                                      >  1:183912/1‑151 (MQ=255)
                       tACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAgc                                                                                                                   >  1:705429/1‑151 (MQ=255)
                                    catTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGc                                                                                                       >  1:748552/1‑150 (MQ=255)
                                        cAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCa                                                                                                   <  1:686261/150‑1 (MQ=255)
                                               cACCCGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGt                                                                                             <  1:403230/149‑1 (MQ=255)
                                                 cccGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGc                                                                                         <  2:5770/151‑1 (MQ=255)
                                                         aCGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATAtt                                                                                 >  2:693644/1‑151 (MQ=255)
                                                            ccATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTccc                                                                              <  1:466475/151‑1 (MQ=255)
                                                                 tGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATGTTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCAcc                                                                         >  2:62242/1‑151 (MQ=255)
                                                                  gTTTCTATCAAACTAAACCGCCCCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCAc                                                                          <  2:772615/149‑1 (MQ=255)
                                                                  gTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCAccg                                                                        >  1:632571/1‑151 (MQ=255)
                                                                      cTATCAAACTAAACCGCCCCGGCAAGAAACGCTCCACCCCCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGCCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCAccgcc                                                                      <  1:537283/149‑1 (MQ=255)
                                                                      cTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCAccgccgc                                                                    >  1:365170/1‑151 (MQ=255)
                                                                        aTCAAACTAACCCGCCCCGGCAAGAAACGCTCCCCCCCCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGGGGCCGTAAATTCCCCCGCCCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCAccgccg                                                                     <  2:740401/148‑1 (MQ=255)
                                                                        aTCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCAccgccgccg                                                                  >  2:490323/1‑151 (MQ=255)
                                                                         tCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCAccgccgccg                                                                  >  2:569782/1‑150 (MQ=255)
                                                                          cAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTg                                                                >  1:453301/1‑151 (MQ=255)
                                                                                aaaCCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGcc                                                          >  1:831731/1‑151 (MQ=255)
                                                                                    cGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGt                                                      >  2:813787/1‑151 (MQ=255)
                                                                                     gCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGt                                                      >  2:652166/1‑150 (MQ=255)
                                                                                      cACCGGCAAGAAACGCTCCCCCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGt                                                      <  2:26317/149‑1 (MQ=255)
                                                                                            cAAGAAACGCTCCCCCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGttt                                               <  2:106202/150‑1 (MQ=255)
                                                                                             aaGAAACGCTCCACCCCCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCa                                             <  2:576170/151‑1 (MQ=255)
                                                                                             aaGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTc                                              >  1:800466/1‑150 (MQ=255)
                                                                                                     cTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCttt                                      <  2:793886/150‑1 (MQ=255)
                                                                                                     cTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTc                                     >  2:745607/1‑151 (MQ=255)
                                                                                                      tCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTc                                     >  2:884900/1‑150 (MQ=255)
                                                                                                       ccaccaCCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACctct                                        <  1:756523/146‑1 (MQ=255)
                                                                                                        caccacCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCAcc                                  <  2:872628/151‑1 (MQ=255)
                                                                                                         accaccGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCAcc                                  >  1:183760/1‑150 (MQ=255)
                                                                                                           caccGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAAc                               >  1:450540/1‑151 (MQ=255)
                                                                                                            accGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACt                              >  1:808952/1‑151 (MQ=255)
                                                                                                               gcgATATTGGTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACtgtg                           <  2:417877/151‑1 (MQ=255)
                                                                                                                       ggTCAGCAGATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAacc                   >  1:448933/1‑151 (MQ=255)
                                                                                                                              agATGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCg             <  1:66571/150‑1 (MQ=255)
                                                                                                                                 tGGCATGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCAc         >  2:651467/1‑151 (MQ=255)
                                                                                                                                     aTGAGGGATGGGCGGCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAc     >  1:546849/1‑151 (MQ=255)
                                                                                                                                       gAGGGCTTGGCGGCCGAAAACTCCCCGGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCCAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTCTCACCAACTGTGCCGCAACCACCTCGGCTCTGAcgc   <  2:709313/151‑1 (MQ=255)
                                                                                                                                         gggATGGGCGGCCGTAAATCCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcgcg  <  2:237062/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
CGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAACTAAACCGCACCGGCAAGAAACGCTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCG  >  NC_000913/3555798‑3556084

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