Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 782,046 T→A 10.0% intergenic (‑201/+49) glyA‑II ← / ← PP_0672 serine hydroxymethyltransferase/sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947782,0460TA10.0% 140.0 / 5.6 50intergenic (‑201/+49)glyA‑II/PP_0672serine hydroxymethyltransferase/sensory box protein
Reads supporting (aligned to +/- strand):  ref base T (22/23);  new base A (3/2);  total (25/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 = 9.96e-01

GCGGCTGGCATATGAAGTGGCGTCAAGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCAC  >  NC_002947/781899‑782181
                                                                                                                                                   |                                                                                                                                       
gcgGCTGGCATATGAAGTGGCGTCAAGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTg                                                                                                                                        >  6:201223/1‑149 (MQ=255)
            tGAAGTGGCGTCAAGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTcgcg                                                                                                                            >  4:212001/1‑149 (MQ=255)
              taGTGGCGTCAAGTGGCAGACAAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGGTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                           <  2:83254/147‑1 (MQ=255)
               aGTGGCGTCAAGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGt                                                                                                                         >  3:305396/1‑149 (MQ=255)
                   gCGTCAAGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAAc                                                                                                                     <  3:192669/149‑1 (MQ=255)
                     gTCAAGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAAccc                                                                                                                   <  5:90545/149‑1 (MQ=255)
                      tCAAGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTcgcg                                                                                                                            >  2:5440/1‑139 (MQ=255)
                      tCAAGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCcg                                                                                                                  >  7:63811/1‑149 (MQ=255)
                      tCAAGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCAAAGCCGgtgccatc                                                                                                                             <  7:27210/138‑6 (MQ=255)
                      tCAAGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCAAAGCCGgtgccatc                                                                                                                             >  8:27210/1‑133 (MQ=255)
                         aGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCgcgc                                                                                                               >  6:264888/1‑149 (MQ=255)
                         aGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCgcgc                                                                                                               >  7:200869/1‑149 (MQ=255)
                             gCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTTTACCcgcgaacac                                                                                                           <  1:166661/149‑6 (MQ=255)
                              cAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACg                                                                                                          <  5:201223/149‑1 (MQ=255)
                                       gggCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCtt                                                                                                  <  8:63811/148‑1 (MQ=255)
                                        ggCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGTTAAACCCGCGCCTACGGGGGCCTtgt                                                                                                <  1:168346/149‑1 (MQ=255)
                                           ttCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTtgtgtg                                                                                             >  6:324995/1‑149 (MQ=255)
                                            tCATCTTGAATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGc                                                                                            >  8:185889/1‑149 (MQ=255)
                                              aTCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCgg                                                                                          <  3:72399/149‑1 (MQ=255)
                                              aTCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCgg                                                                                          <  6:288460/149‑1 (MQ=255)
                                                 ttGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGt                                                                                        >  8:140829/1‑148 (MQ=255)
                                                    cATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGGTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCaa                                                                                     <  3:212001/148‑1 (MQ=255)
                                                        ccGCGCCGTTTCTGGGCGGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGaa                                                                                <  6:172056/149‑1 (MQ=255)
                                                           cgcCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACAt                                                                             <  7:185889/149‑1 (MQ=255)
                                                                           ggcgggggcgggAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGtt                                                              <  8:200869/148‑1 (MQ=255)
                                                                             cgggggcgggAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTg                                                             <  4:257424/147‑1 (MQ=255)
                                                                             cgggggcgggAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTg                                                             >  3:257424/1‑147 (MQ=255)
                                                                                         ttGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAAGATGAACAAGGTCTCGTTGGAGAATTGCGCCTc                                               >  1:135427/1‑149 (MQ=255)
                                                                                           gACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCg                                              >  7:257949/1‑148 (MQ=255)
                                                                                           gACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAAACCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCg                                              >  7:285500/1‑148 (MQ=255)
                                                                                            aCCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGaa                                            <  8:257949/149‑1 (MQ=255)
                                                                                            aCCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGaa                                            >  3:13653/1‑149 (MQ=255)
                                                                                             ccTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAAc                                           <  4:13653/149‑1 (MQ=255)
                                                                                              cTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACt                                          >  5:232843/1‑149 (MQ=255)
                                                                                                       cGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTg                                 >  2:97378/1‑149 (MQ=255)
                                                                                                                  gCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGgccgc                      >  5:302580/1‑149 (MQ=255)
                                                                                                                  gCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGgccgc                      >  3:289891/1‑149 (MQ=255)
                                                                                                                    aGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGTGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGccc                    <  4:289891/149‑1 (MQ=255)
                                                                                                                    aGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGgccgc                      >  7:123596/1‑147 (MQ=255)
                                                                                                                    aGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGgccgc                      <  8:123596/147‑1 (MQ=255)
                                                                                                                            ggcgCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACgg                                                                                                         <  5:317437/56‑1 (MQ=255)
                                                                                                                            ggcgCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACgg                                                                                                         >  6:317437/1‑56 (MQ=255)
                                                                                                                              cgcgGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAg          <  5:324995/149‑1 (MQ=255)
                                                                                                                                cgGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATCGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGcc        <  7:140829/149‑1 (MQ=255)
                                                                                                                                     gcgaacACCGGCAAAGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGgccgcc                     >  6:54111/6‑131 (MQ=39)
                                                                                                                                     gcgaacACCGGCAAAGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGgccgcc                     <  5:54111/126‑1 (MQ=39)
                                                                                                                                     gATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGc    <  1:97378/148‑1 (MQ=255)
                                                                                                                                       gaacACCGGCAAAGCCGGTGTTCGCGGGTAAACCCGCGGCAACCGGGGCCTTGTTGGGGGGTGAATGGAACATGAACAAGGGCGGGTTGGAGAATTTCGGCGCGGACAGATGCGGTGGCATAGGCCGCCCGAACAGATAGCACTGccc  >  2:262151/4‑146 (MQ=255)
                                                                                                                                               gCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGccc       <  2:107611/135‑1 (MQ=255)
                                                                                                                                               gCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGccc       >  1:107611/1‑135 (MQ=255)
                                                                                                                                                   |                                                                                                                                       
GCGGCTGGCATATGAAGTGGCGTCAAGTCGCAGACCAATGGGCTTCATCTTGCATCCCGCGCCGTTTCTGGGCTGGGCGGGGGCGGGAGTTGACCTGGAAGTCCGCTTTTTTATGCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCAC  >  NC_002947/781899‑782181

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