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.C91.7% 69.9 / 0.9 24intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (5/17);  total (7/17)
Fisher's exact test for biased strand distribution p-value = 7.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCA  >  NC_002947/1126498‑1126784
                                                                                                                                                     |                                                                                                                                           
tGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAAAGCGGAGGGTTAGCCCCCAACGCAATAGCGGGGaaa                                                                                                                                              >  8:9619/1‑148 (MQ=255)
 ggAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAc                                                                                                                                            <  2:66266/148‑1 (MQ=255)
               cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                <  6:52948/147‑1 (MQ=255)
                       aaCAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGc                                                                                                                            <  3:112681/143‑1 (MQ=255)
                             tCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCaa                                                                                                                 <  8:141210/148‑1 (MQ=255)
                                          gTGACATGTGGACCTCCTGCGAGGTTGTTGGGTGGAAATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGTGGTATAGCCTCCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg                                                                                                     <  4:151316/147‑1 (MQ=38)
                                          gTGACATGTGGACCTCCTGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCg                                                                                                      <  6:188298/146‑1 (MQ=255)
                                          gTGACATGTGGACCTCCTGCCAGGTGGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg                                                                                                     <  8:188442/147‑1 (MQ=255)
                                             aCATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGCCACCATCGCATTAGCGGGTAAAACCCCTCCTACACAGAGCCCCCCCCCCAAGCCACGTGGCGGGGAg                                                                                                 >  2:19439/1‑148 (MQ=255)
                                                     aCCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCa                                                                                         <  8:38854/148‑1 (MQ=255)
                                                       ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACACCCAATGCGGAGGGATAAGCAACCTCGCAATAGCCGGTAAAACCCCAACAACACAAAGCCCCCCCCCCAAACCCAGTGGCGGGGAGGGCTATCacg                                                                                       >  6:135852/1‑146 (MQ=255)
                                                           tGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGt                                                                                    <  4:63077/147‑1 (MQ=255)
                                                                  ttgttgGGTGTAATTTGAGCCCGTTCTGCGGTGGTTTAACCGCGACTGCGTAGGTATCGCCACAATCGCTTTGGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAATCCAGGTGGCGGGGAGGGATATCAGGCGTCCTGATa                                                                             <  6:18750/147‑1 (MQ=255)
                                                                     ttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGACGCCTGATAGcc                                                                          <  4:95177/147‑2 (MQ=38)
                                                                               aTTGAGCCCTGAAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAACCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGCTAGCTCTCGATCGACg                                                               >  8:150522/1‑148 (MQ=38)
                                                                                    gCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAg                                                          <  7:9619/148‑1 (MQ=255)
                                                                                                 gggTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCCAGCCAGGTGGCGGGGAGGGCTATCAGGCGGCCTGATAGCTCGCGATCGCCGGGCAGGCGCCCACCAGGt                                             >  1:207387/1‑148 (MQ=38)
                                                                                                         ccGCTACTGCGTAGGTCTCGCCCCCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGTCTATCAGGCGTCCTTATAGCTCGCGATCGACGGGAAGGCGCACACCAGGTTGCGGTc                                      <  5:79163/147‑1 (MQ=255)
                                                                                                              aaTGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGaa                                <  7:40870/148‑1 (MQ=255)
                                                                                                                          tcccccccaTCGCATGCGCGGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAc                    <  2:194108/142‑1 (MQ=255)
                                                                                                                             ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                 >  8:7921/1‑148 (MQ=255)
                                                                                                                                  aTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAt                                                                              >  7:168942/1‑82 (MQ=255)
                                                                                                                                   tCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATa                                                                             <  8:168942/82‑1 (MQ=255)
                                                                                                                                       aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc       <  1:119406/148‑1 (MQ=255)
                                                                                                                                          cgcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcgg     >  4:45705/1‑147 (MQ=255)
                                                                                                                                            cgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCGAGGGGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCa  <  5:107974/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
TGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCA  >  NC_002947/1126498‑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: