Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 274,139 T→G 11.9% K144T (AAG→ACG)  PP_0221 ← methionine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947274,1390TG11.9% 77.5 / 5.6 35K144T (AAG→ACG) PP_0221methionine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (15/15);  new base G (2/2);  total (17/18)
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.56e-01

CCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTT  >  NC_002947/273996‑274287
                                                                                                                                               |                                                                                                                                                    
ccTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCGTGGccc                                                                                                                                                >  2:456438/1‑149 (MQ=255)
       cAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcagcagcag                                                                                                                                         <  1:804518/150‑1 (MQ=255)
       cAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcagcagcag                                                                                                                                         <  1:186701/150‑1 (MQ=255)
          gTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATAAACCCTGCTCTGGCgagtagcagcag                                                                                                                                      <  2:634901/150‑1 (MQ=255)
               cggcggcgtAGGCGCGTACCAGTGGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCAc                                                                                                                                 <  2:754123/144‑1 (MQ=255)
                gtcgtcgtAGGCGCGTACCAGTTGCTCCGCCTCGGCCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTCGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACg                                                                                                                                <  1:4359/150‑1 (MQ=255)
                      gtAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCggg                                                                                                                          <  1:340658/150‑3 (MQ=255)
                                  cAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCggg                                                                                                                <  2:162801/148‑1 (MQ=255)
                                  cAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCggg                                                                                                                >  1:162801/1‑148 (MQ=255)
                                       gCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTgg                                                                                                         >  1:185683/1‑150 (MQ=255)
                                                               tCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCCGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAg                                                                                 <  1:517381/150‑1 (MQ=255)
                                                                  gCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGTCGCTGTTGATCGGGTCGTTGGCGTTGGCCACGCTGGCGCCTTCAGGCa                                                                              <  2:346024/150‑1 (MQ=255)
                                                                                             ccTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCg                                                   >  1:76862/1‑150 (MQ=255)
                                                                                                      ccTTGTAGCCGACAACGGGGTGGGGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGc                                          >  2:36141/1‑150 (MQ=255)
                                                                                                           tAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACg                                     >  2:448597/1‑150 (MQ=255)
                                                                                                               ccACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACgttgt                                 >  1:489845/1‑150 (MQ=255)
                                                                                                                 aCGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGttgttg                               <  2:489845/150‑1 (MQ=255)
                                                                                                                     cGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtgatga                           >  1:1764/1‑150 (MQ=255)
                                                                                                                     cGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtgat                             <  2:428454/148‑1 (MQ=255)
                                                                                                                     cGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtgat                             >  1:428454/1‑148 (MQ=255)
                                                                                                                          ttGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGgc                                                                                         >  2:607235/1‑83 (MQ=255)
                                                                                                                          ttGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGgc                                                                                         <  1:607235/83‑1 (MQ=255)
                                                                                                                              gggTGATCAACCCGGCCTTGGACAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGCGTACAGGCCAACGTTGTTGATGATGCCCgggg                  <  1:9563/150‑1 (MQ=255)
                                                                                                                               ggTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGc                 >  2:137048/1‑150 (MQ=255)
                                                                                                                                    tCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAg            >  2:305324/1‑150 (MQ=255)
                                                                                                                                    tCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAg            >  1:446751/1‑150 (MQ=255)
                                                                                                                                     cAACCAGCCCTTGGCCAGCAGCATCAGGACACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGGGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAgg           <  2:369893/150‑1 (MQ=255)
                                                                                                                                     aacagcGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtga                              >  1:155449/6‑131 (MQ=37)
                                                                                                                                     aacagcGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtga                              <  2:155449/126‑1 (MQ=37)
                                                                                                                                        cccGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCa                                                                                  >  1:805867/1‑76 (MQ=255)
                                                                                                                                        cccGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCa                                                                                  <  2:805867/76‑1 (MQ=255)
                                                                                                                                             ccTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGt   <  1:137048/150‑1 (MQ=255)
                                                                                                                                              cTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGt                                                >  1:791273/1‑104 (MQ=255)
                                                                                                                                              cTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGt                                                <  2:791273/104‑1 (MQ=255)
                                                                                                                                              cTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGtt  >  2:167712/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCTGCACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTT  >  NC_002947/273996‑274287

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: