Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,130,140 A→C 10.0% intergenic (‑195/+109) 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,1400AC10.0% 140.2 / 8.2 50intergenic (‑195/+109)gcvH‑I/PP_0990glycine cleavage system protein H/DoxX family protein
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base C (3/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

CTCTTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGC  >  NC_002947/1129992‑1130284
                                                                                                                                                    |                                                                                                                                                    
ctctTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCa                                                                                                                                                      <  1:293496/149‑1 (MQ=255)
     aGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCaa                                                                                                                                                  <  1:329403/148‑1 (MQ=255)
         tAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGATGCTCAAACCCt                                                                                                                                             <  6:135413/149‑1 (MQ=255)
                 ccTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTtgt                                                                                                                                          >  3:144973/1‑144 (MQ=255)
                 ccTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTtgt                                                                                                                                          <  4:144973/144‑1 (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‑‑CGACGGCgaggag                                                                                                                     <  7:20606/149‑1 (MQ=255)
                               cGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGaggagg                                                                                                                    <  2:118951/148‑1 (MQ=255)
                                ggCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCt                                                                                                                  <  6:214966/149‑1 (MQ=255)
                                  ccAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCtgt                                                                                                                <  7:277432/149‑1 (MQ=255)
                                    aaCCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTgtgt                                                                                                              <  4:167306/149‑1 (MQ=255)
                                     aCCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTa                                                                                                             >  4:117933/1‑149 (MQ=255)
                                       cTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCATCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTATATc                                                                                                           >  5:321761/1‑149 (MQ=255)
                                             aaTTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTc                                                                                                     <  2:210096/149‑1 (MQ=255)
                                                            aaTAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCg                                                                                       >  2:320856/1‑148 (MQ=255)
                                                                  ttACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGa                                                                                <  3:103916/149‑1 (MQ=255)
                                                                   tACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGa                                                                                >  1:61706/1‑148 (MQ=255)
                                                                      gTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCTCGCCAGTCCAGGCACCACATGTATCGTTTTCGCTACGTCGTAGCGAAAACGa                                                                                >  5:117829/1‑149 (MQ=18)
                                                                         aaTTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACAtg                                                                          <  2:61706/148‑1 (MQ=255)
                                                                                      gCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTg                                                            >  3:218890/1‑149 (MQ=255)
                                                                                        cGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGc                                                          <  5:1369/149‑1 (MQ=255)
                                                                                          tatttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGa                                                        <  3:117933/149‑1 (MQ=255)
                                                                                           atttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGtttt                                                                                                      <  2:298094/102‑1 (MQ=255)
                                                                                           atttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGtttt                                                                                                      >  1:298094/1‑102 (MQ=255)
                                                                                              tattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGGGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTg                                                    >  1:104572/1‑149 (MQ=255)
                                                                                              tattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTg                                                    >  3:292164/1‑149 (MQ=255)
                                                                                              tattGTGTGGAATGATATCGGTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGt                                                     <  5:26168/148‑1 (MQ=255)
                                                                                                              atCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGt                                     >  4:232683/1‑148 (MQ=255)
                                                                                                               tCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCtt                                         <  8:97090/143‑1 (MQ=255)
                                                                                                               tCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCtt                                         >  7:97090/1‑143 (MQ=255)
                                                                                                                cGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt                                  >  2:281762/1‑149 (MQ=255)
                                                                                                                 gTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtc                                 >  4:156355/1‑149 (MQ=255)
                                                                                                                 gTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtc                                 >  4:151348/1‑149 (MQ=255)
                                                                                                                        gTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACttt                          <  4:218890/149‑1 (MQ=255)
                                                                                                                            ctctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcag                       >  2:89497/1‑148 (MQ=255)
                                                                                                                                 agagGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAgat                  <  1:89497/148‑1 (MQ=255)
                                                                                                                                     gAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc              >  2:173303/1‑148 (MQ=255)
                                                                                                                                      acaaGGGTTTGATGCTCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc              >  8:31967/3‑147 (MQ=255)
                                                                                                                                      acaaGGGTTTGATGCTCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc              >  3:314056/3‑147 (MQ=255)
                                                                                                                                          gggTTTGATGCTCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt                                  >  8:236429/1‑123 (MQ=255)
                                                                                                                                          gggTTTGATGCTCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt                                  <  7:236429/123‑1 (MQ=255)
                                                                                                                                            gTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGaaa      <  1:281762/149‑1 (MQ=255)
                                                                                                                                            gTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGaa       >  5:8036/1‑148 (MQ=255)
                                                                                                                                                gAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCAc                             >  6:69239/1‑122 (MQ=255)
                                                                                                                                                gAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCAc                             <  5:69239/122‑1 (MQ=255)
                                                                                                                                                gAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAAcg   <  6:8036/148‑1 (MQ=255)
                                                                                                                                                    aTCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAAcgc  <  4:3065/145‑1 (MQ=255)
                                                                                                                                                    aTCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAAcgc  >  3:3065/1‑145 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCTTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGC  >  NC_002947/1129992‑1130284

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