Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 344,674 T→C Y167H (TAC→CAC)  yahL → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913344,6740TC98.2% 190.1 / ‑4.6 55Y167H (TAC→CAC) yahLuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (28/26);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

ATATGACCAGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATT  >  NC_000913/344547‑344807
                                                                                                                               |                                                                                                                                     
atatGACCAGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCATATGTATTCTTCTTGAACACCATATGCCt                                                                                                                            <  2:1162060/139‑1 (MQ=255)
               tgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTc                                                                                                             >  2:548460/1‑139 (MQ=255)
               tgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTc                                                                                                             >  2:1688116/1‑139 (MQ=255)
                    ttCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTa                                                                                                         <  2:1693233/138‑1 (MQ=255)
                    ttCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTa                                                                                                         >  1:1693233/1‑138 (MQ=255)
                     tCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTaaa                                                                                                       >  2:289058/1‑139 (MQ=255)
                         gATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGa                                                                                                   <  2:1703033/139‑1 (MQ=255)
                           tatTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                 <  1:289058/139‑1 (MQ=255)
                                aaGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAgtg                                                                                            >  1:1244475/1‑139 (MQ=255)
                                     aaCGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGc                                                                                       >  1:673928/1‑139 (MQ=255)
                                     aaCGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGc                                                                                       <  1:481803/139‑1 (MQ=255)
                                                  gACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATg                                                                          >  1:815161/1‑139 (MQ=255)
                                                   aCGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATgg                                                                         >  2:1588381/1‑139 (MQ=255)
                                                        ttACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTa                                                                                                         <  1:1868195/102‑1 (MQ=255)
                                                        ttACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTa                                                                                                         >  2:1868195/1‑102 (MQ=255)
                                                        ttACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAACGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTcg                                                                    >  1:129591/1‑139 (MQ=255)
                                                              ggCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAAt                                                              <  2:1244475/139‑1 (MQ=255)
                                                                          tAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGa                                                                                                   >  1:1642422/1‑90 (MQ=255)
                                                                          tAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGa                                                                                                   <  2:1642422/90‑1 (MQ=255)
                                                                             tGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAAt                                                              <  2:1005033/124‑1 (MQ=255)
                                                                             tGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAAt                                                              >  1:1005033/1‑124 (MQ=255)
                                                                                  tttATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                 >  1:1010560/1‑84 (MQ=255)
                                                                                  tttATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                 <  2:1010560/84‑1 (MQ=255)
                                                                                  tttATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacag                                          >  1:544937/1‑139 (MQ=255)
                                                                                  tttATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacag                                          >  2:1097613/1‑139 (MQ=255)
                                                                                   ttATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                         <  2:388282/139‑1 (MQ=255)
                                                                                     aTCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcg                                       <  1:548460/139‑1 (MQ=255)
                                                                                       cGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCtt                                                                                                                           <  1:1435326/53‑1 (MQ=255)
                                                                                       cGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCtt                                                                                                                           >  2:1435326/1‑53 (MQ=255)
                                                                                       cGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTATAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                     <  1:1097613/139‑1 (MQ=255)
                                                                                               tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCata                                                                                                                                 >  2:696213/1‑39 (MQ=255)
                                                                                               tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCata                                                                                                                                 <  1:696213/39‑1 (MQ=255)
                                                                                               tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTa                                                                                                         <  2:502573/63‑1 (MQ=255)
                                                                                               tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTa                                                                                                         >  1:502573/1‑63 (MQ=255)
                                                                                               tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtgt                                                                                   <  2:185246/85‑1 (MQ=255)
                                                                                               tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtgt                                                                                   >  1:185246/1‑85 (MQ=255)
                                                                                               tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAAt                                                              <  2:1196458/106‑1 (MQ=255)
                                                                                               tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAAt                                                              >  1:1196458/1‑106 (MQ=255)
                                                                                                 cGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGc                           <  1:129219/139‑1 (MQ=255)
                                                                                                  gAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCa                          <  2:673928/139‑1 (MQ=255)
                                                                                                    aCAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAgg                        >  1:618424/1‑139 (MQ=255)
                                                                                                       gATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTcg                                                                    <  2:1350403/92‑1 (MQ=255)
                                                                                                       gATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGgcgc                     >  2:1053041/1‑139 (MQ=255)
                                                                                                       gATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGAGGAAAAATGGCTTcg                                                                    >  1:1350403/1‑92 (MQ=255)
                                                                                                         ttGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                 <  2:1640/61‑1 (MQ=255)
                                                                                                         ttGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                 >  1:1640/1‑61 (MQ=255)
                                                                                                              aaTGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa              <  2:544937/139‑1 (MQ=255)
                                                                                                              aaTGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa              <  1:1688116/139‑1 (MQ=255)
                                                                                                                 gTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa              <  2:1658180/136‑1 (MQ=255)
                                                                                                                 gTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa              >  1:1658180/1‑136 (MQ=255)
                                                                                                                  tATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa              >  1:156577/1‑135 (MQ=255)
                                                                                                                  tATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa              <  2:156577/135‑1 (MQ=255)
                                                                                                                          ttGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattattatt  >  2:29579/1‑139 (MQ=255)
                                                                                                                            gAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGgaga                                             >  1:136715/1‑94 (MQ=255)
                                                                                                                            gAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGgaga                                             <  2:136715/94‑1 (MQ=255)
                                                                                                                               |                                                                                                                                     
ATATGACCAGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATT  >  NC_000913/344547‑344807

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