Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,547,989 G→A 12.9% A371A (GCC→GCT rnr ← exoribonuclease R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,547,9890GA12.9% 71.6 / 4.0 31A371A (GCC→GCTrnrexoribonuclease R
Reads supporting (aligned to +/- strand):  ref base G (14/13);  new base A (2/2);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

TCCGGCACGACCTCCTTGTACTCTTCGCGCAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTACACCGAGTTACCGCGCTGCTGGGCTTCGACATC  >  NC_002947/5547843‑5548138
                                                                                                                                                  |                                                                                                                                                     
tCCGGCACGACCTCCTTGTACTCTTCGCGCAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCt                                                                                                                                                    >  1:425110/1‑150 (MQ=255)
 ccGGCACGACCTCCTTGTACTCTTCGCTCAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCtt                                                                                                                                                   <  2:65033/150‑1 (MQ=255)
      aCGACCTCCTTGTACTCATCGCGCAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTc                                                                                                                                              >  1:423053/1‑150 (MQ=255)
                      ctTCGCGCAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGa                                                                                                                              >  1:587704/1‑150 (MQ=255)
                           cgcAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATg                                                                                                                         <  2:719918/150‑1 (MQ=255)
                                         cTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGAc                                                                                                           <  1:683028/150‑1 (MQ=255)
                                                    gcgAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATg                                                                                                                                           >  2:186360/1‑107 (MQ=255)
                                                    gcgAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGGGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATg                                                                                                                                           <  1:186360/107‑1 (MQ=255)
                                                      gAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGc                                                                                              <  1:408112/150‑1 (MQ=255)
                                                           cTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGca                                                                                         <  2:423053/150‑1 (MQ=255)
                                                                             cTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATgg                                                                                                                        >  2:708181/1‑101 (MQ=255)
                                                                             cTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATgg                                                                                                                        <  1:708181/101‑1 (MQ=255)
                                                                             cTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCAGCTTTGTTCATGGTCATTTcaca                                                                                                                                >  1:791974/1‑93 (MQ=255)
                                                                             cTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCAGCTTTGTTCATGGTCATTTcaca                                                                                                                                <  2:792407/93‑1 (MQ=255)
                                                                             cTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCAGCTTTGTTCATGGTCATTTcaca                                                                                                                                <  2:791974/93‑1 (MQ=255)
                                                                             cTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCAGCTTTGTTCATGGTCATTTcaca                                                                                                                                >  1:792407/1‑93 (MQ=255)
                                                                              ttgttgTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGctct                                                                      >  2:629663/1‑150 (MQ=255)
                                                                                       gTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGc                                                             >  1:491714/1‑150 (MQ=255)
                                                                                               ggCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGcac                                                     <  2:425110/150‑1 (MQ=255)
                                                                                                     ggAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCAcgcg                                               >  2:792367/1‑150 (MQ=255)
                                                                                                         tGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCg                                           <  2:587704/150‑1 (MQ=255)
                                                                                                             tCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAg                                                                                      <  2:302135/103‑1 (MQ=255)
                                                                                                             tCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAg                                                                                      >  1:302135/1‑103 (MQ=255)
                                                                                                                cGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGt                                    >  1:305812/1‑150 (MQ=255)
                                                                                                                 gCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTa                                   >  1:186630/1‑150 (MQ=255)
                                                                                                                 gCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGATCGCTAAGTTGTTTCAGCGAGCACAGGCCGTTGGACACCTCTTCGGGGCGGATCGGCACCACGCGCTCGGGGAAGTa                                   >  1:574698/1‑150 (MQ=255)
                                                                                                                    ttCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGAc                                                                                                                             >  2:461612/1‑57 (MQ=255)
                                                                                                                    ttCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGAc                                                                                                                             <  1:461612/57‑1 (MQ=255)
                                                                                                                       gTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGTCATGCGGGCTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTACACCGa                             >  1:592281/1‑150 (MQ=255)
                                                                                                                                            cTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTACACCGAGTTACCGCGCTGCTGGGCTTc        <  1:583790/150‑1 (MQ=255)
                                                                                                                                            cTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTACACCGAGTTACCGCGCTGCTGGGCTTc        <  2:305812/150‑1 (MQ=255)
                                                                                                                                                  ggCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTACACCGAGTTACCGCGCTGCTGGGCTTCGACATc  >  1:737218/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCCGGCACGACCTCCTTGTACTCTTCGCGCAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTACACCGAGTTACCGCGCTGCTGGGCTTCGACATC  >  NC_002947/5547843‑5548138

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: