Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,901,085 A→G 7.4% D165G (GAT→GGT)  yqcE → putative MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,901,0850AG 7.4% 166.0 / 3.7 54D165G (GAT→GGT) yqcEputative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/20);  new base G (2/2);  total (32/22)
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.88e-01

TGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAACGAAGAAAAACAGTCATTCCAGCTTAGCGACAT  >  NC_000913/2900940‑2901232
                                                                                                                                                 |                                                                                                                                                   
tgGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGa                                                                                                                                                     <  3:12703/146‑1 (MQ=255)
             ttAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCt                                                                                                                                      >  3:30842/1‑148 (MQ=255)
                  gCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGaaaa                                                                                                                                 >  5:62012/1‑148 (MQ=255)
                   cTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAAc                                                                                                                                >  8:42581/1‑148 (MQ=255)
                     gCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACgg                                                                                                                              >  6:301718/1‑148 (MQ=255)
                         cGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCAt                                                                                                                          >  8:4444/1‑148 (MQ=255)
                                tGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCTGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCTCTTTGCACCGGATGACAGCACCAGCATGAAAACGGTCATTatcatc                                                                                                                   >  7:301922/1‑148 (MQ=255)
                                tGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTatcatc                                                                                                                   >  1:261657/1‑148 (MQ=255)
                                   cGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCtat                                                                                                                <  8:40383/148‑1 (MQ=255)
                                     aTCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAg                                                                                                              >  6:150398/1‑148 (MQ=255)
                                          aGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGGTGACAGCACCAGCCTGaa                                                                                                                                   <  1:112741/122‑1 (MQ=255)
                                          aGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGGTGACAGCACCAGCCTGaa                                                                                                                                   >  2:112741/1‑122 (MQ=255)
                                          aGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGGTGACAGCACCAGCCTGaa                                                                                                                                   <  1:117271/122‑1 (MQ=255)
                                          aGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGGTGACAGCACCAGCCTGaa                                                                                                                                   >  2:117271/1‑122 (MQ=255)
                                                   gggAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCt                                                                                                >  3:132410/1‑148 (MQ=255)
                                                    ggAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCtt                                                                                               <  6:62012/148‑1 (MQ=255)
                                                       aaaTTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTTGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCAt                                                                                                                          <  5:256622/118‑1 (MQ=255)
                                                       aaaTTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCAt                                                                                                                          >  6:256622/1‑118 (MQ=255)
                                                          ttATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTggg                                                                                         <  3:274946/148‑1 (MQ=255)
                                                          ttATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTggg                                                                                         <  4:30842/148‑1 (MQ=255)
                                                                 cTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCtg                                                                                  >  2:302782/1‑148 (MQ=255)
                                                                    gatggaAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAg                                                                                                              <  5:179450/117‑1 (MQ=255)
                                                                    gatggaAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAg                                                                                                              >  6:179450/1‑117 (MQ=255)
                                                                    gatggaAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCtgtg                                                                                <  4:2856/147‑1 (MQ=255)
                                                                     atggaAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCt                                                                              >  6:23941/1‑148 (MQ=255)
                                                                     atggaAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCt                                                                              >  2:206443/1‑148 (MQ=255)
                                                                        gaAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGt                                                                           >  8:185281/1‑148 (MQ=255)
                                                                          aGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCAtt                                                                                                                         >  8:294743/1‑100 (MQ=255)
                                                                          aGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCAtt                                                                                                                         <  7:294743/100‑1 (MQ=255)
                                                                          aGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGttt                                                                         <  7:185281/148‑1 (MQ=255)
                                                                            ggCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGttttt                                                                       <  5:23941/148‑1 (MQ=255)
                                                                             gCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCa                                                                                                                                             >  5:23588/1‑77 (MQ=255)
                                                                             gCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCa                                                                                                                                             <  6:23588/77‑1 (MQ=255)
                                                                                   cggcgTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTc                                                                                                                            <  3:234738/88‑1 (MQ=255)
                                                                                   cggcgTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTc                                                                                                                            >  4:234738/1‑88 (MQ=255)
                                                                                    ggcgTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCg                                                               <  7:42581/148‑1 (MQ=255)
                                                                                       gTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGAt                                                             >  2:14392/1‑147 (MQ=255)
                                                                                       gTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATa                                                            <  1:14392/148‑1 (MQ=255)
                                                                                        tCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATa                                                            >  5:132955/1‑147 (MQ=255)
                                                                                         cGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATaa                                                           >  6:46917/1‑147 (MQ=255)
                                                                                                  gTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTg                                                 >  6:98197/1‑148 (MQ=255)
                                                                                                     gCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTAGATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTgcgc                                              <  5:150398/148‑1 (MQ=255)
                                                                                                            gtgtTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGcc                                        >  6:49383/1‑147 (MQ=255)
                                                                                                            gtgtTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCa                                       <  6:132955/148‑1 (MQ=255)
                                                                                                                ttACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTCAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATaa                                   >  7:275340/1‑148 (MQ=255)
                                                                                                                 tACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAAc                                  >  5:239556/1‑148 (MQ=255)
                                                                                                                  aCCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAAc                                  >  1:218330/1‑147 (MQ=255)
                                                                                                                          ggTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAACGAAGaaaa                          <  5:49383/147‑1 (MQ=255)
                                                                                                                             cTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAACGAAGAAAAACAg                      >  7:82870/1‑148 (MQ=255)
                                                                                                                               tttCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAACGAAGAAAAACAGt                     >  7:116768/1‑147 (MQ=255)
                                                                                                                                ttCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAACGAAGAAAAACAGTc                    <  7:4444/147‑1 (MQ=255)
                                                                                                                                ttCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAACGAAGAAAAACAGTc                    <  5:98197/147‑1 (MQ=255)
                                                                                                                                ttCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGCTTTACATCTTGTAGGGGATACTGTGCTGGTTTATTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAACGAAGAAAAACAGTCa                   >  7:32778/1‑148 (MQ=255)
                                                                                                                                     gCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAACGAAGAAAAACAGTCATTc                <  2:261657/146‑1 (MQ=255)
                                                                                                                                      cTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAACGAAGAAAAACAGTGCTTcccg             >  8:248238/1‑146 (MQ=255)
                                                                                                                                                 aTGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAACGAAGAAAAACAGTCATTCCAGCTTAGCGACAt  >  1:138642/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGATAACAACAACCTGCGCAGTGCCAATAACGAAGAAAAACAGTCATTCCAGCTTAGCGACAT  >  NC_000913/2900940‑2901232

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑