Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,851,227 T→A 14.7% intergenic (+34/+34) ilvD → / ← PP_5129 dihydroxy‑acid dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,851,2270TA14.7% 81.5 / 7.3 34intergenic (+34/+34)ilvD/PP_5129dihydroxy‑acid dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/15);  new base A (3/2);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

AGGTCCGTCCACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGAGCAGCGGCACCTACCCCTCCT  >  NC_002947/5851087‑5851375
                                                                                                                                            |                                                                                                                                                     
aGGTCCGTCCACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCgg                                                                                                                                              <  1:504775/150‑1 (MQ=255)
         cACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCgg                                                                                                                                              >  2:43518/1‑141 (MQ=255)
         cACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCgg                                                                                                                                              <  1:43518/141‑1 (MQ=255)
            gCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGAtt                                                                                                                                 <  2:587552/150‑1 (MQ=255)
            gCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGAtt                                                                                                                                 <  1:62249/150‑1 (MQ=255)
                      caccGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCAGAGCGTAACAAGGCAATACTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGAc                                                                                                                       >  2:30490/1‑150 (MQ=255)
                         cgcgCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGt                                                                                                                    <  1:296404/150‑1 (MQ=255)
                             cTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGt                                                                                                                >  2:353978/1‑150 (MQ=255)
                              tCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTg                                                                                                               >  1:75021/1‑150 (MQ=255)
                                aaGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACAGGCCCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGGGaa                                                                                                             >  1:67332/1‑150 (MQ=255)
                                 aGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAg                                                                                                            >  2:144636/1‑150 (MQ=255)
                                     cTATGCCCTGCTGGCCACCAGTGCCGATAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTCTTTGATTGTTCATTGCCGGTCAGTGAAGCGAg                                                                                                        <  2:885/150‑1 (MQ=255)
                                              gctgGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTggg                                                                                               <  1:76477/150‑1 (MQ=255)
                                                tgGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCGCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTc                                                                                             >  2:164283/1‑150 (MQ=255)
                                                     accaGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGAt                                                                                         <  2:477564/150‑1 (MQ=255)
                                                     accaGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGACTGTGATGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTGTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATc                                                                                        >  1:322841/1‑150 (MQ=255)
                                                                cAAGGGTGCGGTGTGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGAc                                                                             >  2:3621/1‑150 (MQ=255)
                                                                cAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGa                                                                                          <  1:302792/137‑1 (MQ=255)
                                                                cAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGa                                                                                          >  2:302792/1‑137 (MQ=255)
                                                                          gtgcgTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGga                                                                   <  1:144636/150‑1 (MQ=255)
                                                                           tgcgTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGgag                                                                  >  1:80538/1‑150 (MQ=255)
                                                                                aaCAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATACGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCa                                                                                                                  <  1:143216/97‑1 (MQ=255)
                                                                                aaCAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCa                                                                                                                  >  2:143216/1‑97 (MQ=255)
                                                                                        aaTGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGt                                                     <  2:75021/150‑1 (MQ=255)
                                                                                            cTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCt                                                 >  1:459895/1‑150 (MQ=255)
                                                                                              ggAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCg                                               <  2:80538/150‑1 (MQ=255)
                                                                                                     ctctGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGa                                        >  2:258933/1‑150 (MQ=255)
                                                                                                                       gTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGag                      <  1:258933/150‑1 (MQ=255)
                                                                                                                             aaaaaaCCGGCGCCCAGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGAGCAGCgg                >  1:185211/5‑150 (MQ=255)
                                                                                                                                 aaCCGGCGCCCAGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGga                                                                   >  2:208906/1‑95 (MQ=255)
                                                                                                                                 aaCCGGCGCCCAGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGga                                                                   <  1:208906/95‑1 (MQ=255)
                                                                                                                                 aaCCGGCGCCCAGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTc                                                >  2:21115/1‑114 (MQ=255)
                                                                                                                                 aaCCGGCGCCCAGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTc                                                <  1:21115/114‑1 (MQ=255)
                                                                                                                                           cTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGAGCAGCGGCACCTACCcctcct  <  1:164283/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
AGGTCCGTCCACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGAGCAGCGGCACCTACCCCTCCT  >  NC_002947/5851087‑5851375

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