Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,004,984 +CG coding (1516/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,9831.G97.8% 228.0 / ‑3.0 89A505A (GCG→GCCrecQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (50/37);  minor base C (2/0);  total (52/37)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
*NC_0009134,004,9832.C97.8% 228.9 / ‑0.2 89coding (1515/1833 nt)recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (50/37);  total (52/37)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

CGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTCGCG‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCA  >  NC_000913/4004911‑4005042
                                                                            ||                                                           
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                                                  >  2:4400619/1‑76 (MQ=255)
           cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCttt                                                  <  1:2340450/76‑1 (MQ=255)
           cGCAAAATATTGCCCAGCATCCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCttt                                                  <  1:1954907/76‑1 (MQ=255)
            gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                 >  1:3747828/1‑76 (MQ=255)
            gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                 >  2:3822628/1‑76 (MQ=255)
            gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                 >  1:3660345/1‑76 (MQ=255)
            gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                 <  1:1406360/76‑1 (MQ=255)
            gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                 >  1:1891309/1‑76 (MQ=255)
            gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                 <  2:566997/76‑1 (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:2967839/76‑1 (MQ=255)
                  tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                           >  2:2455634/1‑76 (MQ=255)
                  tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                           <  2:1889380/76‑1 (MQ=255)
                  tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                           >  2:2120938/1‑76 (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:298070/1‑76 (MQ=255)
                     tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                        >  2:4577906/1‑76 (MQ=255)
                      gCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccg                                       <  2:3500197/76‑1 (MQ=255)
                       cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                      <  2:2838792/76‑1 (MQ=255)
                       cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                      >  2:3048800/1‑76 (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:368432/1‑76 (MQ=255)
                              ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGt                               >  1:3613339/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                          <  1:4000300/76‑1 (MQ=255)
                                   cccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                          <  2:1586783/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                    >  1:257607/1‑76 (MQ=255)
                                         aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                    <  2:1013549/76‑1 (MQ=255)
                                         aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                    >  2:3644912/1‑76 (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          >  2:2431776/1‑76 (MQ=255)
                                                   ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa          >  1:4847733/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:521360/1‑76 (MQ=255)
                                                           gcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa  <  1:2301068/76‑1 (MQ=255)
                                                                            ||                                                           
CGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTCGCG‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCA  >  NC_000913/4004911‑4005042

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