Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 628,569 A→G N7D (AAT→GAT)  entA → 2,3‑dihydro‑2,3‑dihydroxybenzoate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913628,5690AG95.0% 66.2 / ‑6.2 20N7D (AAT→GAT) entA2,3‑dihydro‑2,3‑dihydroxybenzoate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (5/14);  minor base T (0/1);  total (5/15)
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

AAGCGGCGCTGCGTGAGGTGATCCTGCCGTTGCTGGACGAGTCCGATGAACCGTTCGATGACGACAACCTGATCGACTACGGTCTGGATTCGGTGCGCATGATGGCGCTGGCGGCGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGA  >  NC_000913/628337‑628812
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
aaGCGGCGCTGCGTGAGGTGATCCTGACGTTGCTGGACGAGTCCGATGAACCGTTCGATGACGACAACCTGATCGACTACGGTCTGGATTCGGTGCGCATGATGGCGCTGGCGGCGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGgcg                                                                                                                                                                                                                                    <  1:191792/251‑1 (MQ=255)
                                   gACGAGTCCGAGGCCCCGTTCGATGAACACAACCTGATCGACTACGGTCGGGATTCGGTGCGCATGACGGCGCTGGCGGCGCGCTGGCGCAAAGGGCATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCAGCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGATAAAGGTATCGGCTACGCCAC‑GGCGCTGGCg                                                                                                                                                                                                <  2:203457/251‑1 (MQ=255)
                                          ccGATGAACCGTTCGATGACGACAACCTGATCGACTACGGTCTGGATTCGGTGCGCATGATGGCGCTGGCGGCGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGATGGCGTttgttg                                                                                                                                                                                         >  1:51702/1‑251 (MQ=255)
                                                         aTGACGACAACCTGATCGACTACGGTCTGGATTCGGTGCGCATGATGGCGCTGGCGGCGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAg                                                                                                                                                                          <  1:135948/251‑1 (MQ=255)
                                                         aTGACGACAACCTGATCGACTACGGTCTGGATTCGGAGCGCATGATGGCGCTCGCGACGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGACAAACCCGACCAAAGACGCCTGGTGGAAGCTACTATCCCGCGAGGAGAAATACTGGCTTTCAGCGGTAAGGATGTATGGGCAACCAGCGCAGGTAAAGGTATCGGCTCCGC‑ACAGGCGGTGGCGTTTTGTGAGGCGGGAggcaga                                                                                                                                                                           >  2:270085/1‑245 (MQ=255)
                                                          tGACGACAACCTGATCGACTACGGTCTGGATTCGGTGCGCATGATGGCGCTGGCGGCGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTttgtt                                                                                                                                                                                          >  1:86161/1‑234 (MQ=255)
                                                          tGACGACAACCTGATAGACTACGGTCTGGAGTCGGTGCGCATGATGGCGCTTGCGGCGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTttgtt                                                                                                                                                                                          <  2:86161/234‑1 (MQ=255)
                                                           gacgacAACCTGATCGACTACGGTCTGGATTCGGTGCGCATGATGGCGCTGGCGGCGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGACAAACCCGACCATCGCCGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGAAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGGGGCGGGcgcgcagg                                                                                                                                                                          >  2:210601/1‑245 (MQ=255)
                                                               acAACCTGATCGACTACGGTCTGGATTCGGTGCGCATGATGGCGCTGGCGGCGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGTCGGGAGCGAAAGTTACAg                                                                                                                                                                    <  1:321565/251‑1 (MQ=255)
                                                                  aCCTGATCGACTACGGTCTGGATTCGGTGCGCATGATGGCGCTGGCGGCGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGtt                                                                                                                                                                 <  1:13763/251‑1 (MQ=255)
                                                                              aCGGTCTGGATTCGGTGCGCATGATGGCGCTGGCGGCGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCg                                                                                                                                                      <  1:219801/250‑1 (MQ=255)
                                                                                          cGGTGCGCATGATGGCGCTGGCGGCGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTAGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGc                                                                                                                                         <  1:42292/251‑1 (MQ=255)
                                                                                                                                           tatCGACTTTGTCAAGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGGATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGATCAAGCGTTCACTCAGGGGCAATATCCCTTTGCGAACGAAGAGATGGATGTTGCCGCCGCTGGGCAggg                                                                                        >  1:275349/1‑250 (MQ=255)
                                                                                                                                                        aTGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAgtgtgt                                                                           <  1:348506/251‑1 (MQ=255)
                                                                                                                                                          gCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCa                                                                         <  1:122959/251‑1 (MQ=255)
                                                                                                                                                                    agcccGTCCATCGACGCCTGGTGGAAGCTTCTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTa                                                               <  1:51496/249‑1 (MQ=255)
                                                                                                                                                                             aTCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACg                                                      <  1:178829/251‑1 (MQ=255)
                                                                                                                                                                                          ggAAGCTACTCTCCCGCGAGGTGAAATAATGGCTTTCAGCGGTAAAGATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGAcg                                         <  1:246389/251‑1 (MQ=255)
                                                                                                                                                                                                    cgtccGCTAGGAGAAATAATGGATGTCAGCTGTAAAGATGACTGGGTAACCGGTGCGGGGAAATGTATCGGTTACGGCAC‑TGCGCTGGCGTTTGTTGAGTGGGGAGCGGAAGTTGCAGGGTTTGATCTAGCGTTCACTCAGGATCAATATCCTTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTTTTAGCTTAAACGGAGCGACTGGAtgcggttgtc                                 <  2:313002/246‑7 (MQ=255)
                                                                                                                                                                                                                                 gggTAAATAGGTCTGTGTAACCAGCGCAGGTAAAGGTATCGGCGTCGCCAC‑GGCGATGTCGTTTGTTGGGGCAGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCCCTCAGGAGCAATATCCCTTGGCGACCGAAGTTATGGATGTTGCCGACTCTGCGCAGGTCGCTCAAGTGGGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTTGTCAAGTCGTCTGGAATTTTACGCATGGGCGCGa  <  2:210525/250‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
AAGCGGCGCTGCGTGAGGTGATCCTGCCGTTGCTGGACGAGTCCGATGAACCGTTCGATGACGACAACCTGATCGACTACGGTCTGGATTCGGTGCGCATGATGGCGCTGGCGGCGCGCTGGCGCAAAGTGCATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCAC‑GGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGA  >  NC_000913/628337‑628812

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑