Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,118,147 G→T C59* (TGC→TGA)  metJ ← transcriptional repressor, S‑adenosylmethionine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,118,1470GT95.8% 133.5 / NA 50C59* (TGC→TGA) metJtranscriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base T (20/26);  minor base C (1/1);  total (22/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 = 1.00e+00

TTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATG  >  CP009273/4118001‑4118292
                                                                                                                                                  |                                                                                                                                                 
ttGCCTGAGTATTCCCGCGTCTCCGGGTTAATCCCCCTCTCCCGCATGTCCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTca                                                                                                                                                 <  2:966142/149‑1 (MQ=255)
ttGATTTAGTATCCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTca                                                                                                                                                 <  2:88694/149‑1 (MQ=255)
     ttAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCAcagca                                                                                                                                            >  2:76870/1‑149 (MQ=255)
           ttCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCGGAAACGCTTCTCAcagcag                                                                                                                                           >  1:130177/1‑144 (MQ=255)
               cACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCg                                                                                                                                       <  1:929151/144‑1 (MQ=255)
                acgcccccGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTg                                                                                                                                 <  2:676980/143‑1 (MQ=255)
                          gTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGATCggaagagcgtcg                                                                                                                            >  1:594652/1‑133 (MQ=255)
                          gTTAATCCCCATCGCCCGCATGGTCTCGTTTGCCGCTTCCGGGATTGCGTCGCGGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGGTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGAc                                                                                                                       <  2:323326/149‑1 (MQ=255)
                               tCCCCATCTCACGCATGCTCTCTTTTGCCGCTTCCGGGATTTCGTCGATGCGCTCATTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGTATGCAGAAACGCTTCAAACAGCAGCTCGCTGTTGGTAGCGTGACGCAgg                                                                                                                  <  2:762664/149‑1 (MQ=255)
                                      ctcGCGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgtt                                                                                                                <  1:105894/144‑1 (MQ=255)
                                      cccacgtcttctctcTTTTGCCGCGTCCGGGATTTCGTCGCTGCGCCCTTGACGCAGATCGGCATCATCCGGCAAAGGTGGCCCGGTAAAGGCATGCATAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgtt                                                                                                                <  1:20755/133‑1 (MQ=255)
                                         aCGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgttgtt                                                                                                             >  1:801491/1‑144 (MQ=255)
                                         aCGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGCGGGCCCGGTAAAGGCATGAAGAAACGCTTCCCACAGAAGCTCGCTGTCAGTAGCGTGACGCAGgttgtt                                                                                                             >  1:260252/1‑144 (MQ=255)
                                             aTGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGACACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGCCGCCGGTTGTTCaac                                                                                                         >  1:287125/1‑142 (MQ=255)
                                                 tctcTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                                                                                                     >  1:944792/1‑144 (MQ=255)
                                                 tctcTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                                                                                                     <  1:138407/144‑1 (MQ=255)
                                                              ttCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTc                                                                                        <  1:308946/144‑1 (MQ=255)
                                                                 cGGGATTTCGGCGCCGCGCCCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATc                                                                                     <  1:395758/144‑1 (MQ=255)
                                                                    gATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGt                                                                                  <  1:501504/144‑1 (MQ=255)
                                                                        tcgtcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCCCAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtttt                                                                         >  2:66116/1‑149 (MQ=255)
                                                                            cgCTGCGCTCTTTACGCTGATCGGCATCTTCCGGCAAAGGTTTCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAac                                                                     <  2:351020/149‑1 (MQ=255)
                                                                            cgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAttt                                                                          <  1:600870/144‑1 (MQ=255)
                                                                              cTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAttttt                                                                        <  1:322809/144‑1 (MQ=255)
                                                                                gcgcTCTTTACGCAGATCGGCATACTCCGGCAAAGGTTGCCCGGTCTAATCATCCAGAAACGCTTCTCACAGCAGCTCGCTGTGGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCt                                                                 <  2:237555/149‑1 (MQ=255)
                                                                                   ctctTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTaa                                                              >  2:387682/1‑149 (MQ=255)
                                                                                    tctTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACAcc                                                                  >  1:890439/1‑144 (MQ=255)
                                                                                     ctTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGATATGCACTTTCTTCTACTTACGACGCGTGCGTTCATCGGTGATGATTTTTAACACCTTAaga                                                            <  2:246550/149‑1 (MQ=255)
                                                                                      tttACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTCTTGGTCGCGTGACGCAGGTTGTTCACATGACGACGCATGCGGTCATCGGTGAGGATTGTTAAGACCTTAggag                                                           >  2:216366/1‑149 (MQ=255)
                                                                                          cGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAaga                                                            <  1:512825/144‑1 (MQ=255)
                                                                                           gCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTACAGGCATGCAGAAACGCTTCTCACAGCAGCTCACAGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGAATTTTTAACACCTTAagag                                                           >  1:613774/1‑144 (MQ=255)
                                                                                           gCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATg                                                      >  2:151441/1‑149 (MQ=255)
                                                                                                  ggCACCATCGGGCAAAGGCTGCCCGGGAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCg                                               <  2:69484/149‑1 (MQ=255)
                                                                                                         tcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGtc                                             >  1:235705/1‑143 (MQ=255)
                                                                                                            ggCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCttttt                                     >  2:792074/1‑149 (MQ=255)
                                                                                                                 cggTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCttttt                                     <  1:550058/143‑1 (MQ=255)
                                                                                                                 aGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCttttt                                     <  1:489518/144‑1 (MQ=255)
                                                                                                                  ggTTGCCCGGTAAAGGCGTGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCAGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTg                                    >  1:696593/1‑144 (MQ=255)
                                                                                                                     tGCCCGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTc                            >  2:233032/1‑149 (MQ=255)
                                                                                                                       cccGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTg                               <  1:704398/144‑1 (MQ=255)
                                                                                                                         cGGTAAAGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCATTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGtt                             <  1:59467/144‑1 (MQ=255)
                                                                                                                              aaGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTcttctt                   <  2:829283/149‑1 (MQ=255)
                                                                                                                               aGGCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAccct                       >  1:138685/1‑144 (MQ=255)
                                                                                                                                 gCATGCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtcttc                     >  1:597710/1‑144 (MQ=255)
                                                                                                                                  cATGCAGAAACGCTTCCCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtcttct                    <  1:326695/144‑1 (MQ=255)
                                                                                                                                     gCAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGc                 <  1:77164/144‑1 (MQ=255)
                                                                                                                                      cAGAAACGCTTCTCACAGCAGCTCGCTGTTGGTATCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGcc                <  1:106646/144‑1 (MQ=255)
                                                                                                                                        gAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGCCGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAAGGGAAACCGTAATCTTCTTGACTTGTTCACTCTTCTTGCCGt              >  1:315775/1‑144 (MQ=255)
                                                                                                                                        gAAACGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCa         <  2:874645/149‑1 (MQ=255)
                                                                                                                                           aCGCTTCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGGTTTTTAACACCTTAAGAGGAATGGAAACCGTAATGTTTTTGACTTGTTCACTCTTCTTGCCGTCCTCAGCg      >  2:368248/1‑149 (MQ=255)
                                                                                                                                               ttCTCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGACTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATg  >  2:799262/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATG  >  CP009273/4118001‑4118292

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