Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,670,463 T→C 10.9% intergenic (+52/‑33) 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,4630TC10.9% 143.9 / 7.6 55intergenic (+52/‑33)metF/PP_49785,10‑methylenetetrahydrofolate reductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/24);  new base C (3/3);  total (28/27)
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

TTGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCT  >  NC_002947/5670319‑5670602
                                                                                                                                                |                                                                                                                                           
ttGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTaaa                                                                                                                                          <  5:179989/148‑1 (MQ=255)
ttGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTaaa                                                                                                                                          >  6:179989/1‑148 (MQ=255)
  gCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGc                                                                                                                                        <  4:288398/148‑1 (MQ=255)
   cAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCt                                                                                                                                      >  8:54119/1‑149 (MQ=255)
    aaGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCt                                                                                                                                      >  5:342548/1‑148 (MQ=255)
    aaGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCt                                                                                                                                      >  5:43649/1‑148 (MQ=255)
      ggcggcGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTAt                                                                                                                                   >  1:218495/1‑149 (MQ=255)
      ggcggcGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTAGTCTATGACCAAAGcc                                                                                                                                       <  6:311663/145‑1 (MQ=255)
      ggcggcGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTAGTCTATGACCAAAGcc                                                                                                                                       >  5:311663/1‑145 (MQ=255)
          gcgcgcCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGCGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCAt                                                                                                                               <  6:314207/149‑1 (MQ=255)
               ccAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCg                                                                                                                          <  3:151188/149‑1 (MQ=255)
                  ggCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGAATAAGGCTTTAGTCTATGACCAAAGCCTTAATCAtttcttaag                                                                                                                       <  6:136294/149‑7 (MQ=37)
                          cTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAgg                                                                                                               <  2:91793/149‑1 (MQ=255)
                             cTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTc                                                                                                             <  6:43649/148‑1 (MQ=255)
                                taCCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTtgc                                                                                                         >  2:163786/1‑149 (MQ=255)
                                   cTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTtgcctg                                                                                                      >  5:318485/1‑149 (MQ=255)
                                    ttGAACCAGGCTGAACCGACCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTTGCTAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAATCCTTATTCATTTTCGGGTTAAACAGGTCTtgcctgc                                                                                                     <  7:56587/149‑1 (MQ=255)
                                          cAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATgca                                                                                               <  1:342916/149‑1 (MQ=255)
                                          cAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATg                                                                                                 >  5:66074/1‑147 (MQ=255)
                                            ggCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCtt                                                                                                           <  2:132765/135‑1 (MQ=255)
                                            ggCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCtt                                                                                                           >  1:132765/1‑135 (MQ=255)
                                                      aGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAAcc                                                                                   <  6:318485/149‑1 (MQ=255)
                                                        ccTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCtt                                                                                 >  5:212393/1‑149 (MQ=255)
                                                          tGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTcg                                                                               >  1:176949/1‑149 (MQ=255)
                                                                  tGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCca                                                                         >  3:33131/1‑147 (MQ=255)
                                                                         acCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGta                                                                <  1:184287/149‑1 (MQ=255)
                                                                          ccTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGtat                                                               <  5:120831/149‑1 (MQ=255)
                                                                          ccTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGta                                                                >  6:120831/1‑148 (MQ=255)
                                                                                     cgcgCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCtttt                                                    <  6:317013/149‑1 (MQ=255)
                                                                                       cgcTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTATAGCCTTATTCATTTTCGGGTTAAAGAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTgg                                                  <  2:286360/149‑1 (MQ=255)
                                                                                            aaaTCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCgg                                             >  3:46281/1‑149 (MQ=255)
                                                                                              aTCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGc                                            >  2:141448/1‑148 (MQ=255)
                                                                                              aTCTTTTGGCAAAGTCTGCCTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGc                                            >  4:323120/1‑148 (MQ=255)
                                                                                                   ttGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAAc                                         <  5:161843/146‑1 (MQ=255)
                                                                                                    tGGCAAAGTCTGCCTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTa                                      <  3:323120/148‑1 (MQ=255)
                                                                                                           gTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCCCAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAg                               >  2:227180/1‑148 (MQ=255)
                                                                                                           gTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGACTGCCATGCAGCACACCCAACCTTCGCGCCCAAAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACcagcag                               >  2:227185/1‑148 (MQ=255)
                                                                                                                cTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGACCGGCAACCTACCATCAGGaa                            <  4:66032/146‑1 (MQ=255)
                                                                                                                 ttAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCt                        >  6:52980/1‑149 (MQ=255)
                                                                                                                   aaGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTg                       <  2:176949/148‑1 (MQ=255)
                                                                                                                     gAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCtt                    <  4:46281/149‑1 (MQ=255)
                                                                                                                          gaatAAGGCTTTAGTCTATGACCAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCCGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACGATCAGGAAGCCTGCTTCAGc                >  6:264007/4‑148 (MQ=38)
                                                                                                                           aTTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCACCCTACCATCAGGAAGCCTGCTTCAGCAt              >  8:291094/1‑149 (MQ=255)
                                                                                                                           aTTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCAt              >  6:166512/1‑149 (MQ=255)
                                                                                                                            ttAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCCGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCAGGCTTcctgatg             <  6:342548/149‑7 (MQ=255)
                                                                                                                            ttAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATc             >  5:11368/1‑149 (MQ=255)
                                                                                                                            ttAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATc             <  4:33131/149‑1 (MQ=255)
                                                                                                                              aaGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGt           >  5:128982/1‑149 (MQ=255)
                                                                                                                               aGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGt           <  7:54119/148‑1 (MQ=255)
                                                                                                                                ggCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGt           >  8:298615/1‑147 (MQ=255)
                                                                                                                                ggCTTTAGTCTATGACCAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGaa                            <  2:170493/130‑1 (MQ=255)
                                                                                                                                ggCTTTAGTCTATGACCAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGaa                            >  1:170493/1‑130 (MQ=255)
                                                                                                                                       gTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCt  <  1:141448/149‑1 (MQ=255)
                                                                                                                                            aGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCt                                       >  4:224479/1‑107 (MQ=255)
                                                                                                                                            aGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCt                                       <  3:224479/107‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
TTGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCT  >  NC_002947/5670319‑5670602

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: ‑