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.8% 1603.4 / ‑5.7 474L504R (CTC→CGC) recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (259/214);  minor base A (1/0);  total (260/214)
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.94e-01

Only 100 of 499 total aligned reads displayed.

CGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCT‑CGCG‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004887‑4005060
                                                                                             |                                                                                   
  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgc                                                                                    <  1:16627/93‑3 (MQ=255)
  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgc                                                                                    <  1:918981/93‑3 (MQ=255)
  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgc                                                                                    >  2:3427424/1‑91 (MQ=255)
   gATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcg                                                                                   >  2:7234143/1‑90 (MQ=255)
    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                  >  2:5274259/1‑89 (MQ=255)
       cGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gc                                                                             <  2:6799410/93‑2 (MQ=255)
       cGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gc                                                                             <  2:5060368/93‑2 (MQ=255)
         ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcga                                                                           >  1:5572406/1‑90 (MQ=255)
                aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCctct                                                                      >  1:3548927/1‑93 (MQ=255)
                aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCctct                                                                      <  1:4946365/93‑1 (MQ=255)
                aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCctct                                                                      <  1:744637/93‑1 (MQ=255)
                 ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCtt                                                                     <  1:2233831/93‑1 (MQ=255)
                  tCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCttt                                                                    <  2:464999/93‑1 (MQ=255)
                    aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGc                                                                  <  2:3904953/93‑1 (MQ=255)
                    aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGc                                                                  >  1:6525825/1‑93 (MQ=255)
                    aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGc                                                                  >  1:6267773/1‑93 (MQ=255)
                    aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGc                                                                  >  1:3381653/1‑93 (MQ=255)
                    aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGc                                                                  >  1:1820356/1‑93 (MQ=255)
                      cTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCaa                                                                >  2:3218273/1‑93 (MQ=255)
                       tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAAc                                                               >  1:997873/1‑93 (MQ=255)
                       tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAAc                                                               <  1:1951021/93‑1 (MQ=255)
                       tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAAc                                                               >  2:2522433/1‑93 (MQ=255)
                         ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                                             >  1:2667044/1‑93 (MQ=255)
                          gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtg                                                            >  2:6210420/1‑93 (MQ=255)
                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                                          >  2:5963413/1‑93 (MQ=255)
                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                                          >  1:7045036/1‑93 (MQ=255)
                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                                          >  2:5401958/1‑93 (MQ=255)
                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                                          <  1:5921243/93‑1 (MQ=255)
                             tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccg                                                         >  1:3739119/1‑93 (MQ=255)
                             tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccg                                                         >  1:2754104/1‑93 (MQ=255)
                             tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccg                                                         >  1:2040740/1‑93 (MQ=255)
                             tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAAGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccg                                                         >  1:2927493/1‑93 (MQ=255)
                              ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                                        <  2:6578013/93‑1 (MQ=255)
                              ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                                        <  1:4811017/93‑1 (MQ=255)
                               gTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtg                                                       >  1:5702037/1‑93 (MQ=255)
                               gTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtg                                                       <  2:2910793/93‑1 (MQ=255)
                               gTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtg                                                       >  1:6450009/1‑93 (MQ=255)
                               gTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtg                                                       >  1:3641108/1‑93 (MQ=255)
                               gTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtg                                                       >  2:5476453/1‑93 (MQ=255)
                               gTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtg                                                       >  2:1873559/1‑93 (MQ=255)
                                tGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgc                                                      <  2:4369843/93‑1 (MQ=255)
                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgcc                                                     >  2:2433441/1‑93 (MQ=255)
                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgcc                                                     >  1:3140364/1‑93 (MQ=255)
                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgcc                                                     >  2:404551/1‑93 (MQ=255)
                                  aCGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgccg                                                    <  2:1261228/93‑1 (MQ=255)
                                  aCGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgccg                                                    >  1:3059626/1‑93 (MQ=255)
                                    gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgcg                                                  >  2:6368066/1‑93 (MQ=255)
                                    gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgcg                                                  >  2:1415127/1‑93 (MQ=255)
                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  1:5979996/1‑93 (MQ=255)
                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  2:184230/1‑93 (MQ=255)
                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                <  2:6035656/93‑1 (MQ=255)
                                       aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               >  2:268058/1‑93 (MQ=255)
                                       aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGCGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               <  2:4849654/93‑1 (MQ=255)
                                         atatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCg                                             >  1:6669205/1‑93 (MQ=255)
                                          tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                                            >  2:3561655/1‑93 (MQ=255)
                                          tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                                            <  1:820894/93‑1 (MQ=255)
                                          tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                                            >  2:2231625/1‑93 (MQ=255)
                                           atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                           >  2:5648440/1‑93 (MQ=255)
                                            ttGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                                          >  2:1714211/1‑93 (MQ=255)
                                            ttGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                                          <  2:5200064/93‑1 (MQ=255)
                                             tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgc                                         <  2:1524246/93‑1 (MQ=255)
                                             tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgc                                         <  1:4262456/93‑1 (MQ=255)
                                               cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       >  1:1040997/1‑93 (MQ=255)
                                               cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       >  2:2179459/1‑93 (MQ=255)
                                               cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       >  2:7467052/1‑93 (MQ=255)
                                                ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  2:6937269/1‑93 (MQ=255)
                                                ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  1:5837655/1‑93 (MQ=255)
                                                 cAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTc                                     >  1:6933747/1‑93 (MQ=255)
                                                    cATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaaa                                  <  2:7193393/93‑1 (MQ=255)
                                                     aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 <  2:6719344/93‑1 (MQ=255)
                                                     aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 >  2:27919/1‑93 (MQ=255)
                                                     aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 <  1:1475298/93‑1 (MQ=255)
                                                     aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 <  2:4215258/93‑1 (MQ=255)
                                                      ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                >  2:4518403/1‑93 (MQ=255)
                                                      ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                <  2:7324013/93‑1 (MQ=255)
                                                        cTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGa                              <  2:5561535/93‑1 (MQ=255)
                                                            ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          <  1:5423690/93‑1 (MQ=255)
                                                              tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        >  2:1440436/1‑93 (MQ=255)
                                                              tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        <  2:1197860/93‑1 (MQ=255)
                                                               aCAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGAt                       <  2:5908695/93‑1 (MQ=255)
                                                                cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATg                      >  1:3936010/1‑93 (MQ=255)
                                                                cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATg                      >  1:7112548/1‑93 (MQ=255)
                                                                 aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     >  1:7274068/1‑93 (MQ=255)
                                                                 aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     >  1:5013002/1‑93 (MQ=255)
                                                                 aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     >  2:811792/1‑93 (MQ=255)
                                                                  aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    <  2:3760822/93‑1 (MQ=255)
                                                                   cTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAg                   >  1:521003/1‑93 (MQ=255)
                                                                    tGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGa                  >  1:4526735/1‑93 (MQ=255)
                                                                    tGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGa                  >  1:1605880/1‑93 (MQ=255)
                                                                       cAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               <  1:4067989/93‑1 (MQ=255)
                                                                       cAGAGGCCGCGCGCCCGGTGCTGCGAGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               <  1:2471026/93‑1 (MQ=255)
                                                                       cAGAAGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               >  1:1342211/1‑93 (MQ=255)
                                                                        agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATc              >  1:1752232/1‑93 (MQ=255)
                                                                         gagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             >  2:5741648/1‑93 (MQ=255)
                                                                            gCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  2:3128416/93‑1 (MQ=255)
                                                                            gCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  2:1591176/93‑1 (MQ=255)
                                                                             ccGCGCGCCCGGTGCTGCGCGGCGAATCGTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcg         <  2:876481/93‑1 (MQ=255)
                                                                              cgcgcgCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcgg        <  1:4340033/93‑1 (MQ=255)
                                                                               gcgcgcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       <  2:6502160/93‑1 (MQ=255)
                                                                                    cccGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  >  1:4654191/1‑93 (MQ=255)
                                                                                             |                                                                                   
CGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCT‑CGCG‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004887‑4005060

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