Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,416,622 C→G 100% 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,6220CG100.0% 87.4 / NA 34Q626E (CAG→GAG) ptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (17/17);  total (17/17)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGCAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCC  >  NC_000913/2416489‑2416759
                                                                                                                                     |                                                                                                                                         
tATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTacgacg                                                                                                                                  >  1:409865/1‑143 (MQ=255)
       ccGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGCGTACGACGCTGCGGt                                                                                                                           >  1:127340/1‑143 (MQ=255)
                    tGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGCAAAAAACCACAATGCTGATGATCGACGGCCCGCTGGAGCACGACGCTGCGGTCATGGCCGACgct                                                                                                              >  1:405631/1‑141 (MQ=255)
                     gCTGTCCCACTCCAGCGGTACTGCTGGGGCAGGTAGCGACGTAGAAACAGTTCGCGACGCAACTCGTCGGGCGCAGGAAAGACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTg                                                                                                             <  2:264175/143‑1 (MQ=255)
                           cTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAGGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGCGTACGACGCTGCGGTAATGGCTGACGTTGCGaaag                                                                                                       >  1:524284/1‑142 (MQ=255)
                                  aCCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGCAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGCACGACGCTGCGGTAATGGCTGACGTTGCGAACTCCACAgc                                                                                                >  1:271377/1‑143 (MQ=255)
                                  aCCGGTACTTCGGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAgc                                                                                                <  1:480919/143‑1 (MQ=255)
                                            cTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAActct                                                                                      >  1:308524/1‑143 (MQ=255)
                                                      tAGCGACCGAGAGAGAGTCCGCGAAGCAACGCGTCTGGCGCAGGAAAAACGTCCTGACCTGGTGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAAGGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAg                                                                            <  1:15609/143‑1 (MQ=255)
                                                       agcaccGTAGAAAAAGCTCGCGATGCAACTCGTCTGGCTCAGGAAAAACGTCCTGACCTGATGATCGACAGCCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAgg                                                                           <  2:59166/138‑1 (MQ=255)
                                                       aGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGGTGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAgg                                                                           >  1:125576/1‑143 (MQ=255)
                                                       aGCGACCTATCAGAAGTTCGCGATGCCCCTCGGCTGGGGCAGGATAGACGTCCTGACCTGGTGATACACGGTCTGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAgg                                                                           <  2:59040/143‑1 (MQ=255)
                                                          gACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTGCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGACATCCAAAGCGCCGAACTCTCCGGTTGCAGGTcg                                                                        >  1:177654/1‑143 (MQ=255)
                                                               aGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTa                                                                   >  1:353943/1‑143 (MQ=255)
                                                                     ggTTCGCGAAGCAACTCGCCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCgtgt                                                             <  1:794313/142‑1 (MQ=255)
                                                                      gTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATCATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGtt                                                            >  1:389817/1‑143 (MQ=255)
                                                                       ttCGCGAAGCAACTCTTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTc                                                           <  1:152878/143‑1 (MQ=255)
                                                                       cccccgAAGCAACTCGCCTGGCGCACGCATATCGTCCGGCTCTTATGCTCGACGGTCCGCTGGAGTACGACCCTGCGGTAATCGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTc                                                           <  2:223289/139‑1 (MQ=255)
                                                                           cgAATCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCt                                                       <  1:305018/143‑1 (MQ=255)
                                                                            gAAGCAACTCGTCTGGCGCAGGAAAAACGGCCTGACCTGATGATCGACGGGCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCCCCGAACTCTCCGGTCGCAGGTCCCGCTACCGTGTTCATCtt                                                      >  1:855350/1‑143 (MQ=255)
                                                                               gCAACTGGTCTGGGGCAGGAAAAACGTCCTGACCTGATGACCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTccc                                                   <  1:328597/143‑1 (MQ=255)
                                                                                 agcgcGTCTGGCGCAGGAATCTCGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCgg                                                 <  2:659850/139‑1 (MQ=255)
                                                                                 aaCTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGCAATCCAAAGCGCCGAACTCTCCCGTTGCAGGTCGCGCTACCGTGTCCAGCTTcccac                                                 >  2:128660/1‑141 (MQ=255)
                                                                                           gcgcAGGAAAACCGTCCTGACCTGATGGTCGACGGTCGGCTGGAGTACGCCGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAacac                                       <  2:509268/143‑1 (MQ=255)
                                                                                             gcAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGCACTCTCCGGCTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGCACGCCg                                     >  1:302806/1‑143 (MQ=255)
                                                                                               aGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCCAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGt                                   <  2:868943/143‑1 (MQ=255)
                                                                                                ggAAAAACGTCCTGACCTGATGATCGACGGTCGGGTCGAGCACGACGCGGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTa                                  <  2:180594/143‑1 (MQ=255)
                                                                                                 gAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTaa                                 >  1:143959/1‑143 (MQ=255)
                                                                                                   aaaaCGTCCTGACCGGATGATCCACGGTCCGCTAGAGAGCGAGGTTGCGGGGATGGCTGACGTTGCGAAATCCAAACCGCCGCACTCTCCGGGTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAaca                               <  1:737715/143‑1 (MQ=255)
                                                                                                      aCGGCCTGACCGGATGATCGACGGTCCGCTGGAGTACGGCGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAAcacca                            <  2:385041/143‑1 (MQ=255)
                                                                                                                 tgatgaTCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGCACACCGGTAACACCACCTACAAAGCg                 >  1:207503/1‑143 (MQ=255)
                                                                                                                    gtatCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTAAAAAGCGATa              <  2:590593/141‑1 (MQ=255)
                                                                                                                      atCGACGGTCCGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACa            >  1:47406/1‑143 (MQ=255)
                                                                                                                                cGCTGGAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGcc  >  1:840475/1‑143 (MQ=255)
                                                                                                                                     |                                                                                                                                         
TATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTGACCTGATGATCGACGGTCCGCTGCAGTACGACGCTGCGGTAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCC  >  NC_000913/2416489‑2416759

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