Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,900,327 A→C 24.2% intergenic (‑54/‑265) ygcW ← / → yqcE putative SDR family oxidoreductase/putative MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,900,3270AC24.2% 62.1 / 23.8 33intergenic (‑54/‑265)ygcW/yqcEputative SDR family oxidoreductase/putative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base C (4/4);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

AACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCAC  >  NC_000913/2900186‑2900469
                                                                                                                                             |                                                                                                                                              
aaCCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGaa                                                                                                                                         <  3:76471/149‑1 (MQ=255)
   cGCTCTTCCGATCTGTAAAAATTGCCGGTTTTCCTTTCAAGGAGAAAAAAACCATTGAAAACGCATTGAGAGATTCGATTGACATAGTTAGCCCCAATTTCCCTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAGtt                                                                                                                                      <  4:306841/149‑1 (MQ=255)
   cGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGtt                                                                                                                                      <  3:20935/149‑1 (MQ=255)
                 gTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGtt                                                                                                                                      <  3:880803/135‑1 (MQ=255)
                 gTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGtt                                                                                                                                      >  4:880803/1‑135 (MQ=255)
                                     tttCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGtt                                                                                                                                      >  3:465156/1‑115 (MQ=255)
                                     tttCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGtt                                                                                                                                      <  4:465156/115‑1 (MQ=255)
                                          gggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAggg                                                                                                                                           >  4:353403/1‑105 (MQ=255)
                                          gggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAggg                                                                                                                                           <  3:353403/105‑1 (MQ=255)
                                              gaAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGACAAAAAAAAAGAAAGGTAACTTCCCTCTCGcgag                                                                                                                                             <  1:295416/99‑1 (MQ=255)
                                               aaaaaaTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCgaga                                                                                                                                              <  3:1193351/97‑1 (MQ=255)
                                               aaaaaaTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCgaga                                                                                                                                              >  4:1193351/1‑97 (MQ=255)
                                                    aTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTAc                                                                                                                                    >  4:1018688/1‑102 (MQ=255)
                                                    aTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTAc                                                                                                                                    >  4:794114/1‑102 (MQ=255)
                                                    aTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTAc                                                                                                                                    <  3:794114/102‑1 (MQ=255)
                                                    aTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTAc                                                                                                                                    <  3:1018688/102‑1 (MQ=255)
                                                                                               tCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGt                                                                 >  4:17361/1‑126 (MQ=255)
                                                                                               tCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGt                                                                 <  3:17361/126‑1 (MQ=255)
                                                                                                                             ggTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATc                                                                   <  2:732902/94‑1 (MQ=255)
                                                                                                                             ggTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATc                                                                   >  1:732902/1‑94 (MQ=255)
                                                                                                                             ggTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAgt                                                        <  2:371754/105‑1 (MQ=255)
                                                                                                                             ggTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAgt                                                        >  1:371754/1‑105 (MQ=255)
                                                                                                                              gTAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATc                                                                   <  3:560592/93‑1 (MQ=37)
                                                                                                                              gTAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATc                                                                   >  4:560592/1‑93 (MQ=37)
                                                                                                                                 aCTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGa                                                                     >  1:1216899/1‑88 (MQ=37)
                                                                                                                                 aCTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGa                                                    >  1:1292437/1‑105 (MQ=255)
                                                                                                                                 aCTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGa                                                    <  2:1292437/105‑1 (MQ=255)
                                                                                                                                  cTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTAct       >  1:965748/1‑149 (MQ=255)
                                                                                                                                   ttCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATcaca                                                          <  2:643287/97‑1 (MQ=38)
                                                                                                                                   ttCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAg                       <  1:59625/132‑1 (MQ=255)
                                                                                                                                   ttCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTCGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAg                       >  2:59625/1‑132 (MQ=255)
                                                                                                                                       cccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAAttt                              >  3:716817/1‑121 (MQ=255)
                                                                                                                                       cccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAAttt                              <  4:716817/121‑1 (MQ=255)
                                                                                                                                       cccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTcac  >  3:27116/1‑149 (MQ=255)
                                                                                                                                        ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTa                                            <  2:343527/106‑1 (MQ=255)
                                                                                                                                        ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGAAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTa                                            >  1:343527/1‑106 (MQ=255)
                                                                                                                                             |                                                                                                                                              
AACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCAC  >  NC_000913/2900186‑2900469

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