Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,130,138 G→T 9.8% intergenic (‑193/+111) 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,1380GT 9.8% 152.5 / 6.4 51intergenic (‑193/+111)gcvH‑I/PP_0990glycine cleavage system protein H/DoxX family protein
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base T (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

GCTCTTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACG  >  NC_002947/1129991‑1130283
                                                                                                                                                   |                                                                                                                                                     
gcTCTTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAg                                                                                                                                                       >  1:118951/1‑148 (MQ=255)
 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                                                                                                                                                   >  7:132094/1‑127 (MQ=255)
                         aataatGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATc                                                                                                                                                   <  8:132094/127‑1 (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                                                                                                     >  1:298094/1‑102 (MQ=255)
                                                                                            atttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGtttt                                                                                                     <  2:298094/102‑1 (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                                        >  7:97090/1‑143 (MQ=255)
                                                                                                                tCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCtt                                        <  8:97090/143‑1 (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             >  3:314056/3‑147 (MQ=255)
                                                                                                                                       acaaGGGTTTGATGCTCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc             >  8:31967/3‑147 (MQ=255)
                                                                                                                                           gggTTTGATGCTCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt                                 <  7:236429/123‑1 (MQ=255)
                                                                                                                                           gggTTTGATGCTCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt                                 >  8:236429/1‑123 (MQ=255)
                                                                                                                                             gTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGaaa     <  1:281762/149‑1 (MQ=255)
                                                                                                                                             gTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGaa      >  5:8036/1‑148 (MQ=255)
                                                                                                                                                 gAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCAc                            <  5:69239/122‑1 (MQ=255)
                                                                                                                                                 gAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCAc                            >  6:69239/1‑122 (MQ=255)
                                                                                                                                                 gAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAAcg  <  6:8036/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTCTTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAC‑‑CGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACG  >  NC_002947/1129991‑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: