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 I2 R1 884 257.2 12958910 96.9% 12557183 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,987,586 2 bp→GC coding (814‑815/963 nt) hemC ← porphobilinogen deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,5860CG100.0% 730.5 / NA 214R272A (CGG→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (129/85);  total (129/85)
*NC_0009133,987,5870GC100.0% 724.7 / NA 214R272A (CGG→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (129/85);  total (129/85)

Only 100 of 229 total aligned reads displayed.

AGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCCGGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGG  >  NC_000913/3987499‑3987672
                                                                                       ||                                                                                     
aGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGgcgcgccg                                                                                   <  2:855574/93‑1 (MQ=255)
  ctctTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAc                                                                                 >  2:3278556/1‑93 (MQ=255)
     ttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAg                                                                              >  2:5834634/1‑93 (MQ=255)
     ttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAg                                                                              <  2:3039435/93‑1 (MQ=255)
     ttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAg                                                                              >  2:3095386/1‑93 (MQ=255)
        tGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgcgc                                                                           >  2:3297906/1‑93 (MQ=255)
           cAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacg                                                                        >  2:3792186/1‑93 (MQ=255)
               gAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGc                                                                    >  1:2994225/1‑93 (MQ=255)
               gAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGc                                                                    >  1:3161651/1‑93 (MQ=255)
                aaaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGcc                                                                   <  2:1859226/93‑1 (MQ=255)
                 aaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCa                                                                  >  1:4145598/1‑93 (MQ=255)
                 aaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCa                                                                  >  2:2761787/1‑93 (MQ=255)
                 aaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACAGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCCGCCa                                                                  >  1:3255049/1‑93 (MQ=255)
                  aTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAg                                                                 >  1:5164365/1‑93 (MQ=255)
                   tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGTGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGa                                                                >  1:2754204/1‑93 (MQ=255)
                   tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGa                                                                <  1:1821773/93‑1 (MQ=255)
                   tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGa                                                                <  2:6180384/93‑1 (MQ=255)
                   tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGa                                                                <  2:515806/93‑1 (MQ=255)
                   tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGa                                                                >  1:3224504/1‑93 (MQ=255)
                   tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGGTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAgt                                                                >  2:5388791/1‑92 (MQ=255)
                    ccccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAt                                                               >  1:4008279/1‑93 (MQ=255)
                     cccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAtt                                                              >  2:2046893/1‑93 (MQ=255)
                      ccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAttt                                                             >  2:6387892/1‑93 (MQ=255)
                       cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTc                                                            <  1:1343413/93‑1 (MQ=255)
                       cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTCCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTc                                                            <  2:4527076/93‑1 (MQ=255)
                        aTTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCg                                                           >  2:5349453/1‑93 (MQ=255)
                        aTTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCg                                                           >  1:3769702/1‑93 (MQ=255)
                        aTTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCg                                                           >  1:5735161/1‑93 (MQ=255)
                         tttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGc                                                          <  1:4868267/93‑1 (MQ=255)
                            gTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCAt                                                       >  2:1202362/1‑93 (MQ=255)
                             ttCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATc                                                      >  2:2542872/1‑93 (MQ=255)
                               cGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCaa                                                    <  1:2783969/93‑1 (MQ=255)
                                ggCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAAt                                                   >  2:4138188/1‑93 (MQ=255)
                                 gCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATa                                                  >  1:461435/1‑93 (MQ=255)
                                 gCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATa                                                  >  1:5868021/1‑93 (MQ=255)
                                   aTCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCCGCCAGATTTCGCCATCAATAAg                                                >  2:5093101/1‑93 (MQ=255)
                                    tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGc                                               <  1:6410810/93‑1 (MQ=255)
                                    tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGc                                               >  1:506288/1‑93 (MQ=255)
                                      ttGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTc                                             <  2:1381775/93‑1 (MQ=255)
                                         cggcgCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGc                                          <  2:1508593/93‑1 (MQ=255)
                                           gcgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGt                                        <  1:4272286/93‑1 (MQ=255)
                                           gcgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGt                                        >  2:790485/1‑93 (MQ=255)
                                            cgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTa                                       <  2:2288223/93‑1 (MQ=255)
                                            cgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTa                                       <  2:5643611/93‑1 (MQ=255)
                                             gcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAg                                      >  2:611827/1‑93 (MQ=255)
                                              cACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGc                                     <  1:5302764/93‑1 (MQ=255)
                                               aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGca                                    >  2:5244321/1‑92 (MQ=255)
                                               aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                    <  1:5668758/93‑1 (MQ=255)
                                               aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                    >  1:2351286/1‑93 (MQ=255)
                                               aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                    <  1:2737150/93‑1 (MQ=255)
                                               aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                    >  2:5942487/1‑93 (MQ=255)
                                                 cgcgGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTAc                                  <  1:6394441/93‑1 (MQ=255)
                                                  gcggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTAcc                                 >  2:4391841/1‑93 (MQ=255)
                                                  gcggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTAcc                                 <  2:2368903/93‑1 (MQ=255)
                                                    ggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCaa                               <  1:4826093/93‑1 (MQ=255)
                                                    ggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCaa                               >  1:5251913/1‑93 (MQ=255)
                                                       gTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTg                            <  1:1926761/93‑1 (MQ=255)
                                                        ttCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTgg                           >  1:1560832/1‑93 (MQ=255)
                                                        ttCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTgg                           >  1:5402393/1‑93 (MQ=255)
                                                         tCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGc                          >  2:4153143/1‑93 (MQ=255)
                                                            ccGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAcc                       <  1:6055077/93‑1 (MQ=255)
                                                             cgcgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCt                      <  1:1108688/93‑1 (MQ=255)
                                                              gcgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTg                     >  2:4697614/1‑93 (MQ=255)
                                                               cgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGa                    <  1:151794/93‑1 (MQ=255)
                                                               cgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGa                    >  1:353432/1‑93 (MQ=255)
                                                                gAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGAc                   <  2:959095/93‑1 (MQ=255)
                                                                gAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGAc                   <  1:3860941/93‑1 (MQ=255)
                                                                 aataatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCTCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACa                  >  1:2548746/1‑93 (MQ=255)
                                                                 aataatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACa                  >  2:1618979/1‑93 (MQ=255)
                                                                 aataatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACa                  <  1:952880/93‑1 (MQ=255)
                                                                   taatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATc                <  1:832174/93‑1 (MQ=255)
                                                                   taatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATc                >  1:6329357/1‑93 (MQ=255)
                                                                    aatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATcc               <  1:3067712/93‑1 (MQ=255)
                                                                    aatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATcc               <  1:3595674/93‑1 (MQ=255)
                                                                      tCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGc             >  2:1979956/1‑93 (MQ=255)
                                                                      tCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGc             <  2:4634883/93‑1 (MQ=255)
                                                                      acTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGc             <  1:4188448/92‑1 (MQ=255)
                                                                       cTGCGAACCGGCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGcc            >  1:3932707/1‑93 (MQ=255)
                                                                        tGCGAACCGTCCGGCTCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAAGAAGCTCGGCGTAGCTACCAATGGGCACCTGACATCTGCCt           <  1:3471758/93‑1 (MQ=255)
                                                                         gCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCtt          <  2:3471161/93‑1 (MQ=255)
                                                                          cGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc         <  1:409878/93‑1 (MQ=255)
                                                                          cGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc         <  1:5461227/93‑1 (MQ=255)
                                                                          cGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc         >  2:3716086/1‑93 (MQ=255)
                                                                          cGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc         <  1:2330830/93‑1 (MQ=255)
                                                                          cGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc         >  1:3280961/1‑93 (MQ=255)
                                                                           gAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCg        >  1:1482530/1‑93 (MQ=255)
                                                                           gAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCg        >  2:602607/1‑93 (MQ=255)
                                                                             aCCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgag      >  1:2053052/1‑93 (MQ=255)
                                                                              ccGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga     <  1:474208/93‑1 (MQ=255)
                                                                              ccGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga     >  1:6017929/1‑93 (MQ=255)
                                                                              ccGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga     >  2:2635607/1‑93 (MQ=255)
                                                                               cGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGAc    <  2:6248817/93‑1 (MQ=255)
                                                                               cGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGAc    <  1:5921264/93‑1 (MQ=255)
                                                                               cGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGAc    <  2:2866629/93‑1 (MQ=255)
                                                                                gTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg   >  2:2123072/1‑93 (MQ=255)
                                                                                gTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg   >  1:4388990/1‑93 (MQ=255)
                                                                                gTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg   >  2:4952827/1‑93 (MQ=255)
                                                                                gTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg   >  2:4249508/1‑93 (MQ=255)
                                                                                 gccggcgcgCCGACCAGCGCGCGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACgg  <  1:700684/92‑1 (MQ=255)
                                                                                 gccggcgcgCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACgg  <  1:2885780/92‑1 (MQ=255)
                                                                                       ||                                                                                     
AGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCCGGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGG  >  NC_000913/3987499‑3987672

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCCGGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTA  >  NC_000913/3987498‑3987675
                                                                                        |                                                                                         
TAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCG                                                                                      >  SRR4302534.5607044/1‑94 (MQ=60)
TAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCG                                                                                      >  SRR4302534.2031587/1‑94 (MQ=60)
 AGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGA                                                                                     >  SRR4302534.5974337/1‑94 (MQ=60)
 AGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGA                                                                                     >  SRR4302534.5990761/1‑94 (MQ=60)
   CTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACC                                                                                   >  SRR4302534.2624589/1‑94 (MQ=60)
   CTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACC                                                                                   >  SRR4302534.4753687/1‑94 (MQ=60)
     CTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTTCGAAGCGTCCGGCGCGCCGACCAG                                                                                 >  SRR4302534.14648/1‑94 (MQ=60)
     CTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAG                                                                                 <  SRR4302534.909531/94‑1 (MQ=60)
       TCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCG                                                                               <  SRR4302534.649535/94‑1 (MQ=60)
        CTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGC                                                                              >  SRR4302534.3237513/1‑94 (MQ=60)
         TGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCA                                                                             >  SRR4302534.3798994/1‑94 (MQ=60)
          GCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCAC                                                                            <  SRR4302534.4346751/94‑1 (MQ=60)
           CCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACG                                                                           <  SRR4302534.2811153/94‑1 (MQ=60)
            CAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGC                                                                          >  SRR4302534.3573572/1‑94 (MQ=60)
               CGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCCCGCAGC                                                                       >  SRR4302534.842416/1‑94 (MQ=60)
               CGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGC                                                                       >  SRR4302534.3072467/1‑94 (MQ=60)
               CGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGC                                                                       >  SRR4302534.4239006/1‑94 (MQ=60)
               CGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGC                                                                       >  SRR4302534.3244962/1‑94 (MQ=60)
               CGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGC                                                                       >  SRR4302534.2527924/1‑94 (MQ=60)
               CGAAATCCCCATGTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCCCCGACCAGCGCACGCAGC                                                                       <  SRR4302534.1865235/94‑1 (MQ=60)
                GAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCC                                                                      >  SRR4302534.6578916/1‑94 (MQ=60)
                 AAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCA                                                                     >  SRR4302534.4257482/1‑94 (MQ=60)
                 AAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCA                                                                     >  SRR4302534.3313086/1‑94 (MQ=60)
                 AAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACAGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCCGCCA                                                                     >  SRR4302534.3341743/1‑94 (MQ=60)
                  AATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAG                                                                    <  SRR4302534.2926957/94‑1 (MQ=60)
                  AATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAG                                                                    >  SRR4302534.5303582/1‑94 (MQ=60)
                   ATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGTGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGA                                                                   >  SRR4302534.2826751/1‑94 (MQ=60)
                   ATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGA                                                                   >  SRR4302534.3310152/1‑94 (MQ=60)
                   ATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGA                                                                   >  SRR4302534.1565289/1‑94 (MQ=60)
                    TCCCCATTTGTTCGGCATCTTGCGGCGCCCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCACCGACCAGCGCACGCAGCAAGAT                                                                  >  SRR4302534.4932749/1‑94 (MQ=60)
                    TCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGCCCAGCGCACGCCGCCAGAT                                                                  >  SRR4302534.2807084/1‑94 (MQ=60)
                    TCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAT                                                                  >  SRR4302534.5416064/1‑94 (MQ=60)
                    TCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAT                                                                  <  SRR4302534.1869128/94‑1 (MQ=60)
                    TCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAT                                                                  >  SRR4302534.4115609/1‑94 (MQ=60)
                     CCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATT                                                                 >  SRR4302534.2081072/1‑94 (MQ=60)
                       CCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTC                                                               >  SRR4302534.5072758/1‑94 (MQ=60)
                        CATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCG                                                              <  SRR4302534.1378839/94‑1 (MQ=60)
                        CATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCG                                                              >  SRR4302534.3870919/1‑94 (MQ=60)
                        CATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCG                                                              >  SRR4302534.5890570/1‑94 (MQ=60)
                          TTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCC                                                            <  SRR4302534.5000327/94‑1 (MQ=60)
                              TTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCA                                                        <  SRR4302534.3045332/94‑1 (MQ=60)
                                CGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAAT                                                      <  SRR4302534.2857192/94‑1 (MQ=60)
                                 GGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATTATCTGCGAACCCTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATA                                                     >  SRR4302534.2442926/1‑94 (MQ=60)
                                 GGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATA                                                     >  SRR4302534.6027345/1‑94 (MQ=60)
                                 GGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATA                                                     >  SRR4302534.473091/1‑94 (MQ=60)
                                   CATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAG                                                   <  SRR4302534.5718577/94‑1 (MQ=60)
                                   CATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAG                                                   >  SRR4302534.2877594/1‑94 (MQ=60)
                                    ATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGC                                                  >  SRR4302534.518874/1‑94 (MQ=60)
                                    ATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGC                                                  >  SRR4302534.4392848/1‑94 (MQ=60)
                                    ATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGC                                                  >  SRR4302534.3532987/1‑94 (MQ=60)
                                     TCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCT                                                 <  SRR4302534.6583798/94‑1 (MQ=60)
                                       TTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCG                                               <  SRR4302534.4161215/94‑1 (MQ=60)
                                       TTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCG                                               <  SRR4302534.5317883/94‑1 (MQ=60)
                                       TTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCG                                               >  SRR4302534.6648269/1‑94 (MQ=60)
                                       TTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCG                                               >  SRR4302534.714475/1‑94 (MQ=60)
                                        TGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGG                                              >  SRR4302534.6421845/1‑94 (MQ=60)
                                        TGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGG                                              >  SRR4302534.4178481/1‑94 (MQ=60)
                                        TGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGG                                              <  SRR4302534.3082371/94‑1 (MQ=60)
                                          CGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCG                                            >  SRR4302534.4445896/1‑94 (MQ=60)
                                            GCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTA                                          <  SRR4302534.4387810/94‑1 (MQ=60)
                                             CGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAG                                         >  SRR4302534.32858/1‑94 (MQ=60)
                                             CGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAG                                         >  SRR4302534.4494765/1‑94 (MQ=60)
                                             CGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACACAGCCAGATTTCGCCATCAATAAACTCGGCGTAG                                         >  SRR4302534.2287266/1‑94 (MQ=60)
                                               CACCGCGGCGTTCGCCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCT                                       >  SRR4302534.5289801/1‑94 (MQ=60)
                                               CACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCT                                       <  SRR4302534.5445846/94‑1 (MQ=60)
                                               CACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCT                                       >  SRR4302534.2413693/1‑94 (MQ=60)
                                                ACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTA                                      <  SRR4302534.2809295/94‑1 (MQ=60)
                                                ACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTA                                      <  SRR4302534.5822235/94‑1 (MQ=60)
                                                ACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTA                                      >  SRR4302534.3474816/1‑94 (MQ=60)
                                                  CGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACC                                    <  SRR4302534.6567133/94‑1 (MQ=60)
                                                   GCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCA                                   <  SRR4302534.5591538/94‑1 (MQ=60)
                                                    CGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAA                                  >  SRR4302534.5393824/1‑94 (MQ=60)
                                                    CGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAA                                  >  SRR4302534.495737/1‑94 (MQ=60)
                                                    CGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAA                                  >  SRR4302534.3927009/1‑94 (MQ=60)
                                                     GGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAAT                                 >  SRR4302534.5205360/1‑94 (MQ=60)
                                                     GGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAAT                                 <  SRR4302534.4957057/94‑1 (MQ=60)
                                                      GCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATT                                <  SRR4302534.3651702/94‑1 (MQ=60)
                                                       CGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTG                               >  SRR4302534.1521552/1‑94 (MQ=60)
                                                       CGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTG                               >  SRR4302534.3267787/1‑94 (MQ=60)
                                                       CGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTG                               >  SRR4302534.4702110/1‑94 (MQ=60)
                                                        GTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGG                              >  SRR4302534.5547799/1‑94 (MQ=60)
                                                        GTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGG                              >  SRR4302534.5887435/1‑94 (MQ=60)
                                                        GTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGG                              >  SRR4302534.1292517/1‑94 (MQ=60)
                                                        GTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGG                              <  SRR4302534.1977133/94‑1 (MQ=60)
                                                        GTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGG                              >  SRR4302534.1601688/1‑94 (MQ=60)
                                                        GTTCACCGCGAATAATCTGCAAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTCG                              >  SRR4302534.6403050/1‑94 (MQ=60)
                                                           CACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGcctaactaccaaattggcac                           >  SRR4302534.5458571/1‑74 (MQ=60)
                                                           CACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCCGCGTAGCTACCAATTGGCAC                           >  SRR4302534.6168593/1‑94 (MQ=60)
                                                             CCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCT                         <  SRR4302534.6219737/94‑1 (MQ=60)
                                                             CCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCT                         >  SRR4302534.4021915/1‑94 (MQ=60)
                                                              CGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTG                        >  SRR4302534.6038459/1‑94 (MQ=60)
                                                              CGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTG                        <  SRR4302534.1137596/94‑1 (MQ=60)
                                                               GCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGA                       >  SRR4302534.362476/1‑94 (MQ=60)
                                                                CGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGAC                      <  SRR4302534.155751/94‑1 (MQ=60)
                                                                CGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCACATTTTGCCATCAATAAGATCCGCGTAGCTACCAATTAGCACCTCAC                      >  SRR4302534.2184110/1‑94 (MQ=60)
                                                                 GAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCTCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACA                     >  SRR4302534.2616326/1‑94 (MQ=60)
                                                                 GAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACA                     <  SRR4302534.3964610/94‑1 (MQ=60)
                                                                  AATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACAT                    <  SRR4302534.977142/94‑1 (MQ=60)
                                                                  AATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACAT                    <  SRR4302534.4291749/94‑1 (MQ=60)
                                                                  AATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACAT                    >  SRR4302534.239872/1‑94 (MQ=60)
                                                                   ATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATC                   >  SRR4302534.6500365/1‑94 (MQ=60)
                                                                    TAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCC                  <  SRR4302534.853100/94‑1 (MQ=60)
                                                                     AATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCG                 <  SRR4302534.3692388/94‑1 (MQ=60)
                                                                     AATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCG                 <  SRR4302534.3148281/94‑1 (MQ=60)
                                                                       TCTGCGAACCGGCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCC               >  SRR4302534.4038043/1‑94 (MQ=60)
                                                                       ACTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCC               <  SRR4302534.4301622/94‑1 (MQ=60)
                                                                        CTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCT              <  SRR4302534.1916826/94‑1 (MQ=60)
                                                                        CTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCT              >  SRR4302534.4068879/1‑94 (MQ=60)
                                                                         TGCGAACCGTCCGGCTCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAAGAAGCTCGGCGTAGCTACCAATGGGCACCTGACATCTGCCTT             <  SRR4302534.3565069/94‑1 (MQ=60)
                                                                         TGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTT             >  SRR4302534.4910662/1‑94 (MQ=60)
                                                                         TGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTT             >  SRR4302534.3359324/1‑94 (MQ=60)
                                                                          GCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTC            <  SRR4302534.2827165/94‑1 (MQ=60)
                                                                          GCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTC            >  SRR4302534.3368403/1‑94 (MQ=60)
                                                                           CGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCG           <  SRR4302534.2392664/94‑1 (MQ=60)
                                                                           CGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCG           <  SRR4302534.5608413/94‑1 (MQ=60)
                                                                           CGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCG           <  SRR4302534.3889493/94‑1 (MQ=60)
                                                                           CGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCG           >  SRR4302534.1521526/1‑94 (MQ=60)
                                                                           CGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCG           <  SRR4302534.420358/94‑1 (MQ=60)
                                                                            GAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGA          >  SRR4302534.1567701/1‑94 (MQ=60)
                                                                             AACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAG         >  SRR4302534.2106694/1‑94 (MQ=60)
                                                                              ACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGA        >  SRR4302534.6181443/1‑94 (MQ=60)
                                                                               CCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGAC       <  SRR4302534.6186772/94‑1 (MQ=60)
                                                                               CCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGAC       <  SRR4302534.486143/94‑1 (MQ=60)
                                                                                CGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACG      <  SRR4302534.6082031/94‑1 (MQ=60)
                                                                                   CCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGT   >  SRR4302534.5855824/1‑94 (MQ=60)
                                                                                    CGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTGGGCGTAGCAACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTA  <  SRR4302534.4079111/94‑1 (MQ=60)
                                                                                    CGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTA  >  SRR4302534.5777952/1‑94 (MQ=60)
                                                                                    CGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTA  >  SRR4302534.4162851/1‑94 (MQ=60)
                                                                                        |                                                                                         
TAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCCGGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTA  >  NC_000913/3987498‑3987675

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