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 F2 I204 R1 2 48.7 2644746 96.6% 2554824 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,609,448 A→G D10G (GAC→GGC)  deoB → phosphopentomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,609,4480AG97.3% 228.3 / ‑1.2 73D10G (GAC→GGC) deoBphosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (44/27);  total (45/28)
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.18e-01

GTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTG  >  NZ_CP009273/4609371‑4609532
                                                                             |                                                                                    
gTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCAt                                                                          >  2:97772/1‑90 (MQ=255)
 tataCTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATc                                                                         <  1:97772/90‑1 (MQ=255)
  ataCTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCg                                                                        >  2:428440/1‑90 (MQ=255)
    aCTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgc                                                                      <  2:698503/90‑1 (MQ=255)
    aCTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgc                                                                      <  2:509457/90‑1 (MQ=255)
       gatctgatcATTTTAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  2:752289/1‑90 (MQ=255)
       gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgcgc                                                                    <  2:554117/89‑1 (MQ=255)
       gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgcgc                                                                    >  1:554117/1‑89 (MQ=255)
       gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgc                                                                      <  2:462517/87‑1 (MQ=255)
       gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgc                                                                      >  1:462517/1‑87 (MQ=255)
       gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  1:1052259/1‑90 (MQ=255)
       gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  2:997572/1‑90 (MQ=255)
       gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  2:804278/1‑90 (MQ=255)
       gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  2:1162141/1‑90 (MQ=255)
       gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  1:95411/1‑90 (MQ=255)
       gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  1:426074/1‑90 (MQ=255)
       gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  1:1191485/1‑90 (MQ=255)
       gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  1:122633/1‑90 (MQ=255)
          ctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACa                                                                >  2:474472/1‑90 (MQ=255)
            gatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACaga                                                              >  2:1133219/1‑90 (MQ=255)
            gatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACaga                                                              >  1:780312/1‑90 (MQ=255)
            gatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACaga                                                              >  2:368711/1‑90 (MQ=255)
             atcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACagaa                                                             <  2:495740/90‑1 (MQ=255)
                 tttAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATg                                                         >  1:72407/1‑90 (MQ=255)
                 tttAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATg                                                         >  1:1286023/1‑90 (MQ=255)
                    aaaTTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAg                                                      <  1:908635/90‑1 (MQ=255)
                          gAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCt                                                >  1:390620/1‑90 (MQ=255)
                                        ggAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCggg                                  <  1:1278715/90‑1 (MQ=255)
                                        ggAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCggg                                  <  2:1278387/90‑1 (MQ=255)
                                        ggAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCggg                                  <  2:126436/90‑1 (MQ=255)
                                         gagaACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACagaa                                                             >  1:897588/1‑62 (MQ=255)
                                         gagaACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACagaa                                                             <  2:897588/62‑1 (MQ=255)
                                         gagaACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAGGATGCAGAACGCTTTGGTGACGTCgggg                                 >  1:945991/1‑90 (MQ=255)
                                           gaACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCt                               >  1:1042890/1‑90 (MQ=255)
                                            aaCATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTg                              >  1:910538/1‑90 (MQ=255)
                                              cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgcg                                                                     <  1:7353/49‑1 (MQ=255)
                                              cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgcg                                                                     >  2:7353/1‑49 (MQ=255)
                                              cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCa                                                       <  1:144685/63‑1 (MQ=255)
                                              cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCa                                                       >  2:144685/1‑63 (MQ=255)
                                              cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            <  1:1162141/90‑1 (MQ=255)
                                              cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            >  2:166513/1‑90 (MQ=255)
                                              cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            >  2:497239/1‑90 (MQ=255)
                                              cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            >  1:1253038/1‑90 (MQ=255)
                                              cATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            <  2:1280198/90‑1 (MQ=255)
                                              cATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            >  1:1280198/1‑90 (MQ=255)
                                                tatGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgc                                                                      >  2:1086092/1‑46 (MQ=255)
                                                tatGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgc                                                                      <  1:1086092/46‑1 (MQ=255)
                                                tatGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGacac                          >  1:657221/1‑90 (MQ=255)
                                                 atGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACAcc                         >  2:1270424/1‑90 (MQ=255)
                                                  tGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACAccc                        <  2:72507/90‑1 (MQ=255)
                                                  tGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACAccc                        <  1:166513/90‑1 (MQ=255)
                                                  tGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACAccc                        <  2:159556/90‑1 (MQ=255)
                                                    aaaCGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTg                      <  1:497239/90‑1 (MQ=255)
                                                        gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTggtca                  >  2:212256/1‑87 (MQ=255)
                                                        gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTggggc                  >  2:602667/1‑88 (MQ=255)
                                                        gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTc                  >  2:283989/1‑90 (MQ=255)
                                                        gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTc                  >  2:209564/1‑90 (MQ=255)
                                                        gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTc                  >  2:168117/1‑90 (MQ=255)
                                                        gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTc                  >  1:820878/1‑90 (MQ=255)
                                                        gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTc                  >  1:885436/1‑90 (MQ=255)
                                                          gCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCat                <  2:95411/90‑1 (MQ=255)
                                                          gCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCat                <  2:1088454/90‑1 (MQ=255)
                                                            aTTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCatat              <  2:910652/90‑1 (MQ=255)
                                                                 ttatGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAg         <  1:184144/90‑1 (MQ=255)
                                                                 ttatGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAg         <  1:564264/90‑1 (MQ=255)
                                                                 ttatGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAg         <  2:1191485/90‑1 (MQ=255)
                                                                    tGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAg      >  2:347184/1‑90 (MQ=255)
                                                                      gTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCt    >  1:737807/1‑90 (MQ=255)
                                                                       tGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCtt   >  1:555354/1‑90 (MQ=255)
                                                                       tGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCtt   <  1:428440/90‑1 (MQ=255)
                                                                        gCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            >  2:657655/1‑64 (MQ=255)
                                                                        gCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            <  1:657655/64‑1 (MQ=255)
                                                                        gCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTtg  >  1:587755/1‑90 (MQ=255)
                                                                             |                                                                                    
GTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTG  >  NZ_CP009273/4609371‑4609532

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGC  >  NZ_CP009273/4609350‑4609547
                                                                                                  |                                                                                                   
TCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGC                                                                                                    <  SRR3722086.1253442/100‑1 (MQ=60)
                  ACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCT                                                                                  >  SRR3722086.96365/1‑100 (MQ=60)
                  ACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCT                                                                                  >  SRR3722086.1065801/1‑100 (MQ=60)
                  ACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCT                                                                                  >  SRR3722086.559686/1‑100 (MQ=60)
                  ACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCT                                                                                  >  SRR3722086.467258/1‑100 (MQ=60)
                  ACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCT                                                                                  >  SRR3722086.430496/1‑100 (MQ=60)
                  ACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCT                                                                                  >  SRR3722086.123838/1‑100 (MQ=60)
                  ACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCT                                                                                  >  SRR3722086.1207041/1‑100 (MQ=60)
                      TATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAG                                                                              <  SRR3722086.98742/100‑1 (MQ=60)
                       ATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGA                                                                             >  SRR3722086.789259/1‑100 (MQ=60)
                       ATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACctgtctcttatacacatctccga                                                                             >  SRR3722086.101061/1‑77 (MQ=60)
                         gtctcgtgggctcggagatgtgtataagagacagACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAG                                                                           <  SRR3722086.1100144/66‑1 (MQ=60)
                          tcgtgggctcggagatgtgtataagagacagGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGA                                                                          <  SRR3722086.7405/69‑1 (MQ=60)
                            GATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATG                                                                        >  SRR3722086.73121/1‑100 (MQ=60)
                            GATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATG                                                                        >  SRR3722086.1302846/1‑100 (MQ=60)
                                     ATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCT                                                               >  SRR3722086.394702/1‑100 (MQ=60)
                                        atgtgtataagagacagGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTG                                                            <  SRR3722086.146136/83‑1 (MQ=60)
                                         AAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGG                                                           <  SRR3722086.919947/100‑1 (MQ=60)
                                                    ACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAGGATGCAGAACGCTTTGGTGACGTCGGGG                                                >  SRR3722086.957970/1‑100 (MQ=60)
                                                    ACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACctgtctcttatacacatc                                                >  SRR3722086.908703/1‑82 (MQ=60)
                                                      TGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCT                                              >  SRR3722086.1056321/1‑100 (MQ=60)
                                                       GAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTG                                             >  SRR3722086.921886/1‑100 (MQ=60)
                                                         GTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGAC                                           >  SRR3722086.1269458/1‑100 (MQ=60)
                                                         GTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGAC                                           >  SRR3722086.1296953/1‑100 (MQ=60)
                                                           ACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACAC                                         >  SRR3722086.664480/1‑100 (MQ=60)
                                                             GGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCC                                       <  SRR3722086.1295447/100‑1 (MQ=60)
                                                                   tgtgtataagagacagTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTC                                 <  SRR3722086.664919/84‑1 (MQ=60)
                                                                   CATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTC                                 <  SRR3722086.1177299/100‑1 (MQ=60)
                                                                   CATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTC                                 >  SRR3722086.896347/1‑100 (MQ=60)
                                                                   CATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTC                                 >  SRR3722086.830571/1‑100 (MQ=60)
                                                                       TGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATAT                             <  SRR3722086.168217/100‑1 (MQ=60)
                                                                         AAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCG                           <  SRR3722086.502283/100‑1 (MQ=60)
                                                                                 ATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCT                   >  SRR3722086.746109/1‑100 (MQ=60)
                                                                                  TTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTT                  >  SRR3722086.560933/1‑100 (MQ=60)
                                                                                   TTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTG                 >  SRR3722086.593865/1‑100 (MQ=60)
                                                                                      TTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCtgactcttatacac              >  SRR3722086.1201428/1‑86 (MQ=60)
                                                                                      TTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGC              <  SRR3722086.569950/100‑1 (MQ=60)
                                                                                      TTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGC              >  SRR3722086.1142396/1‑100 (MQ=60)
                                                                                      TTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGC              <  SRR3722086.186039/100‑1 (MQ=60)
                                                                                        ATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCA            >  SRR3722086.220793/1‑100 (MQ=60)
                                                                                           GTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAG         >  SRR3722086.294320/1‑100 (MQ=60)
                                                                                           GTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAG         >  SRR3722086.276060/1‑100 (MQ=60)
                                                                                            TGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGG        <  SRR3722086.432882/100‑1 (MQ=60)
                                                                                                  GCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGC  <  SRR3722086.1147995/100‑1 (MQ=60)
                                                                                                  GCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGC  <  SRR3722086.968939/100‑1 (MQ=60)
                                                                                                  |                                                                                                   
TCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGC  >  NZ_CP009273/4609350‑4609547

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