Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,782,913 A→C 37.8% V250G (GTG→GGG)  gpsA ← glycerol‑3‑phosphate dehydrogenase (NAD+)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,782,9130AC37.8% 21.1 / 7.2 28V250G (GTG→GGG) gpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand):  ref base A (9/8);  new base C (2/9);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCA  >  NC_000913/3782680‑3783150
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
ttAAAGTCATATCTGCCTCGCGGGCGTTTTTTCCGCAATTTAATAATTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGa                                                                                                                                                                                                                                           <  2:882799/238‑1 (MQ=255)
        aaTGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTCTACCCACCGGGTTGGCGGACACGTAAGGACCaaaaacccac                                                                                                                                                                                                                                    >  1:1078126/1‑229 (MQ=255)
                cttggcggtgTTTTTTCCAAAAAATAATTCTTGAAAAATTGCCTGGGTTATTGGAATTTAAACGCCGAAGCGATGCGCCAGTTCGCGTATTTTTTTCGTATTGGGGTAGCCTTCCCCCCCCTGACCAATCTTCACCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCa                                                                                                                                                                                                                          <  2:69072/230‑1 (MQ=255)
                ctcgcggctgTTTTTTCCGATTATAAATCTTTGATAATTTTCCTGGGTTATTGGAATTCAAACCCCGACGGGCTGCGCCAGTTCGCGGACTTTTTTCGTATTGCGGTAGCCTTCCCCCCCCTGCCCAACCTCCTCCTGCGCGCTTGGTACATCCATGCCCTGACCGAGCATCATGCCAAACCGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCa                                                                                                                                                                                                                          <  2:869322/230‑1 (MQ=255)
                                  caaaattaattCTTGGTAAATTTCCCCGGTTTTTGGCATTCCACCGCGGAGGCGATGCGCAGTTTCGCGCACTTCTTTATTATGGCGGTAGCCTGCCACCACCTGACAAATCTTCCCCTGAGGGTGTTGCACAGCCATGCCCTGACCGAGCATCATGCCAAACCGGCGGTTACGCGATTGGTTGCCGGTACAGGTAAGCCCCAGCTCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAg                                                                                                                                                                                                        <  2:212936/233‑1 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCCCGCCTGGTTTGCCGGTAAGGTAACACCCCACAGCCCCAACCCCCCCCAGCCCCACAAAGGGGGCAGGGCGGGAccccg                                                                                                                                                                                           >  1:1077605/1‑236 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACCGGGTGTCGGGACAGGTAAGACCCAAATCGCCAAACCCCCCCATGCCCCTAAAAGGTGCAGGGTCggggcaccg                                                                                                                                                                                           >  1:679730/1‑231 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCCACAGGGTGTCGGTACAGGTAAGCACCACAACGCCCAGCCCCCCCCCGCCCCAAAAAAGGGGAGGGGCGGCACCCaac                                                                                                                                                                                          >  2:510728/1‑237 (MQ=255)
                                                                           tccccgccagggCGATGCGCCGGTTCGCGGCCTTTTTGCGTTTTGCGGTACCTTTCCCCCCCCTGACCAATTCTCCCCTGCGCTCTTTGTACCTCCAGGCCCTGACCGAGCATCAAGCAAAACCGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAGGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                                <  2:646471/228‑1 (MQ=255)
                                                                                  cGAAGCGATGCCCCAGTTCGCGGACTTCTTTCGTATTGCGGTACCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCa                                                                                                                                                         <  1:510728/239‑1 (MQ=255)
                                                                                                                              ccTTCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGACCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGt                                                                                                                                                                                                     <  1:23524/151‑1 (MQ=255)
                                                                                                                                                        cTGCGCGCTTGGTACATCCATGCCCGGACCGAGCATCATGCCAAAACGGGGGTTACGCGACTGGTTGTCGGTACAGGAAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCag                                                                                   <  1:662158/239‑2 (MQ=255)
                                                                                                                                                             cgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTATGCGACTGGGTGTCGGTACCGGTAAGCACCAAAACGCCCAGCCCCCCCATGCCCCCAAAGAAGGGGGGGGCGGCCCCCACCCGCCCCCCCAGAAGCGCCACATCAGCCAGCCCCCCGGGCGTCATCGCCGCACGCGCCCCCGCCCCCAAACCCACACCGCCGGACACCCCCCCCCCACAGGCAaaaaac                                                                              >  2:345522/1‑237 (MQ=255)
                                                                                                                                                                  tGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGAt                                                                                                                                                 <  1:294320/167‑1 (MQ=255)
                                                                                                                                                                         cactGCCGTGACCGGGCAACAGGCCAAAACGGCGGTTCCGCGACTGGTTGTCGGTAACGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACcgcg                                                                  <  1:1072694/236‑1 (MQ=255)
                                                                                                                                                                              cccTGACCGAGCATCATGCCAAAACGGCGGTTTCGCGACTGGGTGTCGGGACAGGTAAGCACCCCATCGCCAAGCCCCGCCATTCCCACAAAGgggggcggg                                                                                                                                                                                                      >  1:854340/1‑97 (MQ=255)
                                                                                                                                                                               ccTGACCGAGCATCATGCCAAAACGGGGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAACCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCa                                                            <  1:246974/239‑1 (MQ=255)
                                                                                                                                                                                    ccgggatacatGCAAAAGGGGGGGTTACGCGACTGGTTGTCGGTACAGGTAGGCCCCAGTTCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg                                                       <  1:244848/231‑1 (MQ=255)
                                                                                                                                                                                      gggcttcagGCCAAAACGGGGGTTCCGGGACTGGTTGCCGGTAAAGTAAAGCCCCAGATCCCAAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCa                                                     <  1:345522/234‑1 (MQ=255)
                                                                                                                                                                                                                  gcaCTGGTTTGTGGGTACGGGTAGCCCCCAGACCCCAAGCCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTaa                           <  1:445159/237‑1 (MQ=255)
                                                                                                                                                                                                                                 aaCAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGaa                                           <  1:366370/207‑1 (MQ=255)
                                                                                                                                                                                                                                   cAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCCATAACGTTTTTCACCCCCCCGCCACGCTCCACGCCCATGCAATCCGGGTTGACTTAAACGCGaaaaacattgcc          >  2:717874/1‑231 (MQ=255)
                                                                                                                                                                                                                                   cAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACCCGGCAACTTTTGcc          >  1:316574/1‑239 (MQ=255)
                                                                                                                                                                                                                                      gTAAGCCCCGAACCGCAAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg         <  1:35910/237‑1 (MQ=255)
                                                                                                                                                                                                                                      gTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCGCTAAACCAGGCAAATTTTTCCGcc       >  1:20391/1‑238 (MQ=255)
                                                                                                                                                                                                                                       gaaGCACCGGATCCCCAGCCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg         <  1:1035229/235‑1 (MQ=255)
                                                                                                                                                                                                                                            cACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGCCATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATg                                             >  1:1055246/1‑196 (MQ=255)
                                                                                                                                                                                                                                            cACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTAACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGCTTGCTGTAAACGCGGAAACTTTTGCCGCAGTgca  >  1:1048000/1‑239 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
GTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCA  >  NC_000913/3782680‑3783150

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