Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,126,647 +C intergenic (‑101/+48) tdcG‑II ← / ← gcvP‑I L‑serine dehydratase/glycine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,126,6451.C90.2% 126.7 / 7.4 41intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base C (17/20);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCA  >  NC_002947/1126500‑1126784
                                                                                                                                                        |                                                                                                                                           
gAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcc                                                                                                                                           <  3:181454/148‑1 (MQ=255)
 aGGGGCCGATGCCAATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccg                                                                                                                                         >  5:334784/1‑147 (MQ=255)
  ggggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAAAGCGGGGGTATAGCCCCCATAGCCATAGCGGGTTAAcccg                                                                                                                                         >  8:213455/1‑146 (MQ=255)
                  ttGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgc                                                                                                                                        <  6:338899/133‑4 (MQ=255)
                  ttGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTAGCGGGTAAAcccgc                                                                                                                                        >  5:338899/1‑130 (MQ=255)
                       cAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCAcc                                                                                                                    <  6:259581/149‑1 (MQ=255)
                       cAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCCCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCAcc                                                                                                                    >  5:367514/1‑149 (MQ=255)
                          gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTTACCCCCGATTCGGAGGTATAGCCACCAAAGCCATAACCGGGGAACCCGCGCCAACACACCGCCCCCCCCCCaa                                                                                                                 >  8:36212/1‑149 (MQ=255)
                          gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCaa                                                                                                                 <  2:12822/149‑1 (MQ=255)
                                     tCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCGTGTAGGCGCGGGTTTACCCCCGAATTCGGAGGTATAGCCACCATCGCAATCGCCGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAACCCAGCTGGc                                                                                                       >  1:87417/1‑149 (MQ=255)
                                            cATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAgg                                                                                                <  6:334784/148‑1 (MQ=255)
                                               gtgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCAAATGCGGAGGTATAGCCACCATCGCATTCGCCGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTa                                                                                            >  5:119681/1‑149 (MQ=255)
                                                tgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAAAGCGGAGGTATAGCCACCATCGCATTAGCGGGGAAACCCCCTCCTAAAAAAAGCCCGCCACCCACACCAGGGGGCGGGGAGGGcaa                                                                                            >  8:274106/1‑146 (MQ=255)
                                                  gACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGTAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCa                                                                                         <  4:122302/149‑1 (MQ=255)
                                                     ctcctGCCCGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc                                                                                      <  4:194802/149‑1 (MQ=255)
                                                        ctGCCAGGTGGTTGGGTGCACATTGAGCCCTGTAGGCGCGGGTTTCCCGGCGAATGCGTAGGTTTAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCg                                                                                     <  8:41704/147‑1 (MQ=37)
                                                          gCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAAACCCCAAAGCGGAGGTATAGACCCCAACGCAATAGCGGGGAAAACCCCGCCAAAACACAGCCCCCCACCCAAACCCGGTGGCCGGGCGGGCGATAACGCGGCCt                                                                                 >  8:172618/1‑149 (MQ=255)
                                                           ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCa                                                                                                                      <  3:36019/111‑1 (MQ=255)
                                                             aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGAATTGGGAGGTATAGCCACCCAAGCATTACCGGGGTAAACCCCGCCGACACAGAGCCCCCCCCCCAACCCCCGTGGCGGGGGGGGCTATAAGGCGTCCGGa                                                                               >  6:336595/1‑148 (MQ=255)
                                                                   ttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                        <  5:336595/149‑1 (MQ=255)
                                                                    tgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                        <  7:274106/148‑1 (MQ=255)
                                                                     gggTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCtctc                                                                       <  7:213455/148‑1 (MQ=255)
                                                                          gACATTGAGCCCTGTAGGCGCGGGTTTTACCCCCAATTGGGAGGGATATCCCCCAAAGCCTTAGCGGGGAAACCCCCGCCAACACACAGCCCCCCCCCCCACCCACGGGGGGGGGCGGGCTATAAGGCGGCCAGATATCGGtcgatcga                                                                 >  6:331028/1‑149 (MQ=255)
                                                                               tGAGCCCTGTAGGCGCGGGTTTAACCGCGCATGCGGAGGGATAGCCCCCATAGCATTAGCGGGTAAAACCCCTCCAACACAAAGCCCCCCCCCCCAACCAGGGGGCGGGGAGGGCTATAAGGCGTCCGGATAGCTCTCGATCGACGGGc                                                            >  2:40945/1‑149 (MQ=255)
                                                                                     cTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                        <  5:348377/131‑1 (MQ=255)
                                                                                     cTGTAGGCGCGGGTTTACCCGCAAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgcgc                                                      <  6:119681/149‑1 (MQ=255)
                                                                                             cgcgGGGTTACCCGGGAATGGGGATGTATAGCCACCATCGCATTCGCGTGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAg                                               <  6:201422/149‑1 (MQ=255)
                                                                                                  ggTTTACCCGGGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTg                                           <  3:200071/149‑1 (MQ=255)
                                                                                                  ggTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTg                                           <  7:172618/149‑1 (MQ=255)
                                                                                                      tACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCgg                                        <  3:335032/148‑1 (MQ=255)
                                                                                                       aCCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc                                      >  5:177259/1‑149 (MQ=255)
                                                                                                        cccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc                                      <  6:177259/148‑1 (MQ=255)
                                                                                                         ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                    >  2:239105/1‑149 (MQ=255)
                                                                                                         ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                    >  2:217213/1‑149 (MQ=255)
                                                                                                         ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                    >  3:38795/1‑149 (MQ=255)
                                                                                                         ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGCTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGAGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                    >  8:51647/1‑149 (MQ=255)
                                                                                                                  cGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACg                            <  6:367514/148‑1 (MQ=255)
                                                                                                                  cGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACg                            <  4:38795/148‑1 (MQ=255)
                                                                                                                       gTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACAGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                       >  5:250495/1‑148 (MQ=255)
                                                                                                                            gCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                 >  5:34854/1‑149 (MQ=255)
                                                                                                                              caccaTCGCATTCGCGGGTAAACCCGCTGCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc               >  1:213201/1‑149 (MQ=255)
                                                                                                                                       cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc       >  7:305124/1‑149 (MQ=255)
                                                                                                                                       cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc       >  1:48254/1‑149 (MQ=255)
                                                                                                                                             gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGTGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCa  >  8:185966/1‑149 (MQ=255)
                                                                                                                                                  gTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                         >  7:19538/1‑106 (MQ=255)
                                                                                                                                                  gTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                         <  8:19538/106‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                           
GAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCA  >  NC_002947/1126500‑1126784

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