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.G99.5% 681.2 / ‑4.3 189coding (1511/1833 nt)recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (97/91);  total (97/92)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 204 total aligned reads displayed.

CGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTCGCG‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004887‑4005060
                                                                                                |                                                                                   
cGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCt                                                                                      <  1:1060373/93‑1 (MQ=255)
cGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCt                                                                                      <  2:2698583/93‑1 (MQ=255)
cGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCt                                                                                      <  2:143774/93‑1 (MQ=255)
 gTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCg‑‑ctgctg                                                                                     >  2:2839820/1‑92 (MQ=255)
  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCg‑‑ctgctgc                                                                                    >  2:2126643/1‑91 (MQ=255)
  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCtgc                                                                                    <  2:4813092/93‑3 (MQ=255)
  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCtgc                                                                                    >  1:295323/1‑91 (MQ=255)
    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCtgcgc                                                                                  <  2:2631111/93‑5 (MQ=255)
        gCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCg‑‑gcg                                                                            >  1:4798325/1‑91 (MQ=255)
        cccAGCTGATTCACCTCGGCCTGGTGACGAAAAATATTGCGCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCg‑‑gcg                                                                            <  2:2570608/92‑3 (MQ=255)
                aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCctct                                                                      >  2:5530869/1‑93 (MQ=255)
                aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCctct                                                                      <  1:609470/93‑1 (MQ=255)
                 ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCtt                                                                     >  1:5648801/1‑93 (MQ=255)
                  tCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCttt                                                                    <  1:3078095/93‑1 (MQ=255)
                  tCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCttt                                                                    <  2:4258183/93‑1 (MQ=255)
                   cACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTg                                                                   <  1:3949835/93‑1 (MQ=255)
                    cccTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGc                                                                  <  1:1881852/92‑1 (MQ=255)
                    cccTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGc                                                                  <  2:1764047/92‑1 (MQ=255)
                       tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAAc                                                               >  2:5626980/1‑93 (MQ=255)
                       tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAAc                                                               >  2:2649514/1‑93 (MQ=255)
                       tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAAc                                                               >  2:647907/1‑93 (MQ=255)
                          gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtg                                                            >  2:4415834/1‑93 (MQ=255)
                          gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtg                                                            <  1:4888766/93‑1 (MQ=255)
                           ccTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgc                                                           <  2:3026203/93‑1 (MQ=255)
                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                                          <  1:2319468/93‑1 (MQ=255)
                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                                          >  2:1945712/1‑93 (MQ=255)
                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                                          <  1:3721942/93‑1 (MQ=255)
                             tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgccg                                                         >  1:4225600/1‑93 (MQ=255)
                             tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgccg                                                         <  1:1366513/93‑1 (MQ=255)
                              ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                                        <  2:3845482/93‑1 (MQ=255)
                              ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                                        <  2:5968971/93‑1 (MQ=255)
                                tGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgc                                                      >  1:1559833/1‑93 (MQ=255)
                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgcc                                                     >  1:737909/1‑93 (MQ=255)
                                     cAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGt                                                 >  1:2814729/1‑93 (MQ=255)
                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  2:1683698/1‑93 (MQ=255)
                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  2:1170528/1‑93 (MQ=255)
                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  2:2370855/1‑93 (MQ=255)
                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                <  1:3630741/93‑1 (MQ=255)
                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  2:4664803/1‑93 (MQ=255)
                                       aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               <  2:3489286/93‑1 (MQ=255)
                                       aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               >  2:3414661/1‑93 (MQ=255)
                                        aaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATc                                              <  2:4431861/93‑1 (MQ=255)
                                         atatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCg                                             <  2:3083277/93‑1 (MQ=255)
                                         atatTCCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCg                                             <  1:2727634/93‑1 (MQ=255)
                                           atTGCCCAGCATTGTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                           <  2:498025/93‑1 (MQ=255)
                                           atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTCTCGTg                                           >  2:6066578/1‑93 (MQ=255)
                                           atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                           <  1:4753517/93‑1 (MQ=255)
                                           atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGAATCGTg                                           >  1:1030131/1‑93 (MQ=255)
                                            ttGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                                          <  2:1037462/93‑1 (MQ=255)
                                            ttGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                                          <  1:5049261/93‑1 (MQ=255)
                                             tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgc                                         <  2:68579/93‑1 (MQ=255)
                                                ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      <  1:3289315/93‑1 (MQ=255)
                                                ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  2:5797120/1‑93 (MQ=255)
                                                  aGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                                    >  2:42412/1‑93 (MQ=255)
                                                  aGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                                    >  2:1363505/1‑93 (MQ=255)
                                                  aGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                                    >  1:4655796/1‑93 (MQ=255)
                                                  aGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                                    >  2:3145927/1‑93 (MQ=255)
                                                   gCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaa                                   >  1:4998299/1‑93 (MQ=255)
                                                    cATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaaa                                  <  1:5977428/93‑1 (MQ=255)
                                                     aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 >  1:5450267/1‑93 (MQ=255)
                                                     aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 >  1:337056/1‑93 (MQ=255)
                                                     aTTCTGCCATACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 <  1:1477040/93‑1 (MQ=255)
                                                        cTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGa                              >  2:4790687/1‑93 (MQ=255)
                                                         tGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaa                             <  2:4809197/93‑1 (MQ=255)
                                                         tGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaa                             >  2:3797509/1‑93 (MQ=255)
                                                           cccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  1:5771225/93‑1 (MQ=255)
                                                            ccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  2:5234938/1‑93 (MQ=255)
                                                            ccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  2:368662/1‑93 (MQ=255)
                                                            ccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  1:892257/1‑93 (MQ=255)
                                                            ccTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  1:2440838/1‑93 (MQ=255)
                                                             cTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGTCGCTCAAACCGAAAGCg                         >  2:201437/1‑93 (MQ=255)
                                                             cTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg                         >  1:2544557/1‑93 (MQ=255)
                                                             cTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg                         <  1:6040366/93‑1 (MQ=255)
                                                              tACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        >  1:1985263/1‑93 (MQ=255)
                                                              tACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        >  1:1510820/1‑93 (MQ=255)
                                                              tACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        <  2:2192508/93‑1 (MQ=255)
                                                               aCAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGAt                       >  1:1690995/1‑93 (MQ=255)
                                                                cAACTGGCAGAGGCCGCGCGCCCGATGCTGCTCGGC‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATg                      <  2:2111542/93‑1 (MQ=255)
                                                                 aaCTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     >  2:5389368/1‑93 (MQ=255)
                                                                 aaCTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     <  2:1693106/93‑1 (MQ=255)
                                                                 aaCTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     <  1:5357166/93‑1 (MQ=255)
                                                                  aCTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    >  1:5131318/1‑93 (MQ=255)
                                                                  aCTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    <  1:1233682/93‑1 (MQ=255)
                                                                  aCTGACAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    <  2:3896197/93‑1 (MQ=255)
                                                                       cAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               <  1:342652/93‑1 (MQ=255)
                                                                       cAGAGGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               <  2:1231476/93‑1 (MQ=255)
                                                                         gagGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             >  2:2973763/1‑93 (MQ=255)
                                                                          agGCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGt            <  1:4044773/93‑1 (MQ=255)
                                                                           ggCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           <  1:3896281/93‑1 (MQ=255)
                                                                           ggCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           <  2:3343673/93‑1 (MQ=255)
                                                                            gCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  1:6071625/93‑1 (MQ=255)
                                                                            gCCGCGCGCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  2:5869253/93‑1 (MQ=255)
                                                                              cgcgcgCCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcgg        >  2:4038532/1‑93 (MQ=255)
                                                                               gcgcgcCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       >  2:315475/1‑93 (MQ=255)
                                                                               gcgcgcCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       <  2:1830231/93‑1 (MQ=255)
                                                                                cgcgcCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggcg      <  1:6141102/93‑1 (MQ=255)
                                                                                cgcgcCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggcg      >  1:4459486/1‑93 (MQ=255)
                                                                                 gcgcCCG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggcgg     <  1:4621932/93‑1 (MQ=255)
                                                                                    cccG‑‑GTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  <  2:3742105/93‑1 (MQ=255)
                                                                                    cccG‑‑GTGCTGCGCGGCGAATCCTATTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  >  2:5425263/1‑93 (MQ=255)
                                                                                                |                                                                                   
CGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCG‑‑GTGCTCGCG‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004887‑4005060

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