Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,087,500 Δ1 bp 10.3% intergenic (+313/+605) PP_1865 → / ← aroH transposase/3‑deoxy‑7‑phosphoheptulonate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,087,5000T.10.3% 114.0 / 7.7 39intergenic (+313/+605)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base T (18/17);  new base . (2/2);  total (20/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 = 5.57e-01

ATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACG  >  NC_002947/2087363‑2087644
                                                                                                                                          |                                                                                                                                                
aTGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGtt                                                                                                                                       >  6:190044/1‑149 (MQ=255)
    ccGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGc                                                                                                                                   >  8:201238/1‑149 (MQ=255)
    ccGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGc                                                                                                                                   <  2:203481/149‑1 (MQ=255)
       tCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCg                                                                                                                                <  7:69410/149‑1 (MQ=255)
                     aaTGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAg                                                                                                                  >  2:114991/1‑149 (MQ=255)
                      aTGCCTCCACTCCGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAgg                                                                                                                 <  1:114991/149‑1 (MQ=255)
                         ccTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGat                                                                                                               >  7:8877/1‑148 (MQ=255)
                          cTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTc                                                                                                                             <  2:197256/133‑1 (MQ=255)
                          cTCCACTCGGCCTCCTGAAGTCGCAATCTGCATCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTc                                                                                                                             >  1:197256/1‑133 (MQ=255)
                              aCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGcc                                                                                                         <  6:255258/149‑1 (MQ=255)
                                              tCGCAATCTGCGCCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGa                                                                                                                                           >  4:165132/1‑99 (MQ=255)
                                              tCGCAATCTGCGCCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGa                                                                                                                                           <  3:165132/99‑1 (MQ=255)
                                                    tCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTgg                                                                                   >  1:212676/1‑149 (MQ=255)
                                                        cGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCTGCATGGCTGCCCTCGGCCTGgt                                                                                  <  7:201238/146‑1 (MQ=255)
                                                                 gACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGgaag                                                                       >  3:102301/1‑148 (MQ=255)
                                                                    tACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACg                                                                   >  1:150178/1‑149 (MQ=255)
                                                                      ccGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGc                                                                  <  2:212676/148‑1 (MQ=255)
                                                                      ccGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGAAATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGcc                                                                 >  7:221364/1‑149 (MQ=255)
                                                                        gcgcgcTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGt                                                               <  3:104514/149‑1 (MQ=255)
                                                                             cTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGttgctc                                                          <  2:150178/149‑1 (MQ=255)
                                                                                aGAACCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGttgctctt                                                        <  8:152524/148‑1 (MQ=255)
                                                                                  aaGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGttgctcttgct                                                     >  5:268596/1‑149 (MQ=255)
                                                                                    gCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGttgctcttgctct                                                   >  2:150056/1‑149 (MQ=255)
                                                                                    gCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCAGttgctcttgctct                                                   >  1:313714/1‑149 (MQ=255)
                                                                                      aaGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGttgctcttgctcttg                                                 <  5:190044/149‑1 (MQ=255)
                                                                                          tCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGTATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGttgctcttgctcttgctc                                              <  8:8877/148‑1 (MQ=255)
                                                                                            aaaaGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAAc                                                                    <  2:276546/124‑1 (MQ=255)
                                                                                                     aaGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCtt                                  <  1:150056/149‑1 (MQ=255)
                                                                                                              cGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGAAATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCt                          <  8:221364/148‑1 (MQ=255)
                                                                                                               ggTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCt                          >  6:119461/1‑147 (MQ=255)
                                                                                                                 tGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAgc                      >  1:74888/1‑149 (MQ=255)
                                                                                                                   aTCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAgcgc                    >  3:233251/1‑149 (MQ=255)
                                                                                                                        ccggcACCGCCG‑ATGCATTCGGCGAGGTTCAGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCgg                >  1:192603/1‑148 (MQ=255)
                                                                                                                            cctcgccgaATGCA‑TCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGt                                                               <  2:298817/94‑1 (MQ=18)
                                                                                                                            cctcgccgaATGCA‑TCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGt                                                               >  1:298817/4‑97 (MQ=18)
                                                                                                                            cctcgccgaATGCA‑TCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCtt                                        <  2:3905/117‑1 (MQ=25)
                                                                                                                            cctcgccgaATGCA‑TCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCtt                                        >  1:3905/4‑120 (MQ=25)
                                                                                                                                   g‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGa    <  2:74888/149‑1 (MQ=255)
                                                                                                                                         aTTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACg  >  7:251818/1‑146 (MQ=255)
                                                                                                                                          ttCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACg  >  4:189427/1‑145 (MQ=255)
                                                                                                                                          |                                                                                                                                                
ATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCG‑ATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACG  >  NC_002947/2087363‑2087644

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