Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,895,005 G→A 8.7% V18V (GTC→GTT)  yegH ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,895,0050GA 8.7% 130.8 / 2.6 46V18V (GTC→GTT) yegHputative transporter
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base A (2/2);  total (24/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 = 4.97e-01

GCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATA  >  NC_002947/3894871‑3895148
                                                                                                                                      |                                                                                                                                               
gCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCaa                                                                                                                                   >  8:63630/1‑149 (MQ=255)
          ccAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGc                                                                                                                         >  2:131238/1‑149 (MQ=255)
           cAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGcc                                                                                                                        <  7:293678/149‑1 (MQ=255)
                   tgatgaGCGCCTGGCTCAGGCCTATGACCCGCGCGCGTTCGCGCTGATGGTGCGGGAACTTGTCGGCAAGGATGGCGATGAACTCCAGGTTCTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTg                                                                                                                <  5:137385/149‑1 (MQ=255)
                      tgaGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTgggg                                                                                                             >  4:61815/1‑149 (MQ=255)
                      tgaGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTgggg                                                                                                             >  4:89197/1‑149 (MQ=255)
                             cAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAg                                                                                                                                >  3:143577/1‑123 (MQ=255)
                             cAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAg                                                                                                                                <  4:143577/123‑1 (MQ=255)
                             cAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGc                                                                                                                         >  4:45027/1‑130 (MQ=255)
                             cAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGc                                                                                                                         <  3:45027/130‑1 (MQ=255)
                              aGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAg                                                                                                     >  6:139577/1‑149 (MQ=255)
                                    aGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccattc                                                                                               <  6:51890/149‑1 (MQ=255)
                                    aGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccatt                                                                                                <  5:134084/148‑1 (MQ=255)
                                      gCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccattcca                                                                                             >  7:272502/1‑149 (MQ=255)
                                        cAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccattccat                                                                                            >  2:276024/1‑148 (MQ=255)
                                           tGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGt                                                                                        <  1:131238/149‑1 (MQ=255)
                                           tGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGt                                                                                        <  3:89197/149‑1 (MQ=255)
                                           tGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGt                                                                                        <  3:61815/149‑1 (MQ=255)
                                           tGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGt                                                                                        <  7:63630/149‑1 (MQ=255)
                                                   gcgcgGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTc                                                                                <  7:298986/149‑1 (MQ=255)
                                                           gcgcTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGa                                                                        <  7:153137/149‑1 (MQ=255)
                                                             gcTGATGGGGCTGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCg                                                                                                                                       <  5:20443/84‑1 (MQ=255)
                                                             gcTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCg                                                                                                                                       >  6:20443/1‑84 (MQ=255)
                                                                  tGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTg                                                                                                                <  6:152078/102‑1 (MQ=255)
                                                                  tGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTg                                                                                                                >  5:152078/1‑102 (MQ=255)
                                                                                 tCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGaa                                                   <  1:276024/148‑1 (MQ=255)
                                                                                     cAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATg                                              >  1:79665/1‑149 (MQ=255)
                                                                                             gCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGaaa                                      <  5:139577/149‑1 (MQ=255)
                                                                                                     cacCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCa                                          >  8:244745/1‑137 (MQ=255)
                                                                                                     cacCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCg                                 >  4:278259/1‑146 (MQ=255)
                                                                                                      acCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGa                             >  7:175761/1‑149 (MQ=255)
                                                                                                        cAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGAc                                                          >  8:185022/1‑118 (MQ=255)
                                                                                                        cAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGAc                                                          <  7:185022/118‑1 (MQ=255)
                                                                                                        cAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACt                           >  3:225256/1‑149 (MQ=255)
                                                                                                         aGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTg                          >  5:229443/1‑149 (MQ=255)
                                                                                                         aGGTTGTCGATACCCAGCACCAGCTCGAGAACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACgggg                                                      >  3:117375/1‑121 (MQ=255)
                                                                                                         aGGTTGTCGATACCCAGCACCAGCTCGAGAACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACgggg                                                      >  3:112703/1‑121 (MQ=255)
                                                                                                         aGGTTGTCGATACCCAGCACCAGCTCGAGAACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACgggg                                                      <  4:117375/121‑1 (MQ=255)
                                                                                                         aGGTTGTCGATACCCAGCACCAGCTCGAGAACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACgggg                                                      <  4:112703/121‑1 (MQ=255)
                                                                                                                      ccAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCcgcg              >  2:272789/1‑148 (MQ=255)
                                                                                                                      ccAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCcgcg              >  7:13687/1‑148 (MQ=255)
                                                                                                                      ccAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCcgcg              >  3:45195/1‑148 (MQ=255)
                                                                                                                      ccAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGACTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCcgcg              >  1:191551/1‑148 (MQ=255)
                                                                                                                        aGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAg           <  1:272789/149‑1 (MQ=255)
                                                                                                                              cgCTCGAGGACCCTACCCGCCAACCCGCCTAGCCCGGCCGTGGGGTCGACACGCCATTCCATTAGTGGGCACTCACGGGGGACACGTCACAAGGACGGGGGAAGCATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTc     <  2:191551/148‑1 (MQ=255)
                                                                                                                                   gAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCata  <  2:79665/147‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                               
GCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATA  >  NC_002947/3894871‑3895148

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: ‑