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

GTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATAT‑AATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCA  >  NC_000913/3782680‑3783150
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
ttAAAGTCATATCTGCCTCGCGGGCGTTTTTTCCGCAATTT‑AATAATTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGa                                                                                                                                                                                                                                           <  2:882799/238‑1 (MQ=255)
        aaTGCTGCCTCGCGCGCGTTTTTTCCGCAATAT‑AATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTCTACCCACCGGG‑TTGGCGG‑ACACGTAA‑GGACCaaaaacccac                                                                                                                                                                                                                                    >  1:1078126/1‑229 (MQ=255)
                cttggcggtgTTTTTTCCAAAAAAT‑AATTCTTGAAAAATTGCCTGGGTTATTGGAATTTAAACGCCGAAGCGATGCGCCAGTTCGCGTATTTTTTTCGTATTGGGGTAGCCTTCCCCCCCCTGACCAATCTTCACCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCa                                                                                                                                                                                                                          <  2:69072/230‑1 (MQ=255)
                ctcgcggctgTTTTTTCCG‑ATTATAAATCTTTGATAATTTTCCTGGGTTATTGGAATTCAAACCCCGACGGGCTGCGCCAGTTCGCGGACTTTTTTCGTATTGCGGTAGCCTTCCCCCCCCTGCCCAACCTCCTCCTGCGCGCTTGGTACATCCATGCCCTGACCGAGCATCATGCCAAACCGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCCCCAGATCGCCAAG‑CCCCGCCa                                                                                                                                                                                                                          <  2:869322/230‑1 (MQ=255)
                                  caaaatt‑aattCTTGGTAAATTTCCCCGGTTTTTGGCATTCCACCGCGGAGGCGATGCGCAGTTTCGCGCACTTCTTTATTATGGCGGTAGCCTGCCACCACCTGACAAATCTTCCCCTGAGGGTGTTGCACAGCCATGCCCTGACCGAGCATCATGCCAAACCGGCGGT‑TACGCGATTGG‑TTGCCGGTACAGGTAA‑GCCCCAGCTCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAg                                                                                                                                                                                                        <  2:212936/233‑1 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TTCCCGCCTGGTTTGCCGGTA‑AGGTAACACCCCACAGCCCCAAC‑CCCCCCCAGCCCCACAAAGGGGGCAGGG‑CGGGAccccg                                                                                                                                                                                           >  1:1077605/1‑236 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCGACCGG‑GTGTCGGGACAGGTAA‑GACCCAAATCGCCAAA‑CCCCCCCATGCCCCTAAAAGGTGCAGGGTCggggcaccg                                                                                                                                                                                           >  1:679730/1‑231 (MQ=255)
                                                 gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCCACAGG‑GTGTCGGTACAGGTAA‑GCACCACAACGCCCAG‑CCCCCCCCCGCCCCAAAAAAGGGGAGGGGCGGCACCCaac                                                                                                                                                                                          >  2:510728/1‑237 (MQ=255)
                                                                           tccccgccagggCGATGCGCCGGTTCGCGGCCTTTTTGCGTTTTGCGGTACCTTTCCCCCCCCTGACCAATTCTCCCCTGCGCTCTTTGTACCTCCAGGCCCTGACCGAGCATCAAGCAAAACCGGCGGT‑TCCGCGACTGG‑TTGTCGGTACAGGTAA‑GCCCCAGATCGCCAGG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                                <  2:646471/228‑1 (MQ=255)
                                                                                  cGAAGCGATGCCCCAGTTCGCGGACTTCTTTCGTATTGCGGTACCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCa                                                                                                                                                         <  1:510728/239‑1 (MQ=255)
                                                                                                                              ccTTCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGACCATCATGCCAAAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGt                                                                                                                                                                                                     <  1:23524/151‑1 (MQ=255)
                                                                                                                                                        cTGCGCGCTTGGTACATCCATGCCCGGACCGAGCATCATGCCAAAACGGGGGT‑TACGCGACTGG‑TTGTCGGTACAGGAAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCag                                                                                   <  1:662158/239‑2 (MQ=255)
                                                                                                                                                             cgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TATGCGACTGG‑GTGTCGGTACCGGTAA‑GCACCAAAACGCCCAG‑CCCCCCCATGCCCCCAAAGAAGGGGGGGGCGGCCCCCACCCGCCCCCCCAGAAGCGCCACATCAGCCAGCCCCCCGGGCGTCATCGCCGCACGCGCCCCCGCCCCCAAACCCACACCGCCGGACACCCCCCCCCCACAGGCAaaaaac                                                                              >  2:345522/1‑237 (MQ=255)
                                                                                                                                                                  tGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TCCGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGAt                                                                                                                                                 <  1:294320/167‑1 (MQ=255)
                                                                                                                                                                         cactGCCGTGACCGGGCAACAGGCCAAAACGGCGGT‑TCCGCGACTGG‑TTGTCGGTAACGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACcgcg                                                                  <  1:1072694/236‑1 (MQ=255)
                                                                                                                                                                              cccTGACCGAGCATCATGCCAAAACGGCGGT‑TTCGCGACTGG‑GTGTCGGGACAGGTAA‑GCACCCCATCGCCAAG‑CCCCGCCATTCCCACAAAGgggggcggg                                                                                                                                                                                                      >  1:854340/1‑97 (MQ=255)
                                                                                                                                                                               ccTGACCGAGCATCATGCCAAAACGGGGGT‑TCCGCGACTGG‑TTGTCGGTACAGGTAA‑GCCCCAGATCGCCAAC‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCa                                                            <  1:246974/239‑1 (MQ=255)
                                                                                                                                                                                    ccgggatacatGCAAAAGGGGGGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAG‑GCCCCAGTTCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg                                                       <  1:244848/231‑1 (MQ=255)
                                                                                                                                                                                      gggcttcagGCCAAAACGGGGGT‑TCCGGGACTGG‑TTGCCGGTAAAGTAAA‑GCCCCAGATCCCAAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCa                                                     <  1:345522/234‑1 (MQ=255)
                                                                                                                                                                                                                  gcaCTGGTTTGTGGGTACGGGTAG‑CCCCCAGA‑CCCCAAGCCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTaa                           <  1:445159/237‑1 (MQ=255)
                                                                                                                                                                                                                                 aaCAGGTAA‑GCCCCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGaa                                           <  1:366370/207‑1 (MQ=255)
                                                                                                                                                                                                                                   cAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCCATAACGTTTTTCACCCCCCCGCCACGCTCCACGCCCATGCAATCCGGGTTGACTTAAACGCGaaaaacattgcc          >  2:717874/1‑231 (MQ=255)
                                                                                                                                                                                                                                   cAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACCCGGCAACTTTTGcc          >  1:316574/1‑239 (MQ=255)
                                                                                                                                                                                                                                      gTAA‑GCCCCGAACCGCAAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg         <  1:35910/237‑1 (MQ=255)
                                                                                                                                                                                                                                      gTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCGCTAAACCAGGCAAATTTTTCCGcc       >  1:20391/1‑238 (MQ=255)
                                                                                                                                                                                                                                       gaa‑GCACCGGATCCCCAGC‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg         <  1:1035229/235‑1 (MQ=255)
                                                                                                                                                                                                                                            cACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGCCATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATg                                             >  1:1055246/1‑196 (MQ=255)
                                                                                                                                                                                                                                            cACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTAACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGCTTGCTGTAAACGCGGAAACTTTTGCCGCAGTgca  >  1:1048000/1‑239 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
GTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATAT‑AATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCA  >  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: ‑