Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,218,358 C→T 9.0% V20M (GTG→ATG)  metR‑I ← homocysteine‑binding DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,218,3580CT 9.0% 117.7 / 3.3 45V20M (GTG→ATG) metR‑Ihomocysteine‑binding DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (17/23);  new base T (2/2);  total (19/26)
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.77e-01

AGCAGACGCAGCCCGGCGCTGGTGAAGCGGATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTGAAGAAACAGGGAATGAATTGAGTTTCTCTCATGATGCGGCATCTGTCGACTGGGCGATCGTCGCAC  >  NC_002947/1218219‑1218505
                                                                                                                                           |                                                                                                                                                   
agcagACGCAGCCCGGCGCTGGTGAAGCGGATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAAGCGTTCGGCGGCCTCCACCAGGCTGTCg                                                                                                                                           >  1:1148373/1‑150 (MQ=255)
  cagACGCAGCCCGGCGCTGGTGAAGCGGATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGc                                                                                                                                         >  1:931236/1‑150 (MQ=255)
       gCAGCCCGGCGCTGGTGAAGCGGATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTcgc                                                                                                                                    >  1:1179451/1‑150 (MQ=255)
            ccGGCGCTGGTGAAGCGGATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAgg                                                                                                                               <  2:623795/150‑1 (MQ=255)
            ccGGCGCTGGTGAAGCGGATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAgg                                                                                                                               <  2:1281282/150‑1 (MQ=255)
                     gTGAAGCGGATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCATCAGGCTGTCGGCTTCGCGCAGGGCATGAAgg                                                                                                                      >  2:1416051/1‑150 (MQ=255)
                     gTGAAGCGGATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCATCAGGCTGTCGGCTTCGCGCAGGGCATGAAgg                                                                                                                      >  2:1416026/1‑150 (MQ=255)
                      tGAAGCGGATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAggg                                                                                                                     <  1:1189503/150‑1 (MQ=255)
                          gCGGATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGtttt                                                                                                                 >  2:1358036/1‑150 (MQ=255)
                             gATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAg                                                                                                              >  2:721783/1‑150 (MQ=255)
                             gATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCCTCGCGCAGGGCATGAAGGGTTTTCAg                                                                                                              <  1:806352/150‑1 (MQ=255)
                             gATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCCTCGCGCAGGGCATGAAGGGTTTTCAg                                                                                                              >  2:806352/1‑150 (MQ=255)
                                 ggCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGa                                                                                                          <  2:1179451/150‑1 (MQ=255)
                                   cTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACg                                                                                                        >  1:1129116/1‑150 (MQ=255)
                                    ttGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCAACAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACgg                                                                                                       <  1:818918/150‑1 (MQ=255)
                                       gTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGAtc                                                                                                    <  1:493941/150‑1 (MQ=255)
                                           tGCGCACGAAAATCGGGAGGCCCAGGCGCTCTGCCAGCTCCTTGAACTGGTGCGACAAAGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCCTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCa                                                                                                <  1:1357861/150‑1 (MQ=255)
                                              gcACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGca                                                                                             <  1:721783/150‑1 (MQ=255)
                                               cACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGcac                                                                                            <  2:931236/150‑1 (MQ=255)
                                                   aaaaTCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGAt                                                                                        <  2:1006611/150‑1 (MQ=255)
                                                         gggAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGAGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTcc                                                                                  <  2:630331/150‑1 (MQ=255)
                                                            aGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCATCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCggg                                                                               <  1:1416026/150‑1 (MQ=255)
                                                            aGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCATCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCggg                                                                               <  1:1416051/150‑1 (MQ=255)
                                                            aGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCggg                                                                               <  1:831011/150‑1 (MQ=255)
                                                              gCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAAGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGgt                                                                             <  2:1148373/150‑1 (MQ=255)
                                                                  aGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGt                                                                                                                    >  2:1207859/1‑107 (MQ=255)
                                                                  aGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGt                                                                                                                    <  1:1207859/107‑1 (MQ=255)
                                                                         ctTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTgaa                                                                  <  2:357807/150‑1 (MQ=255)
                                                                          ttCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCAt                                                                                                                           >  1:859459/1‑92 (MQ=255)
                                                                          ttCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCAt                                                                                                                           <  2:859459/92‑1 (MQ=255)
                                                                          ttCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGgtgt                                                                           >  2:623588/1‑140 (MQ=255)
                                                                          ttCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGgtgt                                                                           <  1:623588/140‑1 (MQ=255)
                                                                          ttCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTgaag                                                                 >  1:365905/1‑150 (MQ=255)
                                                                                 tCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGAt                                                                                        >  1:402989/1‑120 (MQ=255)
                                                                                 tCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGAt                                                                                        <  2:402989/120‑1 (MQ=255)
                                                                                    ttGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAg                                                                                               <  2:783493/110‑1 (MQ=255)
                                                                                    ttGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAg                                                                                               >  1:783493/1‑110 (MQ=255)
                                                                                            gTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTGAAGAAACAGGGAATGAATTGa                                               >  2:365699/1‑150 (MQ=255)
                                                                                            gTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTGAAGAAACAGGGAATGAATTGa                                               <  1:164509/150‑1 (MQ=255)
                                                                                             tGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTGAAGAAACAGGGAATGAATTGAg                                              <  1:1358036/150‑1 (MQ=255)
                                                                                                      gCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTGAAGAAACAGGGAATGAATTGAGTTTCTCTCa                                     >  1:1324153/1‑150 (MQ=255)
                                                                                                                          gCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTgaa                                                                  <  2:1181689/101‑1 (MQ=255)
                                                                                                                          gCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTgaa                                                                  >  1:1181689/1‑101 (MQ=255)
                                                                                                                                       tCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTGAAGAAACAGGGAATGAATTGAGTTTCTCTCatgat                                 <  2:173622/121‑1 (MQ=255)
                                                                                                                                       tCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTGAAGAAACAGGGAATGAATTGAGTTTCTCTCatgat                                 >  1:173622/1‑121 (MQ=255)
                                                                                                                                         caccaGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCg                                                                                 <  2:1001286/71‑1 (MQ=255)
                                                                                                                                         caccaGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCg                                                                                 >  1:1001286/1‑71 (MQ=255)
                                                                                                                                         caccaGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTGAAGAAACAGGGAATGAATTGAGTTTCTCTCATGATGCGGCATCTGTCGACTGGGCGATCGTCGCAc  >  1:781499/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
AGCAGACGCAGCCCGGCGCTGGTGAAGCGGATTGGCTTGGTCTTGCGCACGAAAATCGGGAGGCCCAGGCGCTCTTCCAGCTCCTTGAACTGGTGCGACAACGCCGACTGGGTCAGGTGCAGGCGTTCGGCGGCCTCCACCAGGCTGTCGGCTTCGCGCAGGGCATGAAGGGTTTTCAGGTGACGGATCTCCAGCACAGATGACTCCGGGGTGTGGAATTGAAGAAACAGGGAATGAATTGAGTTTCTCTCATGATGCGGCATCTGTCGACTGGGCGATCGTCGCAC  >  NC_002947/1218219‑1218505

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