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.C92.5% 123.8 / 3.2 40intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (16/21);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAA‑‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGC  >  NC_002947/1126514‑1126782
                                                                                                                                                |                                                                                                                                           
gatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTCACCCG‑‑‑AATGGG‑‑‑‑‑‑GAGGGATAGCCCACCAAGCCATAGCGGGGAAACCCGCGCCAACACagag                                                                                                                           >  4:21231/1‑149 (MQ=255)
gatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCG‑‑‑AAAGCG‑‑‑‑‑‑GAGGTATATCCCCCAAAGCCATAGCGGGTAAACCCCCTCCTACACagag                                                                                                                           >  2:296120/1‑149 (MQ=255)
            gTCGAAGACGCTCAGTGACATGTGCACCTCCTGCCAGGTTGTTGGGTGGACATTCACCCCTGTAGGCGCGGGTTTACCCGCT‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCaa                                                                                                               <  1:149239/149‑1 (MQ=39)
             tCGAAGACGCTCAGGGACATGGGGACCTCCTGCCAGGTTGTTGGGTGGACATTCAGCCCTGTAGGCGCGGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTacac                                                                                                                               <  3:222431/132‑1 (MQ=37)
                           tGACATGTGGACCGCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGa                                                                                                <  3:21231/149‑1 (MQ=255)
                            gACATGTGGACCTCCTGCCTGGTTTTTGGGTGTACATTTAGCCCTGTTGGCGCGGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAA‑‑CGCGCTCCTACACAGAGCCCTCCACCCAAGCGAGGTGGCGGGGa                                                                                                <  8:72234/148‑1 (MQ=37)
                             aCATGTGGCCCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                               <  5:359413/148‑1 (MQ=255)
                             aCATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGTTAAACCCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                               <  1:74505/149‑1 (MQ=255)
                              cATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                             <  1:289688/149‑1 (MQ=255)
                                   ggACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCG‑‑‑AAAGCG‑‑‑‑‑‑GGGGTTAAGGCCCCCACGCAATAGCGGGTAAAA‑‑CCCCACCTAAACAGCACCCCCCCCCCCACCCCGCTGGGGGGGAggggatt                                                                                         >  6:135971/1‑144 (MQ=255)
                                    gACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACCCG‑‑‑AATGCG‑‑‑‑‑‑GAGGGTTAGACACCATCGCATTAGCGGGTAAACCCGCGCCGAAACACAGCCCCCCCCCCAAGCCAGCTGGCGGGGAGGGCTAt                                                                                         >  2:215141/1‑147 (MQ=255)
                                      cctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCCCGCTCTCTCCGCCCCCTTAAGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAgg                                                                                              <  1:296120/149‑1 (MQ=19)
                                             ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTcc                                                                                <  3:249555/147‑1 (MQ=255)
                                                 gttgttGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCCCCATAGCATTAGCGGGTAAACCCGCTCCTACACAGAGCCCGCCCCCCAACCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAg                                                                          >  7:268046/1‑149 (MQ=38)
                                                       gggTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGa                                                                                                <  7:254610/121‑1 (MQ=255)
                                                        ggTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcga                                                                   <  5:135971/149‑1 (MQ=255)
                                                              cATTGAGCCCTGTAGGCGCGGGTTTAACCCCG‑‑‑AATGCG‑‑‑‑‑‑GAGGGATTGGCACCATAGCAATAGCCGGGAAAC‑‑CCCCGCCTAAACAGAGCCCCCCCCCCCACCCACCGGGCCGGGAGGGCTATAAAGCGGGCGGATAGCAGtcgatgcacc                                                             >  4:258151/1‑144 (MQ=255)
                                                               aTTGACCCCTTTAGGCGCGGGTTTACCCGCT‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCTGGCGTCCTGATAGCTCTCGATCGACgg                                                            <  1:192804/149‑1 (MQ=38)
                                                                    gCCCTGTAGGCGCGGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAgg                                                       <  8:268046/149‑1 (MQ=255)
                                                                             gcggGGGTTTACCCGCG‑‑‑ATTGCG‑‑‑‑‑‑GAGTTATACCCACCATCGCTTTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAg                                                                          <  1:316323/121‑1 (MQ=255)
                                                                             gcgcGGGTTTACCCGCT‑‑‑AATGCG‑‑‑‑‑‑TCGGTATCGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCa                                              <  4:27335/149‑1 (MQ=38)
                                                                               gcgGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGa                                                                                                >  4:249042/1‑97 (MQ=255)
                                                                               gcgGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGa                                                                                                <  3:249042/97‑1 (MQ=255)
                                                                               gcgGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAgg                                            >  8:154553/1‑149 (MQ=255)
                                                                               gcgGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAgg                                            >  2:271129/1‑149 (MQ=255)
                                                                               gcgGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAgg                                            >  1:344104/1‑149 (MQ=255)
                                                                                    tttACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                       <  1:215141/149‑1 (MQ=255)
                                                                                         ccGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGGCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                  >  2:284306/1‑149 (MQ=255)
                                                                                                      g‑‑‑‑‑‑gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCg                                                                                   <  1:206016/90‑1 (MQ=255)
                                                                                                      g‑‑‑‑‑‑gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCg                                                                                   >  2:206016/1‑90 (MQ=255)
                                                                                                      g‑‑‑‑‑‑gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTg                                         <  6:155538/132‑1 (MQ=255)
                                                                                                      g‑‑‑‑‑‑gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTg                                         >  5:155538/1‑132 (MQ=255)
                                                                                                                gTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGGACTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                     >  1:295333/1‑148 (MQ=255)
                                                                                                                  ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAc                  >  5:328500/1‑149 (MQ=255)
                                                                                                                   taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc                 >  7:253781/1‑149 (MQ=255)
                                                                                                                       ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg              >  8:177121/1‑149 (MQ=255)
                                                                                                                        caccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc             <  6:328500/149‑1 (MQ=255)
                                                                                                                              cGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcg       >  7:150987/1‑149 (MQ=255)
                                                                                                                                 aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg    <  7:154553/149‑1 (MQ=255)
                                                                                                                                    cgcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc  <  1:271129/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
GATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCG‑‑‑AATGCG‑‑‑‑‑‑GAGGTATAGCCACCATCGCATTCGCGGGTAAA‑‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGC  >  NC_002947/1126514‑1126782

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: