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.C84.1% 123.2 / 17.5 44intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (7/0);  new base C (11/26);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 8.30e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAG  >  NC_002947/1126504‑1126785
                                                                                                                                                 |                                                                                                                                              
ggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgctc                                                                                                                                        <  8:232805/149‑6 (MQ=255)
     atgccgatCTTGAACAGGTCCACGACGCTCAGGCACATGTGGACCTCCTGCCAGGTGGTTGGGTGGACATTGAGCCCGGTAGGCGCGGGTTTACCCGCGAATGCGGAGGAATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTaca                                                                                                                                    <  2:68136/149‑1 (MQ=37)
     atgccgatCTTGAACAGGGCGAAGACGCTCAGTGACATGTGCACCTCCTGCCAGGTTGTTGGGCGCACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTaca                                                                                                                                    <  5:156589/149‑1 (MQ=37)
         cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGAAAGGCGAGGTATAGCCCCCATCGCAATAGCGGGTTAACCCGCGCCTAccca                                                                                                                                  >  6:134931/1‑147 (MQ=255)
         cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCGAAATCGGAGGTATAACCACCATCGCATTAGCGGGTTAACACGCGCCTACACaga                                                                                                                                >  1:349566/1‑149 (MQ=255)
         cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTAACCCCGCAATCGGAGGTATATCCAACCACGCATTAGCCGGGTAAACCCCGCCTAcccaaa                                                                                                                                >  3:381105/1‑147 (MQ=255)
                gAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGCATTCGGGGGTTTATCCACCCACGCATTCGCGGGTTAACCCGCTCCGACACAGAGcccccc                                                                                                                         >  1:290493/1‑149 (MQ=255)
                    aGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCGAATGCCGAGGTATAGCCACCAACGCATTAGCGGGTTAACCCCCGCCTACACACAGCCCGCcccc                                                                                                                      >  1:277011/1‑148 (MQ=255)
                      gTCGAAGACGCTCAGTGACATGGGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCaa                                                                                                                   <  3:24644/149‑1 (MQ=255)
                        cGAAGACGCTCCTTGACATGTGCACCTCCTGCCAGGTTGTGGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGc                                                                                                                <  5:40637/149‑1 (MQ=37)
                               gCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACCCCCATGGGGGGGTTTTGGCCCCCACGCCATAGCGGGTTAACACCCGCCTACACAGATCCCCCCCCCCAAACCAGgggg                                                                                                         >  5:221216/1‑149 (MQ=255)
                                 tCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCCGAGGTATAGCCACCAACGCATTAACCGGTAAACCCCCACCAAAACAAAACCCCCCCCCCAAACCACGTGGCg                                                                                                       >  8:313959/1‑149 (MQ=255)
                                     tGACTTTTGCACCTCCTGCCTGGTTGTTGGGTGTACATTGAGCCCTGTTGGGGGGGTTTTACCCGCTAATGCGTAGTTATAGCCACATTCGCATTCGGGGTTAAACCCGCTCCTACACAGAGCTCGCCACCCAAGCCAGGTGGCgg                                                                                                      <  2:402872/146‑1 (MQ=255)
                                     ggACATGTGGCGCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGa                                                                                                   <  5:127693/148‑1 (MQ=255)
                                       acttgtgTACCTCCTGCGAGGTTGTTGGGTGGACATTGTCCCCTGTAGGCGGGGGTTTACCCGGTATGGCGGAGGTATAGCCACCTTGGCATTCGGGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                                  <  8:60131/145‑1 (MQ=255)
                                            tgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTCCCGCGAATGCGGAGGTTTAGCCACCATAGCATTAGCGGGGAAACCCGCTCCTACACAGCGCCCCCCCCCCCAGCCAGGTGGCGGGGAGGGCTa                                                                                             >  2:39512/1‑148 (MQ=255)
                                             ggACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATc                                                                                           <  7:295228/149‑1 (MQ=255)
                                             ggACCTCCTGCCAGGTTGTGGGGGGGACATTGCCCCCTGTAGGCGCGGGATTACCCGCGAAGGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATc                                                                                           <  5:262448/149‑1 (MQ=38)
                                                cctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGGATAGCCACCATCGCATTAGCGGGTAAAACCGCGCCAACACAGAGCCCCCCCCCCCAGCCCCGTGGCGGGGAGGGCTATCAgg                                                                                        >  6:422435/1‑149 (MQ=255)
                                                 ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAgg                                                                                                            <  1:40768/128‑1 (MQ=255)
                                                    ctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTCACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTc                                                                                    <  8:20844/149‑1 (MQ=255)
                                                              gttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGAGCTATCAGGACGCCTGATAGcc                                                                           <  2:19042/148‑2 (MQ=37)
                                                                   gTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgat                                                                     <  6:221216/149‑1 (MQ=255)
                                                                            gAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGc                                                             <  1:343524/148‑1 (MQ=255)
                                                                            gAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCAGCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGc                                                             <  5:236696/148‑1 (MQ=255)
                                                                              gCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAgg                                                          <  5:422435/149‑1 (MQ=255)
                                                                               cccTGTAGGCGCGGGTTTACCCGCGAATGCCGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgc                                                         <  5:134931/149‑1 (MQ=255)
                                                                                    tAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGca                                                      <  2:290493/147‑1 (MQ=255)
                                                                                     attcgcgGGTAAACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCacac                                                   >  3:303449/4‑149 (MQ=39)
                                                                                      ggCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCacac                                                   <  5:317139/148‑1 (MQ=255)
                                                                                         gcgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAgg                                               >  2:110008/1‑149 (MQ=255)
                                                                                              tttACCCGCGATTGCGGAGGTATCGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                          <  4:360011/149‑1 (MQ=255)
                                                                                              tttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                          <  4:381105/149‑1 (MQ=255)
                                                                                                   ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                     >  2:69773/1‑149 (MQ=255)
                                                                                                    cgcttaTGGGGAGGTATAGCCACCATCGCATTCGGGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                    <  2:349566/144‑1 (MQ=255)
                                                                                                     gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGcc                                                                                                               >  5:166211/1‑73 (MQ=255)
                                                                                                     gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGcc                                                                                                               <  6:166211/73‑1 (MQ=255)
                                                                                                                 gTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                        >  3:392527/1‑148 (MQ=255)
                                                                                                                         accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCAGCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc               >  1:296998/1‑149 (MQ=255)
                                                                                                                           caTCGCATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAACGTCTTGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCg              <  1:346732/149‑1 (MQ=255)
                                                                                                                                 aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGt                                                                                     >  4:263420/1‑71 (MQ=255)
                                                                                                                                 aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGt                                                                                     <  3:263420/71‑1 (MQ=255)
                                                                                                                                  ttCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                        >  4:373364/1‑131 (MQ=255)
                                                                                                                                  ttCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                        <  3:373364/131‑1 (MQ=255)
                                                                                                                                        gggTAAACCCGCTCATACACAGAGCCCGCCACCCAAGCCAGGAGGCGGTGACCGCTATCAGCCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTCCCGTCGCCGAACACGATGTCGACCCGGCCGACCGGCGGCCAg  >  2:105864/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAG  >  NC_002947/1126504‑1126785

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: