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

Only 100 of 218 total aligned reads displayed.

GTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTCGCG‑‑G‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004888‑4005060
                                                                                            |                                                                                      
gTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtg                                                                                        >  2:923177/1‑92 (MQ=255)
gTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtg                                                                                        <  2:3406353/93‑2 (MQ=255)
gTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtg                                                                                        <  2:4352251/93‑2 (MQ=255)
 tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgc                                                                                       >  1:3366189/1‑91 (MQ=255)
   aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                     >  1:3329809/1‑89 (MQ=255)
   aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                     <  1:3097872/93‑5 (MQ=255)
      cGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gc                                                                               <  2:2127759/93‑2 (MQ=255)
       gCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcg                                                                              >  1:4002730/1‑91 (MQ=255)
        ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcga                                                                             <  1:4656267/93‑4 (MQ=255)
        ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcga                                                                             >  2:311149/1‑90 (MQ=255)
        ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcga                                                                             >  2:4103545/1‑90 (MQ=255)
         cAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcgaa                                                                            >  1:985558/1‑89 (MQ=255)
                ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCtt                                                                     >  1:2153598/1‑93 (MQ=255)
                 tCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCttt                                                                    <  1:1831836/93‑1 (MQ=255)
                  cACCTCGGCCTGGTGACGCAAAATATCGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTg                                                                   <  2:3709250/93‑1 (MQ=255)
                   aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGc                                                                  >  2:158430/1‑93 (MQ=255)
                   aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGc                                                                  >  1:2538633/1‑93 (MQ=255)
                     cTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCaa                                                                <  1:2731995/93‑1 (MQ=255)
                     cTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCACGGCG‑‑AATCCTCTTTGCaa                                                                <  2:4578636/93‑1 (MQ=255)
                      tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAAc                                                               >  1:3772551/1‑93 (MQ=255)
                      tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAAc                                                               >  1:3274908/1‑93 (MQ=255)
                        ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACtt                                                             <  2:2665299/93‑1 (MQ=255)
                        ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACtt                                                             <  2:4202152/93‑1 (MQ=255)
                         gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtg                                                            >  1:1516186/1‑93 (MQ=255)
                         gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtg                                                            >  2:2054614/1‑93 (MQ=255)
                         gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtg                                                            >  1:1899526/1‑93 (MQ=255)
                         gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGAG‑‑AATCCTCTTTGAAACTtg                                                            >  2:740478/1‑93 (MQ=255)
                           cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgcc                                                          >  1:4438395/1‑93 (MQ=255)
                           cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgcc                                                          >  1:3828734/1‑93 (MQ=255)
                                 aCGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgccg                                                    <  1:2308603/93‑1 (MQ=255)
                                  cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCcgc                                                   >  2:11/1‑93 (MQ=255)
                                    cAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGt                                                 <  1:1556481/93‑1 (MQ=255)
                                      aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGTCGCGTAt                                               <  1:2318751/93‑1 (MQ=255)
                                      aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               >  2:3777329/1‑93 (MQ=255)
                                       aaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATc                                              >  1:3062488/1‑93 (MQ=255)
                                       aaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATc                                              <  2:2576013/93‑1 (MQ=255)
                                          atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                           <  1:636805/93‑1 (MQ=255)
                                           ttGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                                          <  1:4103545/93‑1 (MQ=255)
                                            tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCATCTTGCCGTGCCGCGTATCGTGgc                                         <  2:790747/93‑1 (MQ=255)
                                            tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgc                                         <  2:2217808/93‑1 (MQ=255)
                                             gCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcg                                        >  2:2009476/1‑93 (MQ=255)
                                              cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       >  1:2101561/1‑93 (MQ=255)
                                               ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  1:1090251/1‑93 (MQ=255)
                                                cAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTc                                     >  2:2799823/1‑93 (MQ=255)
                                                  gCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaa                                   >  1:268601/1‑93 (MQ=255)
                                                   cATTCTGCCATACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaaa                                  <  1:4257405/93‑1 (MQ=255)
                                                    aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCTGCGTATCGTGGCGCTCAAAc                                 <  2:2384219/93‑1 (MQ=255)
                                                    aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 <  2:4322019/93‑1 (MQ=255)
                                                    aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 >  1:285458/1‑93 (MQ=255)
                                                     ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                <  2:3179396/93‑1 (MQ=255)
                                                      tCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCCGCG‑‑AATCCTCTGTGCAACTTGCCGTGCCGCGTATCGTGACGCTCAAAcca                               >  2:38655/1‑92 (MQ=255)
                                                       cTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGa                              >  2:3083271/1‑93 (MQ=255)
                                                       cTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGa                              >  2:3072394/1‑93 (MQ=255)
                                                        tGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaa                             >  2:1446960/1‑93 (MQ=255)
                                                          cccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  1:3536085/93‑1 (MQ=255)
                                                           ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  1:3410768/1‑93 (MQ=255)
                                                           ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  1:293277/1‑93 (MQ=255)
                                                           ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  2:2891754/1‑93 (MQ=255)
                                                           ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          <  1:1112549/93‑1 (MQ=255)
                                                            cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg                         >  2:1352725/1‑93 (MQ=255)
                                                             tACAACTGACAGCGGCCGCGCGCCCGGTGCTGCGCG‑‑GCAAATCCTCTTTGCAACTTGCCGTGCGGCGTATCGTGGCGCTCAAACCGAAAGCGa                        <  2:3353017/93‑1 (MQ=255)
                                                             tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        >  2:169488/1‑93 (MQ=255)
                                                             tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        >  1:705105/1‑93 (MQ=255)
                                                              aCAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGAt                       <  2:3896383/93‑1 (MQ=255)
                                                              aCAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGAt                       <  1:2360563/93‑1 (MQ=255)
                                                              aCAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGAt                       <  2:3448518/93‑1 (MQ=255)
                                                                aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     <  2:587081/93‑1 (MQ=255)
                                                                aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGACG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTAAAACCGAAAGCGATGc                     <  2:2036718/93‑1 (MQ=255)
                                                                aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCG‑‑GCTAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     <  1:2305248/93‑1 (MQ=255)
                                                                 aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    >  1:3468117/1‑93 (MQ=255)
                                                                 aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    <  2:3524250/93‑1 (MQ=255)
                                                                 aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    >  1:3908659/1‑93 (MQ=255)
                                                                   tGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGa                  <  1:3321294/93‑1 (MQ=255)
                                                                   tGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGa                  >  2:869544/1‑93 (MQ=255)
                                                                    gACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaa                 >  2:3588946/1‑93 (MQ=255)
                                                                    gACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaa                 >  2:24557/1‑93 (MQ=255)
                                                                     aCAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaaa                >  2:1505923/1‑93 (MQ=255)
                                                                     aCAGAGGCCGCGCGCCCGGTCCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaaa                <  2:1242535/93‑1 (MQ=255)
                                                                       agagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATc              <  1:4357365/93‑1 (MQ=255)
                                                                       agagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATc              >  1:4230880/1‑93 (MQ=255)
                                                                        gagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATGGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             <  1:868770/93‑1 (MQ=255)
                                                                        gagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             <  1:1972920/93‑1 (MQ=255)
                                                                         agGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGt            <  2:2821227/93‑1 (MQ=255)
                                                                          ggCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           <  1:266656/93‑1 (MQ=255)
                                                                          ggCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           >  1:2592817/1‑93 (MQ=255)
                                                                          ggCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           >  1:1461741/1‑93 (MQ=255)
                                                                          ggCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           <  2:4529505/93‑1 (MQ=255)
                                                                           gCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  2:2517610/93‑1 (MQ=255)
                                                                           gCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          >  2:4069755/1‑93 (MQ=255)
                                                                           gCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  2:1167481/93‑1 (MQ=255)
                                                                           gCCCCGCGCACGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  1:1636489/93‑1 (MQ=255)
                                                                            ccGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcg         >  1:2497850/1‑93 (MQ=255)
                                                                            ccGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcg         <  1:3009913/93‑1 (MQ=255)
                                                                            ccGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcg         <  2:1887946/93‑1 (MQ=255)
                                                                              gcgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       <  2:3006694/93‑1 (MQ=255)
                                                                              gcgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       >  2:467985/1‑93 (MQ=255)
                                                                              gcgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       <  2:1436857/93‑1 (MQ=255)
                                                                              gcgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       >  2:2882959/1‑93 (MQ=255)
                                                                               cgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggcg      <  1:1020940/93‑1 (MQ=255)
                                                                               cgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCACTCAAACCGAAAGCGATGCAGAAATCGTTcggcg      <  1:599894/93‑1 (MQ=255)
                                                                                            |                                                                                      
GTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTCGCG‑‑G‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004888‑4005060

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