Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 117,418 T→G 7.4% K223Q (AAG→CAG)  metQ ← methionine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947117,4180TG 7.4% 159.6 / 3.0 54K223Q (AAG→CAG) metQmethionine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base G (2/2);  total (26/28)
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.54e-01

AGGGGCGGGTTTACCCGCGAAGAATAGAGCTTGATAGGTCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCTGGGTCAGCACACGCGGGA  >  NC_002947/117271‑117564
                                                                                                                                                   |                                                                                                                                                  
aGGGGCGGGTTTACCCGCGAAGAATAGAGCTTGATAGGTCGGATAAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCAGAGTCCttg                                                                                                                                                   >  7:95217/1‑149 (MQ=255)
      gggTTTACCCGCGAAGAATAGAGCTTGATAGGTCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTGGTTGccc                                                                                                                                             >  5:199511/1‑149 (MQ=255)
            aCCCGCGAATAATAGAGCTTGATAGGTCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCcgggcg                                                                                                                                        <  6:267494/148‑1 (MQ=255)
             cccGCGAAGAATAGAGCTTGATAGGTCGGATCAGACCGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCcgggcggg                                                                                                                                      >  7:70727/1‑149 (MQ=255)
              ccGCGAAGAATAGAGCTTGATAGGTCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCcgggcggg                                                                                                                                      >  2:116076/1‑148 (MQ=255)
               cgcgAAGAATAGAGCTTGATAGGTCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCgggcgggc                                                                                                                                     >  1:159402/1‑148 (MQ=255)
                   aagaaTAGAGCTTGATAGGTCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCa                                                                                                                                <  6:124783/149‑1 (MQ=255)
                       aTAGAGCTTGATAGGTCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGCCCGCCCGGACAACAAGGa                                                                                                                             <  8:1391/148‑1 (MQ=255)
                        tAGAGCTTGATAGGTCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTGGTTGCCCGCCCGGACAACAAGGact                                                                                                                           <  6:199511/149‑3 (MQ=255)
                          gagCTTGATAGGTCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGtt                                                                                                                         <  2:159402/149‑1 (MQ=255)
                           agCTTGATAGGTCGTATGAGAACGCCGGAACCACTGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCTAGTCCTGGTTGCCCGCCCGGACAACAAGGact                                                                                                                           <  2:95375/146‑3 (MQ=255)
                           agCTTGATAGGTCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTGGTTGCCCGCCCGGACAACAAGGact                                                                                                                           >  1:95375/1‑144 (MQ=255)
                                   aGGTCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAgg                                                                                                                 <  1:116076/148‑1 (MQ=255)
                                      tCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGc                                                                                                               >  4:237321/1‑147 (MQ=255)
                                        ggATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCa                                                                                                                                <  8:112353/128‑1 (MQ=255)
                                        ggATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCa                                                                                                                                >  7:112353/1‑128 (MQ=255)
                                         gATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTc                                                                                                          >  3:256901/1‑149 (MQ=255)
                                              gAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGc                                                                                                     >  3:247968/1‑149 (MQ=255)
                                              gAACGACGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGAAACCAGGATGTTCACGAAGGCCGAGTCGCTGc                                                                                                     >  8:19476/1‑149 (MQ=255)
                                               aaCGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGcc                                                                                                    >  6:303364/1‑149 (MQ=255)
                                                 cGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCt                                                                                                  >  4:171535/1‑149 (MQ=255)
                                                                gcCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGGGAGTCGCTGCCCTCGATGGCCAGCGCGt                                                                                   >  8:86134/1‑149 (MQ=255)
                                                                 ccTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTc                                                                                  <  5:303364/149‑1 (MQ=255)
                                                                                   atgaACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCtt                                                                <  4:256901/149‑1 (MQ=255)
                                                                                       aCTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCt                                                            >  5:301106/1‑149 (MQ=255)
                                                                                           cTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTc                                                          <  1:133674/147‑1 (MQ=255)
                                                                                           cTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCa                                                        <  4:247968/149‑1 (MQ=255)
                                                                                           cTTCAACTCATGCGAGTGCAGTGCAGCCGCCAGTTTCTTGATGTCGTCCGAGTCCTTGTTGTCCGGGGGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTc                                                          <  8:95217/147‑1 (MQ=255)
                                                                                                ccTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGcgcg                                                   >  1:304339/1‑149 (MQ=255)
                                                                                                 cTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGc                                                             <  3:207042/138‑1 (MQ=255)
                                                                                                 cTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGc                                                             >  4:207042/1‑138 (MQ=255)
                                                                                                    aGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAg                                                <  3:237321/148‑1 (MQ=255)
                                                                                                    aGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAAGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGt                                               >  6:210407/1‑149 (MQ=255)
                                                                                                       cGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTgg                                            <  1:235688/149‑1 (MQ=255)
                                                                                                          gtgcagtgcaGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCggg                                                                        <  2:315580/118‑1 (MQ=255)
                                                                                                          gtgcagtgcaGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCggg                                                                        >  1:315580/1‑118 (MQ=255)
                                                                                                          gtgcagtgcaGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAAGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGgtgt                                         <  5:210407/149‑1 (MQ=255)
                                                                                                             cagtgcagCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGa                                      <  4:79316/149‑1 (MQ=255)
                                                                                                              agtgcagCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGAt                                     <  8:9703/149‑1 (MQ=255)
                                                                                                              agtgcagCCGCCAGGTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGGCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGAt                                     <  8:14263/149‑1 (MQ=255)
                                                                                                                                 tcATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCTgg                  >  3:98128/1‑149 (MQ=255)
                                                                                                                                  cATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCt                    <  5:299066/146‑1 (MQ=255)
                                                                                                                                  cATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCTggg                 >  3:117744/1‑149 (MQ=255)
                                                                                                                                   aTGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCTGGGt                <  8:70727/149‑1 (MQ=255)
                                                                                                                                      tcgtcCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCTGGGTCAg             >  5:269882/1‑149 (MQ=255)
                                                                                                                                        gtcCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTGCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCTGGGTCAGc            <  6:65236/148‑1 (MQ=255)
                                                                                                                                             aGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTC‑TTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCt                        <  6:181833/130‑1 (MQ=255)
                                                                                                                                             aGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTC‑TTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCt                        >  5:181833/1‑130 (MQ=255)
                                                                                                                                              gTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGc                                                                  <  6:148497/88‑1 (MQ=255)
                                                                                                                                              gTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGc                                                                  >  5:148497/1‑88 (MQ=255)
                                                                                                                                                 cTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCTAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCt                    <  4:135001/131‑1 (MQ=255)
                                                                                                                                                 cTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCTAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCt                    >  3:135001/1‑131 (MQ=255)
                                                                                                                                                 cTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCTgg                  >  3:223505/1‑133 (MQ=255)
                                                                                                                                                 cTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCTgg                  <  4:223505/133‑1 (MQ=255)
                                                                                                                                                 cTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCTGGGTCAGCACACGCGGGa  >  3:221026/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGGGGCGGGTTTACCCGCGAAGAATAGAGCTTGATAGGTCGGATCAGAACGCCGGAACCACAGCGCCTTTGTACTTCTCGATGATGAACTGCTTCACCTCAGGCGAGTGCAGTGCAGCCGCCAGTTTCTTCATGTCGTCCGAGTCCTTGTTGTCCGGGCGGGCAACCAGGATGTTCACGTAGGGCGAGTCGCTGCCCTCGATGGCCAGCGCGTCTTTTTCCGGGTTCAGCTTGGCTTCCAGCGCGTAGTTGGTGTTGATCAGGGCGGCATCTACCTGGGTCAGCACACGCGGGA  >  NC_002947/117271‑117564

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