Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 66 46.7 5280542 98.1% 5180211 47.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_017635 2,860,853 C→A R442S (AGG→AGT qseE ← two component system sensor histidine kinase QseE/GlrK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176352,860,8530CA100.0% 310.0 / NA 78R442S (AGG→AGTqseEtwo component system sensor histidine kinase QseE/GlrK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (38/40);  total (38/40)

CGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTT  >  NC_017635/2860810‑2860897
                                           |                                            
cgtcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTgg                                          >  2:413525/1‑48 (MQ=255)
 gtcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACt                                            <  2:1730334/45‑1 (MQ=255)
 gtcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTg                                           <  2:664285/46‑1 (MQ=255)
  tcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGc                                         >  1:1235561/1‑47 (MQ=255)
  tcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCa                                        <  2:144680/48‑1 (MQ=255)
  tcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCa                                        <  1:2126110/48‑1 (MQ=255)
  tcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCa                                        <  1:503900/48‑1 (MQ=255)
  tcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCa                                        <  2:104168/48‑1 (MQ=255)
  tcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCa                                        <  1:189170/48‑1 (MQ=255)
  tcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCa                                        <  1:1761683/48‑1 (MQ=255)
  tcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCa                                        <  2:116772/48‑1 (MQ=255)
  tcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCa                                        <  2:1439274/48‑1 (MQ=255)
  tcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCa                                        <  2:910336/48‑1 (MQ=255)
      ccAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGc                                    <  2:1481338/48‑1 (MQ=255)
       cAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCt                                   >  2:2265256/1‑48 (MQ=255)
        aGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCtt                                  <  2:1747410/48‑1 (MQ=255)
        aGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCtt                                  <  1:2443583/48‑1 (MQ=255)
         gATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTa                                 <  1:1767687/48‑1 (MQ=255)
          aTACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTaa                                >  1:1100381/1‑48 (MQ=255)
          aTACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTaa                                <  2:1663826/48‑1 (MQ=255)
           tACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAAt                               >  1:1025981/1‑48 (MQ=255)
           tACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAAt                               >  1:1385479/1‑48 (MQ=255)
           tACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAAt                               >  2:1309396/1‑48 (MQ=255)
           tACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAAt                               >  1:1737683/1‑48 (MQ=255)
           tACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAAt                               >  1:2597957/1‑48 (MQ=255)
           tACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAAt                               >  1:2547669/1‑48 (MQ=255)
            aCAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATc                              >  1:2018770/1‑48 (MQ=255)
            aCAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATc                              >  1:2142125/1‑48 (MQ=255)
             cAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATc                              <  1:1913166/47‑1 (MQ=255)
              aGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATcc                             >  2:2139962/1‑47 (MQ=255)
              aGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATcc                             >  1:590828/1‑47 (MQ=255)
              aGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATcc                             >  1:87656/1‑47 (MQ=255)
                ttCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATcc                             >  2:737597/1‑45 (MQ=255)
                ttCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATcc                             >  1:85478/1‑45 (MQ=255)
                 tCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGa                         >  2:875533/1‑48 (MQ=255)
                   cccTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGAcc                       >  2:166840/1‑48 (MQ=255)
                       tGCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTg                   >  2:576971/1‑48 (MQ=255)
                        gCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCt                    <  2:1520909/46‑1 (MQ=255)
                         cATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGcc                 <  1:2503024/48‑1 (MQ=255)
                         cATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGcc                 >  2:2397647/1‑48 (MQ=255)
                         cATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGcc                 <  1:1484351/48‑1 (MQ=255)
                         cATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGcc                 <  1:462518/48‑1 (MQ=255)
                         cATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGcc                 <  1:2588223/48‑1 (MQ=255)
                         cATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGcc                 <  2:599832/48‑1 (MQ=255)
                         cATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGc                  <  2:2082365/47‑1 (MQ=255)
                          aTACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGccc                <  1:2309141/48‑1 (MQ=255)
                          aTACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGccc                <  2:1424957/48‑1 (MQ=255)
                          aTACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGccc                <  2:1479997/48‑1 (MQ=255)
                          aTACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGccc                <  1:14593/48‑1 (MQ=255)
                          aTACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGccc                >  1:465234/1‑48 (MQ=255)
                           tACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCt               >  1:855236/1‑48 (MQ=255)
                           tACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCt               >  2:363720/1‑48 (MQ=255)
                            aCGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCtt              >  1:1901355/1‑48 (MQ=255)
                            aCGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCtt              >  1:1561547/1‑48 (MQ=255)
                            aCGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCtt              >  1:540694/1‑48 (MQ=255)
                            aCGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCtt              >  2:2107614/1‑48 (MQ=255)
                             cGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTc             <  2:2146050/48‑1 (MQ=255)
                              ggCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCa            <  1:872159/48‑1 (MQ=255)
                              ggCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCa            <  1:2167407/48‑1 (MQ=255)
                              ggCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCa            <  1:836787/48‑1 (MQ=255)
                                   aTACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcc       <  2:1882975/48‑1 (MQ=255)
                                   aTACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcc       >  1:1489305/1‑48 (MQ=255)
                                   aTACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcc       >  2:1619024/1‑48 (MQ=255)
                                   aTACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcc       >  2:2259413/1‑48 (MQ=255)
                                   aTACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcc       >  2:2364511/1‑48 (MQ=255)
                                    tACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGccc      >  1:1484501/1‑48 (MQ=255)
                                    tACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGccc      >  1:1709013/1‑48 (MQ=255)
                                    tACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGccc      >  2:2542570/1‑48 (MQ=255)
                                     aCAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcccc     >  1:2101762/1‑48 (MQ=255)
                                     aCAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcccc     >  1:17093/1‑48 (MQ=255)
                                      cAATCACTGGCACTGCTTAACCCCAGACCGCTGCCCTTCACCGcccc     <  1:736681/47‑1 (MQ=37)
                                      cAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcccc     <  2:155292/47‑1 (MQ=255)
                                      cAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcccc     <  1:357094/47‑1 (MQ=255)
                                      cAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcccc     <  1:2419317/47‑1 (MQ=255)
                                      cAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcccc     <  1:2510669/47‑1 (MQ=255)
                                       aaTCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCtt   <  2:854080/48‑1 (MQ=255)
                                       aaTCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCtt   <  1:1920077/48‑1 (MQ=255)
                                        aTCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCttt  >  1:1546887/1‑48 (MQ=255)
                                           |                                            
CGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTT  >  NC_017635/2860810‑2860897

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTG  >  NC_017635/2860804‑2860903
                                                 |                                                  
CGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCAC                                                   <  M01420:84:000000000‑ABVEH:1:2112:15201:11305/51‑1 (MQ=60)
 GCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACT                                                  >  M01420:84:000000000‑ABVEH:1:1102:21551:14547/1‑51 (MQ=60)
    CTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTG                                                 <  M01420:84:000000000‑ABVEH:1:1110:19429:16504/49‑1 (MQ=60)
    CTCGTCGACCAGATACAGTTCCCCTTGCATACGACGAATACAATCACTGGC                                               >  M01420:84:000000000‑ABVEH:1:2114:12261:2547/1‑51 (MQ=60)
      CGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCA                                              >  M01420:84:000000000‑ABVEH:1:1111:21017:16491/1‑50 (MQ=60)
       GTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAAT                                            <  M01420:84:000000000‑ABVEH:1:1105:10503:11643/51‑1 (MQ=60)
       GTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAAT                                            <  M01420:84:000000000‑ABVEH:1:2101:19728:27808/51‑1 (MQ=60)
       GTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAAT                                            <  M01420:84:000000000‑ABVEH:1:2105:7580:7282/51‑1 (MQ=60)
       GTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAAT                                            <  M01420:84:000000000‑ABVEH:1:1102:23895:18329/51‑1 (MQ=60)
          GACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCT                                         >  M01420:84:000000000‑ABVEH:1:2111:16887:9112/1‑51 (MQ=60)
             CAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAA                                      <  M01420:84:000000000‑ABVEH:1:2107:25932:24029/51‑1 (MQ=60)
              AGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAAT                                     <  M01420:84:000000000‑ABVEH:1:2102:15880:3383/51‑1 (MQ=60)
              AGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAAT                                     >  M01420:84:000000000‑ABVEH:1:1110:10919:13410/1‑51 (MQ=60)
               GATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATC                                    >  M01420:84:000000000‑ABVEH:1:2109:18603:7809/1‑51 (MQ=60)
               GATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATC                                    >  M01420:84:000000000‑ABVEH:1:1112:23463:22428/1‑51 (MQ=60)
               GATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATC                                    >  M01420:84:000000000‑ABVEH:1:2101:25133:21421/1‑51 (MQ=60)
               GATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATC                                    >  M01420:84:000000000‑ABVEH:1:2110:13181:2966/1‑51 (MQ=60)
               GATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATC                                    >  M01420:84:000000000‑ABVEH:1:2108:27933:20970/1‑51 (MQ=60)
               GATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATC                                    >  M01420:84:000000000‑ABVEH:1:1109:10086:21900/1‑51 (MQ=60)
                ATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCC                                   >  M01420:84:000000000‑ABVEH:1:2105:22409:10702/1‑51 (MQ=60)
                ATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCC                                   >  M01420:84:000000000‑ABVEH:1:2104:15807:8926/1‑51 (MQ=60)
                  ACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCCC                                  >  M01420:84:000000000‑ABVEH:1:1101:6098:21216/1‑50 (MQ=60)
                  ACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCCC                                  <  M01420:84:000000000‑ABVEH:1:2103:28433:10823/50‑1 (MQ=60)
                  ACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCCC                                  >  M01420:84:000000000‑ABVEH:1:1106:9972:5104/1‑50 (MQ=60)
                    AGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCCC                                  >  M01420:84:000000000‑ABVEH:1:1101:15422:20753/1‑48 (MQ=60)
                              GCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCT                     >  M01420:84:000000000‑ABVEH:1:1105:16393:1680/1‑51 (MQ=60)
                              GCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCT                     <  M01420:84:000000000‑ABVEH:1:2108:14643:12569/51‑1 (MQ=60)
                              GCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCT                     <  M01420:84:000000000‑ABVEH:1:1104:20780:27830/51‑1 (MQ=60)
                              GCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCT                     <  M01420:84:000000000‑ABVEH:1:2109:20386:5494/51‑1 (MQ=60)
                              GCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCT                     <  M01420:84:000000000‑ABVEH:1:1113:9401:18306/51‑1 (MQ=60)
                              GCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCC                      <  M01420:84:000000000‑ABVEH:1:2112:23845:9864/50‑1 (MQ=60)
                               CATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTT                    <  M01420:84:000000000‑ABVEH:1:1101:10136:6469/51‑1 (MQ=60)
                               CATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTT                    <  M01420:84:000000000‑ABVEH:1:2106:18006:20561/51‑1 (MQ=60)
                               CATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTT                    >  M01420:84:000000000‑ABVEH:1:1108:8495:11625/1‑51 (MQ=60)
                                ATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTC                   >  M01420:84:000000000‑ABVEH:1:1105:6952:18738/1‑51 (MQ=60)
                                ATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTC                   >  M01420:84:000000000‑ABVEH:1:2103:22014:8380/1‑51 (MQ=60)
                                ATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTC                   >  M01420:84:000000000‑ABVEH:1:1114:9805:10273/1‑51 (MQ=60)
                                   CGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACC                <  M01420:84:000000000‑ABVEH:1:2111:17154:13902/51‑1 (MQ=60)
                                   CGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACC                <  M01420:84:000000000‑ABVEH:1:1108:16330:14883/51‑1 (MQ=60)
                                   CGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACC                <  M01420:84:000000000‑ABVEH:1:2105:5796:15527/51‑1 (MQ=60)
                                   CGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACC                <  M01420:84:000000000‑ABVEH:1:1108:22061:7808/51‑1 (MQ=60)
                                       GAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCC            >  M01420:84:000000000‑ABVEH:1:1113:20023:19243/1‑51 (MQ=60)
                                        AATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCC           >  M01420:84:000000000‑ABVEH:1:2101:16638:15786/1‑51 (MQ=60)
                                        AATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCC           >  M01420:84:000000000‑ABVEH:1:1113:5365:18334/1‑51 (MQ=60)
                                         ATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCT          >  M01420:84:000000000‑ABVEH:1:2104:19901:26190/1‑51 (MQ=60)
                                         ATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCT          >  M01420:84:000000000‑ABVEH:1:2110:23199:20745/1‑51 (MQ=60)
                                         ATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCT          >  M01420:84:000000000‑ABVEH:1:1101:14390:7054/1‑51 (MQ=60)
                                           ATAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTT        >  M01420:84:000000000‑ABVEH:1:2112:4327:14234/1‑51 (MQ=60)
                                           ACAATCACTGGCACTGCTTAACCCCAGACCGCTGCCCTTCACCGCCCCTT         <  M01420:84:000000000‑ABVEH:1:1107:8076:11392/50‑1 (MQ=44)
                                           ACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTT         <  M01420:84:000000000‑ABVEH:1:2108:19803:13983/50‑1 (MQ=60)
                                           ACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTT         <  M01420:84:000000000‑ABVEH:1:2107:2442:18995/50‑1 (MQ=60)
                                           ACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTT         <  M01420:84:000000000‑ABVEH:1:1104:6729:4814/50‑1 (MQ=60)
                                            CAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTT       <  M01420:84:000000000‑ABVEH:1:2103:5937:12181/51‑1 (MQ=60)
                                            CAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTT       >  M01420:84:000000000‑ABVEH:1:1114:18588:7240/1‑51 (MQ=60)
                                                 ACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTG  >  M01420:84:000000000‑ABVEH:1:1102:27274:8304/1‑51 (MQ=60)
                                                 |                                                  
CGCTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTG  >  NC_017635/2860804‑2860903

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