Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,353,249 T→C H201H (CAT→CAC glpA → anaerobic sn‑glycerol‑3‑phosphate dehydrogenase, large FAD/NAD(P)‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,353,2490TC98.2% 185.6 / ‑4.2 54H201H (CAT→CACglpAanaerobic sn‑glycerol‑3‑phosphate dehydrogenase, large FAD/NAD(P)‑binding subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (33/20);  total (34/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

TACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCATGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTCCCGGCGAAAGGATCGCTG  >  NC_000913/2353153‑2353345
                                                                                                |                                                                                                
tACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCa                                                                                               >  2:1742060/1‑100 (MQ=255)
  ccGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCAcc                                                                                             <  2:1081211/100‑1 (MQ=255)
   cGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCt                                                                                            >  1:1229758/1‑100 (MQ=255)
    gCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTg                                                                                           >  2:410617/1‑100 (MQ=255)
    gCTCATGAAGTCACGGGGCCGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTg                                                                                           >  2:759035/1‑100 (MQ=255)
      tCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCATCGGCGAAACTCAGGCCCTTCACGCACCTGTc                                                                                         >  2:467803/1‑100 (MQ=255)
       cATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCg                                                                                        >  2:52650/1‑100 (MQ=255)
         tGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGACGAAACTCAGGCCCTTCACGCACCTGTCGTg                                                                                      >  2:1135771/1‑100 (MQ=255)
           aaGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGt                                                                                    <  1:433043/100‑1 (MQ=255)
            aGTCACGGAGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGtt                                                                                   >  2:1063623/1‑100 (MQ=255)
             gTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTa                                                                                  >  1:1387159/1‑100 (MQ=255)
                    ggCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGc                                                                           <  2:1797544/100‑1 (MQ=255)
                     gCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCt                                                                          >  1:731215/1‑100 (MQ=255)
                       tGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCATCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGTCGCTgg                                                                        <  1:467801/100‑1 (MQ=255)
                        gATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTggg                                                                       >  1:155831/1‑100 (MQ=255)
                          ttCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCTCTTCACTCACCTGTCGTGGTTAATGCCGCTgggg                                                                      <  1:1575328/99‑2 (MQ=255)
                          ttCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGAt                                                                     >  2:654837/1‑100 (MQ=255)
                          ttCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGAt                                                                     <  1:759031/100‑1 (MQ=255)
                                aaGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTgggg                                                               <  2:246290/100‑1 (MQ=255)
                                aaGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTggg                                                                <  2:169941/99‑1 (MQ=255)
                                 aGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTggggg                                                              <  1:1356150/100‑1 (MQ=255)
                                     gcgACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAac                                                          <  1:505818/100‑1 (MQ=255)
                                      ngACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAaca                                                         >  2:66991/2‑100 (MQ=255)
                                      ctaCGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAac                                                          <  2:628577/97‑1 (MQ=255)
                                      cgACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAaca                                                         >  1:863611/1‑100 (MQ=255)
                                      cgACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAaca                                                         >  1:863442/1‑100 (MQ=255)
                                        aCGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAAcaca                                                       <  2:1214150/100‑1 (MQ=255)
                                          ggTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACAtt                                                     <  1:410615/100‑1 (MQ=255)
                                             gtgCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCATGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGcc                                                  >  2:1003502/1‑100 (MQ=255)
                                               gCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGa                                                >  2:953770/1‑100 (MQ=255)
                                                 ggTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAat                                              <  1:560297/100‑1 (MQ=255)
                                                       cGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGAGGGCAACACATTGCCGAATATGCCg                                        >  2:509547/1‑100 (MQ=255)
                                                        gtgtACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATTCCGCTGGGATCTGGGGGCAACACATTGAAGAATATGCCGa                                       >  2:569864/1‑100 (MQ=255)
                                                               tAACCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGAAAAAAAATTCCCAAAAATGccaacctcgca                                >  2:1216260/1‑91 (MQ=255)
                                                                aaCCATCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCa                               >  2:831081/1‑100 (MQ=255)
                                                                    aTCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCg                           >  1:1158151/1‑100 (MQ=255)
                                                                    aTCTCACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCg                           <  1:252058/100‑1 (MQ=255)
                                                                      ctcACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCa                         <  1:415800/100‑1 (MQ=255)
                                                                        cACCGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATg                       >  1:900999/1‑100 (MQ=255)
                                                                          ccGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGtt                     >  2:1253330/1‑100 (MQ=255)
                                                                          ccGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGtt                     >  1:1604172/1‑100 (MQ=255)
                                                                           cGGCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTc                    >  2:1662784/1‑100 (MQ=255)
                                                                             gCGAAACTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGTGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTcc                   <  1:395367/99‑1 (MQ=255)
                                                                                aaaCTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTCCCGGc               <  2:1773369/100‑1 (MQ=255)
                                                                                 aaCTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTCCCGGCg              >  2:986585/1‑100 (MQ=255)
                                                                                 aaCTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTCCCGGCg              >  2:1159667/1‑100 (MQ=255)
                                                                                   cTCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTCCCGGCGaa            >  1:754782/1‑100 (MQ=255)
                                                                                    tCAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTCCCGGCGaaa           >  2:981136/1‑100 (MQ=255)
                                                                                     cAGGCCCTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTCCCGGCGaaa           <  2:450723/99‑1 (MQ=255)
                                                                                       ggCCCTTCACGCACCTGTCGTGGTTAATGACGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTCCGCATTCGCATGTTCCCGGCGAAAGGa        >  1:253600/1‑100 (MQ=255)
                                                                                          ccTTCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTCCCGGCGAAAGGATCg     >  2:824666/1‑100 (MQ=255)
                                                                                             tCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTCCCGGCGAAAGGATCgctg  >  2:1567578/1‑100 (MQ=255)
                                                                                             tCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTCCCGGCGAAAGGATCgctg  >  2:1294882/1‑100 (MQ=255)
                                                                                             tCACGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTCCCGGCGAAAGGATCgctg  <  1:96502/100‑1 (MQ=255)
                                                                                                |                                                                                                
TACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTGTACGTAACCATCTCACCGGCGAAACTCAGGCCCTTCATGCACCTGTCGTGGTTAATGCCGCTGGGATCTGGGGGCAACACATTGCCGAATATGCCGATCTGCGCATTCGCATGTTCCCGGCGAAAGGATCGCTG  >  NC_000913/2353153‑2353345

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