Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,518,153 C→T G112S (GGC→AGC)  aroB ← 3‑dehydroquinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,518,1530CT97.9% 157.4 / NA 47G112S (GGC→AGC) aroB3‑dehydroquinate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/22);  minor base G (1/0);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTTATACTGCTCGCCGTCAGGGAGGAT  >  NC_000913/3518019‑3518285
                                                                                                                                      |                                                                                                                                    
cgcgCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAAGGGAGACCACGCTGATAGCTCGGCGCCGCGAAGCTGGTCa                                                                                                                                 >  2:204267/1‑140 (MQ=255)
cgcgCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAAGGGACATCAGGCTGATACCTCGATGCCGCGAAGCTGGTCa                                                                                                                                 >  2:357211/1‑140 (MQ=255)
cgcgCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCa                                                                                                                                 >  1:127343/1‑140 (MQ=255)
cgcgCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCa                                                                                                                                 >  1:555631/1‑140 (MQ=255)
cgcgCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTGGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCa                                                                                                                                 >  2:327769/1‑140 (MQ=255)
                ttACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAAGGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTAcgcc                                                                                                                 <  1:357211/140‑1 (MQ=255)
                      aGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACgccgccgcc                                                                                                           >  2:734023/1‑140 (MQ=255)
                      aGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTAATAACTCGACGCCGCGAAGCTGGTCTGATCACCCACTACgccgccgcc                                                                                                           >  2:284709/1‑140 (MQ=255)
                       ggggATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCTACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCa                                                                                                          >  2:172104/1‑140 (MQ=255)
                       ggggATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCa                                                                                                          >  2:576572/1‑140 (MQ=255)
                        gggATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCaa                                                                                                         >  2:116891/1‑140 (MQ=255)
                               ttGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGccac                                                                                                  <  2:555631/140‑1 (MQ=255)
                                  aCCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAg                                                                                               >  2:626425/1‑140 (MQ=255)
                                  aCCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAg                                                                                               >  2:29506/1‑140 (MQ=255)
                                  aCCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCGGGTAAGAACGCCCACTACGCCGCCGCCAAGCGCCACCAg                                                                                               >  2:280775/1‑140 (MQ=255)
                                                    aCGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATg                                                                             <  1:29506/140‑1 (MQ=255)
                                                                  cTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCGGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTaaca                                                               >  1:546955/1‑140 (MQ=255)
                                                                    gCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTaacaac                                                             <  2:472602/140‑1 (MQ=255)
                                                                    gCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTaacaac                                                             <  1:734023/140‑1 (MQ=255)
                                                                    gCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTaacaac                                                             <  1:116891/140‑1 (MQ=255)
                                                                    gCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTaacaac                                                             >  1:472602/1‑140 (MQ=255)
                                                                     cGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCCCGCCCATGCTGTTTTTGTAACAACg                                                            >  1:243831/1‑140 (MQ=255)
                                                                       aCAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGcc                                                          >  2:548758/1‑140 (MQ=255)
                                                                       aCAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGcc                                                          <  2:418559/140‑1 (MQ=255)
                                                                              cgtcgtCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGa                                                   <  1:462075/140‑1 (MQ=255)
                                                                              cgtcgtCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGa                                                   <  1:633881/140‑1 (MQ=255)
                                                                                tcgtcgGGACTTGAATGCAACGAACACCGCGGTGACAACTCGCCGCCGAGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACg                                                 <  2:244309/140‑1 (MQ=255)
                                                                                   tcgGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTa                                              <  1:719175/140‑1 (MQ=255)
                                                                                   tcgGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTa                                              >  1:704702/1‑140 (MQ=255)
                                                                                       gACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGa                                          >  1:69374/1‑140 (MQ=255)
                                                                                          ttGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTa                                       >  2:73190/1‑140 (MQ=255)
                                                                                                gaaACGGACACCGCGTTGATAAATCACCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCa                                 <  1:204267/140‑1 (MQ=255)
                                                                                                gaaACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCa                                 <  2:69374/140‑1 (MQ=255)
                                                                                                 aaaCGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAgg                               >  1:138481/1‑140 (MQ=255)
                                                                                                  aaCGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAgg                               >  1:406558/1‑140 (MQ=255)
                                                                                                  aaCGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAgg                               <  2:127343/140‑1 (MQ=255)
                                                                                                  aaCGGACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAgg                               <  1:576572/140‑1 (MQ=255)
                                                                                                     ggACACCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCtt                            <  1:548758/140‑1 (MQ=255)
                                                                                                        cacCGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAGAGAGGGGATCGAGAACAGCCAGGCTTTta                         >  1:237626/1‑140 (MQ=255)
                                                                                                          ccGCGCTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTtata                       <  2:406558/140‑1 (MQ=255)
                                                                                                             cgcTGATAACTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTCTAACAACGACGTGAAGACGGTATCGCGTACAGCCAGGCTTTTATACTg                    >  2:236746/1‑140 (MQ=255)
                                                                                                                    aaCTCGCCGCCGCGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTTATACTGCTCGCCg             <  2:76445/140‑1 (MQ=255)
                                                                                                                      cTCGCCGCCGCGAAGCTGGTCAGATCGCCAACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTTATACTGCTCGCCGTc           <  1:327769/140‑1 (MQ=255)
                                                                                                                         gccgccgcGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTTATACTGCTCGCCGTCAgg        <  1:694463/140‑1 (MQ=255)
                                                                                                                            gccgcGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTTATACTGCTCGCCGTCAgg        >  2:56168/1‑137 (MQ=255)
                                                                                                                            gccgcGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTTATACTGCTCGCCGTCAgg        <  1:56168/137‑1 (MQ=255)
                                                                                                                            gccgcGAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTTATACTGCTCGCCGTCAGggag     >  2:54189/1‑140 (MQ=255)
                                                                                                                               gcgAAGCTGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTTATACTGCTCGCCGTCAGGGAGGAt  <  1:620485/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
CGCGCCAATCATGTTTTTACCGAGGGGATGGTTGACCGCAGTTTTGCCGCCAACGGAGGAATCGACCTGCGACAGTAACGTCGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTTATACTGCTCGCCGTCAGGGAGGAT  >  NC_000913/3518019‑3518285

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