Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,505,552 C→T 7.5% G177R (GGG→AGG)  ydcO ← BenE family inner membrane putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,505,5520CT 7.5% 163.4 / 3.4 53G177R (GGG→AGG) ydcOBenE family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base T (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

TCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCAT  >  NC_000913/1505414‑1505699
                                                                                                                                          |                                                                                                                                                   
tCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGATTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGAGGATCACGATCCCAATAatcat                                                                                                                                            <  3:18968/148‑1 (MQ=255)
 cGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatcatc                                                                                                                                           >  2:33433/1‑148 (MQ=255)
     aCCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATcgcg                                                                                                                                        <  3:238858/147‑1 (MQ=255)
       cAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGc                                                                                                                                      >  1:35982/1‑147 (MQ=255)
          aaaaaGGGGGAGTGCAACGCTCAGGCGGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGCACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAAtt                                                                                                                                  <  3:213814/148‑1 (MQ=255)
             aaGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACCTCAGTTGTGACAACGTCACCTTGCGCGATGACGAGCACGATCCCAATAATCATCGCGGCAATTAcc                                                                                                                               >  4:162285/1‑148 (MQ=255)
                       gaaaCGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGTCGCACGGGTTTAAAGACAACATCATTTGTGACAACGTCACCTTGCGCTATGACTATCCAGATCACAATAATCATAGAGGACATTACCGCATAgaccg                                                                                                                     >  2:193763/3‑144 (MQ=255)
                       gCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGcgcg                                                                                                                     >  6:215045/1‑148 (MQ=255)
                           tgatgaGGCTTTGAGCATACGAATAATCTGGGGTAATATATGTGGGGTGTTCGTGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACTATCACGATCCCAATAATCATCGCGGCAATTAGCGGATAGCtcgg                                                                                                                    <  8:275110/140‑1 (MQ=255)
                            gCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGgcgc                                                                                                                 >  5:96261/1‑147 (MQ=255)
                            gCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCa                                                                                                                <  1:33433/148‑1 (MQ=255)
                                aGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACg                                                                                                             <  6:77146/147‑1 (MQ=255)
                                 gtcTGTGAGCAAACGAAAAATCAGGGGTAATATAAGGGGGGAGAACGGGGTTAAAGACAACATCAGGTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGc                                                                                                           <  5:208300/146‑1 (MQ=255)
                                 ggCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAgcgc                                                                                                                      <  2:298644/137‑1 (MQ=255)
                                 ggCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAgcgc                                                                                                                      >  1:298644/1‑137 (MQ=255)
                                    tgtgAGCAAACGAAACATCAGGGGTTATATAAGTTGGGAGAACGGGTTTAAAGACATCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCtt                                                                                                        <  6:96261/148‑1 (MQ=255)
                                      tgAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCtt                                                                                                        <  2:35982/146‑1 (MQ=255)
                                      tgAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCtt                                                                                                        <  7:136599/146‑1 (MQ=255)
                                      tgAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTgg                                                                                                      <  8:4036/148‑1 (MQ=255)
                                       gAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGgcgc                                                                                                                 <  3:270920/136‑1 (MQ=255)
                                       gAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGgcgc                                                                                                                 >  4:270920/1‑136 (MQ=255)
                                          cAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCTCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgc                                                                                                  <  5:215045/148‑1 (MQ=255)
                                             aCGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgccag                                                                                               >  6:158153/1‑148 (MQ=255)
                                                aaaacTCAGGGGTAATATAAGTGGGGAGAACGGGGTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCTATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgccagccc                                                                                            >  8:317933/1‑147 (MQ=255)
                                                aaaaaTCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgcc                                                                                                 >  8:13940/1‑143 (MQ=255)
                                                aaaaaTCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgcc                                                                                                 <  7:13940/143‑1 (MQ=255)
                                                             aaTATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCAt                                                                                           <  8:79744/136‑1 (MQ=255)
                                                             aaTATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCAt                                                                                           >  7:79744/1‑136 (MQ=255)
                                                               tataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATa                                                                               <  5:164361/146‑1 (MQ=255)
                                                               tataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATa                                                                               >  6:164361/1‑146 (MQ=255)
                                                                ataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgccagcc                                                                                             <  4:269401/131‑1 (MQ=255)
                                                                ataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgccagcc                                                                                             >  3:269401/1‑131 (MQ=255)
                                                                ataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACtt                                                                            >  3:254573/1‑148 (MQ=255)
                                                                ataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACtt                                                                            >  2:257159/1‑148 (MQ=255)
                                                                            aaCGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGt                                                                 >  3:90946/1‑147 (MQ=255)
                                                                            aaCGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTa                                                                <  1:257159/148‑1 (MQ=255)
                                                                             aCGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTa                                                                >  5:53724/1‑147 (MQ=255)
                                                                                  tttAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGa                                                          <  4:90946/148‑1 (MQ=255)
                                                                                   ttAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCTAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGAc                                                         >  7:148673/1‑148 (MQ=255)
                                                                                   ttAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCTAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGAc                                                         >  7:153135/1‑148 (MQ=255)
                                                                                    tAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTaaa                                                              >  4:171327/1‑142 (MQ=255)
                                                                                    tAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTaaa                                                              <  3:171327/142‑1 (MQ=255)
                                                                                          caacaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAg                                                  <  5:158153/148‑1 (MQ=255)
                                                                                                           aCGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGt                                  >  1:180971/1‑147 (MQ=255)
                                                                                                              tCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAAc                              >  2:82199/1‑148 (MQ=255)
                                                                                                                      cgcgATGACGATCACGATCCTAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGt                      <  8:148673/148‑1 (MQ=255)
                                                                                                                      cgcgATGACGATCACGATCCTAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGt                      <  8:153135/148‑1 (MQ=255)
                                                                                                                         gATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGtaa                    >  8:187748/1‑147 (MQ=255)
                                                                                                                                atcacgatcCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATcc            >  6:231931/1‑148 (MQ=255)
                                                                                                                                atcacgatcCCAATAATAATCGCGGCAATTACCGCATAGCGCGTCGCAACGGCCTTTGTGGCCAGCCATACCAGCACCATACTTCCACCCAAAGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATcc            <  1:193763/148‑1 (MQ=255)
                                                                                                                                atcaagataCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTCGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGAAATAAAATcc            >  8:321173/1‑148 (MQ=255)
                                                                                                                                 tcacgatcCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATccc           <  1:82199/148‑1 (MQ=255)
                                                                                                                                     gatcCCAAAAATCAGCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAACCGTAATAAAATCCCGGCa       <  2:24857/148‑1 (MQ=255)
                                                                                                                                          cAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCAt  >  4:289885/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
TCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCAT  >  NC_000913/1505414‑1505699

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