Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 71,354 Δ1 bp 9.1% coding (551/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_00294771,3540C. 9.1% 136.4 / 7.4 44coding (551/2055 nt)glySglycine‑tRNA ligase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (19/21);  new base . (2/2);  total (21/23)
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.95e-01

CCTGCTGCATGGCCAGCTCTGCGGTACGCTTGCTGATCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGCCCAGCGCATGCGCTTGGGAATCGGCAGGTCGTTCA  >  NC_002947/71208‑71501
                                                                                                                                                  |                                                                                                                                                   
ccTGCTGCATGGCCAGCTCTGCGGTACGCTTGCTGATCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCa                                                                                                                                                   >  4:313766/1‑149 (MQ=255)
       cATGGCCAGCTCTGCGGTACGCTTGCTGATCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGAt                                                                                                                                             <  6:185774/148‑1 (MQ=255)
       cATGGCCAGCTCTGCGGTACGCTTGCTGATCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGAt                                                                                                                                             <  3:261001/148‑1 (MQ=255)
       cATGGCCAGCTCTGCGGTACGCTTGCTGATCAGCGCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGAtt                                                                                                                                            <  3:160502/149‑1 (MQ=255)
                 tctGCGGTACGCTTGCTGATCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGG‑CACGGGATTCACGCCCGGcc                                                                                                                                 >  3:162711/1‑149 (MQ=255)
                 tctGCGGTACGCTTGCTGATCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGG‑CACGGGATTCACGCCCGGcc                                                                                                                                 >  3:168164/1‑149 (MQ=255)
                     cGGTACGCTTGCTGATCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGcc                                                                                                                                 <  3:150884/146‑1 (MQ=255)
                       gTACGCTTGCTGATCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTg                                                                                                                            >  7:167232/1‑149 (MQ=255)
                             ttGCTGATCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTg                                                                                                                            <  4:191471/143‑1 (MQ=255)
                             ttGCTGATCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTg                                                                                                                            >  3:191471/1‑143 (MQ=255)
                                   aTCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGG‑CACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGt                                                                                                                <  4:168164/148‑1 (MQ=255)
                                   aTCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGG‑CACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGt                                                                                                                <  4:162711/148‑1 (MQ=255)
                                       gCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATg                                                                                                                  <  3:229125/143‑1 (MQ=255)
                                       gCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATg                                                                                                                  >  4:229125/1‑143 (MQ=255)
                                         tCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGGACGTTTTCCGGATGGTGGGCGCGGTGGCCCCGGGATTCACGCCCGGCCTTCGGCGACAGGATGGTggcatc                                                                                                          >  8:135020/1‑144 (MQ=255)
                                            cggcgCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGac                                                                                                       >  5:131365/1‑149 (MQ=255)
                                            cggcgCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGac                                                                                                       <  7:121294/149‑1 (MQ=255)
                                               cgcTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGaccac                                                                                                    <  8:167232/149‑1 (MQ=255)
                                                  tCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGcc                                                                                                                                 <  2:87065/117‑1 (MQ=255)
                                                  tCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGcc                                                                                                                                 >  1:87065/1‑117 (MQ=255)
                                                      cGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGTTTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCg                                                                                             >  4:246297/1‑149 (MQ=255)
                                                          gTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCca                                                                                         <  3:313766/149‑1 (MQ=255)
                                                           tCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCcag                                                                                        >  1:52014/1‑149 (MQ=255)
                                                                                   aGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCgg                                                                >  6:197835/1‑149 (MQ=255)
                                                                                   aGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCgg                                                                >  6:284450/1‑149 (MQ=255)
                                                                                        ttCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCa                                                           >  6:171413/1‑149 (MQ=255)
                                                                                        ttCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCa                                                           >  1:329679/1‑149 (MQ=255)
                                                                                                               tgatgaGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCAtt                                                                       >  7:321730/1‑114 (MQ=255)
                                                                                                               tgatgaGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCAtt                                                                       <  8:321730/114‑1 (MQ=255)
                                                                                                                gatgaGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGTTTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGccc                                   <  3:246297/149‑1 (MQ=255)
                                                                                                                    aGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGCCCAgc                                <  5:171413/148‑1 (MQ=255)
                                                                                                                      cACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATg                                                                                                                  <  1:150359/64‑1 (MQ=255)
                                                                                                                      cACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATg                                                                                                                  >  2:150359/1‑64 (MQ=255)
                                                                                                                      cACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGCCCAGCGCa                             >  6:232343/1‑149 (MQ=255)
                                                                                                                         gTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGCCCAGCGCATgc                          >  4:249246/1‑149 (MQ=255)
                                                                                                                                 gATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGCCCAGCGCATGCGCTTg                     >  7:97319/1‑146 (MQ=255)
                                                                                                                                 gATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGCCCAGCGCATGCGCTTg                     <  8:97319/146‑1 (MQ=255)
                                                                                                                                    ggtggAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGt                                                                   >  6:186100/1‑97 (MQ=255)
                                                                                                                                    ggtggAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGt                                                                   <  5:186100/97‑1 (MQ=255)
                                                                                                                                      tggAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGCCCAGCGCATGCGCTTGGGAATCgg             <  6:131365/149‑1 (MQ=255)
                                                                                                                                      tggAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGCCCAGCGCATGCGCTTGGGAATCgg             <  5:197835/149‑1 (MQ=255)
                                                                                                                                         aaGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGCCCAGCGCATGCGCTTGGGAATCGGCAg          <  2:329679/149‑1 (MQ=255)
                                                                                                                                               tGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGCCCAGCGCATGCGCTTGGGAATCGGCAGGTCGtt    <  5:284450/149‑1 (MQ=255)
                                                                                                                                                 gCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGCCCAGCGCATGCGCTTGGGAATCGGCAGGTCGTTCa  <  3:249246/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCTGCTGCATGGCCAGCTCTGCGGTACGCTTGCTGATCAGCTCGCGGCGCTCGGCGAAGTCGGCCAGGACGTAGGCCTTGCGCAGGTCTTCGACGTAGTTGGCCGGGGTGGTGATGAGCACGTTTTCCGGATGGTGGAAGCGGTGGCCACGGGATTCACGCCCGGCCTTCTGCGACAGGATGGTGCAATCGACCACCTGATCGCCCAGCAGCATCACCAGCCATTGGGTCGGGCGCACGAACTCTTCACGGCTGGCCGCCCAGCGCATGCGCTTGGGAATCGGCAGGTCGTTCA  >  NC_002947/71208‑71501

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