Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,424,516 C→T *92* (TAG→TAA yjfY ← YhcN family protein, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,424,5160CT97.2% 111.4 / NA 36*92* (TAG→TAAyjfYYhcN family protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/16);  minor base G (1/0);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTC  >  NC_000913/4424391‑4424644
                                                                                                                             |                                                                                                                                
cgcgCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAg                                                                                                                               >  2:469701/1‑129 (MQ=255)
cgcgCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAg                                                                                                                               <  1:469701/129‑1 (MQ=255)
cgcgCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGGTAAatat                                                                                                                     >  1:159065/1‑139 (MQ=255)
cgcgCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGGTAAatat                                                                                                                     <  2:159065/139‑1 (MQ=255)
       cGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCt                                                                                                              >  2:1646021/1‑139 (MQ=255)
                               aaCTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGACAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTcccc                                                                                     >  1:1489938/1‑139 (MQ=255)
                                    gTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTc                                                                                >  1:767498/1‑139 (MQ=255)
                                     tCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATTGCAAATTTTCGGTGGTTATACCTGGTGGTAGCGATAAATATCCGCGCTGGCGTGATATTTCCGTTAATTCCCGGTTcc                                                                               >  1:604969/1‑139 (MQ=255)
                                     tCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTcc                                                                               <  1:1646021/139‑1 (MQ=255)
                                            ccGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCgcgc                                                                                                               >  2:45415/1‑101 (MQ=255)
                                            ccGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCgcgc                                                                                                               <  1:45415/101‑1 (MQ=255)
                                              gCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCcgcg                                                                                                                <  1:61905/98‑1 (MQ=255)
                                              gCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCcgcg                                                                                                                >  2:61905/1‑98 (MQ=255)
                                              gCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGa                                                                      >  2:17739/1‑139 (MQ=255)
                                               cAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGt                                                                                  >  1:58621/1‑126 (MQ=255)
                                               cAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGt                                                                                  <  2:58621/126‑1 (MQ=255)
                                                  ccGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTa                                                                  <  1:918027/139‑1 (MQ=255)
                                                         tCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGgta                                                           <  2:340186/139‑1 (MQ=255)
                                                                                  ggCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGtt                                  >  1:626503/1‑139 (MQ=255)
                                                                                   gcaaagcaaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCcgcg                                                                                                                <  2:1710721/61‑1 (MQ=255)
                                                                                   gcaaagcaaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCcgcg                                                                                                                >  1:1710721/1‑61 (MQ=255)
                                                                                      aagcaaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttctt                              <  2:626503/139‑1 (MQ=255)
                                                                                                tGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGc                    <  1:17739/139‑1 (MQ=255)
                                                                                                 ggCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCcgcg                                                                                                                <  1:740962/47‑1 (MQ=255)
                                                                                                 ggCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCcgcg                                                                                                                >  2:740962/1‑47 (MQ=255)
                                                                                                 ggCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGcc                   >  2:1705303/1‑139 (MQ=255)
                                                                                                  gCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCt                  >  2:1416396/1‑139 (MQ=255)
                                                                                                   cAATGTTTCGTTGGTTATACGTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTg                 >  1:1078805/1‑139 (MQ=255)
                                                                                                   cAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACt                                                                                         <  2:973151/67‑1 (MQ=255)
                                                                                                   cAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACt                                                                                         >  1:973151/1‑67 (MQ=255)
                                                                                                     aTGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc               <  2:604969/139‑1 (MQ=255)
                                                                                                     aTGTTTCGTTGGTTATACCTGCTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc               <  2:1078805/139‑1 (MQ=255)
                                                                                                                gTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                >  2:1696433/1‑79 (MQ=255)
                                                                                                                gTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                <  1:1696433/79‑1 (MQ=255)
                                                                                                                  tataCCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTc  >  2:529504/1‑139 (MQ=255)
                                                                                                                  tataCCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTc  >  1:542306/1‑139 (MQ=255)
                                                                                                                          gtgtTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTg                                            <  1:1049074/89‑1 (MQ=255)
                                                                                                                             |                                                                                                                                
CGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTC  >  NC_000913/4424391‑4424644

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: