Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,185,604 C→G 16.6% L446L (CTC→CTG PP_3681 → putative Helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,185,6040CG16.6% 100.5 / 10.5 41L446L (CTC→CTGPP_3681putative Helicase
Reads supporting (aligned to +/- strand):  ref base C (16/18);  new base G (4/3);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTATGAGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTGGG  >  NC_002947/4185459‑4185743
                                                                                                                                                 |                                                                                                                                           
cTTATGAGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTcagca                                                                                                                                         <  1:1351830/150‑1 (MQ=255)
             gtggAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGAc                                                                                                                               <  2:454055/147‑1 (MQ=255)
             gtggAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGAc                                                                                                                               >  1:454055/1‑147 (MQ=255)
                 aaGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTccc                                                                                                                        >  1:1102835/1‑150 (MQ=255)
                  aGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCg                                                                                                                       >  2:499171/1‑150 (MQ=255)
                   gCTCACAAGCAGCTTCATCCGTGATGTAGGGTATCCAATCCTGACTGCGCTCTCAGGTGGCTACGAGTCGAAACGACGCTGGCGTGAGGAGGGTCGTAAGGGCAGCGATAGACGTCGAGCGCAGCTCACCGGATAAACGCCCATTccctg                                                                                                                      >  2:162528/1‑148 (MQ=255)
                            cagcTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGCTCGAGATGGAGCGCTGCTCATCAGCGAAAAGACCATTCCCGGTGCCCTTGc                                                                                                             <  1:743553/150‑1 (MQ=255)
                                  cATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCt                                                                                                       <  1:179239/150‑1 (MQ=255)
                                   aTCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCtt                                                                                                      >  2:363893/1‑150 (MQ=255)
                                    tCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCttt                                                                                                     >  2:579279/1‑150 (MQ=255)
                                      cGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCa                                                                                                   <  2:1102835/150‑1 (MQ=255)
                                      cGTGATGTAGATTCTCCCATCCTGACTGCGCTCTCAGGTGACTACGCGTCGACCAAACGCAGTTACCAATGCTTTCGTGCGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGATGCCCTTGCCCGGCTTTCa                                                                                                   <  1:878621/150‑1 (MQ=255)
                                        tGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGGGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGa                                                                                                 <  1:882018/150‑1 (MQ=255)
                                               ggTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTccc                                                                                                                        <  1:546761/120‑1 (MQ=255)
                                               ggTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTccc                                                                                                                        >  2:546761/1‑120 (MQ=255)
                                                       aaTCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTcagcag                                                                                                                                        <  2:875575/96‑1 (MQ=255)
                                                       aaTCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTcagcag                                                                                                                                        >  1:875575/1‑96 (MQ=255)
                                                       aaTCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACg                                                                                  >  2:1017241/1‑150 (MQ=255)
                                                         tCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCAtt                                                                                                                           >  1:800952/1‑107 (MQ=255)
                                                         tCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCAtt                                                                                                                           <  2:800952/107‑1 (MQ=255)
                                                               cTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTca                                                                          <  1:579279/150‑1 (MQ=255)
                                                                     tctcAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTGGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATcc                                                                    >  1:1357944/1‑150 (MQ=255)
                                                                         aGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCa                                                                                           <  2:1266178/123‑1 (MQ=255)
                                                                         aGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCa                                                                                           >  1:1266178/1‑123 (MQ=255)
                                                                           gTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTGAGCAGCgc                                                                                                                                     <  2:1274158/79‑2 (MQ=255)
                                                                           gTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTGAGCAGCgc                                                                                                                                     >  1:1274158/1‑78 (MQ=255)
                                                                           gTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACaaaa                                                              <  1:1017241/150‑1 (MQ=255)
                                                                                  cGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAggg                                                       <  1:363893/150‑1 (MQ=255)
                                                                                        gAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGa                                                 >  2:332675/1‑150 (MQ=255)
                                                                                                    gCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGt                                     <  2:1357944/150‑1 (MQ=255)
                                                                                                      aCGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGtca                                   >  1:575387/1‑150 (MQ=255)
                                                                                                       cGAATTCGTTAGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGtcat                                  <  2:640108/150‑1 (MQ=255)
                                                                                                           gTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCAtt                              >  1:542815/1‑150 (MQ=255)
                                                                                                                 cgtGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCggg                        <  1:499171/150‑1 (MQ=255)
                                                                                                                         cAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAAc                <  2:542815/150‑1 (MQ=255)
                                                                                                                          agagATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAg                                                         <  2:1125009/108‑1 (MQ=255)
                                                                                                                          agagATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAg                                                         >  1:1125009/1‑108 (MQ=255)
                                                                                                                                  aGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCttt       >  2:1414761/1‑150 (MQ=255)
                                                                                                                                    tcttttcgcTGCTGAGCAGCGCAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCt     >  1:202435/7‑150 (MQ=255)
                                                                                                                                     cttttcgcTGCTGAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTg    >  2:980261/6‑150 (MQ=255)
                                                                                                                                       gggggctgctGAGCAGCGAAAAGACCATTTCCGGTGCGGTTCGCCTGCCTTCAGACGATCAGGCCGTACGGTTTCTCATACCCCACAAAAGCAAGGGGCGTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTggg  <  2:477552/146‑1 (MQ=255)
                                                                                                                                       gAGCTCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTggg  <  1:162528/150‑1 (MQ=255)
                                                                                                                                        ttcgcTGCTGAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAgg                                                                                         <  2:205882/60‑1 (MQ=255)
                                                                                                                                        ttcgcTGCTGAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAgg                                                                                         >  1:205882/3‑62 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
CTTATGAGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTGGG  >  NC_002947/4185459‑4185743

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: