Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 521 238.2 12299412 96.0% 11807435 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,971,075 2 bp→GC coding (942‑943/1068 nt) rffG → dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,971,0750CG99.6% 903.7 / NA 251C314W (TGC→TGGrffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (105/145);  minor base T (0/1);  total (105/146)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,971,0760GC100.0% 906.8 / NA 251V315L (GTG→CTG) rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (105/146);  total (105/146)

Only 100 of 265 total aligned reads displayed.

TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGCGTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC  >  NC_000913/3970988‑3971162
                                                                                       ||                                                                                      
 gACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGTCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCg                                                                                   <  1:3170413/93‑1 (MQ=255)
 gACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCg                                                                                   >  1:1874804/1‑93 (MQ=255)
   ccTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCa                                                                                 <  1:3425352/93‑1 (MQ=255)
   ccTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTCATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCa                                                                                 >  2:2394430/1‑93 (MQ=255)
    cTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAg                                                                                >  2:5617225/1‑93 (MQ=255)
     tGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAgg                                                                               >  2:46002/1‑93 (MQ=255)
     tGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAgg                                                                               >  2:5786274/1‑93 (MQ=255)
     tGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAgg                                                                               >  1:979374/1‑93 (MQ=255)
      gATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGa                                                                              <  1:196252/93‑1 (MQ=255)
        tCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGaaa                                                                            >  1:5590421/1‑93 (MQ=255)
        tCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAAATTGGCTGGCTGCCGCAGGaac                                                                            >  1:2448453/1‑92 (MQ=255)
         cACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAc                                                                           <  2:4155798/93‑1 (MQ=255)
          aCCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAcc                                                                          <  2:3731343/93‑1 (MQ=255)
           ccTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCt                                                                         <  1:2073705/93‑1 (MQ=255)
            cTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCtt                                                                        >  2:2028905/1‑93 (MQ=255)
            cTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCtt                                                                        <  2:3513308/93‑1 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       <  1:2363050/93‑1 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       >  1:4053925/1‑93 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       <  2:4898749/93‑1 (MQ=255)
                 tCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGaaa                                                                   <  1:5880976/93‑1 (MQ=255)
                 tCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGaaa                                                                   >  1:3378922/1‑93 (MQ=255)
                   gCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGt                                                                 <  1:3974608/93‑1 (MQ=255)
                   gCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGt                                                                 >  1:1106048/1‑93 (MQ=255)
                    cTGACCGTCCGGGGCATGAACTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTg                                                                <  2:4339805/93‑1 (MQ=255)
                     tGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTTCCCGTGAACTTGGCTGGCTGCCGGCGGAAACCTTTGAAAGTgg                                                               >  1:2707488/1‑93 (MQ=255)
                     tGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTgg                                                               <  1:828064/93‑1 (MQ=255)
                     tGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTgg                                                               >  2:3435418/1‑93 (MQ=255)
                      gACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGa                                                              >  1:2013472/1‑93 (MQ=255)
                       gccGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  2:1354213/92‑1 (MQ=255)
                       aCCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  1:618863/93‑1 (MQ=255)
                       aCCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  1:5469066/93‑1 (MQ=255)
                        ccGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAAt                                                            <  1:4174286/93‑1 (MQ=255)
                        ccGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAAt                                                            >  2:3857110/1‑93 (MQ=255)
                          gTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGGATGc                                                          >  1:3598725/1‑93 (MQ=255)
                          gTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGc                                                          <  1:502578/93‑1 (MQ=255)
                           tCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCTGCAGGAAACCTTTGAAAGTGGAATGCg                                                         <  2:3921780/93‑1 (MQ=255)
                            ccGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGt                                                        >  2:1889921/1‑93 (MQ=255)
                            ccGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGt                                                        <  1:1196336/93‑1 (MQ=255)
                             cGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTa                                                       <  2:515545/93‑1 (MQ=255)
                             cGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTa                                                       >  1:1561804/1‑93 (MQ=255)
                               gggCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTaaa                                                     <  1:3248113/93‑1 (MQ=255)
                                ggCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTaaaa                                                    <  1:4629007/93‑1 (MQ=255)
                                 gCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAAc                                                   <  2:2778644/93‑1 (MQ=255)
                                  cATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACg                                                  <  1:5648533/93‑1 (MQ=255)
                                  cATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGACGGAAACCTTTGAAAGTGGAATGCGCaaaagg                                                  >  2:1698602/1‑91 (MQ=255)
                                    tGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGt                                                <  2:3743925/93‑1 (MQ=255)
                                    tGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGt                                                >  2:2763517/1‑93 (MQ=255)
                                     gATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGtt                                               >  2:74270/1‑93 (MQ=255)
                                     gATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCGTTGAAAGTGGAATGCGGAAAACGGtt                                               >  1:5500573/1‑93 (MQ=255)
                                       tCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCa                                             <  2:47413/93‑1 (MQ=255)
                                        cTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAg                                            >  1:2829308/1‑93 (MQ=255)
                                         tGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAACCCTTTGAAAGTGCACTGCGTAAACAGGTTCCGt                                           >  1:2310738/1‑93 (MQ=255)
                                         tGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGt                                           >  1:4807399/1‑93 (MQ=255)
                                         tGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGt                                           <  2:4771448/93‑1 (MQ=255)
                                          gCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTg                                          <  2:1301283/93‑1 (MQ=255)
                                            gTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGt                                        <  1:904990/93‑1 (MQ=255)
                                             ttATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTa                                       <  2:6466/93‑1 (MQ=255)
                                              tATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTAt                                      <  1:2678969/93‑1 (MQ=255)
                                              tATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTAt                                      >  1:1479354/1‑93 (MQ=255)
                                              tATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTAt                                      <  1:5105840/93‑1 (MQ=255)
                                               aTGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATc                                     <  1:4415156/93‑1 (MQ=255)
                                               aTGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATc                                     <  1:4640861/93‑1 (MQ=255)
                                                 gCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTg                                   >  2:5511050/1‑93 (MQ=255)
                                                   cATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTTGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGc                                 >  2:1180638/1‑93 (MQ=255)
                                                   cATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGc                                 >  2:2660842/1‑93 (MQ=255)
                                                   cATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGc                                 >  1:5830430/1‑93 (MQ=255)
                                                    aTTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCt                                >  2:2252962/1‑93 (MQ=255)
                                                    aTTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCt                                >  1:21725/1‑93 (MQ=255)
                                                    aTTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCt                                >  2:5216575/1‑93 (MQ=255)
                                                      tGATGCTTCGAAAATTGCCCGTTAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTaa                              <  1:1704138/93‑1 (MQ=255)
                                                       gATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAAt                             >  2:2936814/1‑93 (MQ=255)
                                                        aTGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTAGCTCATg                            >  2:5403150/1‑93 (MQ=255)
                                                           cTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGaaa                         <  2:3000601/93‑1 (MQ=255)
                                                            ttCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAg                        >  2:5351733/1‑93 (MQ=255)
                                                            ttCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAg                        <  2:2278622/93‑1 (MQ=255)
                                                             tCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGc                       <  1:257105/93‑1 (MQ=255)
                                                               gAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtg                     >  1:5744577/1‑93 (MQ=255)
                                                                aaaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtgg                    <  2:3717246/93‑1 (MQ=255)
                                                                aaaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtgg                    >  1:5210434/1‑93 (MQ=255)
                                                                aaaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtgg                    <  2:245017/93‑1 (MQ=255)
                                                                 aaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggt                   <  2:3193517/93‑1 (MQ=255)
                                                                 aaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggt                   <  2:3599983/93‑1 (MQ=255)
                                                                  aaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggtg                  <  2:2274888/93‑1 (MQ=255)
                                                                  aaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggtg                  >  1:6028191/1‑93 (MQ=255)
                                                                   aTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggtgg                 >  1:114175/1‑93 (MQ=255)
                                                                    ttGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGa                <  2:599936/93‑1 (MQ=255)
                                                                     tGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGaa               <  2:4614730/93‑1 (MQ=255)
                                                                      gCCCGTGAACTTGGCTGTCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAg              <  1:3734448/93‑1 (MQ=255)
                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            <  1:4230413/93‑1 (MQ=255)
                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGGATCTGGCTAATGAAAGCTGGTGGAAGCa            <  1:2990837/93‑1 (MQ=255)
                                                                         cGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAg           <  1:3819375/93‑1 (MQ=255)
                                                                         cGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAg           >  2:952654/1‑93 (MQ=255)
                                                                          gTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          <  1:5932383/93‑1 (MQ=255)
                                                                          gTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          >  1:4237006/1‑93 (MQ=255)
                                                                          gTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          <  2:69928/93‑1 (MQ=255)
                                                                             aaCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTga       >  2:4295349/1‑92 (MQ=255)
                                                                             aaCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGc       >  2:5256613/1‑93 (MQ=255)
                                                                              aCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCa      <  1:1350089/93‑1 (MQ=255)
                                                                                ttGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAgg    >  2:2882194/1‑93 (MQ=255)
                                                                                 tggctggcTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGa   <  1:1235283/93‑1 (MQ=255)
                                                                                       ||                                                                                      
TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGCGTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC  >  NC_000913/3970988‑3971162

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTATCGTGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGCGTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGACG  >  NC_000913/3970982‑3971163
                                                                                             |                                                                                        
CTATCGTGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGG                                                                                          <  SRR4302532.2599152/94‑1 (MQ=60)
     GTGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCC                                                                                     >  SRR4302532.3481032/1‑94 (MQ=60)
     GTGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCC                                                                                     >  SRR4302532.6276873/1‑94 (MQ=60)
      TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCG                                                                                    >  SRR4302532.1936899/1‑94 (MQ=60)
      TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATCCAATTGATGATTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCG                                                                                    >  SRR4302532.6093006/1‑94 (MQ=60)
      TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGAATGATCGGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCG                                                                                    <  SRR4302532.3109178/94‑1 (MQ=60)
       GACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGTCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGC                                                                                   <  SRR4302532.3275713/94‑1 (MQ=60)
       GACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTCAACTTGGCTGGCTGCCGC                                                                                   >  SRR4302532.1022143/1‑94 (MQ=60)
        ACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCA                                                                                  <  SRR4302532.849014/94‑1 (MQ=60)
        ACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCA                                                                                  <  SRR4302532.4265322/94‑1 (MQ=60)
         CCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAG                                                                                 <  SRR4302532.3540107/94‑1 (MQ=60)
          CTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGG                                                                                >  SRR4302532.1012798/1‑94 (MQ=60)
          CTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGG                                                                                >  SRR4302532.5621073/1‑94 (MQ=60)
           TGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGGCGCAGGA                                                                               >  SRR4302532.4829077/1‑94 (MQ=60)
            GATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAA                                                                              >  SRR4302532.4787348/1‑94 (MQ=60)
            GATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAA                                                                              <  SRR4302532.202914/94‑1 (MQ=60)
             ATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAA                                                                             <  SRR4302532.3104237/94‑1 (MQ=60)
             ATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAA                                                                             >  SRR4302532.5784483/1‑94 (MQ=60)
             ATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAAATTGGCTGGCTGCCGCAGGAAC                                                                             >  SRR4302532.2529870/1‑94 (MQ=60)
                ACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCT                                                                          >  SRR4302532.1194762/1‑94 (MQ=60)
                ACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCT                                                                          >  SRR4302532.5781752/1‑94 (MQ=60)
                 CCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTT                                                                         <  SRR4302532.4262093/94‑1 (MQ=60)
                 CCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTT                                                                         <  SRR4302532.6003521/94‑1 (MQ=60)
                 CCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTT                                                                         <  SRR4302532.2142817/94‑1 (MQ=60)
                 CCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTT                                                                         <  SRR4302532.443054/94‑1 (MQ=60)
                  CTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTT                                                                        <  SRR4302532.3645851/94‑1 (MQ=60)
                  CTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTT                                                                        <  SRR4302532.2522959/94‑1 (MQ=60)
                  CTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTT                                                                        >  SRR4302532.4190289/1‑94 (MQ=60)
                   TTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTG                                                                       <  SRR4302532.4175744/94‑1 (MQ=60)
                   TTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTG                                                                       <  SRR4302532.2441871/94‑1 (MQ=60)
                   TTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTG                                                                       <  SRR4302532.2375582/94‑1 (MQ=60)
                   TTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTG                                                                       <  SRR4302532.6021396/94‑1 (MQ=60)
                   TTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTG                                                                       <  SRR4302532.154845/94‑1 (MQ=60)
                   TTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTG                                                                       <  SRR4302532.6090902/94‑1 (MQ=60)
                    TTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGA                                                                      >  SRR4302532.4340040/1‑94 (MQ=60)
                      GTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAA                                                                    >  SRR4302532.3491926/1‑94 (MQ=60)
                       TCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAG                                                                   <  SRR4302532.6084684/94‑1 (MQ=60)
                        agccaagttcacgggaaatttTCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAACTTTGAAAGG                                                                  >  SRR4302532.1638465/22‑94 (MQ=60)
                        CGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGT                                                                  <  SRR4302532.4420343/94‑1 (MQ=60)
                        CGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGT                                                                  >  SRR4302532.1143023/1‑94 (MQ=60)
                        CGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGT                                                                  <  SRR4302532.5267218/94‑1 (MQ=60)
                         GCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTG                                                                 <  SRR4302532.4108483/94‑1 (MQ=60)
                          CTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTTCCCGTGAACTTGGCTGGCTGCCGGCGGAAACCTTTGAAAGTGG                                                                >  SRR4302532.2798533/1‑94 (MQ=60)
                          CTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGG                                                                >  SRR4302532.1286912/1‑94 (MQ=60)
                          CTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGG                                                                <  SRR4302532.2701274/94‑1 (MQ=60)
                          CTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGGGG                                                                >  SRR4302532.5351199/1‑94 (MQ=60)
                           TGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGA                                                               >  SRR4302532.884907/1‑94 (MQ=60)
                           TGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGA                                                               >  SRR4302532.2080623/1‑94 (MQ=60)
                           TGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGA                                                               >  SRR4302532.3675492/1‑94 (MQ=60)
                           TGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGA                                                               <  SRR4302532.856395/94‑1 (MQ=60)
                           TGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGA                                                               >  SRR4302532.2778607/1‑94 (MQ=60)
                             ACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAAT                                                             <  SRR4302532.888098/94‑1 (MQ=60)
                             ACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAAT                                                             <  SRR4302532.640431/94‑1 (MQ=60)
                             ACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAAT                                                             <  SRR4302532.5658675/94‑1 (MQ=60)
                             ACCGTCAGGGTCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCTTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAAT                                                             <  SRR4302532.5941284/94‑1 (MQ=60)
                              CCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATG                                                            <  SRR4302532.5710091/94‑1 (MQ=60)
                              CCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATG                                                            <  SRR4302532.4314287/94‑1 (MQ=60)
                               CGTCCGGGGCATGATCTGGGTGATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGTTGGATGCCGCAGGAAACCTTTGAAAGTGGAATGG                                                           >  SRR4302532.1278785/1‑94 (MQ=60)
                               CGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGGATGC                                                           >  SRR4302532.3719035/1‑94 (MQ=60)
                               CGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGC                                                           <  SRR4302532.1641831/94‑1 (MQ=60)
                               CGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGA                                                           <  SRR4302532.4819956/94‑1 (MQ=60)
                                GTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCG                                                          >  SRR4302532.2433240/1‑94 (MQ=60)
                                GTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCG                                                          >  SRR4302532.3730609/1‑94 (MQ=60)
                                GTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCG                                                          <  SRR4302532.520148/94‑1 (MQ=60)
                                  CCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTA                                                        <  SRR4302532.1236340/94‑1 (MQ=60)
                                  CCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTA                                                        >  SRR4302532.1614627/1‑94 (MQ=60)
                                   CGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAA                                                       >  SRR4302532.1890547/1‑94 (MQ=60)
                                   CGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAA                                                       <  SRR4302532.1144750/94‑1 (MQ=60)
                                   CGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAA                                                       <  SRR4302532.3880279/94‑1 (MQ=60)
                                    GGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAA                                                      >  SRR4302532.2632125/1‑94 (MQ=60)
                                     GGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAA                                                     <  SRR4302532.3355955/94‑1 (MQ=60)
                                      GGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAAC                                                    <  SRR4302532.2863742/94‑1 (MQ=60)
                                      GGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAAC                                                    <  SRR4302532.4786835/94‑1 (MQ=60)
                                      GGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAAC                                                    <  SRR4302532.2104775/94‑1 (MQ=60)
                                      GGCATGATCTGCGTTATGCCATTCATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAAC                                                    <  SRR4302532.4639654/94‑1 (MQ=60)
                                       GCATGATCTGCGTTATGCCATTGATGCTTCGATAATTGCCCGTGATCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACG                                                   <  SRR4302532.6018397/94‑1 (MQ=60)
                                       GCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACG                                                   <  SRR4302532.291078/94‑1 (MQ=60)
                                        CATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGG                                                  >  SRR4302532.2050901/1‑94 (MQ=60)
                                        CATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGG                                                  <  SRR4302532.5845039/94‑1 (MQ=60)
                                          TGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTT                                                >  SRR4302532.6195342/1‑94 (MQ=60)
                                          TGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTT                                                <  SRR4302532.4342018/94‑1 (MQ=60)
                                          GGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCGTTGAAAGTGGAATGCGGAAAACGGTT                                                >  SRR4302532.5691362/1‑94 (MQ=60)
                                            ATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCA                                              <  SRR4302532.3941819/94‑1 (MQ=60)
                                            ATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCA                                              <  SRR4302532.1679777/94‑1 (MQ=60)
                                             TCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAG                                             >  SRR4302532.2923927/1‑94 (MQ=60)
                                              CTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAACCCTTTGAAAGTGCACTGCGTAAAcaggttccgt                                            >  SRR4302532.2387626/1‑84 (MQ=60)
                                              CTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGT                                            >  SRR4302532.4972943/1‑94 (MQ=60)
                                                 CGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGT                                         >  SRR4302532.5786567/1‑94 (MQ=60)
                                                  GTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTA                                        <  SRR4302532.935877/94‑1 (MQ=60)
                                                   TTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTAT                                       >  SRR4302532.1529214/1‑94 (MQ=60)
                                                   TTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTAT                                       <  SRR4302532.2107250/94‑1 (MQ=60)
                                                   TTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTAT                                       <  SRR4302532.3351762/94‑1 (MQ=60)
                                                   TTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTAT                                       <  SRR4302532.1180254/94‑1 (MQ=60)
                                                   TTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCCCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGGAT                                       >  SRR4302532.4973693/1‑94 (MQ=60)
                                                    TATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATC                                      <  SRR4302532.5281372/94‑1 (MQ=60)
                                                    TATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATC                                      <  SRR4302532.2768878/94‑1 (MQ=60)
                                                     ATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCT                                     <  SRR4302532.3958064/94‑1 (MQ=60)
                                                     ATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCT                                     <  SRR4302532.4565550/94‑1 (MQ=60)
                                                     ATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCT                                     <  SRR4302532.4799090/94‑1 (MQ=60)
                                                       GCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGG                                   >  SRR4302532.2368823/1‑94 (MQ=60)
                                                       GCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGG                                   <  SRR4302532.2407188/94‑1 (MQ=60)
                                                        CCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGC                                  >  SRR4302532.6032849/1‑94 (MQ=60)
                                                        CCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGC                                  >  SRR4302532.951432/1‑94 (MQ=60)
                                                         CATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCT                                 >  SRR4302532.22596/1‑94 (MQ=60)
                                                          ATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTA                                >  SRR4302532.3528612/1‑94 (MQ=60)
                                                          ATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTA                                >  SRR4302532.6132438/1‑94 (MQ=60)
                                                            TGATGCTTCGAAAATTGCCCGTTAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAAT                              <  SRR4302532.1761606/94‑1 (MQ=60)
                                                             GATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATG                             >  SRR4302532.123749/1‑94 (MQ=60)
                                                              ATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGA                            >  SRR4302532.1418007/1‑94 (MQ=60)
                                                               TGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAA                           >  SRR4302532.5621188/1‑94 (MQ=60)
                                                                GCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAA                          <  SRR4302532.6052136/94‑1 (MQ=60)
                                                                 CTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAG                         <  SRR4302532.2704439/94‑1 (MQ=60)
                                                                   TCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCT                       <  SRR4302532.266023/94‑1 (MQ=60)
                                                                    CGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTG                      >  SRR4302532.5944560/1‑94 (MQ=60)
                                                                     GAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGG                     >  SRR4302532.5390133/1‑94 (MQ=60)
                                                                      AAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGT                    <  SRR4302532.5282844/94‑1 (MQ=60)
                                                                      AAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGG                    >  SRR4302532.593995/1‑94 (MQ=60)
                                                                       AAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTG                   >  SRR4302532.6236391/1‑94 (MQ=60)
                                                                       AAATTGCCCGGGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTG                   <  SRR4302532.3509526/94‑1 (MQ=60)
                                                                        AATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGG                  <  SRR4302532.2543673/94‑1 (MQ=60)
                                                                        AATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGG                  >  SRR4302532.118094/1‑94 (MQ=60)
                                                                            GCCCGTGAACTTGGCTGTCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGC              <  SRR4302532.3859907/94‑1 (MQ=60)
                                                                              CCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAG            >  SRR4302532.4379160/1‑94 (MQ=60)
                                                                              CCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAG            <  SRR4302532.4372424/94‑1 (MQ=60)
                                                                              CCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGGATCTGGCTAATGAAAGCTGGTGGAAGCAG            <  SRR4302532.3089677/94‑1 (MQ=60)
                                                                               CGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGG           <  SRR4302532.3948317/94‑1 (MQ=60)
                                                                               CGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGG           >  SRR4302532.1441367/1‑94 (MQ=60)
                                                                               CGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGG           >  SRR4302532.4379304/1‑94 (MQ=60)
                                                                               CGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGAGGAAGGCGTAAAACGGTTCAGTGGTATCTGGCTAAGGAAAGCTGGTGGAAGCAGT           >  SRR4302532.117861/1‑94 (MQ=60)
                                                                                GTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGT          <  SRR4302532.2063980/94‑1 (MQ=60)
                                                                                GTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGT          <  SRR4302532.6137733/94‑1 (MQ=60)
                                                                                 TGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTG         >  SRR4302532.4058859/1‑94 (MQ=60)
                                                                                  GAACTTGGCTGGCTGCCGCAGGAAGCCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGC        <  SRR4302532.567868/94‑1 (MQ=60)
                                                                                  GAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGC        >  SRR4302532.2242116/1‑94 (MQ=60)
                                                                                   AACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCA       <  SRR4302532.6226629/94‑1 (MQ=60)
                                                                                    ACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAG      <  SRR4302532.1395229/94‑1 (MQ=60)
                                                                                    ACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAG      <  SRR4302532.792058/94‑1 (MQ=60)
                                                                                    ACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAG      <  SRR4302532.2943427/94‑1 (MQ=60)
                                                                                     CTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGG     <  SRR4302532.1392231/94‑1 (MQ=60)
                                                                                      TTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGA    >  SRR4302532.6118903/1‑94 (MQ=60)
                                                                                      TTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGA    >  SRR4302532.2942077/1‑94 (MQ=60)
                                                                                       TGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC   <  SRR4302532.1276602/94‑1 (MQ=60)
                                                                                       TGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC   <  SRR4302532.4731404/94‑1 (MQ=60)
                                                                                       TGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC   <  SRR4302532.5137562/94‑1 (MQ=60)
                                                                                       TGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC   >  SRR4302532.3189692/1‑94 (MQ=60)
                                                                                        GGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGACG  <  SRR4302532.1816304/94‑1 (MQ=60)
                                                                                        GGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGACG  <  SRR4302532.3285518/94‑1 (MQ=60)
                                                                                             |                                                                                        
CTATCGTGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGCGTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGACG  >  NC_000913/3970982‑3971163

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