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
A2 F5 I176 R1 180 42.9 1710034 96.5% 1650182 119.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 640,333 A→C W43G (TGG→GGG)  rna ← ribonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273640,3330AC100.0% 173.5 / NA 47W43G (TGG→GGG) rnaribonuclease I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (23/24);  total (23/24)

CGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCAGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGG  >  NZ_CP009273/640199‑640464
                                                                                                                                      |                                                                                                                                   
cGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  2:38925/139‑1 (MQ=255)
cGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  1:26847/139‑1 (MQ=255)
cGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  1:653217/139‑1 (MQ=255)
cGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  2:689615/139‑1 (MQ=255)
cGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  2:45417/139‑1 (MQ=255)
cGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  2:517073/139‑1 (MQ=255)
cGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  2:549605/139‑1 (MQ=255)
  aTTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGagag                                                                                                                               >  1:554958/1‑139 (MQ=255)
      aGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGc                                                                                                                            <  2:23519/138‑1 (MQ=255)
          aaTCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGa                                                                                                                       <  1:112268/139‑1 (MQ=255)
           aTCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTCACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGa                                                                                                                       >  2:20739/1‑138 (MQ=255)
                 gCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACTTAGc                                                                                                                 >  1:436095/1‑138 (MQ=255)
                      aGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCa                                                                                                            >  2:463277/1‑138 (MQ=255)
                        aCCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCaaaa                                                                                                         <  1:31196/139‑1 (MQ=255)
                        aCCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCaaaa                                                                                                         >  1:247115/1‑139 (MQ=255)
                            tGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGc                                                                                                                            <  1:447509/116‑1 (MQ=255)
                            tGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGc                                                                                                                            >  2:447509/1‑116 (MQ=255)
                            tGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGc                                                                                                     >  1:363810/1‑139 (MQ=255)
                                  gTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  1:409434/105‑1 (MQ=255)
                                  gTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 >  2:409434/1‑105 (MQ=255)
                                  gTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACt                                                                                               >  2:843550/1‑139 (MQ=255)
                                   tCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTg                                                                                              >  2:453758/1‑139 (MQ=255)
                                      gAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  1:422660/101‑1 (MQ=255)
                                      gAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 >  2:422660/1‑101 (MQ=255)
                                       aaaaTCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACTAGCGATCAAAATCGCCATACTGTTTTg                                                                                         >  1:127351/1‑139 (MQ=255)
                                             aGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGc                                                                                    <  2:510977/139‑1 (MQ=255)
                                                         cGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCg                                                                        >  2:585725/1‑139 (MQ=255)
                                                                gTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCag                                                                 <  2:247115/139‑1 (MQ=255)
                                                                   tGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAggg                                                                                                                             <  2:307172/76‑1 (MQ=255)
                                                                   tGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAggg                                                                                                                             >  1:307172/1‑76 (MQ=255)
                                                                                      gttcgttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCt                                                                              <  2:355444/104‑1 (MQ=255)
                                                                                      gttcgttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCt                                                                              >  1:355444/1‑104 (MQ=255)
                                                                                      gttcgttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCagag                                               <  2:725038/135‑1 (MQ=255)
                                                                                      gttcgttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAAc                                           >  2:96387/1‑139 (MQ=255)
                                                                                       ttcgttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAAcc                                          >  2:754991/1‑139 (MQ=255)
                                                                                                    cGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACgcgg                             >  1:753471/1‑139 (MQ=255)
                                                                                                                   tGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGaa                <  1:463277/137‑1 (MQ=255)
                                                                                                                    ggCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGc             >  1:420748/1‑139 (MQ=255)
                                                                                                                    ggCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGc             >  1:800942/1‑139 (MQ=255)
                                                                                                                     gCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAAt               <  2:440645/136‑1 (MQ=255)
                                                                                                                     gCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCt            <  2:800942/139‑1 (MQ=255)
                                                                                                                           tCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCAt       >  1:451087/1‑138 (MQ=255)
                                                                                                                            ccGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATaa     <  2:362503/139‑1 (MQ=255)
                                                                                                                             cGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAAt    >  1:515837/1‑139 (MQ=255)
                                                                                                                             cGGTTTGCCCGGAGAGGGCCAGGACTAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAAt    <  2:127351/138‑1 (MQ=255)
                                                                                                                             cGGTCTGCTCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTCTCATAAt    <  2:783724/139‑1 (MQ=255)
                                                                                                                               gTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATgg  >  1:505820/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                   
CGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCAGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGG  >  NZ_CP009273/640199‑640464

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGGGCAGCAACCGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCAGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATAC  >  NZ_CP009273/640188‑640483
                                                                                                                                                 |                                                                                                                                                      
gtataagagacaGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGCTTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGA                                                                                                                                                   <  SRR3722130.758493/139‑1 (MQ=60)
  GGCAGCAACCGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGG                                                                                                                                                 >  SRR3722130.563378/1‑151 (MQ=60)
      tttttttttttcaagcagaagacggcatacgatattcgccttagtctcgtgggctcggagatgtgtataagagaCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCC                                                                                                                                              <  SRR3722130.172177/76‑1 (MQ=60)
          CCGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGAC                                                                                                                                         <  SRR3722130.27194/151‑1 (MQ=60)
          CCGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGAC                                                                                                                                         <  SRR3722130.662992/151‑1 (MQ=60)
           gggctcggagatgtgtataagagacagATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGAC                                                                                                                                         <  SRR3722130.429344/123‑1 (MQ=60)
            cggagatgtgtataagagacagGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGAC                                                                                                                                         <  SRR3722130.415926/127‑1 (MQ=60)
               gtataagagacaGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGC                                                                                                                                    <  SRR3722130.454515/139‑1 (MQ=60)
                 AGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACTTAGCG                                                                                                                                   >  SRR3722130.442953/1‑150 (MQ=60)
                    CAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCA                                                                                                                               <  SRR3722130.113849/151‑1 (MQ=60)
                        CCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAAT                                                                                                                           >  SRR3722130.250916/1‑151 (MQ=60)
                            GCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCC                                                                                                                       >  SRR3722130.369509/1‑151 (MQ=60)
                                  GACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTG                                                                                                                 <  SRR3722130.31606/151‑1 (MQ=60)
                                       TGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACTAGCGATCAAAATCGCCATACTGTTTTGC                                                                                                           >  SRR3722130.129187/1‑151 (MQ=60)
                                                                   TCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgc                                                                                >  SRR3722130.311924/1‑99 (MQ=60)
                                                                                      ACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCctgtctcttatacacatctccgagc                                                             >  SRR3722130.360970/1‑126 (MQ=60)
                                                                                                    CGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGC                                               >  SRR3722130.764503/1‑151 (MQ=60)
                                                                                                                    GTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCT                               >  SRR3722130.427399/1‑151 (MQ=60)
                                                                                                                    GTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCT                               >  SRR3722130.812645/1‑151 (MQ=60)
                                                                                                                           CTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATA                         >  SRR3722130.458145/1‑150 (MQ=60)
                                                                                                                             CTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATG                      >  SRR3722130.523749/1‑151 (MQ=60)
                                                                                                                             CTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAA                        <  SRR3722130.470497/149‑1 (MQ=60)
                                                                                                                               GGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGG                    >  SRR3722130.513611/1‑151 (MQ=60)
                                                                                                                                              GCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATAC      >  SRR3722130.333490/1‑150 (MQ=60)
                                                                                                                                                 CGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATAC  <  SRR3722130.460834/151‑1 (MQ=60)
                                                                                                                                                 |                                                                                                                                                      
CGGGCAGCAACCGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCAGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATAC  >  NZ_CP009273/640188‑640483

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