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.1% 769.2 / ‑1.0 226coding (1515/1833 nt)recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (126/97);  total (126/100)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01
*NC_0009134,004,9832.C98.7% 771.9 / 1.8 226coding (1515/1833 nt)recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (126/97);  total (126/100)
Fisher's exact test for biased strand distribution p-value = 8.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 242 total aligned reads displayed.

CCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTCGCG‑‑G‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004893‑4005060
                                                                                            ||                                                                               
ccGCCAGCTGATTCACCTCGGCCTTGCGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTCCTGCGCg‑‑g                                                                                >  1:5322054/1‑93 (MQ=255)
ccGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑g                                                                                <  2:4085213/93‑1 (MQ=255)
 cGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gc                                                                               <  1:5825626/93‑2 (MQ=255)
  gCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGACGCGCGCCCGGTGCTGCGCg‑‑gcg                                                                              >  1:2076544/1‑91 (MQ=255)
   ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcga                                                                             >  2:2635384/1‑90 (MQ=255)
   ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcga                                                                             >  2:2818053/1‑90 (MQ=255)
    cAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcgaa                                                                            >  1:1120973/1‑89 (MQ=255)
          aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCctct                                                                      >  2:2990603/1‑93 (MQ=255)
           ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCtt                                                                     <  2:5799061/93‑1 (MQ=255)
            tCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCttt                                                                    >  1:3408637/1‑93 (MQ=255)
             cACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGACGTCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTg                                                                   >  1:1525301/1‑93 (MQ=255)
               ccTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCa                                                                 >  1:76520/1‑93 (MQ=255)
                  gggCCTGGTGACGCAAAATATTGCGCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACt                                                              <  1:3477166/92‑1 (MQ=255)
                   ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACtt                                                             <  1:1503351/93‑1 (MQ=255)
                   ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACtt                                                             <  2:3224096/93‑1 (MQ=255)
                   ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACtt                                                             <  2:2520177/93‑1 (MQ=255)
                      cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgcc                                                          >  2:2007408/1‑93 (MQ=255)
                      cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgcc                                                          >  2:153339/1‑93 (MQ=255)
                       tGGTGACGCAAAATATTGCCCATCATTCTGCCCTACAACAGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccg                                                         <  1:4169618/93‑1 (MQ=255)
                       tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccg                                                         >  2:193799/1‑93 (MQ=255)
                       tGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccg                                                         <  2:588515/93‑1 (MQ=255)
                        ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgt                                                        >  1:1794270/1‑93 (MQ=255)
                        ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgt                                                        >  1:2201952/1‑93 (MQ=255)
                        ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgt                                                        <  2:4210885/93‑1 (MQ=255)
                        ggTGACGCAAAATATGGCCCATCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgt                                                        <  1:3157119/93‑1 (MQ=255)
                          tGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgc                                                      <  1:6094738/93‑1 (MQ=255)
                          tGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgc                                                      <  2:3050715/93‑1 (MQ=255)
                           gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgcc                                                     >  2:4675464/1‑93 (MQ=255)
                            aCGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgccg                                                    >  1:4502750/1‑93 (MQ=255)
                             cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCcgc                                                   <  1:954626/93‑1 (MQ=255)
                             cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCcgc                                                   <  2:5517915/93‑1 (MQ=255)
                             cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCcgc                                                   <  1:4545938/93‑1 (MQ=255)
                              gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCcgcg                                                  <  1:4057141/93‑1 (MQ=255)
                                 aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               <  2:2303974/93‑1 (MQ=255)
                                  aaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATc                                              >  2:4036507/1‑93 (MQ=255)
                                   atatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCg                                             >  2:3338050/1‑93 (MQ=255)
                                    tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                                            <  1:1706283/93‑1 (MQ=255)
                                     atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                           <  1:5661739/93‑1 (MQ=255)
                                      ttGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                                          <  2:4849035/93‑1 (MQ=255)
                                      ttGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                                          >  2:1453881/1‑93 (MQ=255)
                                       tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgc                                         <  1:4126561/93‑1 (MQ=255)
                                       tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgc                                         <  1:3727594/93‑1 (MQ=255)
                                       tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgc                                         <  2:5154474/93‑1 (MQ=255)
                                        gCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcg                                        >  1:5097392/1‑93 (MQ=255)
                                          ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  1:4308571/1‑93 (MQ=255)
                                          ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  1:3893390/1‑93 (MQ=255)
                                          ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  1:3112196/1‑93 (MQ=255)
                                          ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  1:2583600/1‑93 (MQ=255)
                                           cAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTc                                     >  2:4390949/1‑93 (MQ=255)
                                           cAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTc                                     >  1:1013467/1‑93 (MQ=255)
                                            cgCATTCTCCCCCAGAACTGACAGAGGCCGCGCCCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                                    <  2:1612107/92‑1 (MQ=255)
                                            aGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                                    >  1:685176/1‑93 (MQ=255)
                                            aGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                                    <  1:1246157/93‑1 (MQ=255)
                                            aGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                                    >  2:6018738/1‑93 (MQ=255)
                                             gCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaa                                   <  2:5545321/93‑1 (MQ=255)
                                               aTTCTGCCCTACAACTGACAGAGGGCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 >  1:6226757/1‑93 (MQ=255)
                                               aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 <  1:5592055/93‑1 (MQ=255)
                                                ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                >  2:2156682/1‑93 (MQ=255)
                                                ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                <  1:545846/93‑1 (MQ=255)
                                                 tCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCg                               <  2:775790/93‑1 (MQ=255)
                                                  cTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGa                              >  2:3823796/1‑93 (MQ=255)
                                                   tGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaa                             >  2:3259203/1‑93 (MQ=255)
                                                    gCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCG‑‑GCTAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa                            <  1:1913622/93‑1 (MQ=255)
                                                     cccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  2:640998/93‑1 (MQ=255)
                                                      ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          <  2:2905697/93‑1 (MQ=255)
                                                      ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          <  2:6187260/93‑1 (MQ=255)
                                                       cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg                         <  2:5475407/93‑1 (MQ=255)
                                                        tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        <  2:3732020/93‑1 (MQ=255)
                                                        tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        <  1:2300997/93‑1 (MQ=255)
                                                         aCAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGAt                       >  1:5105000/1‑93 (MQ=255)
                                                           aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     <  1:3588192/93‑1 (MQ=255)
                                                            aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    >  2:2049498/1‑93 (MQ=255)
                                                            aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    >  2:6428107/1‑93 (MQ=255)
                                                            aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    >  2:4486804/1‑93 (MQ=255)
                                                               gACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaa                 >  2:501401/1‑93 (MQ=255)
                                                                aCAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaaa                <  2:3377268/93‑1 (MQ=255)
                                                                aCAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaaa                >  2:860837/1‑93 (MQ=255)
                                                                 cAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               <  2:3524281/93‑1 (MQ=255)
                                                                 cAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               <  2:4994591/93‑1 (MQ=255)
                                                                  agagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATc              <  1:495865/93‑1 (MQ=255)
                                                                   gagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             <  2:5714048/93‑1 (MQ=255)
                                                                   gagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             >  2:1637728/1‑93 (MQ=255)
                                                                   gagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAGCg             >  1:2272424/1‑93 (MQ=255)
                                                                    agGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGt            >  2:1703361/1‑93 (MQ=255)
                                                                    agGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGt            <  1:2679057/93‑1 (MQ=255)
                                                                     ggCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           >  2:4914797/1‑93 (MQ=255)
                                                                     ggCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           <  1:5178181/93‑1 (MQ=255)
                                                                     ggCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           <  2:1285284/93‑1 (MQ=255)
                                                                     ggCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           >  2:1860798/1‑93 (MQ=255)
                                                                      gCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  2:1054313/93‑1 (MQ=255)
                                                                       ccGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcg         <  1:2982148/93‑1 (MQ=255)
                                                                        cgcgcgCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcgg        >  1:536174/1‑93 (MQ=255)
                                                                         gcgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       >  2:5294461/1‑93 (MQ=255)
                                                                         gcgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       >  1:5876755/1‑93 (MQ=255)
                                                                          cgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggcg      <  2:6055417/93‑1 (MQ=255)
                                                                           gcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggcgg     <  2:3060163/93‑1 (MQ=255)
                                                                             gcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCa   >  1:4668474/1‑93 (MQ=255)
                                                                             gcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCa   >  2:5970097/1‑93 (MQ=255)
                                                                              cccGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  >  1:4400735/1‑93 (MQ=255)
                                                                              cccGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGCGGCGCTCAAACCGCAAGCGATGCAGAAATCGTTCGGCGGGaa  >  2:2751905/1‑93 (MQ=255)
                                                                                            ||                                                                               
CCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTCGCG‑‑G‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004893‑4005060

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