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.C94.4% 51.2 / ‑5.8 19intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (6/11);  minor base G (1/0);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  NC_002947/1126512‑1126788
                                                                                                                                        |                                                                                                                                               
ccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                    <  4:17094/148‑1 (MQ=255)
 cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATTGGGAGGTATAGCCACCATCGCATTAGCGGGTAAACCCGCTCCTAcaca                                                                                                                                    >  4:118440/1‑147 (MQ=38)
                       gCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCGAATGCGGAGGTTTAGCCCCCCTAGCATTAGCGGGGAAACCCCCGCCTAAAAACAGCCCCCCACCCCACCCCGGTgg                                                                                                            >  4:37747/1‑149 (MQ=255)
                                cATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                                   <  8:113189/149‑1 (MQ=255)
                                     gccccccctGCCAGGTTGGGGGGGGCACATGGCGCCCTGTCGGCGCGGGTTAACCCGCGAAGGCGGAGGTATCCCCACCAGCGCATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATc                                                                                              <  1:130957/143‑1 (MQ=255)
                                        cctcctGCCAGGTTGTGGGGTGGACATTCAGCCCTGTAGGCGCGGGTTAACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAgg                                                                                           <  7:67050/149‑1 (MQ=38)
                                         ctcctGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc                                                                                          <  4:118373/149‑1 (MQ=255)
                                            ctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGGGCGGTTTTACCCGCGATTGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTc                                                                                       <  6:165235/149‑1 (MQ=38)
                                             tGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTTACCCCCAAAGCGGGGGGGTAGGCCCCAACGCAATAGCGGGGGAAACCCCGCCAACACAGACCCCCCCCCCCCCCCCCCGGGGCGGGGGGGGGGATAAGGCgggc                                                                                      >  6:130321/1‑146 (MQ=255)
                                                                 aTTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg                                                                  <  3:118440/149‑1 (MQ=255)
                                                                        ccTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgcg                                                           <  4:120520/149‑1 (MQ=255)
                                                                                            ccGCGACTGCGGAGGTATAGCCACCATCGCCTTCGCGGGTAAAGCGGCTCGTCCCCGGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTCTCAGGCGGCCTGATAGCTGTCGATCGCCGGGCAGGGGCACCCCCGGTTGCGGTCGc                                        >  4:27382/1‑149 (MQ=255)
                                                                                              gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTacac                                                                                                                                     >  3:120102/1‑54 (MQ=18)
                                                                                              gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTacac                                                                                                                                     <  4:120102/54‑1 (MQ=18)
                                                                                                        aGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGt                            <  8:27564/149‑1 (MQ=255)
                                                                                                                ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGACCACGTTGTCGAACCgg                    >  3:110914/1‑149 (MQ=255)
                                                                                                                          aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg          >  1:19151/1‑149 (MQ=255)
                                                                                                                                 gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGt    <  5:130321/149‑1 (MQ=255)
                                                                                                                                   gTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAc  >  1:66416/1‑149 (MQ=255)
                                                                                                                                        |                                                                                                                                               
CCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  NC_002947/1126512‑1126788

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: