Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,556,033 A→G Y273H (TAC→CAC)  rtcA ← RNA 3'‑terminal phosphate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,0330AG98.4% 216.5 / ‑5.5 63Y273H (TAC→CAC) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (29/33);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACC  >  NC_000913/3555899‑3556153
                                                                                                                                       |                                                                                                                        
cTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGt                                                                                                                      >  2:557056/1‑139 (MQ=255)
cTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGt                                                                                                                      >  2:716378/1‑139 (MQ=255)
cTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCACCGTGCTTGCCAGGTGGCGt                                                                                                                      >  2:314778/1‑139 (MQ=255)
     ccaccGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCAc                                                                                                                 <  1:557056/139‑1 (MQ=255)
     ccaccGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCAc                                                                                                                 <  2:1613197/139‑1 (MQ=255)
           cgATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCttt                                                                                                           >  1:509439/1‑139 (MQ=255)
           cgATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCttt                                                                                                           >  2:129779/1‑139 (MQ=255)
           cgATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCttt                                                                                                           >  2:551895/1‑139 (MQ=255)
                      agcagATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACtgtg                                                                                                <  2:509439/139‑1 (MQ=255)
                             ggCATGAGTGATGGGCGACCGTATATTCCCCCGCGCCCGCCAGCGCCAGCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                         <  1:1014437/139‑1 (MQ=255)
                             ggCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                         <  2:611650/139‑1 (MQ=255)
                             ggCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                         <  1:551895/139‑1 (MQ=255)
                                         gggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                         <  1:1328888/127‑1 (MQ=255)
                                         gggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                         >  2:1328888/1‑127 (MQ=255)
                                          ggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTg                                                                            <  1:1618883/139‑1 (MQ=255)
                                                cGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCt                                                                     >  1:579847/1‑139 (MQ=255)
                                                       tCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCAcc                                                              <  2:579847/139‑1 (MQ=255)
                                                         ccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAAc                                                                                                    >  1:1294818/1‑100 (MQ=255)
                                                         ccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAAc                                                                                                    <  2:1294818/100‑1 (MQ=255)
                                                          cccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               >  2:1620665/1‑88 (MQ=255)
                                                          cccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               <  1:1620665/88‑1 (MQ=255)
                                                          cccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACcga                                                            <  1:716378/139‑1 (MQ=255)
                                                                  cgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGcc                                                                                              >  1:1176528/1‑97 (MQ=255)
                                                                  cgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGcc                                                                                              <  2:1176528/97‑1 (MQ=255)
                                                                      aGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAg                                                 <  2:165204/139‑1 (MQ=255)
                                                                         gcCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                             <  2:695570/139‑1 (MQ=255)
                                                                          ccATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTc                                            <  1:830199/139‑1 (MQ=255)
                                                                                         gTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                         >  2:1329320/1‑79 (MQ=255)
                                                                                         gTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                         <  1:1329320/79‑1 (MQ=255)
                                                                                          ttGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               <  1:1226812/56‑1 (MQ=255)
                                                                                          ttGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               >  2:1226812/1‑56 (MQ=255)
                                                                                          ttGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                      >  1:1267862/1‑81 (MQ=255)
                                                                                          ttGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                      <  2:1267862/81‑1 (MQ=255)
                                                                                                aGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACtgtg                                                                                                >  1:1302670/1‑65 (MQ=255)
                                                                                                aGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACtgtg                                                                                                <  2:1302670/65‑1 (MQ=255)
                                                                                                  cAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                      <  2:476947/73‑1 (MQ=255)
                                                                                                  cAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                      >  1:476947/1‑73 (MQ=255)
                                                                                                        atTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCAc                                                                                                                 >  2:844967/1‑40 (MQ=255)
                                                                                                        atTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCAc                                                                                                                 <  1:844967/40‑1 (MQ=255)
                                                                                                          tCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgaga            >  2:189427/1‑139 (MQ=255)
                                                                                                             ccccACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACg          >  1:501560/1‑139 (MQ=255)
                                                                                                             ccccACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACg          >  1:1057249/1‑139 (MQ=255)
                                                                                                                 aCCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCgg                                                                                 <  1:1382290/64‑1 (MQ=255)
                                                                                                                 aCCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCgg                                                                                 >  2:1382290/1‑64 (MQ=255)
                                                                                                                  ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                      <  2:1004565/58‑1 (MQ=255)
                                                                                                                  ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCt                                                                                    <  2:795819/60‑1 (MQ=255)
                                                                                                                  ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCt                                                                                    >  1:795819/1‑60 (MQ=255)
                                                                                                                  ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTAtt     >  1:1536091/1‑139 (MQ=255)
                                                                                                                  ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTAtt     >  2:1549866/1‑139 (MQ=255)
                                                                                                                  cagccgccgTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                      >  1:1004565/3‑58 (MQ=255)
                                                                                                                     ccgccgTGCTTGCCAGGTGGCGTTTCTCCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  <  2:442260/139‑1 (MQ=255)
                                                                                                                     ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                             >  2:681786/1‑96 (MQ=255)
                                                                                                                     ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                             <  1:681786/96‑1 (MQ=255)
                                                                                                                     ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  <  1:82732/139‑1 (MQ=255)
                                                                                                                     ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  <  1:1549866/139‑1 (MQ=255)
                                                                                                                     ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  <  2:1536091/139‑1 (MQ=255)
                                                                                                                     ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  >  2:82732/1‑139 (MQ=255)
                                                                                                                     ccgccgTGCTTGCCAGGTGGCGTTACACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  >  2:898400/1‑139 (MQ=255)
                                                                                                                     ccaccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  <  1:314778/136‑1 (MQ=255)
                                                                                                                      cgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACcga                                                            >  1:1112747/1‑80 (MQ=255)
                                                                                                                      cgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGAAACCACCTCGGCACTGACGCGCTTTTCACcga                                                            <  2:1112747/80‑1 (MQ=255)
                                                                                                                       gccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCAcc                                                                                                       <  2:1619470/36‑1 (MQ=255)
                                                                                                                       gccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCAcc                                                                                                       >  1:1619470/1‑36 (MQ=255)
                                                                                                                                       |                                                                                                                        
CTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACC  >  NC_000913/3555899‑3556153

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