Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,827,748 T→C 10.5% L187P (CTG→CCG)  recA → recombinase RecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,827,7480TC10.5% 93.3 / 3.5 38L187P (CTG→CCG) recArecombinase RecA
Reads supporting (aligned to +/- strand):  ref base T (16/18);  new base C (2/2);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

GACGTGATCATTGTCGACTCCGTTGCCGCGCTGGTACCGAAGGCCGAGATCGAAGGCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTACCGGTGGTAACGCCCTGAAGTTCTATGCTTCGGTCCGTCTGGACATCCGCCGTACCGGCGCGGTCAAGGAAGGCGACGAAG  >  NC_002947/1827600‑1827894
                                                                                                                                                    |                                                                                                                                                  
gACGTGATCATTGTCGACTCCGTTGCCGCGCTGGTACCGAAGGCCGAGATCGAAGGCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTg                                                                                                                                                   >  2:604234/1‑150 (MQ=255)
      aTCATTGTCGACTCCGTTGCCGGGCTGGTACCGAAGGCCGAGATCGAAGGCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATc                                                                                                                                             <  2:139320/150‑1 (MQ=255)
            gTCGACTCCGTTGCCGCGCTGGTACCGAAGGCCGAGATCGAAGGCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTtcatc                                                                                                                                       <  2:86435/150‑1 (MQ=255)
               gACTCCGTTGCCGCGCTGGTACCGAAGGCCGAGATCGAAGGCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAAc                                                                                                                                    <  1:604234/150‑1 (MQ=255)
                  tCCGTTGCCGCGCTGGTACCGAAGGCCGAGATCGAAGGCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAg                                                                                                                                 >  1:305381/1‑150 (MQ=255)
                          cgcgCTGGTACCGAAGGCCGAGATCGAAGGCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTAt                                                                                                                         >  2:1005442/1‑150 (MQ=255)
                             gcTGGTACCGAAGGCCGAGATCGAAGGCGAGATGGGTGACATGCATGTGGGCATGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGaa                                                                                                                      <  1:660952/150‑1 (MQ=255)
                             gcTGGTACCGAAGGCCGAGATCGAAGGCGAGATGGGTGACATGCATGTGGGCATGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGaa                                                                                                                      <  1:660933/150‑1 (MQ=255)
                                   aCCGAAGGCCGAGATCGAAGGCGAGATGGGTGGCATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCgg                                                                                                                >  2:659597/1‑150 (MQ=255)
                                     cGAAGGCCGAGATCGAAGGCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGgtg                                                                                                              >  2:420783/1‑150 (MQ=255)
                                      gAAGGCCGAGATCGAAGGCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGGAACATCAAAAACGCGAACTGCCTGGTCATCTTCATCAACCAGCTCCGTATGAAGATCGgtgt                                                                                                             >  1:234365/1‑150 (MQ=255)
                                        aGGCCGAGATCGAAGGGGAGATGGGTGACATGCATGTGGGACTACAGTCCCGCCTGATGTCCCAGGCACTGGCTAACTTCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGa                                                                                                           <  1:247764/150‑1 (MQ=255)
                                              agaTCGAAGGCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCg                                                                                                     >  1:940458/1‑150 (MQ=255)
                                                   gAAGGCGAGTTGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTAGGCAGc                                                                                                <  2:305381/150‑1 (MQ=255)
                                                   gAAGGCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGc                                                                                                >  1:541243/1‑150 (MQ=255)
                                                      ggCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCg                                                                                             >  2:90958/1‑150 (MQ=255)
                                                               ggTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACt                                                                                    >  1:866371/1‑150 (MQ=255)
                                                                gTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTa                                                                                   <  1:420783/150‑1 (MQ=255)
                                                                                 cTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGACTGGTCATCTTCATCAACCAGAACCGACTGAAGATTAGCGTGATGTTCGGCAGCCAGGAAACCGCTACCGGTGGAAACGCCCTg                                                                  >  2:346497/1‑150 (MQ=255)
                                                                                       gCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTACCGGTGGTAACGCCCTGAAGTTc                                                            <  1:122022/150‑1 (MQ=255)
                                                                                         ccGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCa                                                                                                                                  <  1:512668/78‑1 (MQ=255)
                                                                                         ccGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCa                                                                                                                                  >  2:512668/1‑78 (MQ=255)
                                                                                         ccGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTACCGGTGGTAACGCCCTGAAGTTCTa                                                          <  2:866371/150‑1 (MQ=255)
                                                                                           gCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGAt                                                                                                                               >  2:549073/1‑79 (MQ=255)
                                                                                           gCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGAt                                                                                                                               <  1:549073/79‑1 (MQ=255)
                                                                                                  gTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTACCGGTGGTAACGCCCTGAAGTTCTATGCTTCGGt                                                 <  2:940458/150‑1 (MQ=255)
                                                                                                      cAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCCGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCa                                                                                                  >  1:896192/1‑97 (MQ=255)
                                                                                                      cAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCCGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCa                                                                                                  >  1:900970/1‑97 (MQ=255)
                                                                                                      cAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCCGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCa                                                                                                  <  2:896192/97‑1 (MQ=255)
                                                                                                      cAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCCGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCa                                                                                                  <  2:900970/97‑1 (MQ=255)
                                                                                                             tGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTa                                                                                   <  2:508721/105‑1 (MQ=255)
                                                                                                             tGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTa                                                                                   >  1:508721/1‑105 (MQ=255)
                                                                                                                gTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTACCGGTGGTAACGCCCTGAAGTTCTa                                                          >  2:219032/1‑127 (MQ=255)
                                                                                                                gTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTACCGGTGGTAACGCCCTGAAGTTCTa                                                          <  1:219032/127‑1 (MQ=255)
                                                                                                                        aCCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTACCGGTGGTAACGCCCTGAAGTTCTATGCTTCGGTCCGTCTGGACATCCGCCGTAcc                           >  2:1033645/1‑150 (MQ=255)
                                                                                                                         ccGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTACCGGTGGTAACGCCCTGAAGTTCTATGCTTCGGTCCGTCTGGACATCCGCCGTACCg                          <  1:90958/150‑1 (MQ=255)
                                                                                                                                         cgcgAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTACCGGTGGTAACGCCCTGAAGTTCTATGCTTCGGTCCGTCTGGACATCCGCCGTACCGGCGCGGTCaaggaagg          <  1:1033645/150‑1 (MQ=255)
                                                                                                                                                 gCCTGGTCATCTTCATCAACCAGATCCGTATGAGGATCGGTGTGGTGTTCGGCAGCCCGGAAACCACTACCGGTGGTAACGCCCTGAAGTTCTATGCTTCGGTCCGTCTGGACATCCGCCGTACCGTCGCGGCCAAGGGAGGCGACGAAg  <  2:910825/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GACGTGATCATTGTCGACTCCGTTGCCGCGCTGGTACCGAAGGCCGAGATCGAAGGCGAGATGGGTGACATGCATGTGGGCCTGCAGGCCCGCCTGATGTCCCAGGCACTGCGTAAGATCACCGGTAACATCAAGAACGCGAACTGCCTGGTCATCTTCATCAACCAGATCCGTATGAAGATCGGTGTGATGTTCGGCAGCCCGGAAACCACTACCGGTGGTAACGCCCTGAAGTTCTATGCTTCGGTCCGTCTGGACATCCGCCGTACCGGCGCGGTCAAGGAAGGCGACGAAG  >  NC_002947/1827600‑1827894

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: