Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,854,423 T→C E144E (GAA→GAG)  ydjF ← putative DNA‑binding transcriptional regulator YdjF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,854,4230TC97.7% 143.9 / NA 43E144E (GAA→GAG) ydjFputative DNA‑binding transcriptional regulator YdjF
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (23/19);  total (24/19)
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

TCCAGCGCGCCGCTGTTAATATCAAGACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGTTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCACTACGGTCCTGTAACAATTTGAGCAATTCCATTACGGT  >  NC_000913/1854304‑1854545
                                                                                                                       |                                                                                                                          
tCCAGCGCGCCGCTGTTAATATCAAGACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCTTTGCAGTGAAAGAGCACCTTTGTTTGGCTCGTCACCCGtcg                                                                                                               >  2:247731/1‑131 (MQ=255)
tCCAGCGCGCCGCTGTTAATATCAAGACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTg                                                                                                               >  1:115335/1‑133 (MQ=255)
  cAGCGCGCCGCTGTTAATATCAAGACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGaa                                                                                                             >  1:127883/1‑133 (MQ=255)
   aGCGCGCCGCTGTTACTATCAAGATCTTCACACCTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGaaa                                                                                                            <  2:162081/133‑1 (MQ=255)
   aGCGCGCCGCTGTTAATATCAAGACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGaaa                                                                                                            <  1:27530/133‑1 (MQ=255)
     cgcgcCGCTGTTAATATCAAGACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTCTTTTTGTCTAGTTCGCCACCCATTGAAACg                                                                                                          >  1:162081/1‑133 (MQ=255)
     cgcgcCGCTGTTAATATCAAGACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACg                                                                                                          >  2:129285/1‑133 (MQ=255)
        gccgcTGTTAATATCAAGACCTTTACAACTCATTACCATGATATCGACATGATTGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACt                                                                                                       >  2:112102/1‑133 (MQ=255)
                   tatCAAGACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTctga                                                                                            <  2:509571/133‑1 (MQ=255)
                         gACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCTTTGAAACGACTTTAATTTCTGACTGAGc                                                                                      >  2:154676/1‑133 (MQ=255)
                         gACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGc                                                                                      >  2:84952/1‑133 (MQ=255)
                         gACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGc                                                                                      >  2:396473/1‑133 (MQ=255)
                         gACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGc                                                                                      >  1:440940/1‑133 (MQ=255)
                                   aCTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGa                                                                            >  1:389007/1‑133 (MQ=255)
                                   aCTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGa                                                                            >  1:333767/1‑133 (MQ=255)
                                    cTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATGTCTGACGGAGCCATAACATGaa                                                                           >  1:130604/1‑133 (MQ=255)
                                    cTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGGAATTCTTCCTTGCAGTGAAAGTGGATTTTTGTTTAGCTCCCCACCCGGTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGaa                                                                           >  1:34486/1‑133 (MQ=255)
                                     tCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAAt                                                                          <  2:127883/133‑1 (MQ=255)
                                      cATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACTTGAATc                                                                         >  2:340994/1‑133 (MQ=255)
                                       aTTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTATAATTTCCGACTGACCCAGAACATGAATCg                                                                        >  1:251738/1‑133 (MQ=255)
                                             aTGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTccg                                                                  <  2:270257/133‑1 (MQ=255)
                                                 tatCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGa                                                              >  1:454434/1‑133 (MQ=255)
                                                  atCGACATGATAGCGTCTGATGATCTCTCTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAg                                                             <  2:152438/133‑1 (MQ=255)
                                                  atCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAg                                                             <  1:84952/133‑1 (MQ=255)
                                                    cGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGtt                                                           <  2:14306/133‑1 (MQ=255)
                                                         tGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTGTGACTGAGACAGCACATGAATCGCTGCCGCTCAGTTGGTTa                                                      >  2:227838/1‑133 (MQ=255)
                                                          gATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAg                                                     <  2:78603/133‑1 (MQ=255)
                                                           aTAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGc                                                    <  2:454434/133‑1 (MQ=255)
                                                                   tgatgaTCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAg                                            <  1:129285/133‑1 (MQ=255)
                                                                    gatgatCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAgg                                           >  1:258084/1‑133 (MQ=255)
                                                                     atgatCTCTTTAGTACTTCTTCCTCGCACTGACAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGc                                          <  2:182391/133‑1 (MQ=255)
                                                                       gatCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCa                                        <  1:396473/133‑1 (MQ=255)
                                                                          ctctTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTCGGTAAGCAACGTCAGGCCACTa                                     >  1:182391/1‑133 (MQ=255)
                                                                            ctTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCACTACg                                   >  1:341542/1‑133 (MQ=255)
                                                                               tAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAAGCCACTACGGTc                                <  2:333767/133‑1 (MQ=255)
                                                                                     ttcttcCTTGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCACTACGGTCCTGTaa                          >  1:296723/1‑133 (MQ=255)
                                                                                       cttcCTTGCAGTGAAAGTGTATGTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCACTACGGTCCTGTAACa                        >  2:313953/1‑133 (MQ=255)
                                                                                             tGCAGTGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCACTACGGTCCTGTAACAATTTGa                  >  2:116526/1‑133 (MQ=255)
                                                                                                  tGAAAGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACTTGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCACTACGGTCCTGTAACAATTTGAGCAAt             <  1:340994/133‑1 (MQ=255)
                                                                                                    aaaGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCACTACGGTCCTGTAACAATTTGAGCAATTc           <  1:468730/133‑1 (MQ=255)
                                                                                                    aaaGTGTATTTTTGTTTAGCTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCACTACGGTCCTGTAACAATTTGAGCAATTc           <  2:341542/133‑1 (MQ=255)
                                                                                                            tttttGTTTAGCTCCCCACCCGTTTAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCACTACGGTCCTGTAACAATTTGAGCAATTCCATTACgg   <  2:173604/133‑1 (MQ=255)
                                                                                                            tctttgtttgACTCCACACCCGTTAATACGACTTTAATTTCCGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCACTACGGTCCTGTAACAATTTGAGCAATTCCATTACgg   <  2:251738/131‑1 (MQ=255)
                                                                                                             ttttGTTTAGCTCCCGACTCGTGGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCACTACGGTCCTGTAACAATTTGAGCAATTCCATTACGGt  <  1:247731/133‑1 (MQ=255)
                                                                                                                       |                                                                                                                          
TCCAGCGCGCCGCTGTTAATATCAAGACCTTTACAACTCATTACCATGATATCGACATGATAGCGCCTGATGATCTCTTTAGTAATTCTTCCTTGCAGTGAAAGTGTATTTTTGTTTAGTTCCCCACCCGTTGAAACGACTTTAATTTCTGACTGAGCCAGAACATGAATCGCTTCCGCCGAGTTGGTTAGCAACGTCAGGCCACTACGGTCCTGTAACAATTTGAGCAATTCCATTACGGT  >  NC_000913/1854304‑1854545

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