Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,922,743 G→A 9.3% A563V (GCG→GTG)  gcvP‑II ← glycine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,922,7430GA 9.3% 116.6 / 2.5 43A563V (GCG→GTG) gcvP‑IIglycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (21/18);  new base A (2/2);  total (23/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 = 8.88e-01

TGCCGACGCCGGGTTGGTGCCGTGGGCCGACGACGGGATCAGGCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATGGGGATCATCTCACTGGTGGCGTTGAGCT  >  NC_002947/5922598‑5922890
                                                                                                                                                 |                                                                                                                                                    
tgCCGACGCCGGGTTGGTGCCGTGGGCCGACGACGGGATCAGGCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCAt                                                                                                                                                   >  2:183908/1‑149 (MQ=255)
  ccGACGCCGGGTTGGTGCCGTGGGCCGACGACGGGATCAGGCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATc                                                                                                                                                  >  6:104785/1‑148 (MQ=255)
        ccGGGTTGGTGCCGTGGGCCGACGACGGGATCAGGCACATCGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCg                                                                                                                                           <  3:350819/149‑1 (MQ=255)
                gTGCCGTGGGCCGACGACGGGATCAGGCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATcgc                                                                                                                                   <  5:104785/149‑1 (MQ=255)
                gTGCCGTGGGCCGACGACGGGATCAGGCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATcgc                                                                                                                                   >  5:185377/1‑149 (MQ=255)
                gTGCCGTGGGCCGACGACGGGATCAGGCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCACATCGAAGCCGGTGATcgc                                                                                                                                   >  3:353558/1‑149 (MQ=255)
                                    gATCAGGCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCACATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGtt                                                                                                               <  4:353558/149‑1 (MQ=255)
                                     aTCAGGCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCCCAGCCAGCTTTCCAGTtc                                                                                                              >  8:337508/1‑149 (MQ=255)
                                       cAGGCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTtcgt                                                                                                            >  3:78473/1‑149 (MQ=255)
                                         ggCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGGCAGCTTTCCAGTtcgtcg                                                                                                          <  4:107004/149‑1 (MQ=255)
                                          gCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTGTCCAGTTCGTCGa                                                                                                         >  4:248906/1‑149 (MQ=255)
                                          gCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCAGCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGa                                                                                                         >  2:99292/1‑149 (MQ=255)
                                          gCACAACGTGCGCACCGGCTGGTGGCGGCTGCAGTGATAGCGCGTAATGGCCTTCAGGCCGGCGTATGCACCCTGAGCACCCGAGTTCGGCTGCATGCAGATCGCATCGAACCCGGTGATCGCGCACAGCCAGCTTTCCAGTTGGTCGa                                                                                                         <  2:111299/149‑1 (MQ=255)
                                           cacaACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGat                                                                                                        <  6:185377/149‑1 (MQ=255)
                                                     gcATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATCGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTg                                                                                              >  7:266489/1‑149 (MQ=255)
                                                                               tAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGc                                                                   <  2:116428/149‑1 (MQ=255)
                                                                                aGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGcc                                                                                                 <  2:222722/119‑1 (MQ=255)
                                                                                aGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGcc                                                                                                 >  1:222722/1‑119 (MQ=255)
                                                                                  cGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGaa                                                                <  7:337508/149‑1 (MQ=255)
                                                                                     gTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACgg                                                             <  1:183908/149‑1 (MQ=255)
                                                                                       aaTGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATcg                                                                                                                                    >  8:130070/1‑77 (MQ=255)
                                                                                       aaTGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATcg                                                                                                                                    <  7:130070/77‑1 (MQ=255)
                                                                                          ggCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGACCGGGtt                                                          >  8:77142/1‑146 (MQ=255)
                                                                                             cATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCt                                                     >  6:114352/1‑149 (MQ=255)
                                                                                                      ggCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCAGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAAc                                            >  1:161361/1‑149 (MQ=255)
                                                                                                         gTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCgg                                         <  2:161361/149‑1 (MQ=255)
                                                                                                          tATTCACCCTGCGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTGTTAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCgg                                         <  1:89868/149‑1 (MQ=255)
                                                                                                          tATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCggg                                        <  4:78473/149‑1 (MQ=255)
                                                                                                               aCCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAgg                                   >  6:121621/1‑149 (MQ=255)
                                                                                                                     aGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCt                                                                                >  5:98929/1‑98 (MQ=255)
                                                                                                                     aGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCt                                                                                <  6:98929/98‑1 (MQ=255)
                                                                                                                     aGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATg                              <  1:297951/148‑1 (MQ=255)
                                                                                                                      gCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATggg                            >  6:24826/1‑149 (MQ=255)
                                                                                                                            gAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATGGGGatcat                      >  2:234022/1‑149 (MQ=255)
                                                                                                                             aggtGGGCTGCATGCAGATCGCATCGAAGCCGGTGACCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATGGGGatcatc                     <  8:148167/146‑1 (MQ=255)
                                                                                                                              gTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATGGGGATCAtct                    <  1:234022/149‑1 (MQ=255)
                                                                                                                                  ggCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATGGGGATCATCTCAc                 <  8:292507/148‑1 (MQ=255)
                                                                                                                                  ggCTGCATGCAGATCACATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATGGGGATCATCTCACt                >  5:75015/1‑149 (MQ=255)
                                                                                                                                       catgcaGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATGGGGATCATCTCACtggtgg           >  1:43257/1‑149 (MQ=255)
                                                                                                                                           caGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCGTTCCAGTTCGGCGATCAGCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGCGGCGAACGGGTGCAGCTGGGCAACGCATGGCGAGGTGAGGGGGATCATCTCACTGGTGGCg         >  4:26005/1‑147 (MQ=255)
                                                                                                                                             gATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATGGGGATCATCTCACTGGTGGCGtt       >  2:116096/1‑147 (MQ=255)
                                                                                                                                             gATCACATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATGGGGATCATCTCACTGGTGGCGTTGa     <  6:75015/149‑1 (MQ=255)
                                                                                                                                                cGCATCGAAGCCGGTGATCACGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATGGGGATCATCTCACTGGTGGCGTTGAGCt  >  5:371562/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGCCGACGCCGGGTTGGTGCCGTGGGCCGACGACGGGATCAGGCACAACGTGCGCATCGGCTGGTGGCGGCTGCAGTGATAGCGGGTAATGGCCATCAGGCCGGCGTATTCACCCTGAGCACCGGAGTTGGGCTGCATGCAGATCGCATCGAAGCCGGTGATCGCGCACAGCCAGCTTTCCAGTTCGTCGATCATCGCCTTG‑TAGCCCGCTGCCTGGGCCGCCGGGGCGAACGGGTGCAGCTGGGCAAAACCGGGCCAGGTGATGGGGATCATCTCACTGGTGGCGTTGAGCT  >  NC_002947/5922598‑5922890

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