Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,813,478 C→A 7.7% A221A (GCC→GCA)  pheS → phenylalanine‑‑tRNA ligase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,813,4780CA 7.7% 161.8 / 3.9 52A221A (GCC→GCA) pheSphenylalanine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (20/28);  new base A (2/2);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

GTGCAGGTGCGGACCATGGAGTCCACCCAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACATTCAGTGCGTGATGTGTTCCGGCAACGGCTG  >  NC_002947/2813332‑2813621
                                                                                                                                                  |                                                                                                                                                
gTGCAGGTGCGGACCATGGAGTCCACCCAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGa                                                                                                                                                <  1:146732/149‑1 (MQ=255)
    aGGTGCGGACCATGGAGTCCACCCAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCt                                                                                                                                             >  1:37407/1‑148 (MQ=255)
       tGCGGACCATGGAGTCCACCCAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAg                                                                                                                                         >  3:77129/1‑149 (MQ=255)
           gACCATGGAGTCCACCCAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCa                                                                                                                                     >  4:12080/1‑149 (MQ=255)
              cATGGAGTCCACCAAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTCCTGATAGATCGCGATATCAACTTCGCCGACCTCAAGGGCCCCa                                                                                                                                  >  3:296019/1‑149 (MQ=255)
                tGGAGTCCACCCAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATc                                                                                                                                <  5:316999/149‑1 (MQ=255)
                tGGAGTCCACCCAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATc                                                                                                                                >  4:348956/1‑149 (MQ=255)
                  gAGTCCACGCAGCCGACCATTCGCATTGTGTGCCCTGGCCGCGTATACCGCTGCGACTCGGATATCACCTACTCGCCGATGTTCCAGGATTTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCAACATCg                                                                                                                               <  5:252394/148‑1 (MQ=255)
                   aGTCCACCCAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCgaa                                                                                                                             <  4:77129/149‑1 (MQ=255)
                   aGTCCACCCAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCgaa                                                                                                                             <  8:333750/149‑1 (MQ=255)
                       cACCCAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAAt                                                                                                                         <  3:348956/149‑1 (MQ=255)
                           cAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCt                                                                                                                     <  2:37407/149‑1 (MQ=255)
                             gccgccCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTgc                                                                                                                   <  3:34276/149‑1 (MQ=255)
                               cgccCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTgcgc                                                                                                                 <  1:67487/149‑1 (MQ=255)
                                        gCATTGTCTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTAttcttc                                                                                                        <  4:75642/149‑1 (MQ=255)
                                          aTTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAg                                                                                                                                         <  5:135442/114‑1 (MQ=255)
                                          aTTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAg                                                                                                                                         >  6:135442/1‑114 (MQ=255)
                                          aTTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATACGCGCAGGAATTCTTC‑ga                                                                                                     <  8:271925/149‑1 (MQ=255)
                                          aTTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGACACCATCGAAGAATACGCGCAGGAATTCTTC‑ga                                                                                                     <  3:91075/149‑1 (MQ=255)
                                           ttGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATACGCGCAGGAATTCTTC‑gat                                                                                                    <  3:12080/149‑2 (MQ=255)
                                                  cccGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑Gagaaagaa                                                                                              <  6:299963/148‑1 (MQ=255)
                                                                   gCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTTATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGc                                                                            <  3:288607/149‑1 (MQ=255)
                                                                   gCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTCTGAGAAAGAACTGGCGGTCCGTTTCCg                                                                             <  1:3033/149‑1 (MQ=255)
                                                                     tGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGcc                                                                           >  6:311714/1‑148 (MQ=255)
                                                                                        aCTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGAc                                                                                                                                               <  6:283835/62‑1 (MQ=255)
                                                                                        aCTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGAc                                                                                                                                               >  5:283835/1‑62 (MQ=255)
                                                                                         cTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCg                                                      <  7:30610/149‑1 (MQ=255)
                                                                                             ccGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCAATCGTTCTTCCCCTTCACCGAgac                                                  >  4:30097/1‑147 (MQ=255)
                                                                                                    tCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTgcgc                                                                                                                 <  6:172559/80‑1 (MQ=255)
                                                                                                    tCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTgcgc                                                                                                                 >  5:172559/1‑80 (MQ=255)
                                                                                                    tCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTg                                             >  7:88869/1‑147 (MQ=255)
                                                                                                      caccaGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAg                                          >  5:307680/1‑148 (MQ=255)
                                                                                                      caccaGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAg                                          <  8:88869/148‑1 (MQ=255)
                                                                                                       accaGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTcc                                                                              >  5:244919/1‑111 (MQ=255)
                                                                                                       accaGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTcc                                                                              <  6:244919/111‑1 (MQ=255)
                                                                                                            gTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGa                                   <  5:356553/149‑1 (MQ=255)
                                                                                                               gAAGGCCTGCTGATCGATGGCGATATCAACTTCGCAGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACAt                                >  6:221644/1‑149 (MQ=255)
                                                                                                               gAAGGCCTGCTGATCGATAGCGATATCAACTTCGCAGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACAt                                >  6:226809/1‑149 (MQ=255)
                                                                                                                    ccTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACATTCAGt                           >  4:333372/1‑149 (MQ=255)
                                                                                                                        ctgATTGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCACATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCTACATTCAGTGCGt                       >  6:78524/1‑149 (MQ=255)
                                                                                                                        ctgATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACATTCAGTGCGt                       <  5:311714/149‑1 (MQ=255)
                                                                                                                           atcgatcgCGATATCAACTTCGCAGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACATTCAGTGCGTGAt                    <  5:221644/149‑1 (MQ=255)
                                                                                                                           atcgatcgCGATATCAACTTCGCAGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACATTCAGTGCGTGAt                    <  5:226809/149‑1 (MQ=255)
                                                                                                                                  gcgATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACATTCAGTGCGTGATGTGTTcc             <  1:148095/149‑1 (MQ=255)
                                                                                                                                   cgATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACAtt                               >  2:68870/1‑130 (MQ=255)
                                                                                                                                   cgATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACAtt                               <  1:68870/130‑1 (MQ=255)
                                                                                                                                    gATATCATCTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCTCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCGGTCTGCAGAAGTCGACATTCAGTGCGTGATGTGTTCCgg           <  3:30097/149‑1 (MQ=255)
                                                                                                                                    gATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACATTCAGTGCGTGATGTGTTCCgg           >  6:79136/1‑149 (MQ=255)
                                                                                                                                      tatCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACATTCAGTGCGTGATGTGTTCCGGc          >  8:97954/1‑148 (MQ=255)
                                                                                                                                          aaCTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACATTCAGTGCGTGATGTGTTCCGGCAACg      <  7:97954/148‑1 (MQ=255)
                                                                                                                                            cTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACATTCAGTGCGTGATGTGTTCCGGCAACGGCt   >  7:293591/1‑149 (MQ=255)
                                                                                                                                             ttCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTACCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACATTCAGTGCGTGATGTGTTCCGGCAACGGCTg  >  4:49799/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTGCAGGTGCGGACCATGGAGTCCACCCAGCCGCCCATTCGCATTGTGTGCCCGGGCCGCGTATACCGCTGCGACTCGGATATCACCCACTCGCCGATGTTCCACCAGGTCGAAGGCCTGCTGATCGATCGCGATATCAACTTCGCCGACCTCAAGGGCACCATCGAAGAATTCCTGCGCGTATTCTTC‑GAGAAAGAACTGGCGGTCCGTTTCCGCCCATCGTTCTTCCCCTTCACCGAGCCGTCTGCAGAAGTCGACATTCAGTGCGTGATGTGTTCCGGCAACGGCTG  >  NC_002947/2813332‑2813621

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: ‑