Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,667,961 C→T 14.8% A307A (GCC→GCT ffh → signal recognition particle GTPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,667,9610CT14.8% 56.8 / 3.8 27A307A (GCC→GCTffhsignal recognition particle GTPase subunit
Reads supporting (aligned to +/- strand):  ref base C (9/14);  new base T (2/2);  total (11/16)
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.14e-01

GTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGACAAGCTGC  >  NC_002947/1667824‑1668103
                                                                                                                                         |                                                                                                                                              
gTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCGCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATc                                                                                                                                    <  1:179834/150‑1 (MQ=255)
       cAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaagg                                                                                                                             <  1:239599/150‑1 (MQ=255)
               ggCAACCCGATCAAGTTCATCGGTAGGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGCCGACGTGCTCAGGCTGAGCGAGCAGGCCGAGCAGACCATCGACaaggccaaggcc                                                                                                                     <  2:12033/150‑1 (MQ=255)
                    aCCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACaa                                                                                                                <  2:433019/150‑1 (MQ=255)
                       gATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGc                                                                                                              <  1:194901/149‑1 (MQ=255)
                       gATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGc                                                                                                              >  2:194901/1‑149 (MQ=255)
                           aaGTTCATCGGTATGGGAGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGCTCCGGCAGGCCGTGGTGCACAGCGACAAGGGCGATGCCGACAAGCTgggc                                                                                                         >  2:313562/1‑148 (MQ=255)
                             gTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCaa                                                                                                       <  1:403777/150‑1 (MQ=255)
                                   gtttgccggTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGCGCCCACCCAGCACACCATTGACAAAGACAACGCTTACACGCGGGCCAAGAAACt                                                                                                 >  2:67677/8‑150 (MQ=255)
                                    ggTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGc                                                                                                          >  2:74508/1‑140 (MQ=255)
                                    ggTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGc                                                                                                          <  1:74508/140‑1 (MQ=255)
                                    ggTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTg                                                                                                <  1:259326/150‑1 (MQ=255)
                                               gaagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTgaa                                                                                              <  1:448455/141‑1 (MQ=255)
                                               gaagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTgaa                                                                                              >  2:448455/1‑141 (MQ=255)
                                                aagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggcaag                                                                                    <  2:478295/150‑1 (MQ=255)
                                                      gAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTc                                                                              <  2:533560/150‑1 (MQ=255)
                                                                         aCCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCcgcg                                                           >  1:229036/1‑150 (MQ=255)
                                                                             cGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCATGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGCAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCa                                                       <  1:245753/150‑1 (MQ=255)
                                                                                          tCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGCGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATg                                          <  1:313562/150‑1 (MQ=255)
                                                                                                             gCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTcg                       >  1:129630/1‑150 (MQ=255)
                                                                                                                 tctGCTCGGCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTcggcg                    <  1:523297/147‑1 (MQ=255)
                                                                                                                 tctGCTCGGCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTcggcg                    >  2:523297/3‑149 (MQ=255)
                                                                                                                      ttggCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagg                                                                                         <  1:267892/72‑1 (MQ=37)
                                                                                                                      ttggCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagg                                                                                         >  2:267892/4‑75 (MQ=37)
                                                                                                                          ccTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGa          >  1:474424/1‑150 (MQ=255)
                                                                                                                            tGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGACa        >  2:205419/1‑150 (MQ=255)
                                                                                                                                  aGCAGGCCGGGCGGACCATCGACAAGGCCAAGGCCGACAATCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGACAAGCTGc  <  2:34803/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                              
GTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGACAAGCTGC  >  NC_002947/1667824‑1668103

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: