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 1624 250.6 12620048 98.0% 12367647 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,969,103 2 bp→GC coding (232‑233/1263 nt) wecC → UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,969,1030CG100.0% 762.1 / NA 241R78A (CGC→GCC) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (103/138);  total (103/138)
*NC_0009133,969,1040GC99.6% 716.5 / NA 228R78A (CGC→GCC) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (92/134);  total (94/134)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 247 total aligned reads displayed.

ATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATCGCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCG  >  NC_000913/3969016‑3969189
                                                                                        ||                                                                                     
atCGTCGAACCTGATTTGGCGAGCGTAGTAAAAACTGCCGTCGCAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGcctgg                                                                                   >  1:595690/1‑93 (MQ=255)
   gtcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGcctggctg                                                                                >  1:454781/1‑93 (MQ=255)
   gtcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGcctggctg                                                                                >  1:1532831/1‑93 (MQ=255)
    tcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGcctggctgt                                                                               >  1:421252/1‑92 (MQ=255)
    tcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGa                                                                               >  1:2411042/1‑93 (MQ=255)
      gAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAtt                                                                             >  2:893286/1‑93 (MQ=255)
      gAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAtt                                                                             >  1:294372/1‑93 (MQ=255)
       aaCCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTg                                                                            >  2:3191146/1‑93 (MQ=255)
        aCCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGc                                                                           <  1:700626/93‑1 (MQ=255)
         ccTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GCGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCt                                                                          >  2:1851994/1‑93 (MQ=255)
         ccTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCt                                                                          >  1:6280688/1‑93 (MQ=255)
         ccTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCt                                                                          >  2:2543315/1‑93 (MQ=255)
         ccTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCt                                                                          >  2:4102349/1‑93 (MQ=255)
          cTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTg                                                                         >  1:1623253/1‑93 (MQ=255)
           tGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGt                                                                        >  1:101022/1‑93 (MQ=255)
           tGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGt                                                                        <  2:515917/93‑1 (MQ=255)
            gATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTa                                                                       >  1:4070955/1‑93 (MQ=255)
             aTTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAc                                                                      >  1:1876544/1‑93 (MQ=255)
              tttGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAcc                                                                     >  1:6077885/1‑93 (MQ=255)
                tGGCGAGTGTAGTAAAACCTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCa                                                                   <  2:1746401/93‑1 (MQ=255)
                tGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCa                                                                   >  2:2789007/1‑93 (MQ=255)
                 ggCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAc                                                                  >  2:558904/1‑93 (MQ=255)
                    gAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGcc                                                               >  1:4529691/1‑93 (MQ=255)
                    gAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGcc                                                               <  2:531598/93‑1 (MQ=255)
                     aGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCg                                                              >  2:3134942/1‑93 (MQ=255)
                       tgtAGTAAAAACTGCCGTAGAAGGCGTTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGtt                                                            <  1:2235898/93‑1 (MQ=255)
                       tgtAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGtt                                                            >  1:1886519/1‑93 (MQ=255)
                        gtagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGttt                                                           >  1:304865/1‑93 (MQ=255)
                         tagaaAAACCTGCGGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTa                                                          <  2:5750688/93‑1 (MQ=255)
                         cagtaAAACCTGCCGCAGAAGTCCGTTACTCACGACC‑GAGCACGACGCCAGTTGAAGCTGATGCCTGGCTGATTGCTCTACCCACGCCGTTTa                                                          <  2:2613024/92‑1 (MQ=255)
                          agaaaaaaCGGCGGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTaa                                                         <  2:223701/90‑1 (MQ=255)
                           gtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAg                                                        <  2:2112384/93‑1 (MQ=255)
                           gtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGACTGGCTGATTGCTGTACCCACGCCGTTTAAg                                                        >  2:6263287/1‑93 (MQ=255)
                           gaaaaaaCTGCCGAAGAGGGCGGTTTTTTCCGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAg                                                        <  2:1385257/91‑1 (MQ=255)
                            taAAAACTGCCGTAGAAGGCGGTTTTTTACGAAC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAgg                                                       >  1:5484762/1‑93 (MQ=255)
                             aaaaaCTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAggg                                                      >  1:5394142/1‑93 (MQ=255)
                             aaaaaCTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAggg                                                      <  1:406707/93‑1 (MQ=255)
                              aaaaCTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGc                                                     <  1:4817707/93‑1 (MQ=255)
                                accTGCGGTAGAGGGCGTTTTTTAACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGa                                                   <  2:3393888/91‑1 (MQ=255)
                                aaCTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGa                                                   <  2:4262016/93‑1 (MQ=255)
                                 acagCCGTAGAAGGTGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGAt                                                  <  1:5587286/90‑1 (MQ=255)
                                 aCTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGAt                                                  <  2:1676204/93‑1 (MQ=255)
                                  ctgctgTAGAATTCGGTTTTTCGCGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATc                                                 <  2:5791488/93‑1 (MQ=255)
                                  cTGCCGTAGAGGGCGTTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATc                                                 <  2:5420381/93‑1 (MQ=255)
                                  cTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATc                                                 >  1:3891166/1‑93 (MQ=255)
                                     ccGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAGCCACGCCGTTTAAGGGCGATCATg                                              >  2:5507328/1‑93 (MQ=255)
                                       gTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAg                                            >  2:2057553/1‑93 (MQ=255)
                                        tAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGc                                           >  2:605403/1‑93 (MQ=255)
                                         aGAAGGCGTTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGcc                                          <  2:2488476/93‑1 (MQ=255)
                                          gAAGTCGTTTTTTTACGAGC‑GGGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                                         <  1:4671838/93‑1 (MQ=255)
                                          gAAGGCGTTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                                         <  2:2340068/93‑1 (MQ=255)
                                          gAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                                         <  1:5197215/93‑1 (MQ=255)
                                          gAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                                         <  1:1439000/93‑1 (MQ=255)
                                            aGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGa                                       <  1:5043958/93‑1 (MQ=255)
                                            aGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGa                                       <  1:1177475/93‑1 (MQ=255)
                                             ggCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                                      >  1:5774497/1‑93 (MQ=255)
                                              gCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGata                                     >  2:3973437/1‑93 (MQ=255)
                                               gcgTTTTTTAGGCGCTGATCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGata                                     <  1:3284790/91‑1 (MQ=255)
                                               cGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGatat                                    >  2:6286858/1‑93 (MQ=255)
                                                ggTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATg                                   <  2:5069564/93‑1 (MQ=255)
                                                 gTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGa                                  <  1:2857940/93‑1 (MQ=255)
                                                   tttttACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAcc                                >  1:912060/1‑93 (MQ=255)
                                                    ttttACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAcca                               >  2:4544896/1‑92 (MQ=255)
                                                     tttACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTa                              <  2:3915091/93‑1 (MQ=255)
                                                      ttACGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTAc                             <  1:1361903/93‑1 (MQ=255)
                                                      taaCGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTAc                             <  2:3508881/91‑1 (MQ=255)
                                                        aCGAGC‑GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGt                           <  2:949222/93‑1 (MQ=255)
                                                         ggagc‑gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGTCCTACGtt                          <  2:2835680/92‑1 (MQ=255)
                                                           agc‑gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGa                        <  2:2728706/93‑1 (MQ=255)
                                                           agc‑gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGa                        >  1:3649183/1‑93 (MQ=255)
                                                            gc‑gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGaa                       <  1:3855459/93‑1 (MQ=255)
                                                             c‑gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAAt                      <  2:1081001/93‑1 (MQ=255)
                                                             c‑gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAAt                      <  2:404724/93‑1 (MQ=255)
                                                               gcgcACGACGGCAGTTGAAGCGGATGCCTGGCTGATTGCTGTGCCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATc                     <  1:5109218/91‑1 (MQ=255)
                                                               gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATc                     <  1:2707082/93‑1 (MQ=255)
                                                               gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATc                     <  2:5725955/93‑1 (MQ=255)
                                                                agcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCg                    >  1:4305192/1‑93 (MQ=255)
                                                                  cACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgc                  >  2:3831237/1‑93 (MQ=255)
                                                                  cACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgc                  <  1:4853293/93‑1 (MQ=255)
                                                                  cACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCAATCATGAGCCAGATATGACCTACGTTGAATCGgc                  <  1:4001241/93‑1 (MQ=255)
                                                                   acgacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgct                 >  2:5878593/1‑93 (MQ=255)
                                                                    cgacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgctg                >  2:1040259/1‑93 (MQ=255)
                                                                     gacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgctgc               <  2:4395448/93‑1 (MQ=255)
                                                                       cgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctc             <  1:3258406/93‑1 (MQ=255)
                                                                         ccAGTTGAAGCGGATGCCTGGCGGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcgc           <  1:4231308/93‑1 (MQ=255)
                                                                           aGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcgctc         <  2:434334/93‑1 (MQ=255)
                                                                            gTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTcc        >  1:4425059/1‑93 (MQ=255)
                                                                             ttGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCa       >  2:1332417/1‑93 (MQ=255)
                                                                              tGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCAt      <  2:4619116/93‑1 (MQ=255)
                                                                              tGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCAt      <  1:2186038/93‑1 (MQ=255)
                                                                              tGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCAt      <  1:2215366/93‑1 (MQ=255)
                                                                              tGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCAt      <  1:4772169/93‑1 (MQ=255)
                                                                               gAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCAtt     <  2:494710/93‑1 (MQ=255)
                                                                                aaGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTg    >  1:1324542/1‑93 (MQ=255)
                                                                                aaGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTg    <  1:1123409/93‑1 (MQ=255)
                                                                                 aGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTgc   <  2:2092446/93‑1 (MQ=255)
                                                                                 aGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTgc   <  2:5997346/93‑1 (MQ=255)
                                                                                 aGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTgc   <  2:3344211/93‑1 (MQ=255)
                                                                                 aGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTgc   >  1:2074558/1‑93 (MQ=255)
                                                                                 aGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTgc   >  1:1405474/1‑93 (MQ=255)
                                                                                        ||                                                                                     
ATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGC‑GAGCACGACGCCAGTTGAAGCGGATCGCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCG  >  NC_000913/3969016‑3969189

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATCGCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCA  >  NC_000913/3969010‑3969192
                                                                                             |                                                                                         
ATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATG                                                                                           <  SRR4302533.6260621/94‑1 (MQ=60)
ATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATG                                                                                           >  SRR4302533.2633935/1‑94 (MQ=60)
     TATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGG                                                                                      >  SRR4302533.832820/1‑94 (MQ=60)
     TATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGG                                                                                      <  SRR4302533.5242902/94‑1 (MQ=60)
     TATCGTCGAACCTGATTTGGCGAGCGTAGTAAAAACTGCCGTCGCAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGG                                                                                      >  SRR4302533.599203/1‑94 (MQ=60)
      ATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGC                                                                                     >  SRR4302533.3858135/1‑94 (MQ=60)
       TCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCT                                                                                    >  SRR4302533.5526949/1‑94 (MQ=60)
       TCGTCGAACCTGATTTGGCGAGTGAAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCT                                                                                    <  SRR4302533.303227/94‑1 (MQ=60)
        CGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTG                                                                                   <  SRR4302533.5982765/94‑1 (MQ=60)
        CGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTG                                                                                   >  SRR4302533.1542038/1‑94 (MQ=60)
        CGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTG                                                                                   >  SRR4302533.457498/1‑94 (MQ=60)
        CGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTG                                                                                   <  SRR4302533.1931929/94‑1 (MQ=60)
         GTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGA                                                                                  >  SRR4302533.2425934/1‑94 (MQ=60)
           CGACCCTGATTGGGCGAGTGTAGTAAAACCTGCGGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATT                                                                                <  SRR4302533.4567767/94‑1 (MQ=60)
           CGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATT                                                                                >  SRR4302533.296169/1‑94 (MQ=60)
           CGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATT                                                                                <  SRR4302533.1629374/94‑1 (MQ=60)
           CGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATT                                                                                <  SRR4302533.2919519/94‑1 (MQ=60)
           CGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCCGATT                                                                                >  SRR4302533.460966/1‑94 (MQ=60)
           CGAACCTGATTTGGCGAGTGTAGTAAAAACGGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATT                                                                                <  SRR4302533.2919520/94‑1 (MQ=60)
             AACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGC                                                                              >  SRR4302533.1454498/1‑94 (MQ=60)
              cctgttttgggggaatgtaaaaaaacCTGCCGTAGAGGGCGTTTTTTAACGAGCGGGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCT                                                                             <  SRR4302533.3257414/68‑1 (MQ=40)
              ACCTTATTTGGCGAGTGTAGTAAAACCTGCCGTAGAAGGCGTTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCT                                                                             <  SRR4302533.176538/94‑1 (MQ=60)
              ACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCT                                                                             >  SRR4302533.2974219/1‑94 (MQ=60)
              ACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCT                                                                             >  SRR4302533.6318176/1‑94 (MQ=60)
              ACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCT                                                                             <  SRR4302533.704773/94‑1 (MQ=60)
              ACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCT                                                                             >  SRR4302533.2220623/1‑94 (MQ=60)
               CCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTG                                                                            >  SRR4302533.1633084/1‑94 (MQ=60)
               CCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTG                                                                            >  SRR4302533.1328697/1‑94 (MQ=60)
                CTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGT                                                                           >  SRR4302533.101650/1‑94 (MQ=60)
                CTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGT                                                                           <  SRR4302533.6001896/94‑1 (MQ=60)
                 TGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTA                                                                          >  SRR4302533.4095737/1‑94 (MQ=60)
                 TGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTA                                                                          >  SRR4302533.846320/1‑94 (MQ=60)
                  GATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAC                                                                         <  SRR4302533.2100525/94‑1 (MQ=60)
                  GATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAC                                                                         >  SRR4302533.1887982/1‑94 (MQ=60)
                   ATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACC                                                                        >  SRR4302533.6114645/1‑94 (MQ=60)
                   ATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACC                                                                        <  SRR4302533.6264244/94‑1 (MQ=60)
                     TTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCA                                                                      >  SRR4302533.3895748/1‑94 (MQ=60)
                      TGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAC                                                                     <  SRR4302533.1389893/94‑1 (MQ=60)
                       GGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACG                                                                    >  SRR4302533.672655/1‑94 (MQ=60)
                        GCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGC                                                                   <  SRR4302533.4138515/94‑1 (MQ=60)
                         CGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCC                                                                  >  SRR4302533.4557121/1‑94 (MQ=60)
                          GAGTGTAGTAAAACCTGCGGTAGAAGGCGTTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCG                                                                 <  SRR4302533.5433410/94‑1 (MQ=60)
                          GAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGACAGTTGAAGCGGAGGTCTCGCTGATTGCTGTTTCCACGCCG                                                                 >  SRR4302533.2321490/1‑94 (MQ=60)
                           AGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAGCCACACCGT                                                                >  SRR4302533.1134921/1‑94 (MQ=60)
                            GTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTT                                                               <  SRR4302533.3176266/94‑1 (MQ=60)
                            GTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTT                                                               >  SRR4302533.1898040/1‑94 (MQ=60)
                             TGTAGTAAAAACTGCCGTAGAAGGCGTTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTT                                                              <  SRR4302533.2249721/94‑1 (MQ=60)
                             TGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTT                                                              >  SRR4302533.306729/1‑94 (MQ=60)
                             TGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTT                                                              >  SRR4302533.3101969/1‑94 (MQ=60)
                              GTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTA                                                             >  SRR4302533.6243089/1‑94 (MQ=60)
                                AGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAG                                                           >  SRR4302533.4740232/1‑94 (MQ=60)
                                 GTAAAAACTGCCGTAGAAGGCGTTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGG                                                          <  SRR4302533.3627946/94‑1 (MQ=60)
                                 GTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAACGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGG                                                          >  SRR4302533.5517819/1‑94 (MQ=60)
                                  TAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGG                                                         <  SRR4302533.4306492/94‑1 (MQ=60)
                                  TAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGG                                                         <  SRR4302533.284481/94‑1 (MQ=60)
                                  TAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGG                                                         >  SRR4302533.5426554/1‑94 (MQ=60)
                                   AAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGC                                                        >  SRR4302533.4254277/1‑94 (MQ=60)
                                   AAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGC                                                        <  SRR4302533.409151/94‑1 (MQ=60)
                                    AAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCG                                                       <  SRR4302533.2032590/94‑1 (MQ=60)
                                    AAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCG                                                       <  SRR4302533.4846908/94‑1 (MQ=60)
                                    AAAACTGCCGTAGAAGGCGGGTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAATGGCG                                                       >  SRR4302533.4741045/1‑94 (MQ=60)
                                       ACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATC                                                    >  SRR4302533.3914908/1‑94 (MQ=60)
                                       ACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATC                                                    <  SRR4302533.2716259/94‑1 (MQ=60)
                                       ACAGCCGTAGAAGGTGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATC                                                    <  SRR4302533.5621042/94‑1 (MQ=60)
                                           CCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGA                                                <  SRR4302533.4038199/94‑1 (MQ=60)
                                            CGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTACAGGCGATCATGAG                                               >  SRR4302533.439735/1‑94 (MQ=60)
                                            CGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAG                                               >  SRR4302533.1297033/1‑94 (MQ=60)
                                            CGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAG                                               >  SRR4302533.4457688/1‑94 (MQ=60)
                                            CGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAG                                               >  SRR4302533.1894059/1‑94 (MQ=60)
                                             GTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGC                                              >  SRR4302533.4093937/1‑94 (MQ=60)
                                             GTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGC                                              >  SRR4302533.5110221/1‑94 (MQ=60)
                                               AGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCA                                            <  SRR4302533.1745471/94‑1 (MQ=60)
                                                GAAGTCGTTTTTTTACGAGCGGGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAG                                           <  SRR4302533.4700146/94‑1 (MQ=60)
                                                GAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAG                                           <  SRR4302533.1447684/94‑1 (MQ=60)
                                                GAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAG                                           <  SRR4302533.5228557/94‑1 (MQ=60)
                                                  AGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGAT                                         <  SRR4302533.5074481/94‑1 (MQ=60)
                                                  AGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGAT                                         >  SRR4302533.5809369/1‑94 (MQ=60)
                                                  AGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGAT                                         <  SRR4302533.1184540/94‑1 (MQ=60)
                                                    gcgttttttaggcgctgatCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATAT                                       <  SRR4302533.3305065/75‑1 (MQ=60)
                                                    GCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATAT                                       <  SRR4302533.1366170/94‑1 (MQ=60)
                                                    GCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATAT                                       <  SRR4302533.395226/94‑1 (MQ=60)
                                                       GTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAC                                    <  SRR4302533.2875524/94‑1 (MQ=60)
                                                       GTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAC                                    >  SRR4302533.4928543/1‑94 (MQ=60)
                                                        TTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACC                                   >  SRR4302533.917410/1‑94 (MQ=60)
                                                          TTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTA                                 >  SRR4302533.4801955/1‑94 (MQ=60)
                                                           tataccgccgacCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTAC                                <  SRR4302533.4317047/82‑1 (MQ=60)
                                                           TTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTAC                                <  SRR4302533.5566590/94‑1 (MQ=60)
                                                            TTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACG                               <  SRR4302533.1370082/94‑1 (MQ=60)
                                                               CGAGGGATCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTG                            <  SRR4302533.5378271/94‑1 (MQ=60)
                                                               CGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTG                            <  SRR4302533.1097982/94‑1 (MQ=60)
                                                                GAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGA                           >  SRR4302533.3671621/1‑94 (MQ=60)
                                                                GAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGA                           >  SRR4302533.3907391/1‑94 (MQ=60)
                                                                  GCGAGCACGACGCCAGTTGAAGCTGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAAT                         <  SRR4302533.5120226/94‑1 (MQ=60)
                                                                  GCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAAT                         <  SRR4302533.4280200/94‑1 (MQ=60)
                                                                  GCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAAT                         <  SRR4302533.3879001/94‑1 (MQ=60)
                                                                   CGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATC                        <  SRR4302533.235364/94‑1 (MQ=60)
                                                                    GCGCACGACGGCAGTTGAAGCGGATGCCTGGCTGATTGCTGTGCCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCG                       <  SRR4302533.5140091/94‑1 (MQ=60)
                                                                    GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCG                       <  SRR4302533.2246107/94‑1 (MQ=60)
                                                                    GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCG                       >  SRR4302533.4331095/1‑94 (MQ=60)
                                                                    GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCG                       >  SRR4302533.2222469/1‑94 (MQ=60)
                                                                    GAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCG                       <  SRR4302533.2723690/94‑1 (MQ=60)
                                                                     AGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGG                      <  SRR4302533.4045642/94‑1 (MQ=60)
                                                                     AGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGG                      <  SRR4302533.5159951/94‑1 (MQ=60)
                                                                       gtCGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCT                    <  SRR4302533.4882704/92‑1 (MQ=60)
                                                                       CACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCT                    <  SRR4302533.4882705/94‑1 (MQ=60)
                                                                       CACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCAATCATGAGCCAGATATGACCTACGTTGAATCGGCT                    <  SRR4302533.4025605/94‑1 (MQ=60)
                                                                        ACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTCTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGCCCTACGTTGAATCGGCTG                   <  SRR4302533.5874225/94‑1 (MQ=60)
                                                                          GACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCT                 <  SRR4302533.3285446/94‑1 (MQ=60)
                                                                            CGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCG               <  SRR4302533.1373119/94‑1 (MQ=60)
                                                                            CGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCG               >  SRR4302533.3485415/1‑94 (MQ=60)
                                                                            CGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCG               <  SRR4302533.3278503/94‑1 (MQ=60)
                                                                             GCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGC              <  SRR4302533.5632038/94‑1 (MQ=60)
                                                                              CCAGTTGAAGCGGATGCCTGGCGGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCT             <  SRR4302533.4256884/94‑1 (MQ=60)
                                                                               CAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTC            >  SRR4302533.1426426/1‑94 (MQ=60)
                                                                                AGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCC           >  SRR4302533.4451807/1‑94 (MQ=60)
                                                                                 GTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCA          <  SRR4302533.4479580/94‑1 (MQ=60)
                                                                                   TGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATT        <  SRR4302533.2229040/94‑1 (MQ=60)
                                                                                   TGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATT        <  SRR4302533.2199512/94‑1 (MQ=60)
                                                                                   TGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATT        <  SRR4302533.4801089/94‑1 (MQ=60)
                                                                                    GAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTG       >  SRR4302533.1332486/1‑94 (MQ=60)
                                                                                     AAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGC      <  SRR4302533.417227/94‑1 (MQ=60)
                                                                                     AAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGC      >  SRR4302533.2087332/1‑94 (MQ=60)
                                                                                     AAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGC      <  SRR4302533.1490427/94‑1 (MQ=60)
                                                                                     AAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGC      >  SRR4302533.1413977/1‑94 (MQ=60)
                                                                                     AAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGC      <  SRR4302533.1130074/94‑1 (MQ=60)
                                                                                      AGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCG     <  SRR4302533.5418949/94‑1 (MQ=60)
                                                                                       GCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGC    <  SRR4302533.4865365/94‑1 (MQ=60)
                                                                                         GGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCA  >  SRR4302533.3136446/1‑94 (MQ=60)
                                                                                             |                                                                                         
ATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATCGCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCA  >  NC_000913/3969010‑3969192

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