Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 70,755:1 +G 8.3% coding (1150/2055 nt) glyS ← glycine‑tRNA ligase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294770,7551.G 8.3% 183.2 / 9.4 48coding (1150/2055 nt)glySglycine‑tRNA ligase subunit beta
Reads supporting (aligned to +/- strand):  ref base . (24/20);  new base G (2/2);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCGGCATGTACTGCTCGTTCAGGGCCAGGGCTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCACGTTCTTGAGGCGCTCGTTGAAGGTTTCCAGTGG  >  NC_002947/70608‑70899
                                                                                                                                                    |                                                                                                                                                
cgcgcgGCATGTACTGCTCGTTCAGGGCCAGGGCTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCa                                                                                                                                                 >  2:729306/1‑149 (MQ=255)
cgcgcgGCATGTACTGCTCGTTCAGGGCCAGGGCTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCa                                                                                                                                                 <  1:729306/149‑1 (MQ=255)
              tGCTCGTTCAGGGCCAGGGCTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCGACGCCCAGCGCGCTg                                                                                                                                   >  1:453612/1‑150 (MQ=255)
                cTCGTTCAGGGCCAGGGCTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGc                                                                                                                                >  2:829779/1‑150 (MQ=255)
                cTCGTTCAGGGCCAGGGCTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGc                                                                                                                                >  2:829976/1‑150 (MQ=255)
                cTCGTTCAGGGCCAGGGCTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGc                                                                                                                                >  1:408696/1‑150 (MQ=255)
                    ttCAGGTCCAGGACTACGTCTTCCGGTTCACCGTCGGTGGGCGTGTAGTAGTAAACGGCAACGCCCTGCGTCTCAGGGAATTCGCGGACCATCGCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGGCAGCGCGCTGGtcgtcg                                                                                                                            >  2:961888/1‑150 (MQ=255)
                     tCAGGGCCAGGGCTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCGACGCCCAGCGCGCTGGCCGtgc                                                                                                                            <  2:453612/150‑3 (MQ=255)
                      cagtgccaggtcTACGTCTCCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATGCGCCGACCATCTCGGTGGCCANGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTcgcc                                                                                                                          <  1:623187/150‑1 (MQ=255)
                        gggccagggcTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTcgccgc                                                                                                                        >  2:1106757/1‑150 (MQ=255)
                          gccagggcTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGc                                                                                                                                >  1:540187/1‑140 (MQ=255)
                          gccagggcTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGc                                                                                                                                <  2:540187/140‑1 (MQ=255)
                                gcTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGc                                                                                                                >  2:912754/1‑150 (MQ=255)
                                                 aCCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGACCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGACCAATGAAgccggc                                                                                               >  2:431498/1‑150 (MQ=255)
                                                              gcgTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCGACGCCCAGCGCGCTGGCCGt                                                                                                                              >  1:245386/1‑107 (MQ=255)
                                                              gcgTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCGACGCCCAGCGCGCTGGCCGt                                                                                                                              <  2:245386/107‑1 (MQ=255)
                                                              gcgTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAaca                                                                                  >  2:211370/1‑150 (MQ=255)
                                                                gtagtagtaACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGGCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACg                                                                                >  2:362104/1‑150 (MQ=255)
                                                                        aaCCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCt                                                                        >  1:902945/1‑150 (MQ=255)
                                                                         aCCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCt                                                                                                                                    <  1:601561/89‑1 (MQ=255)
                                                                         aCCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCt                                                                                                                                    >  2:601561/1‑89 (MQ=255)
                                                                            ggCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCCCTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAACAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTc                                                                    <  1:1149976/150‑1 (MQ=255)
                                                                            ggCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTc                                                                    <  1:829779/150‑1 (MQ=255)
                                                                            ggCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTc                                                                    <  1:912754/150‑1 (MQ=255)
                                                                            ggCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTc                                                                    <  1:829976/150‑1 (MQ=255)
                                                                                   cccTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTAAACg                                                             <  2:1121672/150‑1 (MQ=255)
                                                                                    ccTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACgg                                                            >  1:820017/1‑150 (MQ=255)
                                                                                      tGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGa                                                                                                             <  2:115869/99‑1 (MQ=255)
                                                                                      tGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGa                                                                                                             >  1:115869/1‑99 (MQ=255)
                                                                                          tttCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGa                                                      <  2:902945/150‑1 (MQ=255)
                                                                                           ttCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCg                                                                            <  2:1067294/127‑1 (MQ=255)
                                                                                           ttCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCg                                                                            >  1:1067294/1‑127 (MQ=255)
                                                                                              aGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCg                                                                            <  1:959379/124‑1 (MQ=255)
                                                                                              aGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCg                                                                            >  2:959379/1‑124 (MQ=255)
                                                                                                  aaTTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGcc                                              >  2:969578/1‑150 (MQ=255)
                                                                                                    ttCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTg                                            <  2:408696/150‑1 (MQ=255)
                                                                                                       gCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGaa                                         <  2:486528/150‑1 (MQ=255)
                                                                                                            ccATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAAcacca                                    >  2:174021/1‑150 (MQ=255)
                                                                                                                   ggTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGACGGCGGCCCGCAAGGCCAAGCGTTATTCCTGGTCGGGGACGGTGCCGTGCTGAGCCTGGAACACCACGTTCtt                             >  1:871799/1‑150 (MQ=255)
                                                                                                                    gTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCACGTTCTTg                            <  1:211370/150‑1 (MQ=255)
                                                                                                                       gCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCGCGTTCTTGAgg                         <  1:969578/150‑1 (MQ=255)
                                                                                                                           ggTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCACGTTCTTGAGGCGCt                     >  2:725232/1‑150 (MQ=255)
                                                                                                                            gTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCACGTTCTTGAGGCGCTc                    >  2:448083/1‑150 (MQ=255)
                                                                                                                                   ttCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCACGTTCTTGAGGCGCTCGTTGAAg             >  1:365243/1‑150 (MQ=255)
                                                                                                                                      gACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCACGTTCTTGAGGCGCTCGTTGaa              <  2:1205247/146‑1 (MQ=255)
                                                                                                                                      gACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCACGTTCTTGAGGCGCTCGTTGaa              >  1:1205247/1‑146 (MQ=255)
                                                                                                                                      caCAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCAAATTTTCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCACGTTCTTGAGGCGCTCGTTGAAGGtt          <  1:316905/149‑1 (MQ=255)
                                                                                                                                       aCAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCACGTTCTTGAGGCGCTCGTTGAAGGttt         >  1:786690/1‑150 (MQ=255)
                                                                                                                                            agGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCACGTTCTTGAGGCGCTCGTTGAAGGTTTCCAGt    <  1:638536/150‑1 (MQ=255)
                                                                                                                                              gCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCACGTTCTTGAGGCGCTCGTTGAAGGTTTCCAGTgg  <  1:1106757/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGCGCGGCATGTACTGCTCGTTCAGGGCCAGGGCTACGTCTTCCGGCTCACCGTCGTTGAGCGCGTAGTAGTAACCGGCAACACCCTGCATTTCAGGGAATTCGCCGACCATCTCGGTGGCCAGGTCGCACTTCGACAGCAGGCCGGCACGGCCAGCGCGCTGGGCGTCGCCGCCGATCAGCGGAGCAATGAAGGCGGCCAGCCTGGAAACACGTTCGGCCTTGTCGTACACGGTGCCGAGCTGAGCCTGGAACACCACGTTCTTGAGGCGCTCGTTGAAGGTTTCCAGTGG  >  NC_002947/70608‑70899

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