Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,111,939 T→G 6.3% T92T (ACA→ACC)  yiiQ ← DUF1454 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,111,9390TG 6.3% 228.9 / 7.1 63T92T (ACA→ACC) yiiQDUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base G (2/2);  total (32/31)
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.66e-01

TTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGCAGGCTGGGATTCTGGCTGTTAAAGTCTTCTCGAA  >  NC_000913/4111794‑4112069
                                                                                                                                                 |                                                                                                                                  
ttttGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACAGCGCTGGAGCGCGGt                                                                                                                                  >  1:767868/1‑148 (MQ=255)
 tttGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTa                                                                                                                                 >  3:205167/1‑148 (MQ=255)
 tttGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACAGCGCTGGAGCGCGGTa                                                                                                                                 <  2:767868/148‑1 (MQ=255)
         ttAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAg                                                                                                                             <  1:264741/144‑1 (MQ=255)
         ttAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAg                                                                                                                             >  2:264741/1‑144 (MQ=255)
                              tCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGa                                                                                                                               >  4:448419/1‑121 (MQ=255)
                              tCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGa                                                                                                                               <  3:448419/121‑1 (MQ=255)
                                                ccTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGAtt                                                                                  >  4:356143/1‑148 (MQ=255)
                                                ccTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGAtt                                                                                  <  3:356143/148‑1 (MQ=255)
                                                       tttcgctttcgcGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTtctc                                                                                                                >  4:959899/1‑111 (MQ=255)
                                                       tttcgctttcgcGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGAtctc                                                                                                                <  3:959899/111‑1 (MQ=255)
                                                                  cgGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCa                                                                                                                           >  2:692313/1‑89 (MQ=255)
                                                                  cgGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCa                                                                                                                           <  1:692313/89‑1 (MQ=255)
                                                                  cgGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGActgct                                                                >  3:909049/1‑148 (MQ=255)
                                                                      tttttGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCg                                                            <  4:909049/148‑1 (MQ=255)
                                                                       ttttGCTCTGGCCGCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAAtttt                                                                                                       <  3:532500/104‑1 (MQ=255)
                                                                       ttttGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAAtttt                                                                                                       >  4:532500/1‑104 (MQ=255)
                                                                        tttGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGAt                                                          <  4:205167/148‑1 (MQ=255)
                                                                              ctGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTg                                                                                 <  3:869804/119‑1 (MQ=255)
                                                                              ctGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTg                                                                                 >  4:869804/1‑119 (MQ=255)
                                                                                    ccTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTtctc                                                                                                                >  3:463977/1‑82 (MQ=255)
                                                                                    ccTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTtctc                                                                                                                <  4:463977/82‑1 (MQ=255)
                                                                                      tGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTg                                                                                                                                   <  3:890277/61‑1 (MQ=255)
                                                                                      tGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTg                                                                                                                                   >  4:890277/1‑61 (MQ=255)
                                                                                         aTGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGtt                                                                                                                   <  1:1071000/74‑1 (MQ=255)
                                                                                         aTGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGtt                                                                                                                   >  2:1071000/1‑74 (MQ=255)
                                                                                              tAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGc                                    >  3:210245/1‑148 (MQ=255)
                                                                                                  cAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTgcagca                                                                                             <  4:1047156/87‑1 (MQ=255)
                                                                                                  cAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTgcagca                                                                                             >  3:1047156/1‑87 (MQ=255)
                                                                                                      tCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTGATTTTACTTGCAGCACGAGTGAGATTGGCtt                                                                             <  4:576027/99‑1 (MQ=255)
                                                                                                      tCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTGATTTTACTTGCAGCACGAGTGAGATTGGCtt                                                                             >  3:576027/1‑99 (MQ=255)
                                                                                                         tttGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCggg                                                                      >  4:491639/1‑103 (MQ=255)
                                                                                                         tttGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCggg                                                                      <  3:491639/103‑1 (MQ=255)
                                                                                                         tttGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTc                                              >  1:415151/1‑127 (MQ=255)
                                                                                                         tttGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTc                                              <  2:415151/127‑1 (MQ=255)
                                                                                                         tttGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGCAGGCTGGGAtt                         >  3:687591/1‑148 (MQ=255)
                                                                                                              aTGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTaa                                                                                                           <  3:222903/61‑1 (MQ=255)
                                                                                                              aTGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTaa                                                                                                           >  4:222903/1‑61 (MQ=255)
                                                                                                               tGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGActgc                                                                 >  3:1049025/1‑102 (MQ=255)
                                                                                                               tGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGActgc                                                                 <  4:1049025/102‑1 (MQ=255)
                                                                                                               tGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGCAGGCTGGGATTCTgg                     <  2:819052/146‑1 (MQ=255)
                                                                                                               tGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGCAGGCTGGGATTCTgg                     >  1:819052/1‑146 (MQ=255)
                                                                                                                   tttGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGc                                                                                                                            <  2:78758/39‑1 (MQ=255)
                                                                                                                   tttGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGc                                                                                                                            >  1:78758/1‑39 (MQ=255)
                                                                                                                   tttGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGCAGGCTGGGATTCTGGCTGTTa               >  3:426771/1‑148 (MQ=255)
                                                                                                                    ttGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCg                                                            <  3:227563/102‑1 (MQ=255)
                                                                                                                    ttGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCg                                                            >  4:227563/1‑102 (MQ=255)
                                                                                                                         ttttAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGaaa                                                 <  4:49372/108‑1 (MQ=255)
                                                                                                                         ttttAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGaaa                                                 >  3:49372/1‑108 (MQ=255)
                                                                                                                           ttAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGCAGGCTGGGATTCTGGCTGTTAAAGtcttc       <  4:687591/148‑1 (MQ=255)
                                                                                                                                gTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGCAGGCTGGGATTCTGGCTGTTAAAGTCTTCTCGaa  <  4:210245/148‑1 (MQ=255)
                                                                                                                                 tACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAg                                        <  4:264528/109‑1 (MQ=255)
                                                                                                                                 tACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAg                                        >  3:264528/1‑109 (MQ=255)
                                                                                                                                   ccGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCggg                                                                      <  4:837326/77‑1 (MQ=255)
                                                                                                                                   ccGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCggg                                                                      >  3:837326/1‑77 (MQ=255)
                                                                                                                                   ccGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGCAGGCTGGGATTCTGGCTGTTaa              <  4:444929/133‑1 (MQ=255)
                                                                                                                                   ccGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGCAGGCTGGGATTCTGGCTGTTaa              >  3:444929/1‑133 (MQ=255)
                                                                                                                                      cgcTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGtt                                           >  1:485591/1‑101 (MQ=255)
                                                                                                                                      cgcTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGtt                                           <  2:485591/101‑1 (MQ=255)
                                                                                                                                          ccAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGCAGGCTGGGATTCTGGCTGTTAAAGt           >  1:330734/1‑129 (MQ=255)
                                                                                                                                          ccAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGCAGGCTGGGATTCTGGCTGTTAAAGt           <  2:330734/129‑1 (MQ=255)
                                                                                                                                           gaGCGCGGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGt                                                              <  4:157967/76‑1 (MQ=255)
                                                                                                                                           gaGCGCGGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGt                                                              >  3:157967/2‑77 (MQ=255)
                                                                                                                                                 |                                                                                                                                  
TTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAGGTCATTTGAATGCTTTTGATTTTTAAGGTACCGCGCTCCAGCGCTGTAGAAGCATACAAGTTCTCATTAATTTTACTTGCAGCACGAGTGAGATTGGCTTTGTCGGGACTGCTGTCGATGGCACGAAATTCGTTCAGTGGCAGGCTGGGATTCTGGCTGTTAAAGTCTTCTCGAA  >  NC_000913/4111794‑4112069

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