Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,518,075 A→G S138P (TCC→CCC)  aroB ← 3‑dehydroquinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,518,0750AG97.9% 159.9 / NA 47S138P (TCC→CCC) aroB3‑dehydroquinate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (25/21);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACGCTAACTCACGCGGGGGAAGCGTTTTCAGACAGTCGAGATCCACCACCACTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCAT  >  NC_000913/3517946‑3518209
                                                                                                                                 |                                                                                                                                      
ccGACGCTAACTCACCCGGGGGAAGCGTTTTCAGACACTCGAGATCCACCACCACTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCTGCCAACGGGGGAATCGAcc                                                                                                                              <  1:328038/140‑1 (MQ=255)
                   gggAAGCGTTTTCAGACAGTCGAGATCCACCACCACTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGtcgtcg                                                                                                           >  2:48856/1‑140 (MQ=255)
                    ggAAGCGTTTTCAGACAGTCGAGATCCACCACCACTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCgg                                                                                                          <  2:51492/140‑1 (MQ=255)
                    gaaaGCGTTTTCAGACAGTCGAGATCCACCACCACTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCgg                                                                                                          <  2:191166/138‑1 (MQ=255)
                     gAAGCGTTTTCAGACAGTCGAGATCCACCACCACTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCggg                                                                                                         >  1:96281/1‑140 (MQ=255)
                       aGCGTTTTCAGACAGTCGAGATCCACCACCACTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGAc                                                                                                       <  2:124924/140‑1 (MQ=255)
                                      tCGAGATCCACCACCACTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGac                                                                                        >  2:826309/1‑140 (MQ=255)
                                      tCGAGATCCACCACCACTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGac                                                                                        >  2:1611/1‑140 (MQ=255)
                                      tCGAGATCCACCACCACTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCAGCGACAGGAACGTCATCGGGAGTTGAATGAAACGGac                                                                                        >  1:219757/1‑140 (MQ=255)
                                          gaTCCACCACCACTGAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACcgc                                                                                   <  1:422656/140‑1 (MQ=255)
                                           aTCCACCACCACTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACcgc                                                                                   <  2:35802/140‑1 (MQ=255)
                                                ccaccacTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGa                                                                              >  2:451419/1‑140 (MQ=255)
                                                    cacTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAAc                                                                          >  2:534837/1‑140 (MQ=255)
                                                        gAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTcgc                                                                      <  1:756737/140‑1 (MQ=255)
                                                        gAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTcgc                                                                      >  1:200073/1‑140 (MQ=255)
                                                         aaGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTcgcc                                                                     >  1:44160/1‑140 (MQ=255)
                                                           gCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTcgccgc                                                                   >  1:424753/1‑140 (MQ=255)
                                                                  ggTAGAACGCGCCAATCTTGTTCTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAg                                                            <  2:287773/140‑1 (MQ=255)
                                                                   gTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGc                                                           >  2:827569/1‑140 (MQ=255)
                                                                   gTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGc                                                           >  1:802196/1‑140 (MQ=255)
                                                                   gTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGc                                                           >  2:683591/1‑140 (MQ=255)
                                                                   gTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATAAAACGGAAACCGCGTTGATAAGTCGCCGCCGCGAAGc                                                           >  1:228296/1‑140 (MQ=255)
                                                                     aGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCg                                                         >  2:443346/1‑140 (MQ=255)
                                                                        aCGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTc                                                      <  1:443346/140‑1 (MQ=255)
                                                                         cgcgCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCa                                                     >  2:475361/1‑140 (MQ=255)
                                                                           cgcCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGa                                                   <  2:18626/140‑1 (MQ=255)
                                                                           cgcCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGa                                                   <  2:424753/140‑1 (MQ=255)
                                                                            gcCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGAt                                                  >  2:716613/1‑140 (MQ=255)
                                                                            gcCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGAt                                                  >  2:738398/1‑140 (MQ=255)
                                                                            gcCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGAt                                                  >  2:300881/1‑140 (MQ=255)
                                                                             ccAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGAGACCGCTCAGATAACTCGCCGCGACGAAGCCCGTCAGATc                                                 >  2:374343/1‑140 (MQ=255)
                                                                                    tGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACt                                          <  1:451419/140‑1 (MQ=255)
                                                                                              gAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTAcgccgccgc                                >  1:190627/1‑140 (MQ=255)
                                                                                              gAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTAcgccgccgc                                >  1:151153/1‑140 (MQ=255)
                                                                                               aGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACgccgccgcc                               >  1:487222/1‑140 (MQ=255)
                                                                                                  ggATGGTTAACCGCACATTTGCCGCCAACGGAGGAATCGTCCTGCGACCGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAg                            <  2:404498/140‑1 (MQ=255)
                                                                                                      ggTTGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGcc                        <  2:151153/140‑1 (MQ=255)
                                                                                                        ttGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGccac                      <  1:827569/140‑1 (MQ=255)
                                                                                                        ttGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGccac                      <  1:716613/140‑1 (MQ=255)
                                                                                                        ttGACCGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGccac                      <  1:475361/140‑1 (MQ=255)
                                                                                                             cGCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCg                 <  2:190627/140‑1 (MQ=255)
                                                                                                              gCAGTTTTGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCgt                <  1:300881/140‑1 (MQ=255)
                                                                                                              gCAGTTTTGCCGCCAACAGGGGCATCGACCTCCGATAGTCACGTCGTCTGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCgt                <  2:228296/140‑1 (MQ=255)
                                                                                                                  ttttGCCGCCAACGGGGGAATCGACGTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCgtagta            <  1:312822/140‑1 (MQ=255)
                                                                                                                  ttttGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGccacca                    >  1:725884/1‑132 (MQ=255)
                                                                                                                  ttttGCCGCCAACGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGccacca                    <  2:725884/132‑1 (MQ=255)
                                                                                                                            aaCGGGGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCAt  >  1:643040/1‑140 (MQ=255)
                                                                                                                                 |                                                                                                                                      
CCGACGCTAACTCACGCGGGGGAAGCGTTTTCAGACAGTCGAGATCCACCACCACTGAAGCAGGTTGGTAGAACGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCAT  >  NC_000913/3517946‑3518209

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