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 F3 I160 R1 220 104.1 6051714 95.9% 5803593 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,180,586 (A)7→6 intergenic (+333/+77) yjaZ → / ← thiH heat‑shock protein C/2‑iminoacetate synthase ThiH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,180,5800A.100.0% 277.9 / NA 62intergenic (+327/+83)yjaZ/thiHheat‑shock protein C/2‑iminoacetate synthase ThiH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (21/41);  total (21/41)

TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCT  >  NZ_CP009273/4180492‑4180663
                                                                                        |                                                                                   
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                    <  1:116550/90‑1 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                    >  2:975743/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                    >  2:739825/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                    >  2:600095/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                    >  2:2502579/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                    >  2:2419158/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                    >  1:400527/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                    >  1:2703218/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                    >  1:2538833/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                    >  1:2495561/1‑90 (MQ=255)
 tGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaa                                                                                   >  2:1609108/1‑90 (MQ=255)
 tGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaa                                                                                   >  2:183252/1‑90 (MQ=255)
  gCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaa                                                                                  <  1:2448973/90‑1 (MQ=255)
  gCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaa                                                                                  <  1:2374219/90‑1 (MQ=255)
  gCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaa                                                                                  <  2:2105525/90‑1 (MQ=255)
         ggAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACaa                                                                          <  1:2730416/90‑1 (MQ=255)
         ggAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACaa                                                                          <  1:1526477/90‑1 (MQ=255)
         ggAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACaa                                                                          <  1:1636599/90‑1 (MQ=255)
         ggAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACaa                                                                          <  2:2445447/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                       <  1:2098976/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                       <  1:2039123/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                       <  1:1811633/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                       <  1:2857603/90‑1 (MQ=255)
             tttttAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCt                                                                      <  1:975743/90‑1 (MQ=255)
             tttttAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCt                                                                      <  1:278077/90‑1 (MQ=255)
              ttttAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCtt                                                                     <  1:1415729/90‑1 (MQ=255)
              ttttAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCtt                                                                     <  1:2650604/90‑1 (MQ=255)
                   aGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTcaac                                                                <  1:1484653/90‑1 (MQ=255)
                      aaGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCaacaac                                                             >  2:2339296/1‑90 (MQ=255)
                            gggCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                       >  2:1062306/1‑74 (MQ=255)
                            gggCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                       <  1:1062306/74‑1 (MQ=255)
                            gggCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGc                                                        >  2:2544592/1‑89 (MQ=255)
                            gggCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGc                                                        <  1:2544592/89‑1 (MQ=255)
                                 aaaCGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGttttt                                                  <  1:913183/90‑1 (MQ=255)
                                  aaCGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCtt                                                                     <  1:2663592/70‑1 (MQ=255)
                                  aaCGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCtt                                                                     >  2:2663592/1‑70 (MQ=255)
                                  aaCGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTc                                                 >  1:2171687/1‑90 (MQ=255)
                                   aCGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCa                                                <  1:509948/90‑1 (MQ=255)
                                      aGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCt                                             >  1:139616/1‑90 (MQ=255)
                                      aGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCt                                             >  2:1667001/1‑90 (MQ=255)
                                      aGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCt                                             <  1:600095/90‑1 (MQ=255)
                                      aGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCt                                             <  2:1440719/90‑1 (MQ=255)
                                                 tttCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGc                                  >  1:1680361/1‑90 (MQ=255)
                                                   tCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc                                <  2:2553627/90‑1 (MQ=255)
                                                   tCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc                                <  2:18614/90‑1 (MQ=255)
                                                   tCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc                                <  2:1516279/90‑1 (MQ=255)
                                                   tCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc                                <  2:1260127/90‑1 (MQ=255)
                                                   tCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc                                <  1:739825/90‑1 (MQ=255)
                                                   tCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc                                <  1:417215/90‑1 (MQ=255)
                                                         tGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATa                          >  2:1016953/1‑90 (MQ=255)
                                                         tGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATa                          >  2:359214/1‑90 (MQ=255)
                                                         tGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATa                          >  1:1597701/1‑90 (MQ=255)
                                                               acaTTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTa                    >  2:1687185/1‑90 (MQ=255)
                                                                caTTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTAc                   <  1:68660/90‑1 (MQ=255)
                                                                caTTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTAc                   <  1:2502579/90‑1 (MQ=255)
                                                                caTTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTAc                   <  2:2171687/90‑1 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt                 <  1:1894434/90‑1 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt                 <  1:2339296/90‑1 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt                 <  1:2325485/90‑1 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTAc                   <  2:1912428/86‑1 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTAc                   >  1:1912428/1‑86 (MQ=255)
                                                                      tgctgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCtt             >  2:276169/1‑90 (MQ=255)
                                                                       gctgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTa            >  1:1206961/1‑90 (MQ=255)
                                                                         tgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTTCa          >  2:2549128/1‑90 (MQ=255)
                                                                         tgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACa          >  2:438984/1‑90 (MQ=255)
                                                                          gACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTTTTACAt         >  1:213220/1‑90 (MQ=255)
                                                                          gACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAt         <  1:276169/90‑1 (MQ=255)
                                                                          gACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAt         <  1:2559090/90‑1 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc       >  2:1509564/1‑90 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc       <  2:295773/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg     <  2:295309/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg     <  2:451786/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg     <  1:2938118/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg     <  1:2720415/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg     <  1:1703962/90‑1 (MQ=255)
                                                                                 cTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGtct  >  1:921221/1‑90 (MQ=255)
                                                                                 cTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGtct  >  1:865736/1‑90 (MQ=255)
                                                                                        |                                                                                   
TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCT  >  NZ_CP009273/4180492‑4180663

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAGGCGCGTC  >  NZ_CP009273/4180487‑4180685
                                                                                                   |                                                                                                   
tataagagacagACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAC                                                                                                     <  SRR3722030.2484904/88‑1 (MQ=60)
  ACTTTGCATTACGGAATTTTTAAGAAAGGTAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACAC                                                                                                  <  SRR3722030.30598/100‑1 (MQ=60)
  ACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACAC                                                                                                  <  SRR3722030.2084821/100‑1 (MQ=60)
  ACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACAC                                                                                                  <  SRR3722030.2928475/100‑1 (MQ=60)
  ACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACAC                                                                                                  <  SRR3722030.136882/100‑1 (MQ=60)
     TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAAC                                                                                               <  SRR3722030.118246/100‑1 (MQ=60)
       GCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGC                                                                                             <  SRR3722030.2417391/100‑1 (MQ=60)
       GCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGC                                                                                             <  SRR3722030.2493688/100‑1 (MQ=60)
              GGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAAC                                                                                      <  SRR3722030.1663166/100‑1 (MQ=60)
              GGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAAC                                                                                      <  SRR3722030.2780834/100‑1 (MQ=60)
              GGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAAC                                                                                      <  SRR3722030.1551105/100‑1 (MQ=60)
                 agacagAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAAC                                                                                   <  SRR3722030.1078256/94‑1 (MQ=60)
                 ATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAAC                                                                                   <  SRR3722030.2910380/100‑1 (MQ=60)
                 ATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAAC                                                                                   <  SRR3722030.2136374/100‑1 (MQ=60)
                 ATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAAC                                                                                   <  SRR3722030.1842225/100‑1 (MQ=60)
                 ATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAAC                                                                                   <  SRR3722030.2075213/100‑1 (MQ=60)
                  TTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACC                                                                                  <  SRR3722030.282109/100‑1 (MQ=60)
                  TTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACC                                                                                  <  SRR3722030.990597/100‑1 (MQ=60)
                   taagagacagGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCT                                                                                 <  SRR3722030.2712669/90‑1 (MQ=60)
                   TTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCT                                                                                 <  SRR3722030.1438076/100‑1 (MQ=60)
                   TTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCT                                                                                 <  SRR3722030.2699411/100‑1 (MQ=60)
                        AGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGT                                                                            <  SRR3722030.1508418/100‑1 (MQ=60)
                             GGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTC                                                                       >  SRR3722030.2210624/1‑100 (MQ=60)
                                AGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATC                                                                    <  SRR3722030.2591199/100‑1 (MQ=60)
                                 GGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCT                                                                   >  SRR3722030.141627/1‑100 (MQ=60)
                                      AAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGT                                                              <  SRR3722030.927179/100‑1 (MQ=60)
                                        ACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAG                                                            <  SRR3722030.517763/100‑1 (MQ=60)
                                           AGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGG                                                         <  SRR3722030.609391/100‑1 (MQ=60)
                                            GGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGC                                                        >  SRR3722030.1707646/1‑100 (MQ=60)
                                                    CTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATA                                                >  SRR3722030.1623618/1‑100 (MQ=60)
                                                        cgtcagatgtgtataagagacagGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATT                                            <  SRR3722030.1981479/77‑1 (MQ=60)
                                                        TCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATT                                            <  SRR3722030.423514/100‑1 (MQ=60)
                                                        TCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATT                                            <  SRR3722030.751231/100‑1 (MQ=60)
                                                               GGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTC                                     >  SRR3722030.1945553/1‑100 (MQ=60)
                                                                  CAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTA                                  >  SRR3722030.1224798/1‑100 (MQ=60)
                                                                     CATTCCTGCTGACGGCACTTCCATAAAAAACACAACGCTTTCAACAACCTgtctcttatacacatctgacgctgccgacgacatggccagtgtagatctc                               >  SRR3722030.962061/1‑50 (MQ=54)
                                                                     CATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACctgtctcttatacacatctgacgctgccgacgacatggccagtgtagatctcggt                               >  SRR3722030.1532980/1‑45 (MQ=60)
                                                                     CATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACctgtctcttatacacatctgacgctgccgacgacatggccagtgtagatctcggt                               >  SRR3722030.446320/1‑45 (MQ=60)
                                                                     CATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTTTTACAT                               >  SRR3722030.216268/1‑100 (MQ=60)
                                                                     CATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAT                               <  SRR3722030.2548379/100‑1 (MQ=60)
                                                                     CATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAT                               <  SRR3722030.69707/100‑1 (MQ=60)
                                                                       TTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTC                             <  SRR3722030.2367692/100‑1 (MQ=60)
                                                                       TTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTC                             <  SRR3722030.1927113/100‑1 (MQ=60)
                                                                       TTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTC                             <  SRR3722030.2381798/100‑1 (MQ=60)
                                                                            GCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTTCATTCCGTCT                        >  SRR3722030.935332/1‑100 (MQ=60)
                                                                            GCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCT                        >  SRR3722030.879028/1‑100 (MQ=60)
                                                                               GACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCAT                     <  SRR3722030.280176/100‑1 (MQ=60)
                                                                               GACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCAT                     <  SRR3722030.2606003/100‑1 (MQ=60)
                                                                                 CGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCctgtctcttatacacatctgacgctgccgacgaca                   >  SRR3722030.1176834/1‑65 (MQ=60)
                                                                                   GCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTC                 <  SRR3722030.1731736/100‑1 (MQ=60)
                                                                                   GCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTC                 <  SRR3722030.2992269/100‑1 (MQ=60)
                                                                                   GCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTC                 <  SRR3722030.2770678/100‑1 (MQ=60)
                                                                                    gtgtataagagacagACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTC                 <  SRR3722030.2316215/85‑1 (MQ=60)
                                                                                        ACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACCTTCTTACATTCCGTGTCATTGTCTTTGC            >  SRR3722030.2662678/1‑100 (MQ=60)
                                                                                         CCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCG           <  SRR3722030.2536042/100‑1 (MQ=60)
                                                                                         CCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCG           <  SRR3722030.2665070/100‑1 (MQ=60)
                                                                                         CCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCG           <  SRR3722030.2630164/100‑1 (MQ=60)
                                                                                            atAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAG         >  SRR3722030.2942488/3‑100 (MQ=60)
                                                                                            atAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAG         >  SRR3722030.1973487/3‑100 (MQ=60)
                                                                                             TAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATGGTCTTTGCGAGG        <  SRR3722030.2306029/100‑1 (MQ=60)
                                                                                              AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTAAGTTCTTACATTCCGTCTCATAGTCctgtctctta       >  SRR3722030.1162341/1‑90 (MQ=60)
                                                                                                   ACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAGGCGCGTC  >  SRR3722030.2867202/1‑100 (MQ=60)
                                                                                                   |                                                                                                   
AGACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAGGCGCGTC  >  NZ_CP009273/4180487‑4180685

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