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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

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.94.8% 214.4 / 8.5 58intergenic (+327/+83)yjaZ/thiHheat‑shock protein C/2‑iminoacetate synthase ThiH
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base . (20/35);  total (21/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCG  >  NZ_CP009273/4180492‑4180660
                                                                                        |                                                                                
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  1:1004676/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  2:811736/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  2:1499589/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  2:798188/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  1:1352980/1‑90 (MQ=255)
  gCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaa                                                                               <  1:761728/90‑1 (MQ=255)
  gCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaa                                                                               <  1:1662/90‑1 (MQ=255)
         ggAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACaa                                                                       <  2:49123/90‑1 (MQ=255)
         ggAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACaa                                                                       <  2:1296105/90‑1 (MQ=255)
         ggAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACaa                                                                       <  1:1407858/90‑1 (MQ=255)
          gAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAAc                                                                      >  2:1220392/1‑90 (MQ=255)
          gAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAAc                                                                      >  2:1366278/1‑90 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  1:1420100/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:1235863/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:872700/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  1:760975/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:1396252/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:1518928/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:72931/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:49963/90‑1 (MQ=255)
              ttttAAGAAAGGCAGGGGGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCtt                                                                  <  1:40516/90‑1 (MQ=255)
              ttttAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCtt                                                                  <  2:425258/90‑1 (MQ=255)
               tttAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCttt                                                                 >  2:661698/1‑90 (MQ=255)
                   aGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTcaac                                                             <  1:469940/90‑1 (MQ=255)
                        ggCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCGCTACCATAAAAAACACAACGCTTTCAACAACCt                                                        <  2:310104/90‑1 (MQ=255)
                        ggCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCt                                                        >  1:1273462/1‑90 (MQ=255)
                              gCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGtt                                                  <  1:1133263/90‑1 (MQ=255)
                                   aCGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCa                                             >  2:17840/1‑90 (MQ=255)
                                      aGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCt                                          >  2:990021/1‑90 (MQ=255)
                                              gCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGc                                  <  1:1578503/90‑1 (MQ=255)
                                                   tCGCCCTGGCCAACCTCCCTGCCGACGGCACTACCATAAAAAACACACCGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc                             <  2:1174703/90‑1 (MQ=255)
                                                   tCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACACCCTTGCGTTTTTCATCTTCCGTAGCGGCTc                             <  1:990021/90‑1 (MQ=255)
                                                   tCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc                             <  1:1366278/90‑1 (MQ=255)
                                                        cTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:640549/46‑1 (MQ=255)
                                                        cTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    >  1:640549/1‑46 (MQ=255)
                                                         tGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGc                                  >  1:303273/1‑79 (MQ=255)
                                                         tGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGc                                  <  2:303273/79‑1 (MQ=255)
                                                               acaTTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTa                 >  1:1036805/1‑90 (MQ=255)
                                                               acaTTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTa                 >  1:325151/1‑90 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAGCCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt              <  2:1064395/90‑1 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt              <  2:294209/90‑1 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt              <  1:1499589/90‑1 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAAACACAACGCTTTcaac                                                             >  1:1351025/1‑42 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAAACACAACGCTTTcaac                                                             <  2:1351025/42‑1 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTc            >  2:206942/1‑90 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTc            >  2:1283952/1‑90 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCCGGATAAATTTACGTTc            >  1:729552/1‑90 (MQ=255)
                                                                      tgctgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCtt          <  2:95376/90‑1 (MQ=255)
                                                                       gctgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTa         >  1:1201301/1‑90 (MQ=255)
                                                                         tgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACa       >  2:220467/1‑90 (MQ=255)
                                                                         tgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACa       >  2:611106/1‑90 (MQ=255)
                                                                          gACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAt      >  2:1567987/1‑90 (MQ=255)
                                                                          gACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAt      <  2:203548/90‑1 (MQ=255)
                                                                           aCGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAtt     >  2:718647/1‑90 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATTTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc    <  1:942384/90‑1 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc    <  1:718647/90‑1 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc    <  1:766667/90‑1 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc    >  1:619462/1‑90 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc    >  2:1464541/1‑90 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  1:234680/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  1:295907/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  1:940751/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  2:156199/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  1:111852/90‑1 (MQ=255)
                                                                               cACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  2:1352980/90‑1 (MQ=255)
                                                                                        |                                                                                
TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCG  >  NZ_CP009273/4180492‑4180660

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 31 ≤ 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
ACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAGGCGCGTC  >  NZ_CP009273/4180489‑4180685
                                                                                                 |                                                                                                   
ACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACAC                                                                                                  <  SRR3722078.96170/100‑1 (MQ=60)
     GCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGC                                                                                             <  SRR3722078.1678/100‑1 (MQ=60)
     GCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGC                                                                                             <  SRR3722078.771519/100‑1 (MQ=60)
            GGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAAC                                                                                      <  SRR3722078.1430697/100‑1 (MQ=60)
               ATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAAC                                                                                   <  SRR3722078.1443161/100‑1 (MQ=60)
               ATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAAC                                                                                   <  SRR3722078.770755/100‑1 (MQ=60)
                 TTTTAAGAAAGGCAGGGGGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACCACCT                                                                                 <  SRR3722078.40996/100‑1 (MQ=60)
                 TTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCT                                                                                 >  SRR3722078.1293683/1‑100 (MQ=60)
                      AGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGT                                                                            <  SRR3722078.475527/100‑1 (MQ=60)
                                 GCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTC                                                                 <  SRR3722078.1150621/100‑1 (MQ=60)
                                                 GCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTgtctcttatacacatctgacgctgccgacgat                                                 >  SRR3722078.647919/1‑68 (MQ=60)
                                                 GCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAAT                                                 <  SRR3722078.1604193/100‑1 (MQ=60)
                                                  CTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATC                                                >  SRR3722078.306973/1‑100 (MQ=60)
                                                      agcgtcagatgtgtataagagacagCGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATT                                            <  SRR3722078.398159/75‑1 (MQ=60)
                                                      TCGCCCTGGACAACATTCCTGCTGACGGCCCTACCATAAAAAACACAACGCTTTCAACACCCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATT                                            <  SRR3722078.1004393/100‑1 (MQ=60)
                                                      TCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATT                                            <  SRR3722078.1388261/100‑1 (MQ=60)
                                                        GCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTA                                          >  SRR3722078.1052240/1‑100 (MQ=60)
                                                        GCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTA                                          >  SRR3722078.329139/1‑100 (MQ=60)
                                                             GGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTctgtctcttatacacatctgacgctgccgacgatggcg                                      >  SRR3722078.1372722/1‑62 (MQ=60)
                                                             GGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCCGGATAAATTTACGTTC                                     >  SRR3722078.738668/1‑100 (MQ=60)
                                                                CAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTA                                  >  SRR3722078.1220145/1‑100 (MQ=60)
                                                                     TTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTC                             <  SRR3722078.1524042/100‑1 (MQ=60)
                                                                     TTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTC                             >  SRR3722078.626600/1‑100 (MQ=60)
                                                                               CGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATTTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAG                   <  SRR3722078.955503/100‑1 (MQ=60)
                                                                               CGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAG                   >  SRR3722078.1117905/1‑100 (MQ=60)
                                                                               CGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAG                   <  SRR3722078.776552/100‑1 (MQ=60)
                                                                               CGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAG                   <  SRR3722078.727526/100‑1 (MQ=60)
                                                                                 GCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTC                 <  SRR3722078.237538/100‑1 (MQ=60)
                                                                                 GCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTC                 <  SRR3722078.953819/100‑1 (MQ=60)
                                                                                 GCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTC                 <  SRR3722078.299517/100‑1 (MQ=60)
                                                                                 GCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTC                 <  SRR3722078.113173/100‑1 (MQ=60)
                                                                                  taagagacagAAAAAACACAACACTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTC                 <  SRR3722078.123003/90‑1 (MQ=60)
                                                                                      ACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGC            >  SRR3722078.544859/1‑100 (MQ=60)
                                                                                      ACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGGCTCATAGTCTTTGC            >  SRR3722078.1320868/1‑100 (MQ=60)
                                                                                       CCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCG           <  SRR3722078.1024603/100‑1 (MQ=60)
                                                                                       CCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCG           <  SRR3722078.1018694/100‑1 (MQ=60)
                                                                                                 ACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAGGCGCGTC  >  SRR3722078.917928/1‑100 (MQ=60)
                                                                                                 |                                                                                                   
ACTTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAGGCGCGTC  >  NZ_CP009273/4180489‑4180685

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