Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 829,599 G→T M358I (ATG→ATT dinG → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273829,5990GT95.1% 124.5 / ‑3.9 42M358I (ATG→ATTdinGATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (9/30);  total (11/31)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACA  >  CP009273/829534‑829660
                                                                 |                                                             
gCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGAttct                                                             >  1:4193767/1‑68 (MQ=255)
 cGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGAttctg                                                            >  1:5783896/1‑68 (MQ=255)
  gATGGGCGAACTGCCTGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGAttctgc                                                           >  1:2538664/1‑68 (MQ=255)
  gATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGTCAAAACTCACCGAGAttctgc                                                           >  1:5174984/1‑68 (MQ=255)
  gATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGAttctgc                                                           >  1:4805592/1‑68 (MQ=255)
  gATGGGCGAACTGCCAGATGAAGTGATGGAGATCTTCCAGCGGCTGACAAAACTCACCGAGAttctgc                                                           >  1:5876014/1‑68 (MQ=255)
      ggCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTgg                                                       <  1:144963/68‑1 (MQ=255)
       gCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGc                                                      <  1:2255725/68‑1 (MQ=255)
       gCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGc                                                      <  1:4087886/68‑1 (MQ=255)
        cGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGcc                                                     >  1:1780765/1‑68 (MQ=255)
           aCTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCAGCGTGGCCTgg                                                  <  1:3615486/68‑1 (MQ=255)
                cAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAg                                             <  1:2270429/68‑1 (MQ=255)
                cAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAg                                             <  1:5714721/68‑1 (MQ=255)
                cAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAg                                             <  1:2343232/68‑1 (MQ=255)
                  gATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGtt                                           >  1:2193988/1‑68 (MQ=255)
                     gAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTAtt                                        >  1:4174053/1‑68 (MQ=255)
                     gAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTCCGTGGCCTGGCGGAGTTAtt                                        <  1:1302288/68‑1 (MQ=255)
                       aGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTc                                      <  1:4544230/68‑1 (MQ=255)
                       aGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTc                                      <  1:5968017/68‑1 (MQ=255)
                         tGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCtt                                    >  1:3696597/1‑68 (MQ=255)
                             ggAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACg                                <  1:1962589/68‑1 (MQ=255)
                             ggAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACg                                <  1:5607990/68‑1 (MQ=255)
                                 aTCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGAttt                            >  1:271696/1‑68 (MQ=255)
                                    tGCCAGCGGCTGGCTAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAg                         >  1:6098103/1‑68 (MQ=255)
                                    tGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAg                         <  1:3065256/68‑1 (MQ=255)
                                     gCCAGCGGCTGGCAATACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTAGATCTTAACGATTGAAGt                        >  1:2852158/1‑68 (MQ=255)
                                     gCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGt                        >  1:3169576/1‑68 (MQ=255)
                                     gCCAGCGGCTGGCAAAACTCACCGAGATACTGCCTGGCCTGGCGGAGTTATTTCTTAACGATTTAAga                        >  1:850899/1‑67 (MQ=39)
                                      ccAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTg                       >  1:672213/1‑68 (MQ=255)
                                      ccAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTg                       <  1:1875419/68‑1 (MQ=255)
                                       cAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTga                      <  1:5785109/68‑1 (MQ=255)
                                         gCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTgaga                    <  1:3133632/68‑1 (MQ=255)
                                          cGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGaa                   >  1:5899873/1‑68 (MQ=255)
                                          cGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCGGGCGTAGTTATTTCTAAACGATTTAAGTGAGaa                   <  1:5232483/68‑1 (MQ=255)
                                            gCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGaaaa                 <  1:4901244/68‑1 (MQ=255)
                                                   aaaCTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAg          <  1:1103691/68‑1 (MQ=255)
                                                    aaCTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGc         <  1:5721763/68‑1 (MQ=255)
                                                    aaCTCACCGAGATGCTGCGTGGCCTGGCGGAGTTAATTCTTCACGATTTAAGTGAGAAAACCGGCAGc         <  1:5515760/68‑1 (MQ=255)
                                                     aCTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGcc        <  1:4436049/68‑1 (MQ=255)
                                                     aCTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGcc        <  1:6098240/68‑1 (MQ=255)
                                                       tCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCAt      <  1:1879677/68‑1 (MQ=255)
                                                       tCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCAt      <  1:4537668/68‑1 (MQ=255)
                                                          ccGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGAc   <  1:3572789/68‑1 (MQ=255)
                                                          ccGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGAc   <  1:4226906/68‑1 (MQ=255)
                                                          ccGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGAc   <  1:5200114/68‑1 (MQ=255)
                                                           cGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACa  <  1:4025567/68‑1 (MQ=255)
                                                           cGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCAATGACa  <  1:842177/68‑1 (MQ=255)
                                                           cGAGATTCTGCGGGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACa  <  1:4989191/68‑1 (MQ=255)
                                                                 |                                                             
GCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACA  >  CP009273/829534‑829660

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