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.G99.5% 569.2 / NA 188A505A (GCG→GCArecQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base G (93/93);  minor base A (0/1);  total (93/95)
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.99e-01
*NC_0009134,004,9832.C99.5% 575.8 / ‑4.0 189coding (1515/1833 nt)recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (93/95);  total (93/96)
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.41e-01

Only 100 of 210 total aligned reads displayed.

CCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTCGC‑‑G‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004893‑4005060
                                                                                                ||                                                                             
ccGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑g‑‑g                                                                              <  1:2720483/93‑1 (MQ=255)
ccGCCAGATGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑g‑‑g                                                                              <  1:2898594/93‑1 (MQ=255)
   ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑g‑‑gcga                                                                           >  2:311149/1‑90 (MQ=255)
   ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑g‑‑gcga                                                                           >  1:3649080/1‑90 (MQ=255)
    cAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑g‑‑gcgaa                                                                          >  2:2888600/1‑89 (MQ=255)
    cAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑g‑‑gcgaa                                                                          >  1:985558/1‑89 (MQ=255)
          aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCctct                                                                      <  1:2685500/93‑1 (MQ=255)
           ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCtt                                                                     >  1:2153598/1‑93 (MQ=255)
              aCCTCGGCCTGGTGTCGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTCCTGCCC‑‑GGCGAATCCTCTTTGc                                                                  >  2:3555238/1‑93 (MQ=255)
              aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGc                                                                  >  1:2538633/1‑93 (MQ=255)
               ccTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCa                                                                 <  1:603750/93‑1 (MQ=255)
                cTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCaa                                                                <  1:1874533/93‑1 (MQ=255)
                cTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCAC‑‑GGCGAATCCTCTTTGCaa                                                                <  2:4578636/93‑1 (MQ=255)
                 tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGG‑‑GAATCCTCTTTGCAAc                                                               >  2:4110046/1‑93 (MQ=255)
                 tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAAc                                                               >  1:3274908/1‑93 (MQ=255)
                 tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAAc                                                               <  1:3807126/93‑1 (MQ=255)
                   ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACtt                                                             <  2:2665299/93‑1 (MQ=255)
                    gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTtg                                                            >  2:2054614/1‑93 (MQ=255)
                    gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTtg                                                            >  1:1516186/1‑93 (MQ=255)
                     ccTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTtgc                                                           <  2:2250776/93‑1 (MQ=255)
                      cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTtgcc                                                          >  1:2314699/1‑93 (MQ=255)
                      cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTtgcc                                                          <  2:4337371/93‑1 (MQ=255)
                       tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTtgccg                                                         >  1:1732575/1‑93 (MQ=255)
                       tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTtgccg                                                         <  2:2899863/93‑1 (MQ=255)
                       tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTtgccg                                                         >  1:3066609/1‑93 (MQ=255)
                        ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTtgccgt                                                        >  1:44171/1‑93 (MQ=255)
                         gTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTtgccgtg                                                       >  2:917263/1‑93 (MQ=255)
                             cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCcgc                                                   >  2:11/1‑93 (MQ=255)
                              gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCcgcg                                                  >  2:3206024/1‑93 (MQ=255)
                              gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCcgcg                                                  >  1:397388/1‑93 (MQ=255)
                               cAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGt                                                 <  1:1556481/93‑1 (MQ=255)
                                aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  1:1622143/1‑93 (MQ=255)
                                 aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               >  2:3777329/1‑93 (MQ=255)
                                  aaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATc                                              <  2:2576013/93‑1 (MQ=255)
                                  aaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATc                                              >  2:1577305/1‑93 (MQ=255)
                                     atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                           >  2:1794647/1‑93 (MQ=255)
                                     atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                           <  1:636805/93‑1 (MQ=255)
                                      ttGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                                          >  2:2520985/1‑93 (MQ=255)
                                      ttGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGAGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                                          <  1:2512785/93‑1 (MQ=255)
                                         cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       >  1:2101561/1‑93 (MQ=255)
                                          ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  1:1090251/1‑93 (MQ=255)
                                           cAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTc                                     >  1:4667810/1‑93 (MQ=255)
                                           cAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTc                                     >  2:3447522/1‑93 (MQ=255)
                                              cATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaaa                                  <  2:804879/93‑1 (MQ=255)
                                               aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 <  2:3650070/93‑1 (MQ=255)
                                               aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 >  1:2308255/1‑93 (MQ=255)
                                                ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                <  2:3179396/93‑1 (MQ=255)
                                                 tCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCC‑‑GCGAATCCTCTGTGCAACTTGCCGTGCCGCGTATCGTGACGCTCAAAcca                               >  2:38655/1‑92 (MQ=255)
                                                 tCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCg                               >  2:1983702/1‑93 (MQ=255)
                                                  cTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGa                              >  2:3072394/1‑93 (MQ=255)
                                                  cTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGa                              >  2:3083271/1‑93 (MQ=255)
                                                  cTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGa                              >  1:2492350/1‑93 (MQ=255)
                                                   tGCCCTACACCGGACAGAGGCCGCGCTCCCGGTGCTGCGC‑‑GACGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaa                             <  2:2501729/93‑1 (MQ=255)
                                                   tGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaa                             >  2:1446960/1‑93 (MQ=255)
                                                    gCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTTCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa                            >  1:4301599/1‑93 (MQ=255)
                                                    gCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa                            >  2:2801709/1‑93 (MQ=255)
                                                     cccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  1:3536085/93‑1 (MQ=255)
                                                      ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          <  2:1603966/93‑1 (MQ=255)
                                                      ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          <  2:3710452/93‑1 (MQ=255)
                                                      ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          <  1:1112549/93‑1 (MQ=255)
                                                      ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  1:3410768/1‑93 (MQ=255)
                                                      ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  1:293277/1‑93 (MQ=255)
                                                       cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg                         >  2:1352725/1‑93 (MQ=255)
                                                       cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg                         >  2:2100698/1‑93 (MQ=255)
                                                       cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg                         <  2:238665/93‑1 (MQ=255)
                                                        tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        >  1:705105/1‑93 (MQ=255)
                                                        tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        <  1:1931694/93‑1 (MQ=255)
                                                         acacCTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGAt                       <  2:4231789/93‑1 (MQ=255)
                                                         aCAACTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGAt                       <  2:3896383/93‑1 (MQ=255)
                                                          cAACTGACAGAGGCCGCGCGCCCGGTGCTGAGC‑‑GGCGAATACTCTTTGCAACTTGACGTGCCGCGCATCGTGGCGCTCAAACCGAAAGCGATg                      >  2:2872550/1‑93 (MQ=255)
                                                           aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GACGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTAAAACCGAAAGCGATGc                     <  2:2036718/93‑1 (MQ=255)
                                                            aCTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCGGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    >  2:561852/1‑93 (MQ=255)
                                                            aCTGACAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    >  1:3468117/1‑93 (MQ=255)
                                                                aCAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaaa                <  1:2762237/93‑1 (MQ=255)
                                                                aCAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaaa                >  2:1505923/1‑93 (MQ=255)
                                                                aCAGAGGCCGCGCGCCCGGTCCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaaa                <  2:1242535/93‑1 (MQ=255)
                                                                 cAGAGGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               <  1:3322027/93‑1 (MQ=255)
                                                                  agagGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATc              >  1:4230880/1‑93 (MQ=255)
                                                                   gagGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATGGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             <  1:868770/93‑1 (MQ=255)
                                                                   gagGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             <  2:2200322/93‑1 (MQ=255)
                                                                   gagGCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             <  2:761356/93‑1 (MQ=255)
                                                                     ggCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           <  2:4529505/93‑1 (MQ=255)
                                                                     ggCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           >  1:1461741/1‑93 (MQ=255)
                                                                     ggCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           >  1:2592817/1‑93 (MQ=255)
                                                                      gCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  2:1167481/93‑1 (MQ=255)
                                                                      gCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  1:2360333/93‑1 (MQ=255)
                                                                      gCCGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          >  2:4069755/1‑93 (MQ=255)
                                                                       ccGCGCGCCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcg         <  1:3009913/93‑1 (MQ=255)
                                                                         gggcgcCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       <  1:2302815/91‑1 (MQ=255)
                                                                         gcgcgcCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       <  2:3006694/93‑1 (MQ=255)
                                                                         gcgcgcCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       <  2:634105/93‑1 (MQ=255)
                                                                         gcgcgcCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       <  2:1436857/93‑1 (MQ=255)
                                                                         gcgcgcCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       >  1:4187246/1‑93 (MQ=255)
                                                                         gcgcgcCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       >  2:467985/1‑93 (MQ=255)
                                                                          cgcgcCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggcg      <  1:1020940/93‑1 (MQ=255)
                                                                            cgcCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCggcggc    >  2:3520779/1‑93 (MQ=255)
                                                                            cgcCCGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCggcggc    <  1:3839777/93‑1 (MQ=255)
                                                                             ccccGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCa   <  1:3318750/92‑1 (MQ=255)
                                                                              cccGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  >  1:4052102/1‑93 (MQ=255)
                                                                              cccGGTGCTGCGC‑‑GGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  >  1:4500301/1‑93 (MQ=255)
                                                                                                ||                                                                             
CCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTCGC‑‑G‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004893‑4005060

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