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 4,008,473 +G intergenic (+40/‑193) 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,4721.G98.8% 570.8 / ‑0.9 167intergenic (+39/‑194)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (80/84);  total (82/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 169 total aligned reads displayed.

GAACTGG‑CTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGG‑‑AAA‑G‑CTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑GAAGTGATTGGTACTAATGCCGACGACGCGGTGC  >  NC_000913/4008398‑4008558
                                                                                   |                                                                                        
aaaCTGGTCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACTGGATGACCGACGTCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATg                                                                              <  2:5137555/92‑1 (MQ=255)
 actggg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATggg                                                                            >  1:2061326/5‑93 (MQ=255)
  ctggg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGc                                                                           >  2:3258484/4‑93 (MQ=255)
  ctggg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGc                                                                           <  2:414326/90‑1 (MQ=255)
   tggg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCa                                                                          <  1:4044542/91‑1 (MQ=255)
    ggg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAg                                                                         >  1:4217619/2‑93 (MQ=255)
     gg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGt                                                                        <  2:2856839/93‑1 (MQ=255)
      g‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTg                                                                       <  2:3894347/93‑1 (MQ=255)
        cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTgc                                                                      <  2:6266855/93‑2 (MQ=255)
        cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGCAAG‑G‑CTGCATTATGGGCAGTgc                                                                      >  1:2241773/1‑92 (MQ=255)
         tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTgcg                                                                     <  1:4886730/93‑3 (MQ=255)
         tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTgcg                                                                     >  1:1832504/1‑91 (MQ=255)
         tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTgcg                                                                     >  1:2160175/1‑91 (MQ=255)
         tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTgcg                                                                     <  2:4395037/93‑3 (MQ=255)
             gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCgcacc                                                                 >  1:5509161/3‑93 (MQ=255)
             gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCgcacc                                                                 >  1:2591349/3‑93 (MQ=255)
              ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCgcacca                                                                >  1:1682782/2‑93 (MQ=255)
                tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGc                                                              <  2:443092/93‑1 (MQ=255)
                tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGc                                                              >  1:2156490/1‑93 (MQ=255)
                 gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCg                                                             <  2:5212621/93‑1 (MQ=255)
                 gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCg                                                             >  2:2799496/1‑93 (MQ=255)
                   aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTc                                                           <  2:4358748/93‑1 (MQ=255)
                   aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTc                                                           <  1:3669804/93‑1 (MQ=255)
                   aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTc                                                           <  2:1387586/93‑1 (MQ=255)
                   aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTc                                                           <  2:6133189/93‑1 (MQ=255)
                    ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCt                                                          >  1:1624860/1‑93 (MQ=255)
                    ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCt                                                          <  2:1178316/93‑1 (MQ=255)
                                 tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                                             <  1:4074959/89‑1 (MQ=255)
                                   tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                           <  1:3102204/91‑1 (MQ=255)
                                   tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                           <  2:2888756/91‑1 (MQ=255)
                                    ggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                                          >  1:2995218/2‑93 (MQ=255)
                                    ggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                                          <  2:3794421/92‑1 (MQ=255)
                                    ggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                                          >  2:2321542/2‑93 (MQ=255)
                                     gtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                                         <  2:1706613/93‑1 (MQ=255)
                                     gtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                                         <  1:5936064/93‑1 (MQ=255)
                                      tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGGGCACCAGCGTCTGAAAGACCTTCATCCCgg                                        >  2:6450103/1‑92 (MQ=255)
                                       gACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg                                       <  1:4789934/93‑1 (MQ=255)
                                       gACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg                                       >  2:234961/1‑93 (MQ=255)
                                       gACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg                                       >  2:5144414/1‑93 (MQ=255)
                                       gACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg                                       >  2:5892163/1‑93 (MQ=255)
                                        aCGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAgt                                      <  2:2117644/93‑2 (MQ=255)
                                         cGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt                                     <  2:2663374/93‑1 (MQ=255)
                                         cGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt                                     <  2:975712/93‑1 (MQ=255)
                                          gggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑g                                   <  1:5626457/93‑1 (MQ=255)
                                          gggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑g                                   <  2:6035405/93‑1 (MQ=255)
                                          gggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑g                                   <  2:3208939/93‑1 (MQ=255)
                                           ggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGGCCTTCATCCCGAGCT‑gg                                  >  2:1857379/1‑92 (MQ=255)
                                           ggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑gg                                  >  1:1546250/1‑92 (MQ=255)
                                            gATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑gga                                 >  1:402963/1‑91 (MQ=255)
                                                aaCGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                              >  2:6162777/1‑93 (MQ=255)
                                                aaCGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                              <  2:4770873/93‑1 (MQ=255)
                                                aaCGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                              >  2:4176227/1‑93 (MQ=255)
                                                 acgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                             >  1:1306340/1‑93 (MQ=255)
                                                 acgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                             <  1:2412012/93‑1 (MQ=255)
                                                   gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                           <  1:4446015/93‑1 (MQ=255)
                                                   gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                           >  1:3019064/1‑93 (MQ=255)
                                                   gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                           <  2:4338389/93‑1 (MQ=255)
                                                   gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                           >  2:1748516/1‑93 (MQ=255)
                                                    acgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                          >  2:5834659/1‑93 (MQ=255)
                                                    acgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                          >  1:4063721/1‑93 (MQ=255)
                                                     cgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTgg                         <  2:3583813/93‑1 (MQ=255)
                                                     cgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTgg                         <  1:5366461/93‑1 (MQ=255)
                                                      gCCGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                        <  2:5885363/93‑1 (MQ=255)
                                                       ccGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                       <  1:5983823/93‑1 (MQ=255)
                                                        cGAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc                      <  1:4649788/93‑1 (MQ=255)
                                                         gAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTCTGGGCAATCCGCCCCAGCGTCTGAAAGACCTTCATCTCGAGCTGGAAGTGATTGGTACt                     >  1:5318218/1‑93 (MQ=255)
                                                         gAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                     >  1:6048341/1‑93 (MQ=255)
                                                         gAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                     <  2:3161757/93‑1 (MQ=255)
                                                         gAAATGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                     >  2:4869803/1‑93 (MQ=255)
                                                          aaaTGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTa                    >  2:3471036/1‑93 (MQ=255)
                                                          aaaTGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTa                    >  1:3352200/1‑93 (MQ=255)
                                                          aaaTGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTa                    <  1:237725/93‑1 (MQ=255)
                                                           aaTGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                   >  2:2833463/1‑93 (MQ=255)
                                                           aaTGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                   >  1:3853885/1‑93 (MQ=255)
                                                           aaTGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                   >  1:3449673/1‑93 (MQ=255)
                                                           aaTGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                   <  2:1260781/93‑1 (MQ=255)
                                                            aTGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAt                  >  1:52960/1‑93 (MQ=255)
                                                            aTGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAt                  <  2:4115529/93‑1 (MQ=255)
                                                            aTGCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAt                  <  2:2304196/93‑1 (MQ=255)
                                                              gCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc                >  1:1587885/1‑93 (MQ=255)
                                                              gCTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc                <  1:1050943/93‑1 (MQ=255)
                                                               cTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc               <  1:4084567/93‑1 (MQ=255)
                                                               cTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc               <  2:2232111/93‑1 (MQ=255)
                                                               cTGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc               >  1:4874260/1‑93 (MQ=255)
                                                                tGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg              <  2:2879411/93‑1 (MQ=255)
                                                                tGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg              >  1:1407707/1‑93 (MQ=255)
                                                                tGTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg              <  1:6470595/93‑1 (MQ=255)
                                                                 gTCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga             >  1:2757252/1‑93 (MQ=255)
                                                                  tCGATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac            >  2:678886/1‑93 (MQ=255)
                                                                    gATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga          <  1:279741/93‑1 (MQ=255)
                                                                    gATGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga          <  2:5894190/93‑1 (MQ=255)
                                                                     aTGGCGGGG‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac         >  2:6196961/1‑93 (MQ=255)
                                                                       ggcgggg‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgc       >  1:2807597/5‑93 (MQ=255)
                                                                       ggcgggg‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgc       <  2:3351338/89‑1 (MQ=255)
                                                                          gggg‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt    >  2:5789390/2‑93 (MQ=255)
                                                                          gggg‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt    >  1:2879800/2‑93 (MQ=255)
                                                                          gggg‑‑AAA‑GCCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt    <  2:5624508/92‑1 (MQ=255)
                                                                           ggg‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTg   >  1:5004402/1‑93 (MQ=255)
                                                                           ggg‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTg   <  2:2113744/93‑1 (MQ=255)
                                                                            gg‑‑AAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGc  <  1:5765469/93‑1 (MQ=255)
                                                                                   |                                                                                        
GAACTGG‑CTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGG‑‑AAA‑G‑CTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑GAAGTGATTGGTACTAATGCCGACGACGCGGTGC  >  NC_000913/4008398‑4008558

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
TGGAAGGCTGGCGAACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGG‑‑AAA‑GCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑GAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTA  >  NC_000913/4008386‑4008565
                                                                                                |                                                                                             
gcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAA‑G                                                                                               <  SRR4302534.5495819/76‑1 (MQ=53)
       cgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCAT                                                                                         >  SRR4302534.93156/12‑94 (MQ=49)
        gaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATT                                                                                        >  SRR4302534.3152769/11‑94 (MQ=49)
        gaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATT                                                                                        >  SRR4302534.1179196/11‑94 (MQ=49)
        gaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATT                                                                                        >  SRR4302534.6188388/11‑94 (MQ=49)
         aagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTA                                                                                       >  SRR4302534.135347/10‑94 (MQ=48)
          agaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTAT                                                                                      >  SRR4302534.1117497/9‑94 (MQ=48)
          agaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTAT                                                                                      >  SRR4302534.5535424/9‑94 (MQ=48)
            aaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGG                                                                                    <  SRR4302534.3347287/88‑1 (MQ=51)
             AACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGG                                                                                   >  SRR4302534.2115228/1‑94 (MQ=50)
               ctgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCA                                                                                 >  SRR4302534.2941981/4‑94 (MQ=57)
                tgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAG                                                                                >  SRR4302534.4331355/3‑94 (MQ=58)
                tgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAG                                                                                <  SRR4302534.4153009/92‑1 (MQ=58)
                  GGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTG                                                                              <  SRR4302534.9300/94‑1 (MQ=60)
                   GCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGCAAG‑GCTGCATTATGGGCAGTGC                                                                             >  SRR4302534.2301180/1‑94 (MQ=60)
                    CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGcg                                                                            >  SRR4302534.1880159/1‑92 (MQ=60)
                    CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGcg                                                                            >  SRR4302534.2834767/1‑92 (MQ=60)
                    CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGcg                                                                            >  SRR4302534.2217226/1‑92 (MQ=60)
                    CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGcg                                                                            >  SRR4302534.1663414/1‑92 (MQ=60)
                    CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGcg                                                                            <  SRR4302534.1349270/94‑3 (MQ=60)
                     TACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGcgc                                                                           <  SRR4302534.5019250/94‑4 (MQ=60)
                     TACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGcgc                                                                           >  SRR4302534.5254671/1‑91 (MQ=60)
                     TACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGcgc                                                                           >  SRR4302534.3830216/1‑91 (MQ=60)
                      ACAGCCTGAAGGATTGTATTGCGTTTGGTGAGGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGcgca                                                                          >  SRR4302534.140335/1‑90 (MQ=50)
                       CAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGcgcac                                                                         >  SRR4302534.5463609/1‑89 (MQ=60)
                        agcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACC                                                                        >  SRR4302534.5657772/4‑94 (MQ=60)
                        agcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACC                                                                        >  SRR4302534.2659947/4‑94 (MQ=60)
                         gcgTGAAGGATTGTATTGCGTTGGGTGACGGAATGAACGACGCCGAAATGCTGTCAATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCA                                                                       <  SRR4302534.1699882/91‑1 (MQ=37)
                         gcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCA                                                                       >  SRR4302534.1726551/3‑94 (MQ=60)
                           CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGcgcaccaga                                                                     <  SRR4302534.2274463/94‑10 (MQ=60)
                           CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGC                                                                     >  SRR4302534.2213419/1‑94 (MQ=60)
                               AGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCT                                                                 >  SRR4302534.1667254/1‑94 (MQ=57)
                               AGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCT                                                                 <  SRR4302534.4875070/94‑1 (MQ=57)
                               AGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCT                                                                 <  SRR4302534.3768462/94‑1 (MQ=57)
                                GGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTG                                                                >  SRR4302534.3928917/1‑94 (MQ=60)
                                GGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTG                                                                <  SRR4302534.4932693/94‑1 (MQ=60)
                                GGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTG                                                                <  SRR4302534.5863326/94‑1 (MQ=60)
                                   TTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAA                                                             >  SRR4302534.779915/1‑94 (MQ=60)
                                     gtattgcgtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGA                                                           >  SRR4302534.6032154/13‑94 (MQ=60)
                                     gtattgcgtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGA                                                           >  SRR4302534.969933/13‑94 (MQ=60)
                                        TTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                        <  SRR4302534.4827236/94‑1 (MQ=60)
                                        TTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                        <  SRR4302534.5009941/94‑1 (MQ=60)
                                        TTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                        <  SRR4302534.1155258/94‑1 (MQ=60)
                                         TGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTT                                                       <  SRR4302534.4798702/94‑1 (MQ=60)
                                         TGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTT                                                       <  SRR4302534.4973369/94‑1 (MQ=60)
                                          GCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTC                                                      <  SRR4302534.2778805/94‑1 (MQ=60)
                                           CGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCA                                                     >  SRR4302534.646002/1‑94 (MQ=60)
                                           CGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCA                                                     <  SRR4302534.6172232/94‑1 (MQ=60)
                                             tttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATC                                                   <  SRR4302534.4184323/90‑1 (MQ=60)
                                             tttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATC                                                   <  SRR4302534.5526803/90‑1 (MQ=60)
                                             TG‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCC                                                 <  SRR4302534.3183809/94‑1 (MQ=60)
                                             TG‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCC                                                 >  SRR4302534.3073533/1‑94 (MQ=60)
                                                 GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGA                                               <  SRR4302534.6097252/94‑1 (MQ=60)
                                                  TGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAG                                              >  SRR4302534.3004236/1‑94 (MQ=60)
                                                  TGACGAGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAG                                              <  SRR4302534.2377994/94‑1 (MQ=60)
                                                   GACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGC                                             <  SRR4302534.4919994/94‑1 (MQ=60)
                                                    ACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT                                            >  SRR4302534.6396086/1‑94 (MQ=60)
                                                     CGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑G                                          >  SRR4302534.5120675/1‑94 (MQ=60)
                                                     CGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑G                                          >  SRR4302534.4137848/1‑94 (MQ=60)
                                                      GGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑GG                                         >  SRR4302534.1586773/1‑94 (MQ=60)
                                                      GGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑GG                                         <  SRR4302534.5778724/94‑1 (MQ=60)
                                                       GGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑GGA                                        >  SRR4302534.413263/1‑94 (MQ=60)
                                                         ATGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAG                                       <  SRR4302534.6072789/94‑1 (MQ=60)
                                                          TGAACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGT                                      >  SRR4302534.4101236/1‑94 (MQ=60)
                                                            AACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGA                                    >  SRR4302534.1340804/1‑94 (MQ=60)
                                                             ACGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAT                                   <  SRR4302534.2476150/94‑1 (MQ=60)
                                                              CGACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATT                                  >  SRR4302534.3098112/1‑94 (MQ=60)
                                                               GACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTG                                 >  SRR4302534.4172744/1‑94 (MQ=60)
                                                               GACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTG                                 <  SRR4302534.4566139/94‑1 (MQ=60)
                                                                ACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GGTGCAATATGGGCAGTGCGCACCAGCGTCTGAAACACCTTCATCCCGAGgtggaagtcaatgg                                >  SRR4302534.1717783/1‑80 (MQ=39)
                                                                ACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGG                                <  SRR4302534.215487/94‑1 (MQ=60)
                                                                ACGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGG                                >  SRR4302534.3816131/1‑94 (MQ=60)
                                                                 CGCCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGT                               <  SRR4302534.5511102/94‑1 (MQ=60)
                                                                   CCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAC                             >  SRR4302534.4448491/1‑94 (MQ=60)
                                                                   CCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAC                             <  SRR4302534.6146366/94‑1 (MQ=60)
                                                                   CCGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAC                             <  SRR4302534.4159739/94‑1 (MQ=60)
                                                                    CGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTCTGGGCAATCCGCCCCAGCGTCTGAAAGACCTTCATCTCGAGCTGGAAGTGATTGGTACT                            >  SRR4302534.5461671/1‑94 (MQ=26)
                                                                    CGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACT                            >  SRR4302534.6212755/1‑94 (MQ=60)
                                                                    CGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACT                            >  SRR4302534.3439663/1‑94 (MQ=60)
                                                                    CGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACT                            <  SRR4302534.4776016/94‑1 (MQ=60)
                                                                    CGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACT                            >  SRR4302534.2022668/1‑94 (MQ=60)
                                                                    CGAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGCGCGCACTAGCGTCTGAAAGACCTTTATCCCCAGCTGGTAGTGCTTGGTACT                            >  SRR4302534.1665812/1‑94 (MQ=23)
                                                                     GAAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTA                           >  SRR4302534.3441984/1‑94 (MQ=60)
                                                                      AAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAA                          >  SRR4302534.3957355/1‑94 (MQ=60)
                                                                      AAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAA                          <  SRR4302534.243944/94‑1 (MQ=60)
                                                                      AAATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAA                          >  SRR4302534.3542398/1‑94 (MQ=60)
                                                                       AATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAT                         >  SRR4302534.4181944/1‑94 (MQ=60)
                                                                       AATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAT                         >  SRR4302534.54380/1‑94 (MQ=60)
                                                                        ATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATG                        >  SRR4302534.4258082/1‑94 (MQ=60)
                                                                        ATGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATG                        >  SRR4302534.3121531/1‑94 (MQ=60)
                                                                         TGCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGC                       >  SRR4302534.1629423/1‑94 (MQ=60)
                                                                          GCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCC                      >  SRR4302534.5006458/1‑94 (MQ=60)
                                                                          GCTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCC                      <  SRR4302534.1078154/94‑1 (MQ=60)
                                                                           CTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCG                     <  SRR4302534.4194240/94‑1 (MQ=60)
                                                                           CTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCG                     <  SRR4302534.5468422/94‑1 (MQ=60)
                                                                           CTGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCG                     >  SRR4302534.1444730/1‑94 (MQ=60)
                                                                            TGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGA                    <  SRR4302534.6644884/94‑1 (MQ=60)
                                                                            TGTCGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGA                    >  SRR4302534.2829878/1‑94 (MQ=60)
                                                                               CGATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGGGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTCATGCCGACGA                 >  SRR4302534.1486913/1‑94 (MQ=57)
                                                                                GATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAC                <  SRR4302534.6129864/94‑1 (MQ=60)
                                                                                GATGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAC                <  SRR4302534.286975/94‑1 (MQ=60)
                                                                                 ATGGC‑GGGTAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGAGCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACG               <  SRR4302534.3699943/94‑1 (MQ=60)
                                                                                  TGGC‑GGGGAAAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGC              >  SRR4302534.2881431/1‑94 (MQ=60)
                                                                                    GG‑GGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGT           <  SRR4302534.2772061/94‑1 (MQ=60)
                                                                                      gcggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGG            >  SRR4302534.5204516/11‑94 (MQ=60)
                                                                                       cggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGT           >  SRR4302534.2955382/10‑94 (MQ=60)
                                                                                       GGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGC         >  SRR4302534.2583954/1‑94 (MQ=60)
                                                                                        ggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTG          >  SRR4302534.5139687/9‑94 (MQ=60)
                                                                                        ggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑Gggggagagtgtttctaataaacgcgacgagcgg         >  SRR4302534.5255105/9‑61 (MQ=60)
                                                                                        GG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCC        <  SRR4302534.5921882/94‑1 (MQ=60)
                                                                                            AAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATTCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCT     >  SRR4302534.3353768/1‑94 (MQ=60)
                                                                                            AAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCA     >  SRR4302534.1749025/1‑94 (MQ=60)
                                                                                            AAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCA     >  SRR4302534.4997272/1‑94 (MQ=60)
                                                                                            AAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCA     >  SRR4302534.748747/1‑94 (MQ=60)
                                                                                            AAG‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCA     >  SRR4302534.2407346/1‑94 (MQ=60)
                                                                                            AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGC      >  SRR4302534.3660162/1‑94 (MQ=60)
                                                                                            AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGC      >  SRR4302534.2864792/1‑94 (MQ=60)
                                                                                            AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGC      >  SRR4302534.4512792/1‑94 (MQ=60)
                                                                                              G‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATT   >  SRR4302534.2303568/1‑94 (MQ=60)
                                                                                              G‑GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATT   <  SRR4302534.3074145/94‑1 (MQ=60)
                                                                                                GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTA  <  SRR4302534.3220387/94‑1 (MQ=60)
                                                                                                |                                                                                             
TGGAAGGCTGGCGAACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGG‑‑AAA‑GCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑GAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTA  >  NC_000913/4008386‑4008565

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