Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,670,459 A→T 7.4% intergenic (+48/‑37) metF → / → PP_4978 5,10‑methylenetetrahydrofolate reductase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,670,4590AT 7.4% 159.2 / 3.3 54intergenic (+48/‑37)metF/PP_49785,10‑methylenetetrahydrofolate reductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/28);  new base T (2/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

CAGCTGTTGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATC  >  NC_002947/5670313‑5670591
                                                                                                                                                  |                                                                                                                                    
cAGCTGTTGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGa                                                                                                                                    <  7:229490/149‑1 (MQ=255)
 aGCTGTTGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTTGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGAc                                                                                                                                   <  2:120569/149‑1 (MQ=255)
 aGCTGTTGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGAc                                                                                                                                   >  6:305660/1‑149 (MQ=255)
      ttGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAg                                                                                                                              >  1:14657/1‑149 (MQ=255)
        gCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCCTTGGTCATAGACTAAAGcc                                                                                                                            >  4:98722/1‑149 (MQ=255)
         cAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGcc                                                                                                                            <  2:14657/148‑1 (MQ=255)
         cAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGAATAAGGCTTTAGTCTATGACCAAAGcc                                                                                                                            <  6:107725/148‑1 (MQ=255)
          aaGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCt                                                                                                                           >  1:199773/1‑148 (MQ=255)
             gcggcgCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTAtt                                                                                                                       <  5:305660/149‑1 (MQ=255)
               ggagcgcCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGACTGTCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCa                                                                                                                     <  4:201387/146‑1 (MQ=255)
                   cgcCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGAGCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCAtttt                                                                                                                 <  2:222260/149‑1 (MQ=255)
                    gcCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTc                                                                                                                >  5:12593/1‑149 (MQ=255)
                      cAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTGAAGCCTTATTCATTTTCgg                                                                                                              <  8:210279/149‑1 (MQ=255)
                      cAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCgg                                                                                                              <  4:26941/149‑1 (MQ=255)
                       aGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCCTTGGTCATAGACTAAAGCCTTATTCATTTTCggg                                                                                                             <  3:98722/149‑1 (MQ=255)
                         gCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGtt                                                                                                           <  4:275394/149‑1 (MQ=255)
                                   cTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCt                                                                                                 >  8:65085/1‑149 (MQ=255)
                                    tataCCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCtt                                                                                                <  3:252568/149‑1 (MQ=255)
                                     ataCCTTCAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGATTTGGTCATAGACTAAAGCCTTATTCATTTTCTGGTTAAACAGTTCTtg                                                                                               >  4:69219/1‑149 (MQ=255)
                                         cTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTtgcctg                                                                                           >  6:249845/1‑149 (MQ=255)
                                                cAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATgca                                                                                    <  2:130928/149‑1 (MQ=255)
                                                    cTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATgcagca                                                                                 >  4:223563/1‑148 (MQ=255)
                                                        aCCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACAcc                                                                             >  6:93866/1‑148 (MQ=255)
                                                          cGAGCCTGGCGATCTGGAACAACCTGAAGCTGGCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCaa                                                                          <  5:62965/149‑1 (MQ=255)
                                                               cTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTc                                                                     <  3:288995/149‑1 (MQ=255)
                                                                 ggCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAAc                                                                         <  8:125617/143‑1 (MQ=255)
                                                                 ggCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAAc                                                                         >  7:125617/1‑143 (MQ=255)
                                                                  gCGATCTGGAACAACCTGAAGCTGCCTCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTcgc                                                                   <  4:254633/148‑1 (MQ=255)
                                                                  gCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTcgc                                                                   <  7:65085/148‑1 (MQ=255)
                                                                       cTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTtgcctg                                                                                           <  2:50171/119‑1 (MQ=255)
                                                                       cTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTtgcctg                                                                                           >  1:50171/1‑119 (MQ=255)
                                                                          ggacaacCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCca                                                              <  6:332393/143‑1 (MQ=255)
                                                                          ggacaacCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCca                                                              >  5:332393/3‑145 (MQ=255)
                                                                                ccTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGtat                                                    <  3:69219/149‑1 (MQ=255)
                                                                                    aaGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATAcc                                                 <  6:12593/148‑1 (MQ=255)
                                                                                              gcTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGc                                      >  4:63466/1‑149 (MQ=255)
                                                                                                  aaaTCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCg                                   <  3:63466/148‑1 (MQ=255)
                                                                                                  aaaTCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCg                                   <  3:223563/148‑1 (MQ=255)
                                                                                                       ttttGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAAcc                             >  4:195047/1‑149 (MQ=255)
                                                                                                         ttGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAAc                              <  2:199773/146‑1 (MQ=255)
                                                                                                            gCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCa                        >  3:32658/1‑149 (MQ=255)
                                                                                                            gCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCa                        >  2:267485/1‑149 (MQ=255)
                                                                                                                      cTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGcc              >  2:202/1‑149 (MQ=255)
                                                                                                                       ttAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTgg                                       >  6:251080/1‑123 (MQ=255)
                                                                                                                       ttAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTgg                                       <  5:251080/123‑1 (MQ=255)
                                                                                                                         cgaaaaTGAATAAGGCTTTAGTCTATGACCAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGc           >  5:93860/4‑149 (MQ=38)
                                                                                                                              aTGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCg                                   >  6:119593/1‑120 (MQ=255)
                                                                                                                              aTGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCg                                   <  5:119593/120‑1 (MQ=255)
                                                                                                                                      ggCTTTAGTCTATGACCAAAGCCTTAATCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATc  <  6:172451/145‑1 (MQ=255)
                                                                                                                                      ggCTTTAGTCTATGACCAAAGCCTTAATCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATc  >  5:172451/1‑145 (MQ=255)
                                                                                                                                             gTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTg            >  1:72096/1‑128 (MQ=255)
                                                                                                                                             gTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTg            <  2:72096/128‑1 (MQ=255)
                                                                                                                                                aTAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCaa                               <  1:89549/106‑1 (MQ=255)
                                                                                                                                                aTAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCaa                               >  2:89549/1‑106 (MQ=255)
                                                                                                                                                  |                                                                                                                                    
CAGCTGTTGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATC  >  NC_002947/5670313‑5670591

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