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 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,491 2 bp→CG intergenic (+58/‑174) b4555 → / → yigM hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,4910GC98.8% 526.7 / NA 169intergenic (+58/‑175)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base C (70/97);  total (72/97)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,008,4920CG100.0% 523.1 / NA 165intergenic (+59/‑174)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (68/97);  total (68/97)

Only 100 of 174 total aligned reads displayed.

AGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAA‑‑GCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAC  >  NC_000913/4008409‑4008576
                                                                                           ||                                                                                      
gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCgcacc                                                                                    >  1:4332377/3‑93 (MQ=255)
gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCgcacc                                                                                    >  2:4246988/3‑93 (MQ=255)
gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCgcacc                                                                                    >  2:4074519/3‑93 (MQ=255)
 ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCgcaccg                                                                                   >  1:750969/2‑92 (MQ=255)
  cTGAAGGATTGTATTGCGTTTGGTGACGGGGTGAACGGCGCCGTAATGCTGTCGAAGGCGGGG‑‑AAA‑‑GCCTGCATTATGGGCAGGGGGTACCAg                                                                                  >  2:1076668/1‑93 (MQ=255)
  cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAg                                                                                  >  1:1192872/1‑93 (MQ=255)
  cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAg                                                                                  <  1:3892728/93‑1 (MQ=255)
  cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAg                                                                                  <  1:960598/93‑1 (MQ=255)
   tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                 >  1:336362/1‑93 (MQ=255)
   tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                 <  1:2756240/93‑1 (MQ=255)
    gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCg                                                                                >  2:1980376/1‑93 (MQ=255)
    gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGTAAG‑‑GCTGCATTATGGGCAGTGCGCACCAGCg                                                                                >  2:3586240/1‑93 (MQ=255)
     aaGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGt                                                                               <  2:1341619/93‑1 (MQ=255)
      aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                                              >  2:2104761/1‑93 (MQ=255)
      aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                                              >  2:1623503/1‑93 (MQ=255)
       ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt                                                                             >  2:399184/1‑93 (MQ=255)
       ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt                                                                             <  1:4531157/93‑1 (MQ=255)
       ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt                                                                             >  2:4164986/1‑93 (MQ=255)
                    tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                                                                >  2:1121649/5‑93 (MQ=255)
                      tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                                              <  2:3075612/91‑1 (MQ=255)
                      tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                                              <  2:2683642/91‑1 (MQ=255)
                      tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                                              >  2:1871860/3‑93 (MQ=255)
                         tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                                           >  1:81779/1‑93 (MQ=255)
                         tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                                           <  1:1791295/93‑1 (MQ=255)
                           aCGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc                                                        <  2:110176/93‑1 (MQ=255)
                             gggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTg                                                     <  1:1271192/93‑1 (MQ=255)
                              ggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgg                                                    >  1:1064987/1‑92 (MQ=255)
                              ggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgg                                                    <  2:1252333/93‑2 (MQ=255)
                               gATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgga                                                   <  1:3208493/93‑3 (MQ=255)
                               gATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgga                                                   <  1:1244152/93‑3 (MQ=255)
                                aTGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                                  <  1:2648648/93‑4 (MQ=255)
                                aTGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                                  <  2:439256/93‑4 (MQ=255)
                                 tGAACGACGCGGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                                                 <  1:3075847/93‑5 (MQ=255)
                                 tGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                                                 <  1:3060309/93‑5 (MQ=255)
                                   aaCGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                                                <  1:3223860/93‑1 (MQ=255)
                                    acgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                                               <  2:4116374/93‑1 (MQ=255)
                                    acgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                                               >  2:2467002/1‑93 (MQ=255)
                                     cgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                                              <  2:3949962/93‑1 (MQ=255)
                                      gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                                             >  2:2223405/1‑93 (MQ=255)
                                      gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                                             <  1:4441255/93‑1 (MQ=255)
                                       acgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                                            >  1:910991/1‑93 (MQ=255)
                                       acgCCGAAATGCTCTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                                            <  1:2678236/93‑1 (MQ=255)
                                         gCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                                          <  2:507341/93‑1 (MQ=255)
                                           cGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc                                        >  1:1390274/1‑93 (MQ=255)
                                           cGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc                                        >  1:4173572/1‑93 (MQ=255)
                                            gAAATGCTGTCGCTGGCGGG‑‑AAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCCCGCAGTGATTGGTACt                                       >  2:841908/1‑93 (MQ=255)
                                            gAAATGCTGGCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGAACTTCATCCCGAGCTGGAAGGGATTGGTACt                                       >  1:2236980/1‑93 (MQ=255)
                                              aaTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                     <  1:3768821/93‑1 (MQ=255)
                                              aaTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                     <  2:2085270/93‑1 (MQ=255)
                                               aTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAt                                    >  2:3008086/1‑93 (MQ=255)
                                                tGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                                   >  1:2107200/1‑93 (MQ=255)
                                                 gCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc                                  >  1:3024996/1‑93 (MQ=255)
                                                 gCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc                                  >  2:3400004/1‑93 (MQ=255)
                                                 gCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc                                  <  2:2580693/93‑1 (MQ=255)
                                                  cTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc                                 <  1:1463743/93‑1 (MQ=255)
                                                    gTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGTAGTGATTGGTACTAATGCcga                               >  1:920091/1‑93 (MQ=255)
                                                    gTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                               <  1:3254248/93‑1 (MQ=255)
                                                    gTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                               <  2:1568400/93‑1 (MQ=255)
                                                       gATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga                            <  2:2638556/93‑1 (MQ=255)
                                                        aTGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGGGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac                           >  2:430081/1‑93 (MQ=255)
                                                        aTGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCCCCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac                           <  1:680626/93‑1 (MQ=255)
                                                        aTGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac                           <  2:2255970/93‑1 (MQ=255)
                                                        aTGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac                           <  1:3757164/93‑1 (MQ=255)
                                                           gcgggg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg                        >  1:1452411/4‑93 (MQ=255)
                                                           gcgggg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg                        >  1:4103495/4‑93 (MQ=255)
                                                           gcgggg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg                        <  2:1336706/90‑1 (MQ=255)
                                                            cgggg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCgg                       <  2:4654945/91‑1 (MQ=255)
                                                             gggg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCggg                      >  1:338655/2‑92 (MQ=255)
                                                             gggg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt                      <  1:4461489/92‑1 (MQ=255)
                                                             gggg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt                      <  1:3579281/92‑1 (MQ=255)
                                                             gggg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt                      <  1:3202892/92‑1 (MQ=255)
                                                               gg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGc                    >  2:1678308/1‑93 (MQ=255)
                                                                    aaaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCg                  >  1:3156432/5‑93 (MQ=255)
                                                                     aaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGc                 >  2:54098/4‑93 (MQ=255)
                                                                      aggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCa                >  2:4294051/3‑93 (MQ=255)
                                                                       ggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCAt               <  1:1817108/92‑1 (MQ=255)
                                                                       ggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCAt               <  1:2349192/92‑1 (MQ=255)
                                                                        gCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCAtt              <  1:2504309/93‑1 (MQ=255)
                                                                          cTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTa             <  1:3807552/93‑1 (MQ=255)
                                                                           tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAt            >  1:1058156/1‑93 (MQ=255)
                                                                            gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc           >  2:3526253/1‑93 (MQ=255)
                                                                            gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc           >  1:620171/1‑93 (MQ=255)
                                                                            gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc           <  1:1236797/93‑1 (MQ=255)
                                                                            gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc           <  2:4100270/93‑1 (MQ=255)
                                                                            gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc           >  1:4114126/1‑93 (MQ=255)
                                                                             cATTAGGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt          <  2:3728277/93‑1 (MQ=255)
                                                                               ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc        <  1:2484749/93‑1 (MQ=255)
                                                                               ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc        <  1:2201822/93‑1 (MQ=255)
                                                                               ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc        >  2:3920307/1‑93 (MQ=255)
                                                                               ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc        <  1:2283480/93‑1 (MQ=255)
                                                                               ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc        <  2:2170902/93‑1 (MQ=255)
                                                                               ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc        >  2:1487235/1‑93 (MQ=255)
                                                                                tATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCg       <  2:342322/93‑1 (MQ=255)
                                                                                  tGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTa     >  2:857525/1‑93 (MQ=255)
                                                                                  tGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCACATTATCTGCGTa     <  2:3531500/93‑1 (MQ=255)
                                                                                   gggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaa    >  1:2370385/1‑93 (MQ=255)
                                                                                    ggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaaa   >  2:1733177/1‑93 (MQ=255)
                                                                                    ggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaaa   <  1:3503039/93‑1 (MQ=255)
                                                                                    ggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaaa   <  1:710534/93‑1 (MQ=255)
                                                                                     gCAGTGCGCACCAGCGGCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAc  <  2:199647/93‑1 (MQ=255)
                                                                                           ||                                                                                      
AGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAA‑‑GCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAC  >  NC_000913/4008409‑4008576

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ 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
GCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCT  >  NC_000913/4008404‑4008579
                                                                                              |                                                                                         
GCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGC                                                                                           >  SRR4302535.2050421/1‑94 (MQ=60)
     agcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACC                                                                                      >  SRR4302535.4362005/4‑94 (MQ=60)
      CCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAG                                                                                    >  SRR4302535.1201421/1‑94 (MQ=60)
        CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC                                                                                   <  SRR4302535.967423/94‑1 (MQ=60)
        CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC                                                                                   <  SRR4302535.3919415/94‑1 (MQ=60)
        CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC                                                                                   >  SRR4302535.1255636/1‑94 (MQ=60)
        CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC                                                                                   >  SRR4302535.338688/1‑94 (MQ=60)
         TGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCG                                                                                  <  SRR4302535.2775887/94‑1 (MQ=57)
            AGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCT                                                                               <  SRR4302535.2825932/94‑1 (MQ=57)
             GGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTG                                                                              <  SRR4302535.4562173/94‑1 (MQ=60)
              GATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGA                                                                             >  SRR4302535.57736/1‑94 (MQ=60)
               ATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAA                                                                            <  SRR4302535.912953/94‑1 (MQ=60)
               ATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAA                                                                            <  SRR4302535.3356104/94‑1 (MQ=60)
                 tgtattgcgtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAG                                                                          >  SRR4302535.2401960/14‑94 (MQ=60)
                 tgtattgcgtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAG                                                                          >  SRR4302535.2795015/14‑94 (MQ=60)
                 tgtattgcgtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAG                                                                          <  SRR4302535.1707369/81‑1 (MQ=60)
                  gtattgcgtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGA                                                                         <  SRR4302535.2467044/82‑1 (MQ=60)
                  TATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGAC                                                                        <  SRR4302535.582540/94‑1 (MQ=60)
                    ATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACC                                                                       >  SRR4302535.635963/1‑94 (MQ=60)
                    ATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACC                                                                       >  SRR4302535.237125/1‑94 (MQ=60)
                    ATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACC                                                                       >  SRR4302535.2381732/1‑94 (MQ=60)
                     gtgcgtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                                      <  SRR4302535.2084021/85‑1 (MQ=60)
                     TTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                                      >  SRR4302535.1034242/1‑94 (MQ=60)
                     TTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                                      >  SRR4302535.2252083/1‑94 (MQ=60)
                     TTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                                      <  SRR4302535.2345907/94‑1 (MQ=60)
                     TTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                                      <  SRR4302535.2084022/94‑1 (MQ=60)
                      TGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTT                                                                     <  SRR4302535.992706/94‑1 (MQ=60)
                      TGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTT                                                                     <  SRR4302535.941555/94‑1 (MQ=60)
                        CGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCA                                                                   <  SRR4302535.3945934/94‑1 (MQ=60)
                        CGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCA                                                                   >  SRR4302535.1778076/1‑94 (MQ=60)
                         gtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAT                                                                  <  SRR4302535.1862662/89‑1 (MQ=60)
                          TG‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCC                                                               <  SRR4302535.1460901/94‑1 (MQ=60)
                              GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGA                                                             >  SRR4302535.82362/1‑94 (MQ=60)
                               TGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAG                                                            <  SRR4302535.1804248/94‑1 (MQ=60)
                               GGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAG                                                            <  SRR4302535.2249525/94‑1 (MQ=60)
                                   GGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGG                                                       >  SRR4302535.4193544/1‑94 (MQ=60)
                                   GGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGG                                                       <  SRR4302535.1280275/94‑1 (MQ=60)
                                   GGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGG                                                       >  SRR4302535.1072564/1‑94 (MQ=60)
                                     GATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGgaa                                                     <  SRR4302535.3230716/94‑4 (MQ=60)
                                     GATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGgaa                                                     <  SRR4302535.1253032/94‑4 (MQ=60)
                                      ATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAG                                                     <  SRR4302535.2667620/94‑1 (MQ=60)
                                      ATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAG                                                     >  SRR4302535.3196930/1‑94 (MQ=60)
                                       TGAACGACGCGGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGT                                                    <  SRR4302535.3097308/94‑1 (MQ=60)
                                       TGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGT                                                    >  SRR4302535.985089/1‑94 (MQ=60)
                                       TGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGT                                                    <  SRR4302535.3081714/94‑1 (MQ=60)
                                         AACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGA                                                  <  SRR4302535.3246209/94‑1 (MQ=60)
                                          ACGACGCCGAGATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAT                                                 <  SRR4302535.454719/94‑1 (MQ=60)
                                          ACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAT                                                 <  SRR4302535.1126090/94‑1 (MQ=60)
                                            GACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTG                                               >  SRR4302535.917468/1‑94 (MQ=60)
                                            GACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTG                                               <  SRR4302535.4471707/94‑1 (MQ=60)
                                             ACGCCGAAATGCTCTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGG                                              <  SRR4302535.2697397/94‑1 (MQ=60)
                                                CCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAC                                           >  SRR4302535.1400195/1‑94 (MQ=60)
                                                CCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAC                                           >  SRR4302535.4202139/1‑94 (MQ=60)
                                                 CGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACT                                          <  SRR4302535.2293395/94‑1 (MQ=60)
                                                  GAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTA                                         >  SRR4302535.1820182/1‑94 (MQ=60)
                                                   AAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAA                                        >  SRR4302535.2344464/1‑94 (MQ=60)
                                                   AAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAA                                        >  SRR4302535.3415029/1‑94 (MQ=60)
                                                    AATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAT                                       <  SRR4302535.897072/94‑1 (MQ=60)
                                                    AATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAT                                       <  SRR4302535.3058682/94‑1 (MQ=60)
                                                    AATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAT                                       >  SRR4302535.2651406/1‑94 (MQ=60)
                                                    AATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAT                                       <  SRR4302535.3794835/94‑1 (MQ=60)
                                                     ATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATG                                      >  SRR4302535.2122399/1‑94 (MQ=60)
                                                      TGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGC                                     <  SRR4302535.359317/94‑1 (MQ=60)
                                                      TGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGC                                     >  SRR4302535.3046173/1‑94 (MQ=60)
                                                       GCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGGGCACCAGCGTCTGAAAGACCTTCATACCGAGCTGGAAGGGATTGGTACTTATGCC                                    >  SRR4302535.3150404/1‑94 (MQ=57)
                                                        CTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCG                                   <  SRR4302535.1474228/94‑1 (MQ=60)
                                                        CTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCG                                   >  SRR4302535.1702494/1‑94 (MQ=60)
                                                         TGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGTAGTGATTGGTACTAATGCCGA                                  >  SRR4302535.926616/1‑94 (MQ=60)
                                                          GTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGAC                                 <  SRR4302535.3276858/94‑1 (MQ=60)
                                                          GTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGAC                                 <  SRR4302535.2070629/94‑1 (MQ=60)
                                                           TCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACG                                <  SRR4302535.4645279/94‑1 (MQ=60)
                                                            CGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGA                               <  SRR4302535.380924/94‑1 (MQ=60)
                                                             GATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAC                              <  SRR4302535.2520935/94‑1 (MQ=60)
                                                              ATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCCCCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACG                             <  SRR4302535.685392/94‑1 (MQ=60)
                                                              ATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACG                             <  SRR4302535.3783146/94‑1 (MQ=60)
                                                              ATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACG                             >  SRR4302535.811241/1‑94 (MQ=60)
                                                                GGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCG                           >  SRR4302535.4131587/1‑94 (MQ=60)
                                                                GGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCG                           >  SRR4302535.1462804/1‑94 (MQ=60)
                                                                GGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCG                           <  SRR4302535.3247612/94‑1 (MQ=60)
                                                                  gcggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGG                          <  SRR4302535.1316714/84‑1 (MQ=60)
                                                                   cggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGT                         <  SRR4302535.801094/85‑1 (MQ=60)
                                                                   cggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGT                         >  SRR4302535.3756306/10‑94 (MQ=60)
                                                                   cggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGT                         <  SRR4302535.1236952/85‑1 (MQ=60)
                                                                   cggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGG                         >  SRR4302535.341000/10‑94 (MQ=60)
                                                                    ggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTG                        <  SRR4302535.3225069/86‑1 (MQ=60)
                                                                    ggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTG                        <  SRR4302535.3604091/86‑1 (MQ=60)
                                                                    ggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTG                        <  SRR4302535.4492088/86‑1 (MQ=60)
                                                                    GG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCC                      >  SRR4302535.3849761/1‑94 (MQ=60)
                                                                     G‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCG                     >  SRR4302535.3178349/1‑94 (MQ=60)
                                                                        AAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCA                   >  SRR4302535.3997677/1‑94 (MQ=60)
                                                                        AAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCA                   >  SRR4302535.941217/1‑94 (MQ=60)
                                                                          GGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATT                 <  SRR4302535.1830240/94‑1 (MQ=60)
                                                                          GGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATT                 <  SRR4302535.2366038/94‑1 (MQ=60)
                                                                            GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTA                <  SRR4302535.2522363/94‑1 (MQ=60)
                                                                             CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAT               >  SRR4302535.1065684/1‑94 (MQ=60)
                                                                             CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAT               <  SRR4302535.3833812/94‑1 (MQ=60)
                                                                              TGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCCTTATC              >  SRR4302535.2736397/1‑94 (MQ=60)
                                                                              TGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATC              >  SRR4302535.624483/1‑94 (MQ=60)
                                                                              TGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATC              >  SRR4302535.4142286/1‑94 (MQ=60)
                                                                               GCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCT             <  SRR4302535.4302279/94‑1 (MQ=60)
                                                                               GCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCT             <  SRR4302535.1245631/94‑1 (MQ=60)
                                                                               GCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCT             <  SRR4302535.745094/94‑1 (MQ=60)
                                                                                 ATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGC           >  SRR4302535.3893893/1‑94 (MQ=60)
                                                                                  TTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCG          <  SRR4302535.2502661/94‑1 (MQ=60)
                                                                                  TTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCG          <  SRR4302535.2299796/94‑1 (MQ=60)
                                                                                  TTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCG          <  SRR4302535.2217515/94‑1 (MQ=60)
                                                                                     TGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAA       <  SRR4302535.3239698/94‑1 (MQ=60)
                                                                                     TGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAA       >  SRR4302535.665689/1‑94 (MQ=60)
                                                                                     TGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAA       >  SRR4302535.2387404/1‑94 (MQ=60)
                                                                                      GGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAA      >  SRR4302535.3313448/1‑94 (MQ=60)
                                                                                       GGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAC     <  SRR4302535.3527387/94‑1 (MQ=60)
                                                                                       GGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAC     <  SRR4302535.715520/94‑1 (MQ=60)
                                                                                        GCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACT    <  SRR4302535.2423258/94‑1 (MQ=60)
                                                                                        GCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACT    <  SRR4302535.4400185/94‑1 (MQ=60)
                                                                                          AGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCT  <  SRR4302535.4028998/94‑1 (MQ=60)
                                                                                          AGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCT  >  SRR4302535.1720175/1‑94 (MQ=60)
                                                                                              |                                                                                         
GCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCT  >  NC_000913/4008404‑4008579

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