Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,022,462 G→A 10.8% V183V (GTC→GTT PP_4427 ← amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,022,4620GA10.8% 97.3 / 3.0 37V183V (GTC→GTTPP_4427amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (16/17);  new base A (2/2);  total (18/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.19e-01

AAGAGCTGTATTGCATAGTTCCAGTCGAAGAAATTCATCACAGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGATGTGTTTTTCAACAGCTCGATCAGCAGATTAGTTCCAGGCCCCATCGCTGCGAG  >  NC_002947/5022317‑5022599
                                                                                                                                                 |                                                                                                                                         
aaGAGCTGTATTGCATAGTTCCAGTCGAAGAAATTCATCACAGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGAcc                                                                                                                                        <  8:6870/149‑1 (MQ=255)
     cTGTATTGCATAGTTCCAGTCGAAGAAATTCATCACAGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCGGCGCTAGGGCGAAGACCtt                                                                                                                                      >  6:66982/1‑146 (MQ=255)
            gCATAGTTCCAGTCGAAGAAATTCATCACAGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGtctc                                                                                                                             >  4:122102/1‑148 (MQ=255)
            gCATAGTTCCAGTCGAAGAAATTCATCACAGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGtctc                                                                                                                             >  3:171723/1‑148 (MQ=255)
            gCATAGTTCCAGTCGAAGAAATTCATCACAGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGtctc                                                                                                                             <  2:148260/148‑1 (MQ=255)
                          gaagaaATTCATCACAGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATc                                                                                                              >  1:105918/1‑149 (MQ=255)
                                aTTCATCACAGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTg                                                                                                        >  1:145219/1‑149 (MQ=255)
                                    atcaCAGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGt                                                                                                                                <  1:130459/121‑1 (MQ=255)
                                    atcaCAGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGt                                                                                                                                >  2:130459/1‑121 (MQ=255)
                                       acaGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGcc                                                                                                 >  6:19086/1‑149 (MQ=255)
                                        caGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCtt                                                                                                                                      <  2:725/111‑1 (MQ=255)
                                        caGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCtt                                                                                                                                      >  1:725/1‑111 (MQ=255)
                                           gCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGa                                                                                             <  3:122102/149‑1 (MQ=255)
                                               ccTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCt                                                                                         >  1:151338/1‑149 (MQ=255)
                                               ccTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCt                                                                                         >  1:245798/1‑149 (MQ=255)
                                                 ttGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTc                                                                                        >  8:30561/1‑148 (MQ=255)
                                                   ggACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCa                                                                                       <  2:145219/147‑1 (MQ=255)
                                                    gACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCTATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCt                                                                                                                                       <  8:14745/98‑1 (MQ=255)
                                                    gACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCAATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCt                                                                                                                                       >  7:14745/1‑98 (MQ=255)
                                                            cACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGc                                                                            <  2:31339/149‑1 (MQ=255)
                                                            cACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGc                                                                            <  4:171723/149‑1 (MQ=255)
                                                              cGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGt                                                                          <  5:107889/149‑1 (MQ=255)
                                                               gggCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCGGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGtt                                                                         <  5:66982/149‑1 (MQ=255)
                                                                       tGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGa                                                                 <  1:93309/149‑1 (MQ=255)
                                                                               tAAACGCCGCATTGCCAAGTTAATCGCTTGGTCCATTACGAAGTACATCAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAg                                                           <  4:42454/147‑1 (MQ=255)
                                                                                 aaCGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGa                                                        <  2:245798/148‑1 (MQ=255)
                                                                                          ttGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGATGTGTTTTTc                                              <  7:191312/149‑1 (MQ=255)
                                                                                             ccAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAAACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGATGTGtt                                                  <  6:57332/142‑1 (MQ=255)
                                                                                             ccAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAAACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGATGTGtt                                                  <  6:53528/142‑1 (MQ=255)
                                                                                             ccAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAAACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGATGTGtt                                                  >  5:57332/1‑142 (MQ=255)
                                                                                             ccAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAAACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGATGTGtt                                                  >  5:53528/1‑142 (MQ=255)
                                                                                                      tCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGATGTGTTTTTCAACAGCTCGATc                                  <  7:30561/149‑1 (MQ=255)
                                                                                                          ttGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGATGTGTTTTTCAACAGCTCGATcag                                >  3:175994/1‑147 (MQ=255)
                                                                                                                                  cagcGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGATGTGTTTTTCAACAGCTCGATCAGCAGATTAGTTCCAGGCCCCATCGCTg      >  5:3814/1‑149 (MQ=255)
                                                                                                                                     cgcTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGATGTGTTTTTCAACAGCTCGATcag                                <  4:5188/120‑1 (MQ=255)
                                                                                                                                     cgcTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGATGTGTTTTTCAACAGCTCGATcag                                >  3:5188/1‑120 (MQ=255)
                                                                                                                                      gcTAGGGCGAAGACCTTGAGAGACTCGAATGTAGCCTGATCGAGTTCGCGCGCCCTGAAGCTCAAATCTGCAAGCGTAATAAGGGAAACAAGGGATGTGTTTTTCAACAGCTCGATCAGGCGATTAGTTGCAGGCCCCATCGCAGCGAg  >  6:74964/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
AAGAGCTGTATTGCATAGTTCCAGTCGAAGAAATTCATCACAGGCTCCCTTGGACACGACCACGGGCCAGATGTCGCTCTAAACGCCGCATTGCCAAGTTAATCGCTTGGGCCATTACGAAGTACATGAGCAGCGCTAGGGCGAAGACCTTGAGAGTCTCGAATGTAGCCTGATCGAGTTGGCGCGCCCTGAAGCTCAAATCTGCAAGCGTTATAAGGGAAACAAGGGATGTGTTTTTCAACAGCTCGATCAGCAGATTAGTTCCAGGCCCCATCGCTGCGAG  >  NC_002947/5022317‑5022599

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