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.C88.2% 94.9 / 8.0 34intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base C (12/18);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 3.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTA  >  NC_002947/1126500‑1126787
                                                                                                                                                     |                                                                                                                                               
gatgggCCTCTGCCGAGCTTGAACAGGTCTAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAc                                                                                                                                                <  4:206109/144‑1 (MQ=39)
gAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGATGGTGGCTATACCTCCGCATTCGCGGGTAAAcc                                                                                                                                               <  8:62479/148‑1 (MQ=25)
     gccgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTCCCGCGAATGCGGGGGTATAGCCCCCATCGCATTAGCGGGTAAAcccgcgcc                                                                                                                                         >  6:286519/1‑143 (MQ=255)
     gccgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGCTTTAAACGCCAATGCGGAGGGTTAGCCACCCTAGCATTAGCAGGTAAAcacggaa                                                                                                                                          >  5:88212/1‑142 (MQ=255)
          tgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGCCCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTacac                                                                                                                                    <  1:35206/149‑1 (MQ=255)
                            cGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGc                                                                                                                  <  3:296955/149‑1 (MQ=255)
                                     tCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCCAATGCGGAGGTATAGCCACCATCGCATTAGCGGGTTAAACCGCTCCTACACAGAGCCCCCCCCCCCAGCCAGGTGGCg                                                                                                         >  6:243160/1‑149 (MQ=255)
                                        gTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTCCCGCGCATGGGGAGGGTTTGGGACCATCGCCTTAGCGGGGTAACCCGCTCCGACCCAGAGCCCGCCCCCCAAGCCAgggggggggg                                                                                                      >  7:282712/1‑149 (MQ=255)
                                           aCAGGTGCACCGCCGCCCCGGTGGTTGGGCGCACAGTCCGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATCGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGa                                                                                                     <  4:174221/147‑1 (MQ=21)
                                            cATGGGGCCCTCCTGCCAGGGTGTTGGGTGCACATTGAGCCCTGTAGGCGCGGGATTACCCGCGAATGCGGAGGTATACCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                                  <  4:28986/149‑1 (MQ=37)
                                                tgGACCACCTGCCAGGTTGGTGGGTGGACATTGCGCCCTGTAGGCGCGGGGTTACCCGCGAATGCGGAGGTATAGCCACCAGCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTAt                                                                                              <  3:257445/149‑1 (MQ=38)
                                                                 tgttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGACGCCTGATAGcc                                                                             <  5:129023/149‑2 (MQ=38)
                                                                  gttgGGTGGACATTGTGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGACGCCTGATAGcc                                                                             <  1:180644/148‑2 (MQ=37)
                                                                     gggTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCCAAATCGGAGGTATAGCCAACATAGCATTATCGGGTAAAACCGCAAATACACAGAACCCGCCCCCCAAACCAAGAGGCGGGGAGGGCAATAACGCGTGCTGATaacagac                                                                          >  4:228191/1‑142 (MQ=255)
                                                                        tGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGGGGAGGGTTAGGCACCCTCGCATTAGCGGGTTAACCCGCTCCTACACAGAGCCCGCCACCCCAGCCAGGTGGCGGGGGGGGCTTTCAGGCGTGCTGATTGctg                                                                            >  5:82221/1‑142 (MQ=255)
                                                                          gACATTGAGCCGTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGTCTATCAGGCGTCCTGATAGCTCtcgatcg                                                                     <  6:88212/148‑1 (MQ=255)
                                                                                          ggCGCGGTTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCacac                                                     <  7:108434/148‑1 (MQ=255)
                                                                                                ggTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgc                                                           <  6:274673/136‑1 (MQ=255)
                                                                                                ggTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgc                                                           >  5:274673/1‑136 (MQ=255)
                                                                                                 gTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCtctc                                                                          <  7:76027/120‑1 (MQ=255)
                                                                                                 gTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCtctc                                                                          >  8:76027/1‑120 (MQ=255)
                                                                                                      cccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCg                                        >  5:101617/1‑149 (MQ=255)
                                                                                                       ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGt                                          <  4:284773/146‑1 (MQ=255)
                                                                                                       ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGt                                          >  3:284773/1‑146 (MQ=255)
                                                                                                                cGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc                                         >  6:291567/1‑138 (MQ=255)
                                                                                                                cGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc                                         <  5:291567/138‑1 (MQ=255)
                                                                                                                     gTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCCCACCAGGTTGCGGTCGCCGAACACGTTGTCg                         >  2:218937/1‑149 (MQ=255)
                                                                                                                     gTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                          <  3:228191/148‑1 (MQ=255)
                                                                                                                      tataGCCACCATCGCATTCGCGGGTAAACCCGCTCGCCACACTGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGCATATCAGGCGACGTGATAGATGTCGAACGACGGGCAGGCGCACCCCAGGTTGCGGTCGCCGAACACGTTGCCg                         >  1:206829/1‑149 (MQ=255)
                                                                                                                         aGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGttc                                              >  6:320642/1‑123 (MQ=255)
                                                                                                                          gCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                    >  6:82182/1‑149 (MQ=255)
                                                                                                                           ccaccaTCGCATTCGCGGGTAAACCCGCTCCCACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGCCGACCCgg                   <  2:206829/149‑1 (MQ=255)
                                                                                                                               caTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGa               <  6:101617/149‑1 (MQ=255)
                                                                                                                                         gcgGGTAAACCCGCTCCTACACAGAGCCCGCCAcct                                                                                                                      <  6:256008/36‑2 (MQ=14)
                                                                                                                                         gcgGGTAAACCCGCTCCTACACAGAGCCCGCCAcct                                                                                                                      >  5:256008/1‑35 (MQ=14)
                                                                                                                                            gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg    >  4:47942/1‑148 (MQ=255)
                                                                                                                                            gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg    >  7:163226/1‑148 (MQ=255)
                                                                                                                                             ggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTa  <  7:259548/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTA  >  NC_002947/1126500‑1126787

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: