Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,925,640 Δ1 bp 9.3% coding (561/1083 nt) gcvT‑II ← aminomethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,925,6400G. 9.3% 128.2 / 8.5 43coding (561/1083 nt)gcvT‑IIaminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base . (2/2);  total (22/21)
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.21e-01

AGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCGGGCACGCGGCCCCTGGATGGCAAGGATGGCCAGATCAGG  >  NC_002947/5925500‑5925786
                                                                                                                                             |                                                                                                                                                  
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                                                                                                                                <  1:415509/112‑1 (MQ=255)
                                                 ttGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCa                                                                                                                                >  2:415509/1‑112 (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                                                                                          <  1:742686/118‑1 (MQ=255)
                                                                                 aagaagGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTcgc                                                                                          >  2:742686/1‑118 (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)
                                                                                                                                          ataGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCGGGCACGCGGCCCCTGGATGGCAAGGATGGCCAGATCAgg  <  2:423403/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
AGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGCCTGATCGCCCGGGAAGATGATCTCCAGGCCGTCCTCACCGGTATAGCCGGTGCGGGCAATGAACCAGTCGCCGTCGGCAACGCCCTCGAACGGGCGCAGTTCGCGGATGAGTGCAGCGCGCGCGGGGCTGAGCAGGGCCGCCACCTTTTCGCGGGCACGCGGCCCCTGGATGGCAAGGATGGCCAGATCAGG  >  NC_002947/5925500‑5925786

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: