Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,821,953 T→A 12.5% L173Q (CTG→CAG)  metW → methionine biosynthesis protein MetW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,821,9530TA12.5% 108.3 / 7.9 40L173Q (CTG→CAG) metWmethionine biosynthesis protein MetW
Reads supporting (aligned to +/- strand):  ref base T (14/21);  new base A (2/3);  total (16/24)
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.72e-01

CCGCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGC  >  NC_002947/5821807‑5822093
                                                                                                                                                  |                                                                                                                                            
ccGCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGc                                                                                                                                           <  2:550508/150‑1 (MQ=255)
    tGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTc                                                                                                                                       <  1:301967/150‑1 (MQ=255)
      gTACCTGGTGACCAAAGGCCGCATGCCGGTATCTGACTTCATGCCGTATGCCTGTTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGa                                                                                                                                     <  2:15496/150‑1 (MQ=255)
       tACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGAc                                                                                                                                    <  1:239313/150‑1 (MQ=255)
          cgggCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACACCCTGCCAAGGTACTCGACCGCCTGGCGGTCGACCAc                                                                                                                                 <  1:328367/148‑1 (MQ=255)
          cTGGCGACCAAAGGCCGTATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTCCACCTTCCCCGACTTCGAATAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGagtac                                                                                                                                 <  2:57295/150‑5 (MQ=255)
           tGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACt                                                                                                                                >  1:470971/1‑150 (MQ=255)
            ggCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACtt                                                                                                                               <  2:386247/150‑1 (MQ=255)
            ggCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACtt                                                                                                                               >  1:519949/1‑150 (MQ=255)
                aCCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCAc                                                                                                                           <  2:519949/150‑1 (MQ=255)
                aCCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCAc                                                                                                                           <  2:394889/150‑1 (MQ=255)
                aCCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCAc                                                                                                                           <  2:270033/150‑1 (MQ=255)
                 ccAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCAcc                                                                                                                          >  1:155298/1‑150 (MQ=255)
                            aTGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtgg                                                                                                               <  2:387039/150‑1 (MQ=255)
                              gCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtgggg                                                                                                             <  2:356173/150‑1 (MQ=255)
                                           ttCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGa                                                                                                                                     >  2:74060/1‑113 (MQ=255)
                                           ttCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGa                                                                                                                                     <  1:74060/113‑1 (MQ=255)
                                              aTGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCt                                                                                             >  2:412129/1‑150 (MQ=255)
                                                            gTACAACACGCCGAACATCCACTTCTGCACCTTCGTCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTga                                                                               <  1:470797/150‑1 (MQ=255)
                                                                     gCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCAt                                                                      <  2:365166/150‑1 (MQ=255)
                                                                               cACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCACAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTc                                                            >  2:262000/1‑150 (MQ=255)
                                                                                aCTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGCACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTg                                                                                                                              >  1:219096/1‑83 (MQ=255)
                                                                                aCTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTg                                                                                                                              <  2:219096/83‑1 (MQ=255)
                                                                                aCTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTa                                                                                            >  2:161429/1‑117 (MQ=255)
                                                                                aCTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTa                                                                                            <  1:161429/117‑1 (MQ=255)
                                                                                              gCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATg                                                                                                 >  1:40427/1‑98 (MQ=255)
                                                                                              gCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATg                                                                                                 <  2:40427/98‑1 (MQ=255)
                                                                                                       gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCAt                                    >  2:519304/1‑150 (MQ=255)
                                                                                                        aagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATc                                   >  2:481603/1‑150 (MQ=255)
                                                                                                                gtgCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTcgcg                           <  1:412129/150‑1 (MQ=255)
                                                                                                                    aCGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCt                       >  1:282707/1‑150 (MQ=255)
                                                                                                                      gAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGc                                                <  1:134154/123‑1 (MQ=255)
                                                                                                                      gAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGc                                                >  2:134154/1‑123 (MQ=255)
                                                                                                                      gAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGa                     >  1:350404/1‑150 (MQ=255)
                                                                                                                       aaCGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGAc                    <  2:496158/150‑1 (MQ=255)
                                                                                                                                      tACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAAt     <  2:282707/150‑1 (MQ=255)
                                                                                                                                      tACTCGACCGCCTGGCGGTCGACCACTTGCACCACAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAAt     <  1:262000/150‑1 (MQ=255)
                                                                                                                                         tCGACCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGGATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGc  >  2:556499/1‑150 (MQ=255)
                                                                                                                                               gCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCAt                                    <  1:255083/110‑1 (MQ=255)
                                                                                                                                               gCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCAt                                    >  2:255083/1‑110 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CCGCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGC  >  NC_002947/5821807‑5822093

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: