Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 556,978 C→T 8.1% R9R (CGC→CGT rpsN → 30S ribosomal protein S14

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947556,9780CT 8.1% 133.2 / 3.2 50R9R (CGC→CGTrpsN30S ribosomal protein S14
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base T (2/2);  total (27/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 = 3.95e-01

GATCGATGCTCTGCGCGGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTCG  >  NC_002947/556833‑557125
                                                                                                                                                 |                                                                                                                                                   
gATCGATGCTCTGCGCGGTTTGGACATCACCCTGACCACCAGTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAg                                                                                                                                                  <  5:233561/149‑1 (MQ=255)
             cgcgGTTTGGACATCACGCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTc                                                                                                                                        <  3:159843/146‑1 (MQ=255)
             cgcgGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAg                                                                                                                                      <  1:269559/148‑1 (MQ=255)
                       aCATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGTGAGCTGACGCTTCAGCTCGCGGtttt                                                                                                                           <  2:80669/149‑4 (MQ=255)
                          tCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCaaa                                                                                                                        <  4:254188/149‑1 (MQ=255)
                            aCCCTGACCACCACTGCTCTTTCTGACGACGAAGGCCGCGCTATGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTTGAGTAGGAAATTGGCCAAGAATATCATGTAAAACCGCGATCTTAAGCGTCAGCTCACGGTAGAAAAGt                                                                                                                      <  2:240117/149‑1 (MQ=255)
                            aCCCTGACCACCACTGCTCGTTAGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACAGCGAGATGAAGCGTCAGCTCACGGTAGCaaact                                                                                                                      >  1:240117/1‑147 (MQ=255)
                                  accaccACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCa                                                                                                                                       <  5:221756/126‑1 (MQ=255)
                                  accaccACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCa                                                                                                                                       >  6:221756/1‑126 (MQ=255)
                                   ccaccaCTGCTCGTTCTGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTa                                                                                                                <  7:216923/148‑1 (MQ=255)
                                    caccacTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCt                                                                                                                 >  1:171224/1‑146 (MQ=255)
                                         cTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCg                                                                                                         <  3:209480/149‑1 (MQ=255)
                                               gTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTgct                                                                                                     >  1:1420/1‑147 (MQ=255)
                                                          gaAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGa                                                                                        >  3:193065/1‑149 (MQ=255)
                                                               ccGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATc                                                                                   >  6:134298/1‑149 (MQ=255)
                                                                cgcgCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGATGAAGCGTCAGCTGACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGACCTGAAAGCGACCATCg                                                                                  >  7:182594/1‑149 (MQ=255)
                                                                    ctctGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCaa                                                                              <  4:193065/149‑1 (MQ=255)
                                                                           gcGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCCGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAAc                                                                       >  4:129690/1‑149 (MQ=255)
                                                                            cGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTGACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTAACCCTGAACg                                                                      >  8:201137/1‑149 (MQ=255)
                                                                                  ttCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCtctc                                                                >  3:128502/1‑149 (MQ=255)
                                                                                           ccGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGt                                                                                                        >  7:173882/1‑100 (MQ=255)
                                                                                           ccGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGt                                                                                                        <  8:173882/100‑1 (MQ=255)
                                                                                            cGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCg                                                      >  7:7637/1‑149 (MQ=255)
                                                                                            cGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCg                                                      >  8:44926/1‑149 (MQ=255)
                                                                                              ttCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGtt                                                    >  2:156702/1‑149 (MQ=255)
                                                                                                ccGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGttt                                                   <  1:156702/148‑1 (MQ=255)
                                                                                                  gCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGc                                                <  4:128502/149‑1 (MQ=255)
                                                                                                     aCTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGAGAGCTGAAGCGTCAGATCACGGTAGCAAAGTTCGCTAAGAAGATTGCTGAGCTGAAAGAGACCATCCTCAACATGAACGCCTCTCCAGAAGAGCGTTTCGATg                                              >  7:224433/1‑148 (MQ=255)
                                                                                                                tAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGctc                                  >  5:201237/1‑149 (MQ=255)
                                                                                                                  ggAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGc                                                                     <  8:100099/112‑1 (MQ=255)
                                                                                                                  ggAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGc                                                                     >  7:100099/1‑112 (MQ=255)
                                                                                                                  ggAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGc                                             >  3:226231/1‑136 (MQ=255)
                                                                                                                  ggAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGc                                             <  4:226231/136‑1 (MQ=255)
                                                                                                                  ggAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGctct                                 <  2:1420/148‑1 (MQ=255)
                                                                                                                  ggAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTg                                >  8:245873/1‑149 (MQ=255)
                                                                                                                    aaaaTGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCa                              <  2:171224/149‑1 (MQ=255)
                                                                                                                      aaTGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTg                                        <  6:199851/137‑1 (MQ=255)
                                                                                                                      aaTGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTg                                        >  5:199851/1‑137 (MQ=255)
                                                                                                                          tccAAGAAGAGCATGAAAAACCTCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAagca                        <  5:134298/148‑1 (MQ=255)
                                                                                                                            cAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAagcagc                      >  1:75132/1‑149 (MQ=255)
                                                                                                                               taagaGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGCACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAACCAGCCAc                   <  7:44926/148‑1 (MQ=255)
                                                                                                                                aagaGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCAcg                  >  7:66161/1‑149 (MQ=255)
                                                                                                                                  gagCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCg                                                 >  6:85229/1‑116 (MQ=255)
                                                                                                                                  gagCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCg                                                 <  5:85229/116‑1 (MQ=255)
                                                                                                                                     cATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACg             >  8:183823/1‑149 (MQ=255)
                                                                                                                                              ccGTGAGCTGACGCTTCAGCTCGCGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGAGCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCACTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCt    >  2:65144/1‑149 (MQ=255)
                                                                                                                                              ccGTGAGCTGACGCTTCAGCTCGCGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGttt                                                   <  1:110458/102‑1 (MQ=39)
                                                                                                                                              ccGTGAGCTGACGCTTCAGCTCGCGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGttt                                                   >  2:110458/1‑102 (MQ=39)
                                                                                                                                              ccGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAg         <  2:6698/144‑1 (MQ=255)
                                                                                                                                              ccGCGAGCTGAAGCGTCAGCTCACGGGAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAg         >  1:6698/1‑144 (MQ=255)
                                                                                                                                                gcgAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTcg  >  7:50077/1‑149 (MQ=255)
                                                                                                                                                gcgAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTcg  <  7:245873/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GATCGATGCTCTGCGCGGTTTGGACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTCG  >  NC_002947/556833‑557125

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