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,5990GT97.6% 128.0 / ‑3.6 42M358I (ATG→ATTdinGATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (11/30);  total (11/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

GATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATT  >  CP009273/829536‑829662
                                                               |                                                               
gATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGAttctgc                                                             >  1:4132442/1‑68 (MQ=255)
    ggCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTgg                                                         <  1:3544185/68‑1 (MQ=255)
     gCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGc                                                        <  1:223213/68‑1 (MQ=255)
     gCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGc                                                        <  1:2774658/68‑1 (MQ=255)
     gCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGc                                                        <  1:1680578/68‑1 (MQ=255)
       gAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCt                                                      <  1:2857245/68‑1 (MQ=255)
         aCTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTgg                                                    <  1:901246/68‑1 (MQ=255)
          cTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGc                                                   >  1:4134237/1‑68 (MQ=255)
           tGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCg                                                  >  1:1125796/1‑68 (MQ=255)
              cAGATTAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAg                                               <  1:1785121/68‑1 (MQ=255)
              cAGATGAAGTGCTGGAGATCTGCCAGCGGCTGTCAAAACTCACCGAGATTCTGCGTGGTCTGGCGGAg                                               <  1:1936102/68‑1 (MQ=255)
              cAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCTCCGAGATTCTGCGTGGCCTGGCGGAg                                               <  1:186800/68‑1 (MQ=255)
              cAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAg                                               <  1:2060420/68‑1 (MQ=255)
              cAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAg                                               <  1:2519884/68‑1 (MQ=255)
              cAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAg                                               <  1:451672/68‑1 (MQ=255)
                   gAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTAtt                                          >  1:3956254/1‑68 (MQ=255)
                       tGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCtt                                      >  1:2157999/1‑68 (MQ=255)
                           ggAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACg                                  >  1:3885021/1‑68 (MQ=255)
                           ggAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACg                                  >  1:3715598/1‑68 (MQ=255)
                             agaTCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGAt                                <  1:1656178/68‑1 (MQ=255)
                             agaTCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGAt                                <  1:880030/68‑1 (MQ=255)
                                tCTGCCAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTa                             >  1:4016177/1‑68 (MQ=255)
                                    ccAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTg                         >  1:4360634/1‑68 (MQ=255)
                                    ccAGCGGCTGGCAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTACTTAACGATTTAAGTg                         >  1:758756/1‑68 (MQ=255)
                                    ccAGCGGCGGGCAAAACTCAGCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTg                         <  1:2530866/68‑1 (MQ=255)
                                               cAAAACTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGc              <  1:1635040/68‑1 (MQ=255)
                                                aaaaCTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCa             <  1:2816753/68‑1 (MQ=255)
                                                  aaCTCACCGAGAATCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGc           <  1:3258879/68‑1 (MQ=255)
                                                   aCTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGcc          <  1:2370394/68‑1 (MQ=255)
                                                   aCTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGcc          <  1:4047070/68‑1 (MQ=255)
                                                   aCTCACCGAGATTCAGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGcc          <  1:4529087/68‑1 (MQ=255)
                                                    cTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCa         >  1:3345239/1‑68 (MQ=255)
                                                    cTCACCGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCa         >  1:3962033/1‑68 (MQ=255)
                                                        ccGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTTAGTGAGAAAACCGGCAGCCATGAc     <  1:1209247/68‑1 (MQ=255)
                                                        ccGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTATGTGAGAAAACCGGCAGCCATGAc     <  1:83831/68‑1 (MQ=255)
                                                        ccGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGAc     <  1:231403/68‑1 (MQ=255)
                                                        ccGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGAc     <  1:2041783/68‑1 (MQ=255)
                                                         cGAGATTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACa    <  1:3691899/68‑1 (MQ=255)
                                                         cGAGATTCTGAGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACa    <  1:2992535/68‑1 (MQ=255)
                                                          gagaTTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACAt   <  1:2239330/68‑1 (MQ=255)
                                                           tgaTTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACAtt  <  1:3177499/67‑1 (MQ=255)
                                                           agaTTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACAtt  <  1:306300/68‑1 (MQ=255)
                                                           agaTTCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACAtt  <  1:146676/68‑1 (MQ=255)
                                                               |                                                               
GATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATT  >  CP009273/829536‑829662

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: