Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,445,152 G→A R137C (CGC→TGC)  yfcA ← TauE/TSUP family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,445,1520GA98.0% 183.6 / NA 50R137C (CGC→TGC) yfcATauE/TSUP family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (22/27);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAAAGGTCATGGC  >  NC_000913/2445006‑2445299
                                                                                                                                                  |                                                                                                                                                   
gTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTgg                                                                                                                                                 >  1:438695/1‑151 (MQ=255)
gTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTgg                                                                                                                                                 >  2:589898/1‑151 (MQ=255)
             ttGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCttcttc                                                                                                                                      <  2:315331/149‑1 (MQ=255)
                   aGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTcccc                                                                                                                                   >  2:524811/1‑146 (MQ=255)
                     gTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCctt                                                                                                                            >  1:308320/1‑151 (MQ=255)
                        aaaTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTgg                                                                                                                          <  2:216109/150‑1 (MQ=255)
                          aTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCAt                                                                                                                       <  2:90861/151‑1 (MQ=255)
                          aTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCCAGCTTTAGCa                                                                                                                        >  2:670057/1‑150 (MQ=255)
                             cGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAg                                                                                                                    <  1:289636/151‑1 (MQ=255)
                               cacaGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCa                                                                                                                  <  1:494242/151‑1 (MQ=255)
                                     gTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGaaa                                                                                                             >  2:202962/1‑150 (MQ=255)
                                      tAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGACGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCCAGACTTTGCATCAGCAAGaaa                                                                                                             >  2:462323/1‑149 (MQ=255)
                                       aaCGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAAt                                                                                                            >  2:549020/1‑149 (MQ=255)
                                          gAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAgg                                                                                                       <  2:324578/151‑1 (MQ=255)
                                          gAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAgg                                                                                                       <  1:727116/151‑1 (MQ=255)
                                                 aGCCCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATAc                                                                                                <  1:3039/151‑1 (MQ=255)
                                                  gcgcGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATa                                                                                                 >  2:822220/1‑149 (MQ=255)
                                                     cgTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAAt                                                                                            <  1:160001/151‑1 (MQ=255)
                                                      gTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAAtt                                                                                           <  2:669372/151‑1 (MQ=255)
                                                      gTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAAtt                                                                                           <  1:542952/151‑1 (MQ=255)
                                                        aaaaCGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTAc                                                                                         <  2:777123/151‑1 (MQ=255)
                                                                 cGGCATCCGGGACAGAGAATCCATCATAAAATCCGAGGCAGCCACGGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGGTGACTTCTTCGCCTAGCTTTGGCATCAGCAAAAAATATAGGCCGATGCAAAATACCAGAATg                                                                                  <  1:351171/149‑1 (MQ=255)
                                                                             cAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCg                                                                    <  2:603155/151‑1 (MQ=255)
                                                                              aaaGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGcc                                                                     >  2:620305/1‑149 (MQ=255)
                                                                              aaaGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGt                                                                   >  2:273090/1‑151 (MQ=255)
                                                                                 gAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAg                                                                <  1:341274/151‑1 (MQ=255)
                                                                                  aaTCCATCATAAAACCCGACGCAGCCACCGGCTATTAACGCTAAGGGTAAACCGTACAACCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGa                                                               <  1:729637/151‑1 (MQ=255)
                                                                                    tCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGAc                                                              <  2:10152/150‑1 (MQ=255)
                                                                                     ccATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACAt                                                            <  2:133673/151‑1 (MQ=255)
                                                                                      catcatAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATc                                                           <  2:757866/151‑1 (MQ=255)
                                                                                            aaaaaCCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACAACCGCCTg                                                     >  1:738269/1‑151 (MQ=255)
                                                                                             aaaaCCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTg                                                     >  1:308991/1‑150 (MQ=255)
                                                                                               aaCCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAAc                                                  <  1:502740/151‑1 (MQ=255)
                                                                                                aCCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGGCTTTTTCACCTAGCATTGGCATCAGCAAGAAATAAAGGGCGCTAAAAATTACCAGAATTGGCAAAATCTGCCGTAAGCCATCTGCCTGAACg                                                 >  1:236478/1‑151 (MQ=255)
                                                                                                 ccGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCCTCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGt                                                >  1:196969/1‑151 (MQ=255)
                                                                                                 ccGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGt                                                >  2:553457/1‑151 (MQ=255)
                                                                                                  cGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGt                                                >  1:155287/1‑150 (MQ=255)
                                                                                                     cGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAACTTTGGCCTCAGCCAGAAATAGAGGCCGATACAAATTACCAGAATGGGCCAAATCTGCCGTAAGACATCCGCCTGCACGTAtt                                             >  2:359131/1‑150 (MQ=255)
                                                                                                           cACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACcag                                      <  1:331613/151‑1 (MQ=255)
                                                                                                              cGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACcagc                                     >  2:322768/1‑149 (MQ=255)
                                                                                                                     tAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGcgccg                            >  1:755374/1‑151 (MQ=255)
                                                                                                                     tAACGCGAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGcgcc                             >  1:347381/1‑150 (MQ=255)
                                                                                                                           gAAGGGTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCAtt                      <  2:854444/151‑1 (MQ=255)
                                                                                                                                gTAAACCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTTGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGcc                 >  2:488965/1‑151 (MQ=255)
                                                                                                                                   aaCCGTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGAc              <  1:856255/151‑1 (MQ=255)
                                                                                                                                       gTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAAAg          <  2:513586/151‑1 (MQ=255)
                                                                                                                                       gTACATCCGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAAAg          <  2:299924/151‑1 (MQ=255)
                                                                                                                                             ccGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTGCCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAGAGGTCATg    <  2:777014/151‑1 (MQ=255)
                                                                                                                                             ccGGCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAAAGGTCATg    <  2:516228/151‑1 (MQ=255)
                                                                                                                                               ggCACTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAAAGGTCATGGc  <  2:4395/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACAAATTACCAGAATGGGCAAAATCTGCCGTAAGACATCAGCCTGAACGTATTGCACCAGCAGTGCGCCGCTCATTGAGCCGACAAAGGTCATGGC  >  NC_000913/2445006‑2445299

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