Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,185,604 C→G 9.2% 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,6040CG 9.2% 128.9 / 5.6 44L446L (CTC→CTG) PP_3681putative Helicase
Reads supporting (aligned to +/- strand):  ref base C (20/19);  new base G (2/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

TGAGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTGGGGCAAGCTCG  >  NC_002947/4185463‑4185752
                                                                                                                                             |                                                                                                                                                    
tGAGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGa                                                                                                                                              <  1:823549/150‑1 (MQ=255)
    gACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAg                                                                                                                                          >  2:872702/1‑150 (MQ=255)
      cTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGAc                                                                                                                                        <  2:313758/150‑1 (MQ=255)
       tggtggAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGAcc                                                                                                                                       <  2:142837/150‑1 (MQ=255)
          tggAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTGAGCAGCgct                                                                                                                                             >  2:254713/1‑139 (MQ=255)
          tggAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTGAGCAGCgct                                                                                                                                             <  1:254713/141‑3 (MQ=255)
                   acaAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGcc                                                                                                                           >  1:1164439/1‑150 (MQ=255)
                     aaGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGTAGCGCTGCTCATCAGCGAAAATACCATTCCCGGTGCCCt                                                                                                                         >  2:197413/1‑150 (MQ=255)
                          gcTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGccc                                                                                                                    >  2:92299/1‑150 (MQ=255)
                           cTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCg                                                                                                                   <  1:92299/150‑1 (MQ=255)
                                tCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTcag                                                                                                                                                    <  2:671664/112‑1 (MQ=255)
                                tCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTcag                                                                                                                                                    >  1:671664/1‑112 (MQ=255)
                                 ccGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGaa                                                                                                                                             <  1:173910/118‑1 (MQ=255)
                                 ccGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGaa                                                                                                                                             >  2:173910/1‑118 (MQ=255)
                                          aGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCt                                                                                                                <  2:622101/138‑1 (MQ=255)
                                          aGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCt                                                                                                                >  1:622101/1‑138 (MQ=255)
                                            gTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAgg                                                                                                  >  1:470495/1‑150 (MQ=255)
                                            gTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAgg                                                                                                  >  1:393497/1‑150 (MQ=255)
                                                tCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTGAGCAGCgc                                                                                                                                              <  1:780115/102‑2 (MQ=255)
                                                tCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTGAGCAGCgc                                                                                                                                              >  2:780115/1‑101 (MQ=255)
                                                tCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGTCCATATATGTCCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGt                                                                                              >  2:31227/1‑150 (MQ=255)
                                                  cAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTAc                                                                                            <  1:872702/150‑1 (MQ=255)
                                                  cAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTAc                                                                                            >  2:280992/1‑150 (MQ=255)
                                                     tCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCt                                                                                         >  2:887396/1‑150 (MQ=255)
                                                        tGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTtc                                                                                      <  2:940128/150‑1 (MQ=255)
                                                        tGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTtc                                                                                      <  1:862296/150‑1 (MQ=255)
                                                              cgcTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAAGTCAGAGATCGCGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACaa                                                                                <  2:870884/150‑1 (MQ=255)
                                                                 tctcAGGTGACTACGAGTCGAAACTACGATTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATcc                                                                             <  1:887396/150‑1 (MQ=255)
                                                                  ctcAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCca                                                                            <  1:280992/150‑1 (MQ=255)
                                                                   tcAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCcac                                                                           <  1:197413/150‑1 (MQ=255)
                                                                       gTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACaaaa                                                                       >  1:610351/1‑150 (MQ=255)
                                                                        tGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAg                                                                      >  1:427500/1‑150 (MQ=255)
                                                                              cGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAggg                                                                <  2:393497/150‑1 (MQ=255)
                                                                              cGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAggg                                                                >  2:398741/1‑150 (MQ=255)
                                                                               gAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAgggg                                                               >  1:195325/1‑150 (MQ=255)
                                                                                      aaCTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGt                                                        <  2:610351/150‑1 (MQ=255)
                                                                                             cTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCacac                                                 >  1:186542/1‑150 (MQ=255)
                                                                                               tGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACgg                                               <  2:186542/150‑1 (MQ=255)
                                                                                                   cGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGtcat                                           >  1:155568/1‑150 (MQ=255)
                                                                                                           ttcgtGAGGTCAGAGATGGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGCTTCACACGGTCATCATTCTCg                                   >  2:618168/1‑150 (MQ=255)
                                                                                                                gAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTc                              <  1:398741/150‑1 (MQ=255)
                                                                                                                    tCAGCGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGTCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGaaa                          >  1:398491/1‑150 (MQ=255)
                                                                                                                    tCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGaaa                          <  2:639367/150‑1 (MQ=255)
                                                                                                                         aaTCGACTTAGAGCGCTGCTAAGCAGCAAAAAGACCATTCCCGGTGCCCTTGCCCGGCATTCCAACGATCAGTCCGTACGCTTTCTCACAATCCACAAAGGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCaa                     <  2:427500/149‑1 (MQ=255)
                                                                                                                                ttGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCGAGCTTTCt              <  2:155568/150‑1 (MQ=255)
                                                                                                                                            tCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTGGGGCAAGCTCg  >  2:190131/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
TGAGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTGGGGCAAGCTCG  >  NC_002947/4185463‑4185752

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑