Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,130,139 C→G 9.8% intergenic (‑194/+110) gcvH‑I ← / ← PP_0990 glycine cleavage system protein H/DoxX family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,130,1390CG 9.8% 153.2 / 6.6 51intergenic (‑194/+110)gcvH‑I/PP_0990glycine cleavage system protein H/DoxX family protein
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base G (3/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

CTCTTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACG  >  NC_002947/1129992‑1130283
                                                                                                                                                   |                                                                                                                                                    
ctctTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCa                                                                                                                                                     <  1:293496/149‑1 (MQ=255)
     aGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCaa                                                                                                                                                 <  1:329403/148‑1 (MQ=255)
         tAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGATGCTCAAACCCt                                                                                                                                            <  6:135413/149‑1 (MQ=255)
                 ccTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTtgt                                                                                                                                         <  4:144973/144‑1 (MQ=255)
                 ccTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTtgt                                                                                                                                         >  3:144973/1‑144 (MQ=255)
                 ccTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAc‑‑c                                                                                                                                  >  8:277432/1‑149 (MQ=255)
                 ccTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAc‑‑c                                                                                                                                  >  8:20606/1‑149 (MQ=255)
                   ttAGCAATTATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCa                                                                                                                                                  >  1:213743/1‑133 (MQ=255)
                   ttAGCAATTATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCa                                                                                                                                                  <  2:213743/133‑1 (MQ=255)
                        aataatGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATc                                                                                                                                                   <  8:132094/127‑1 (MQ=255)
                        aataatGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATc                                                                                                                                                   >  7:132094/1‑127 (MQ=255)
                             tGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCgagg‑ag                                                                                                                    <  7:20606/149‑1 (MQ=255)
                               cGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGagg‑agg                                                                                                                   <  2:118951/148‑1 (MQ=255)
                                ggCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCt                                                                                                                 <  6:214966/149‑1 (MQ=255)
                                  ccAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCtgt                                                                                                               <  7:277432/149‑1 (MQ=255)
                                    aaCCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTgtgt                                                                                                             <  4:167306/149‑1 (MQ=255)
                                     aCCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTa                                                                                                            >  4:117933/1‑149 (MQ=255)
                                       cTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCATCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTATATc                                                                                                          >  5:321761/1‑149 (MQ=255)
                                             aaTTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTc                                                                                                    <  2:210096/149‑1 (MQ=255)
                                                            aaTAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCg                                                                                      >  2:320856/1‑148 (MQ=255)
                                                                  ttACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGa                                                                               <  3:103916/149‑1 (MQ=255)
                                                                   tACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGa                                                                               >  1:61706/1‑148 (MQ=255)
                                                                      gTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCTCGCCAGTCCAGGCACCACATGTATCGTTTTCGCTACGTCGTAGCGAAAACGa                                                                               >  5:117829/1‑149 (MQ=18)
                                                                         aaTTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACAtg                                                                         <  2:61706/148‑1 (MQ=255)
                                                                                      gCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTg                                                           >  3:218890/1‑149 (MQ=255)
                                                                                        cGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGc                                                         <  5:1369/149‑1 (MQ=255)
                                                                                          tatttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGa                                                       <  3:117933/149‑1 (MQ=255)
                                                                                           atttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGtttt                                                                                                     >  1:298094/1‑102 (MQ=255)
                                                                                           atttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGtttt                                                                                                     <  2:298094/102‑1 (MQ=255)
                                                                                              tattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGGGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTg                                                   >  1:104572/1‑149 (MQ=255)
                                                                                              tattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTg                                                   >  3:292164/1‑149 (MQ=255)
                                                                                              tattGTGTGGAATGATATCGGTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGt                                                    <  5:26168/148‑1 (MQ=255)
                                                                                                              atCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGt                                    >  4:232683/1‑148 (MQ=255)
                                                                                                               tCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCtt                                        >  7:97090/1‑143 (MQ=255)
                                                                                                               tCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCtt                                        <  8:97090/143‑1 (MQ=255)
                                                                                                                cGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt                                 >  2:281762/1‑149 (MQ=255)
                                                                                                                 gTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtc                                >  4:156355/1‑149 (MQ=255)
                                                                                                                 gTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtc                                >  4:151348/1‑149 (MQ=255)
                                                                                                                        gTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACttt                         <  4:218890/149‑1 (MQ=255)
                                                                                                                            ctctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcag                      >  2:89497/1‑148 (MQ=255)
                                                                                                                                 agagGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAgat                 <  1:89497/148‑1 (MQ=255)
                                                                                                                                     gAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc             >  2:173303/1‑148 (MQ=255)
                                                                                                                                      acaaGGGTTTGATGCTCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc             >  3:314056/3‑147 (MQ=255)
                                                                                                                                      acaaGGGTTTGATGCTCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc             >  8:31967/3‑147 (MQ=255)
                                                                                                                                          gggTTTGATGCTCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt                                 <  7:236429/123‑1 (MQ=255)
                                                                                                                                          gggTTTGATGCTCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt                                 >  8:236429/1‑123 (MQ=255)
                                                                                                                                            gTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGaaa     <  1:281762/149‑1 (MQ=255)
                                                                                                                                            gTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGaa      >  5:8036/1‑148 (MQ=255)
                                                                                                                                                gAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCAc                            <  5:69239/122‑1 (MQ=255)
                                                                                                                                                gAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCAc                            >  6:69239/1‑122 (MQ=255)
                                                                                                                                                gAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAAcg  <  6:8036/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTCTTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGAC‑GGCGAGG‑AGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACG  >  NC_002947/1129992‑1130283

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑