Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,004,979 +G coding (1511/1833 nt) recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,9791.G97.8% 238.4 / ‑0.9 92coding (1511/1833 nt)recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (52/38);  total (52/40)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTCGCG‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCA  >  NC_000913/4004904‑4005042
                                                                             |                                                                  
ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTTCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCt                                                                     >  2:1879397/1‑76 (MQ=255)
   aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcg                                                                  >  1:1090881/1‑73 (MQ=255)
   aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcg                                                                  >  1:451850/1‑73 (MQ=255)
       cGGGCTGGTGACGCAAAATATTGCCCAGCATGCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gc                                                           <  2:360751/76‑2 (MQ=255)
       cGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gc                                                           <  2:3431158/76‑2 (MQ=255)
        ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcg                                                          <  1:3877596/76‑3 (MQ=255)
                gcggCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCctct                                                    <  2:251432/73‑1 (MQ=255)
                gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCctct                                                    >  1:3993385/1‑76 (MQ=255)
                 aCGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCtt                                                   <  2:4299271/76‑1 (MQ=255)
                  cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCttt                                                  <  2:879175/76‑1 (MQ=255)
                  cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCttt                                                  <  1:2340450/76‑1 (MQ=255)
                  cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCttt                                                  >  2:4400619/1‑76 (MQ=255)
                  cGCAAAATATTGCCCAGCATCCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCttt                                                  <  1:1954907/76‑1 (MQ=255)
                   gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                 >  1:3747828/1‑76 (MQ=255)
                   gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                 >  1:1891309/1‑76 (MQ=255)
                   gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                 >  1:3660345/1‑76 (MQ=255)
                   gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                 <  1:1406360/76‑1 (MQ=255)
                   gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                 <  2:566997/76‑1 (MQ=255)
                   gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                 >  2:3822628/1‑76 (MQ=255)
                    cAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGc                                                >  1:1786731/1‑76 (MQ=255)
                     agaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCa                                               >  1:1898600/3‑76 (MQ=255)
                     aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCa                                               >  1:2191795/1‑76 (MQ=255)
                        atgttGCGCAGCATTCTGCGCTCCCACGGACGGAGGCCGCGCGCCCGGGGCTGCGCGGCGAATCCTCTTTGCAACt                                            <  2:3994437/73‑1 (MQ=255)
                        atatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACt                                            <  2:4482721/76‑1 (MQ=255)
                         tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                           >  2:2455634/1‑76 (MQ=255)
                         tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                           >  2:2120938/1‑76 (MQ=255)
                         tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                           <  2:1889380/76‑1 (MQ=255)
                         tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                           <  2:2967839/76‑1 (MQ=255)
                         tatTGCCCAGCATTCGGCCCTACAACTGACAGAGGCCGCGCGACCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                           >  1:3182386/1‑76 (MQ=255)
                          atggCCCAGCATTCTGCCCTACAACTGGCAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtg                                          <  2:249620/73‑1 (MQ=255)
                          atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtg                                          >  1:4149781/1‑76 (MQ=255)
                            tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                        >  2:4577906/1‑76 (MQ=255)
                            tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                        >  2:298070/1‑76 (MQ=255)
                             gCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccg                                       <  2:3500197/76‑1 (MQ=255)
                              cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                      >  2:3048800/1‑76 (MQ=255)
                              cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                      <  2:2838792/76‑1 (MQ=255)
                              cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                      >  2:3052208/1‑76 (MQ=255)
                                cAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgc                                    >  1:1069943/1‑76 (MQ=255)
                                 aGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgcc                                   <  1:2379308/76‑1 (MQ=255)
                                   cATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgcccc                                 >  2:3532228/1‑74 (MQ=255)
                                   cATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgc                                 <  2:2836003/76‑1 (MQ=255)
                                   cATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgc                                 >  1:1643138/1‑76 (MQ=255)
                                    aTTCTGCCCTACAGCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgcg                                <  2:3720050/76‑1 (MQ=255)
                                    aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgcg                                >  2:2052186/1‑76 (MQ=255)
                                    aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgcg                                >  1:848132/1‑76 (MQ=255)
                                     ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGt                               >  1:3613339/1‑76 (MQ=255)
                                     ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGt                               >  1:368432/1‑76 (MQ=255)
                                      tCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGtc                              >  1:1499711/1‑75 (MQ=255)
                                      tCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                              >  1:79260/1‑76 (MQ=255)
                                        tGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCCGCGAATGCTCTTTGCAACTTGCCGTGCCGCGtctc                            >  1:4269513/1‑76 (MQ=255)
                                          cccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                          <  2:1586783/76‑1 (MQ=255)
                                          cccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                          <  1:4000300/76‑1 (MQ=255)
                                           ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                         <  2:2619070/76‑1 (MQ=255)
                                           ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                         <  2:2150769/76‑1 (MQ=255)
                                            cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTCTCGTgg                        >  2:1421115/1‑76 (MQ=255)
                                            cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                        <  1:3436026/76‑1 (MQ=255)
                                            cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                        >  1:4518850/1‑76 (MQ=255)
                                             tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgc                       <  1:298070/76‑1 (MQ=255)
                                              aCAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCgcggcg                      >  2:1134173/1‑76 (MQ=255)
                                               cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                     >  2:889597/1‑76 (MQ=255)
                                                aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                    >  2:3644912/1‑76 (MQ=255)
                                                aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                    >  1:257607/1‑76 (MQ=255)
                                                aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                    <  2:1013549/76‑1 (MQ=255)
                                                 aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACCTGCCGTGCCGCGTATCGTGCCGCTc                   >  2:612722/1‑76 (MQ=255)
                                                  cTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                  >  1:2900668/1‑76 (MQ=255)
                                                   tGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaa                 >  2:2482646/1‑76 (MQ=255)
                                                    gagagagGCCGCGCGTCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGTCGCGTATCGTGGCGCTCaaa                <  2:4053771/73‑1 (MQ=255)
                                                     gcggagGCCGCGCGCCCGGTGCGGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc               <  2:4246228/73‑1 (MQ=255)
                                                     aCAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc               <  1:300051/76‑1 (MQ=255)
                                                      tagagGCCCCGCGCCCGCTGCTGCGCCGCGAATCCTCTTTGCAACTTGCCGCACCGCGTATCGCGGCGCTCAAAcc              >  2:3812501/2‑76 (MQ=255)
                                                      cAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc              <  2:1891308/76‑1 (MQ=255)
                                                         agGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaa           <  1:3768866/76‑1 (MQ=255)
                                                         agGCCGCGAGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCATACCGaa           >  1:756549/1‑76 (MQ=255)
                                                          ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCACACCGaaa          >  2:4081858/1‑76 (MQ=255)
                                                          ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa          >  1:4847733/1‑76 (MQ=255)
                                                          ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa          >  2:2431776/1‑76 (MQ=255)
                                                          ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGNATCGTGGCGCTCAAACCGaaa          <  1:1879398/76‑1 (MQ=255)
                                                          cgCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa          <  1:2408859/75‑1 (MQ=255)
                                                           gCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg         >  1:1703592/1‑76 (MQ=255)
                                                           gCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg         >  1:2873375/1‑76 (MQ=255)
                                                            ccGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGCAAGc        >  1:332696/1‑76 (MQ=255)
                                                            ccGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc        <  2:2348738/76‑1 (MQ=255)
                                                             cgcgcgCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg       >  1:1066024/1‑76 (MQ=255)
                                                             cgcgcgCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg       >  1:85647/1‑76 (MQ=255)
                                                              gcgcgcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa      >  1:4848016/1‑76 (MQ=255)
                                                              gcgcgcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa      >  1:2530357/1‑76 (MQ=255)
                                                               cgcgcCCGGGGCGTCGCGGCGAATCCTCTTTGCAACTTGCCGTGCGGCGTATCGTGGCGCTCAAACCGAAAGCGAt     <  2:2377989/76‑1 (MQ=255)
                                                                gcgcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATg    >  1:4446691/1‑76 (MQ=255)
                                                                 cgcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc   <  1:4679961/76‑1 (MQ=255)
                                                                 cgcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCCAAAGCGATGc   >  2:4477923/1‑76 (MQ=255)
                                                                  gcCCGGTGCTGCGCGGCGAATCCTCTTTGCAGCTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa  <  1:2890481/76‑1 (MQ=255)
                                                                  gcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCg             >  2:4569291/1‑65 (MQ=255)
                                                                  gcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCg             <  1:4569291/65‑1 (MQ=255)
                                                                  gcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa  <  1:2301068/76‑1 (MQ=255)
                                                                  gcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa  >  1:521360/1‑76 (MQ=255)
                                                                             |                                                                  
TTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTCGCG‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCA  >  NC_000913/4004904‑4005042

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