Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,207,035 G→T 8.2% C247* (TGC→TGA)  PP_3691 ← DNA helicase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,207,0350GT 8.2% 153.4 / 4.1 49C247* (TGC→TGA) PP_3691DNA helicase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base T (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

CCGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCATAGTCGAGCTTCATCCGTGCGACGGCGCTCCAAACCGC  >  NC_002947/4206891‑4207181
                                                                                                                                                |                                                                                                                                                  
ccGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACAGCGAGTGAATACGATGACGGGGTTGAGCCCGATTTGCTACAGGGCAGCGCAg                                                                                                                                                 >  8:75973/1‑148 (MQ=255)
  gCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCTCAGaa                                                                                                                                               >  3:367349/1‑148 (MQ=255)
     gCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAAcag                                                                                                                                            >  7:220880/1‑148 (MQ=255)
      cTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAAcagca                                                                                                                                          <  5:160322/149‑1 (MQ=255)
         cACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAgg                                                                                                                                        <  1:114507/148‑1 (MQ=255)
           cccGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTg                                                                                                                                     <  6:270475/149‑1 (MQ=255)
            ccGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTga                                                                                                                                    <  1:303171/149‑1 (MQ=255)
                     aaCCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATg                                                                                                                            >  5:150861/1‑148 (MQ=255)
                                 atcatcGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGTGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGc                                                                                                                    <  4:143213/144‑1 (MQ=255)
                                 atcatcGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGc                                                                                                                    >  3:143213/1‑144 (MQ=255)
                                 atcatcGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAAc                                                                                                               >  1:255730/1‑149 (MQ=255)
                                         ttAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCTCAGAACAGCAgg                                                                                                                                        >  8:353605/1‑116 (MQ=255)
                                         ttAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCTCAGAACAGCAgg                                                                                                                                        <  7:353605/116‑1 (MQ=255)
                                         ttAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCTCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCa                                                                                                       <  4:367349/149‑1 (MQ=255)
                                                aGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCa                                                                                                                          >  3:324996/1‑123 (MQ=255)
                                                aGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCa                                                                                                                          <  4:324996/123‑1 (MQ=255)
                                                 gTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAg                                                                                                                         <  7:320211/123‑1 (MQ=255)
                                                 gTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAg                                                                                                                         >  8:320211/1‑123 (MQ=255)
                                                       aaCTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTg                                                                                          <  6:150861/148‑1 (MQ=255)
                                                            ctctGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTAc                                                                                    >  6:228047/1‑149 (MQ=255)
                                                                     cATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCttt                                                                            <  2:255730/148‑1 (MQ=255)
                                                                        ccACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGAc                                                                        >  3:242749/1‑149 (MQ=255)
                                                                         cacaGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCtttt                                                                           >  6:382637/1‑145 (MQ=255)
                                                                         cacaGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCtttt                                                                           <  5:382637/145‑1 (MQ=255)
                                                                         cacaGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGAcc                                                                       <  5:228047/149‑1 (MQ=255)
                                                                          acaGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCt                                                                      >  2:212514/1‑149 (MQ=255)
                                                                                cTGTGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGt                                                                >  7:295446/1‑149 (MQ=255)
                                                                                 tGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAATCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGt                                                                >  3:207775/1‑148 (MQ=255)
                                                                                     aaaGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCaaa                                                           <  7:75973/149‑1 (MQ=255)
                                                                                               cgcgAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCtggctggc                                                 >  2:284885/1‑149 (MQ=255)
                                                                                                    gTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCaa                                            <  8:220880/149‑1 (MQ=255)
                                                                                                      gAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAATCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAAt                                           <  4:207775/148‑1 (MQ=255)
                                                                                                      gAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAAt                                           <  4:242749/148‑1 (MQ=255)
                                                                                                       aaTACGATGACCGGGTTGAGCACGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGc                                         <  5:101410/149‑1 (MQ=255)
                                                                                                        aTACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCa                                        >  5:44868/1‑149 (MQ=255)
                                                                                                        aTACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCa                                        <  1:284885/149‑1 (MQ=255)
                                                                                                        aTACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCATGGGCTAAGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCa                                        >  7:292884/1‑149 (MQ=255)
                                                                                                            gATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAgctggct                                                    <  7:106031/133‑1 (MQ=255)
                                                                                                            gATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAgctggct                                                    >  8:106031/1‑133 (MQ=255)
                                                                                                            gATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCATAGt                                    >  3:297235/1‑149 (MQ=255)
                                                                                                                  cGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCtt                                                                             <  2:354057/102‑1 (MQ=255)
                                                                                                                  cGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCtt                                                                             >  1:354057/1‑102 (MQ=255)
                                                                                                                  cGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATg                                          >  7:10649/1‑137 (MQ=255)
                                                                                                                  cGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATg                                          <  8:10649/137‑1 (MQ=255)
                                                                                                                  cGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCATAGTCGAGc                               >  4:339717/1‑148 (MQ=255)
                                                                                                                  cGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCATAGTCGAGc                               <  3:339717/148‑1 (MQ=255)
                                                                                                                    ggTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCATAGTCGAGCTTc                            <  8:295446/149‑1 (MQ=255)
                                                                                                                                     tCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCATAGTCGAGCTTCATCCGTg                     <  1:188516/139‑1 (MQ=255)
                                                                                                                                     tCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCATAGTCGAGCTTCATCCGTg                     >  2:188516/1‑139 (MQ=255)
                                                                                                                                              gcgcAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCATAGTCGAGCTTCATCCGTGCGACGGCGCTCCAAACCGc  <  6:44868/149‑1 (MQ=255)
                                                                                                                                               cgcAGCACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCATAGTCGCGCTTCATCCGTGCGACGGCGCTCCAAACCg   >  3:92820/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCGCAGCTTCACCCGCTTGCGAACCGCCGTAACATCATCGATTAGAGCAGTCGTGAACTCCTCTGCCTTCATCCACAGTCCTGCGAAAGCGTGACCGCGAGTGAATACGATGACCGGGTTGAGCCCGATTTGCTCCAGGGCAGCGCAGAACAGCAGGGTGAGATCCATGCAGGTAGCCAAACCAGTATCAGCAATCTGGCTTGGGCTACGCACCTTTTGACCTGACTGTTCAAAGCTGGCTGGCGGCAATGCATAGTCGAGCTTCATCCGTGCGACGGCGCTCCAAACCGC  >  NC_002947/4206891‑4207181

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