Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 834,352 A→T Q428L (CAG→CTG)  dinG → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913834,3520AT98.4% 213.6 / ‑4.0 61Q428L (CAG→CTG) dinGATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (26/34);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TGTTCGAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCAC  >  NC_000913/834253‑834442
                                                                                                   |                                                                                          
tGTTCGAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCa                                                                                            <  2:1366962/100‑1 (MQ=255)
   tCGAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGgctgct                                                                                         >  2:1248338/1‑100 (MQ=255)
      aGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTaca                                                                                      <  1:1628995/100‑1 (MQ=255)
        gcgcAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACAcc                                                                                    <  1:718858/100‑1 (MQ=255)
            aaagcaaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTg                                                                                <  2:776782/100‑1 (MQ=255)
            aaagcaaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTg                                                                                >  1:862871/1‑100 (MQ=255)
             aagcaaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTgg                                                                               <  2:355862/100‑1 (MQ=255)
             aagcaaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTgg                                                                               >  1:1235838/1‑100 (MQ=255)
              agcaaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGt                                                                              >  1:1248944/1‑100 (MQ=255)
                caaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGttt                                                                            >  2:1639813/1‑100 (MQ=255)
                 aaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTTCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACTCCTATGGTTTc                                                                           <  1:1026461/100‑1 (MQ=255)
                   aCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCAc                                                                         <  2:994481/100‑1 (MQ=255)
                   aCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCAc                                                                         <  2:172689/100‑1 (MQ=255)
                    ctctGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACt                                                                        >  2:334278/1‑100 (MQ=255)
                    ctctGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACt                                                                        >  2:353157/1‑100 (MQ=255)
                     tctGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTg                                                                       >  1:1875303/1‑100 (MQ=255)
                        ggcggcTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGt                                                                    <  2:510951/100‑1 (MQ=255)
                           ggctggctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTggg                                                                 <  1:915521/100‑1 (MQ=255)
                           ggctggctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTggg                                                                 <  1:999464/100‑1 (MQ=255)
                             ctggctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGaa                                                               <  2:726002/100‑1 (MQ=255)
                               ggctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATa                                                             <  2:136850/100‑1 (MQ=255)
                               ggctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATa                                                             <  2:890932/100‑1 (MQ=255)
                               ggctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATa                                                             >  2:714242/1‑100 (MQ=255)
                                 ctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACg                                                           <  1:1742990/100‑1 (MQ=255)
                                 ctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACg                                                           >  2:1686295/1‑100 (MQ=255)
                                   tCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACgtg                                                         <  1:919684/100‑1 (MQ=255)
                                     gCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTc                                                       <  1:1105370/100‑1 (MQ=255)
                                     gCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTc                                                       <  2:474798/100‑1 (MQ=255)
                                     gCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTc                                                       >  1:1140718/1‑100 (MQ=255)
                                            cAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATc                                                >  2:448050/1‑100 (MQ=255)
                                            cAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATc                                                <  2:1346579/100‑1 (MQ=255)
                                            cAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGCTTCACTGCGTGGGAATACGTGTCAGCGAt                                                 >  2:221610/1‑99 (MQ=255)
                                               tCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCGCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGc                                             >  1:872127/1‑100 (MQ=255)
                                                cTTCCGGTGCACCGGGGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCt                                            <  1:1148065/100‑1 (MQ=255)
                                                  tCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTgg                                          >  2:1293178/1‑100 (MQ=255)
                                                     ggTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGaaa                                       >  1:1136574/1‑100 (MQ=255)
                                                      gTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAg                                      <  2:1633702/100‑1 (MQ=255)
                                                             ggTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGctgctg                               <  2:1166504/100‑1 (MQ=255)
                                                             ggTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGctgctg                               >  2:1917474/1‑100 (MQ=255)
                                                               tGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCtgtg                             >  2:1501174/1‑100 (MQ=255)
                                                                gACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTgg                            <  2:928985/100‑1 (MQ=255)
                                                                 aCCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGgc                           <  2:1929625/100‑1 (MQ=255)
                                                                  ccAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGgcg                          <  1:1293171/100‑1 (MQ=255)
                                                                   cAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGgcgc                         <  1:1391425/100‑1 (MQ=255)
                                                                     aaTGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAg                       <  2:1466934/100‑1 (MQ=255)
                                                                     aaTGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAg                       <  1:1721018/100‑1 (MQ=255)
                                                                       tGGGCGACGCGGGAAGAGCGCGAAGGGCTACTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTa                     <  2:1309675/100‑1 (MQ=255)
                                                                          gCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTc                  <  2:1421161/100‑1 (MQ=255)
                                                                           cGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTcc                 <  1:1248334/100‑1 (MQ=255)
                                                                           cGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTcc                 >  2:1935183/1‑100 (MQ=255)
                                                                             acgtgtgAAGAGCGCGAATGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGc               <  2:1056897/94‑1 (MQ=255)
                                                                             aCGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGc               <  2:434766/100‑1 (MQ=255)
                                                                             aCGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGc               <  2:1586780/100‑1 (MQ=255)
                                                                                cgatctgagCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGcaca            <  1:1927831/94‑1 (MQ=255)
                                                                                cgGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGcaca            <  2:1942694/100‑1 (MQ=255)
                                                                                cgGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGcaca            >  2:1716501/1‑100 (MQ=255)
                                                                                     agagCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACattatt       >  1:95930/1‑100 (MQ=255)
                                                                                      gagCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGcacaagatcg      >  2:1927844/1‑94 (MQ=255)
                                                                                      gagCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTg      >  1:345310/1‑100 (MQ=255)
                                                                                      gagCGCGAACGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTg      >  2:93853/1‑100 (MQ=255)
                                                                                         cgcgAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCa   >  2:1684160/1‑100 (MQ=255)
                                                                                         cgcgAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCa   >  1:945099/1‑100 (MQ=255)
                                                                                          gcgAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCAc  >  2:1030763/1‑100 (MQ=255)
                                                                                                   |                                                                                          
TGTTCGAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCAC  >  NC_000913/834253‑834442

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