Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,925,635:1 +C 9.5% coding (566/1083 nt) gcvT‑II ← aminomethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,925,6351.C 9.5% 152.9 / 11.7 42coding (566/1083 nt)gcvT‑IIaminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base . (20/18);  new base C (2/2);  total (22/20)
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.45e-01

CTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCGGGCACGCGGCCCCTGGATGGCAAGGATGGCCAGA  >  NC_002947/5925490‑5925781
                                                                                                                                                  |                                                                                                                                                  
cTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGTGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGtata                                                                                                                                                >  2:182494/1‑150 (MQ=255)
          aGGTTCATGCCGGCTTCCAGCCGCAGGGTACCCCGGGCGCCCAGGCCGCTTGGCGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGACGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCgg                                                                                                                                      <  1:714836/150‑1 (MQ=255)
           ggTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCggg                                                                                                                                     <  2:251488/150‑1 (MQ=255)
            gTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGCTATACCGGTGAgga                                                                                                                                     >  1:65938/1‑148 (MQ=255)
            gTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGCTATACCGGTGAgga                                                                                                                                     <  2:65938/149‑2 (MQ=255)
                   ccGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAAc                                                                                                                             >  2:492723/1‑150 (MQ=255)
                    cGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAAcc                                                                                                                            >  1:114957/1‑150 (MQ=255)
                    cGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAAcc                                                                                                                            >  2:379210/1‑150 (MQ=255)
                         tCCAGGGGCAGGGTATCGCGGCCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGGCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCg                                                                                                                       <  2:287518/150‑1 (MQ=255)
                              gcgcAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTc                                                                                                                  <  1:245818/150‑1 (MQ=255)
                                  aGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCa                                                                                                              >  1:431123/1‑150 (MQ=255)
                                              gcgcCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAAc                                                                                                  <  2:279983/150‑1 (MQ=255)
                                                  ccAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGgc                                                                                              >  2:35892/1‑150 (MQ=255)
                                                    aGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGgcgc                                                                                            >  2:607798/1‑150 (MQ=255)
                                                         gCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGCTATACCGGTgaggac                                                                                                                                    <  1:525593/105‑6 (MQ=255)
                                                         gCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGCTATACCGGTgaggac                                                                                                                                    >  2:525593/1‑100 (MQ=255)
                                                           ttGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGCCGACGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTcgc                                                                                     >  1:7039/1‑150 (MQ=255)
                                                           ttGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCa                                                                                                                           >  2:415509/1‑112 (MQ=255)
                                                           ttGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCa                                                                                                                           <  1:415509/112‑1 (MQ=255)
                                                           ttGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTcgc                                                                                     <  1:379210/150‑1 (MQ=255)
                                                               ggCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGAt                                                                                 <  2:141216/150‑1 (MQ=255)
                                                                 cGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGa                                                                               <  2:505414/150‑1 (MQ=255)
                                                                    tACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTg                                                                            >  2:66089/1‑150 (MQ=255)
                                                                                ccAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCgg                                                                <  2:114957/150‑1 (MQ=255)
                                                                                           aagaagGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTcgc                                                                                     >  2:742686/1‑118 (MQ=255)
                                                                                           aagaagGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTcgc                                                                                     <  1:742686/118‑1 (MQ=255)
                                                                                           aagaagGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAgg                                                     <  2:173661/150‑1 (MQ=255)
                                                                                            agaagGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGc                                                                                                                                        <  1:311929/66‑1 (MQ=255)
                                                                                            agaagGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGc                                                                                                                                        >  2:311929/1‑66 (MQ=255)
                                                                                                  ccACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGCGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCa                                              >  1:423403/1‑150 (MQ=255)
                                                                                                    aCCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCAcc                                            >  2:92092/1‑150 (MQ=255)
                                                                                                    aCCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCAcc                                            >  1:148264/1‑150 (MQ=255)
                                                                                                             tCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCgg                                   <  1:66089/150‑1 (MQ=255)
                                                                                                             tCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCgg                                   <  1:607798/150‑1 (MQ=255)
                                                                                                                        gatgatCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCGGGCACGCGGccc                        >  1:320743/1‑150 (MQ=255)
                                                                                                                                 cAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGa                                                          >  1:707705/1‑107 (MQ=255)
                                                                                                                                 cAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGa                                                          <  2:707705/107‑1 (MQ=255)
                                                                                                                                 cAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCGGGCACGCGGCCCCTGGATGGc               >  2:577606/1‑150 (MQ=255)
                                                                                                                                  aGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCGGGCACGCGGCCCCTGGATGGCa              <  2:745847/150‑1 (MQ=255)
                                                                                                                                    gCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCGGGCACGCGGCCCCTGGATGGCAAg            <  1:500052/150‑1 (MQ=255)
                                                                                                                                      cGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCGGGCACGCGGCCCCTGGATGGCAAGGa          >  2:485286/1‑150 (MQ=255)
                                                                                                                                             aCCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCGGGCACGCGGCCCCTGGATGGCAAGGATGGCCAg   >  2:319234/1‑150 (MQ=255)
                                                                                                                                              tcGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCGGGCACGCGGCCCCTGGATGGCAAGGATGGCCAGa  >  2:639993/2‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCGGGCACGCGGCCCCTGGATGGCAAGGATGGCCAGA  >  NC_002947/5925490‑5925781

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