Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,416,622 C→G Q626E (CAG→GAG)  pta → phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,416,6220CG97.1% 84.4 / NA 35Q626E (CAG→GAG) ptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (16/18);  minor base T (0/1);  total (16/19)
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.98e-01

TCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGCAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCG  >  NC_000913/2416491‑2416752
                                                                                                                                   |                                                                                                                                   
tCGAACCACGCGTTGCTTTTCTTGCCTACTCCTCCTTTACTTCTTTTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCt                                                                                                                          <  1:363344/143‑1 (MQ=255)
  gAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGc                                                                                                                        >  1:703571/1‑143 (MQ=255)
             tGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGACAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGCAATGGCTg                                                                                                            >  1:827960/1‑143 (MQ=255)
              gCTATGTTCTCCTACTCCACGAGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGa                                                                                                           <  2:140433/143‑1 (MQ=255)
              gCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTAGCGAAGCAACTCATCTGGCGCAGGAAACACGTCCTGACCTGATGATAGACCGTCCGAGGGAGTACGACGCTGCGGCAGTGGCCGa                                                                                                           >  2:741324/1‑143 (MQ=255)
                  tGCTCTCCTTCTCCTCCGTTATTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGtt                                                                                                       <  1:172028/143‑1 (MQ=255)
                   gCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTg                                                                                                      >  2:168053/1‑143 (MQ=255)
                       tcCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGCCGGTTTGTTGGACCCCGACGCATCGGTAATGGCTGACGTTGCGaa                                                                                                  >  1:782734/1‑143 (MQ=255)
                       tcCTACTCCACCGGTACTTCTGGTGCAGGCAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCACGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTCATGGCTGACGTCGCGaa                                                                                                  >  2:493409/1‑143 (MQ=255)
                          tACTCCACCTGTAATTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAATCAACTCTTCTGGCTCAGGAAAAACGTCCTGACCTGATCATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATc                                                                                               <  1:239830/143‑1 (MQ=255)
                              ccACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGCAGCAACTCGTCTGGCGCCGGAACACCGTCCTGCCCTCATGATCGACGGTCCGCTGGAGTACGAAAAAAAGGAAATAGCAGACGTTACGAAATACaaa                                                                                           >  1:5452/1‑143 (MQ=255)
                               cACCTGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTGCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAg                                                                                          <  2:501081/143‑1 (MQ=255)
                                aCCGGTACTTCTGGTGCTGTTAGCGACGTATATATATTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAgc                                                                                         <  1:150524/143‑1 (MQ=255)
                                     tACTTCTGGTGCAGGTTTCTACGTATAATATTTTCTTGAATCAACTCGTCTGGCGCAGGAAAGACGTCCGGACCTGAGGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGa                                                                                    <  1:166705/143‑1 (MQ=255)
                                             gTGCAGGTAGCGACGGAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCg                                                                            <  2:120625/143‑1 (MQ=255)
                                                 aGGTAGCGACGTAGAAAAAGTTCGCTAAGCAACCCGACTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGCCGTTGCGCAATCCAAAACGCCGAACGCCCCGGtac                                                                        >  2:494216/1‑141 (MQ=255)
                                                   gTAGCGTCGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCa                                                                      <  1:508996/143‑1 (MQ=255)
                                                   gTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAGAACGTCCTGACCTGATGATCGACGGTCCGATGGAGTAAGACGATGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCa                                                                      >  1:189201/1‑143 (MQ=255)
                                                         aCGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAGATCCAAAGCGCCGAACTCTCCGGTTGCAGGTcgc                                                                >  1:140079/1‑143 (MQ=255)
                                                          cGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGACAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTATGGTGGTGAGCACAACCCCTGTGTTGCGAAATCCCAAGCCCCGAACTCCCCGCTTGCAGGCcgca                                                               >  1:368178/1‑142 (MQ=255)
                                                                aaaaGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCg                                                         >  2:365061/1‑143 (MQ=255)
                                                                         cgAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCACAGCGCCGAACTCTCCGGCTGCAGGTCCTGCTACCGTGTTCGTCt                                                >  2:300313/1‑143 (MQ=255)
                                                                           aaGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCAGCTGGAGTACGACCCAACGGTAATGGCTGACGTTGCGAAATCCAAACCGCCGAACTCTCCGGTTGCAACTCGCACTACCGAGTTCatcatc                                              >  2:402523/1‑143 (MQ=255)
                                                                            aGCACCTCGTCTTGCGCAGGAAAATCGTCCTGACCTGATGTTCGATGGTCCGCTGGAGTACTACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTcc                                             <  1:218237/143‑1 (MQ=255)
                                                                            aGCAACTCTTCTGGTGCAGGAAAAACGTCATGACCTGATGATCGACGGTCTGCTGGAGTACGACGCTGCGGGAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTcc                                             <  1:537466/143‑1 (MQ=255)
                                                                                aCTCGTCTGGCGCTGTAAAAACGTCCTGACCTGATGATCGACGGTCTGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGa                                         <  1:772543/143‑1 (MQ=255)
                                                                                 cTCGTCTGGTGAAGGAAAAACGTCCTGACCTGATGATCGATGTTCCGCTGTAGTACTACGCTGCGGTAATTGCTGACCTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAt                                        <  1:441845/143‑1 (MQ=255)
                                                                                  tCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGTTGGAGTATGACGTTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATc                                       <  1:748205/143‑1 (MQ=255)
                                                                                   cGTCTGGCGCATGTATAATGTCTTGTTCTGATGAGCGACGGTCAGGTGGAGTACGACGCGGCGGTAGTGGCTGACGTTGCGAAATCCAAAGCGGCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCt                                      <  1:304861/143‑1 (MQ=255)
                                                                                           gccggAAAACCGTCCTGAGCTGATGAGCGGCGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGGACTCTCCGGTTGCAGTTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCg                              <  2:769604/140‑1 (MQ=255)
                                                                                                aaaaaCGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAac                         <  1:480804/143‑1 (MQ=255)
                                                                                                aaaaaCGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAac                         <  2:703571/143‑1 (MQ=255)
                                                                                                      tttctgccctgaTGACCGACGGTCCGCTGGAGTAGGTCGCGGCGGTCATGGCTGACGTTGCGAAATCCAGAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACaccacc                   <  2:218016/140‑1 (MQ=255)
                                                                                                         ctgacctgaTGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAAGCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCTTGTTCATCTTCCCCGATCTGAACACCGGTAACCCCCCCtcc                >  2:565495/1‑141 (MQ=255)
                                                                                                                       gACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCAGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAAAACCGGTAACACCACCTACAAAGCGGTACAGCg  >  1:187146/1‑143 (MQ=255)
                                                                                                                                   |                                                                                                                                   
TCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGCAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCG  >  NC_000913/2416491‑2416752

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