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 I184 R1 278 29.6 1166432 97.3% 1134938 124.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,026,705 A→C V57G (GTG→GGG)  yccA ← FtsH protease modulator YccA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,026,7050AC95.1% 136.2 / 1.0 41coding (170/660 nt)yccAFtsH protease modulator YccA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (20/19);  minor base . (1/1);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

TCCGGCAGACAGATAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAGTGAT  >  NZ_CP009273/1026572‑1026839
                                                                                                                                     |                                                                                                                                      
tCCGGCAGACAGATAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGc                                                                                                                                   >  1:394314/1‑139 (MQ=255)
                   ttCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCg                                                                                                                   >  1:9179/1‑136 (MQ=255)
                   ttCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCg                                                                                                                   >  1:558605/1‑136 (MQ=255)
                        aaTAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGgc                                                                                                           <  2:209267/139‑1 (MQ=255)
                        aaTAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGgc                                                                                                           <  2:31081/139‑1 (MQ=255)
                             gTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCAt                                                                                                      >  1:133624/1‑139 (MQ=255)
                                  aGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCa                                                                                                 >  1:321888/1‑139 (MQ=255)
                                  aGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCa                                                                                                 <  2:283666/139‑1 (MQ=255)
                                   ggATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCAc                                                                                                <  2:133624/139‑1 (MQ=255)
                                             cAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCAt                                                                                                      >  2:61076/1‑123 (MQ=255)
                                             cAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCAt                                                                                                      <  1:61076/123‑1 (MQ=255)
                                               gAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCa                                                                                                 <  2:62723/126‑1 (MQ=255)
                                               gAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCa                                                                                                 >  1:62723/1‑126 (MQ=255)
                                                                  tGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAg                                                                 >  2:454798/1‑139 (MQ=255)
                                                                  tGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAg                                                                 >  1:180814/1‑139 (MQ=255)
                                                                       agaTAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGc                                                               >  1:156589/1‑136 (MQ=255)
                                                                       agaTAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGc                                                               <  2:156589/136‑1 (MQ=255)
                                                                                 ggTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCt                                                   >  1:200977/1‑138 (MQ=255)
                                                                                   tCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCa                                                             >  1:139496/1‑126 (MQ=255)
                                                                                   tCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCAGCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCa                                                             <  2:139496/126‑1 (MQ=255)
                                                                                   tCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCAt                                                                                                      >  1:32681/1‑84 (MQ=255)
                                                                                   tCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCAt                                                                                                      <  2:32681/84‑1 (MQ=255)
                                                                                           ttCGCCGTTTTATAGGTCAGGGACATCAAACCATACATACCCCCCAGCGGCAGAATCAGACGCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAAt                                        <  2:94767/139‑1 (MQ=255)
                                                                                           ttCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAAt                                        <  2:383444/139‑1 (MQ=255)
                                                                                           ttCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAAt                                        <  2:321888/139‑1 (MQ=255)
                                                                                                cGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTa                                   <  2:200977/139‑1 (MQ=255)
                                                                                                 gTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCa                                                                                                 >  1:453027/1‑76 (MQ=255)
                                                                                                 gTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCa                                                                                                 <  2:453027/76‑1 (MQ=255)
                                                                                                     tataGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACg                              >  1:140931/1‑139 (MQ=255)
                                                                                                          gTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCgg                                                                                                                  >  1:259432/1‑50 (MQ=255)
                                                                                                          gTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCgg                                                                                                                  <  2:259432/50‑1 (MQ=255)
                                                                                                            cAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGGTCAGCAGAACATAGGTATTACGCAGCAcc                       >  2:100737/1‑139 (MQ=255)
                                                                                                             aGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCt                      >  2:180759/1‑139 (MQ=255)
                                                                                                                aaCATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTAt                   <  1:560280/139‑1 (MQ=255)
                                                                                                                         ccATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACTTATGGGTGCTAAg         <  2:266207/139‑1 (MQ=255)
                                                                                                                          catacataCCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAg         <  2:78646/139‑1 (MQ=255)
                                                                                                                            tacatacCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAg         >  1:391318/1‑137 (MQ=255)
                                                                                                                             acatacCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTa                                   <  2:310998/110‑1 (MQ=255)
                                                                                                                             acatacCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTa                                   >  1:310998/1‑110 (MQ=255)
                                                                                                                              catacCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAg      <  1:63942/138‑1 (MQ=255)
                                                                                                                                 acCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAGTGAt  >  1:244703/1‑139 (MQ=255)
                                                                                                                                  ccccccAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAGTGAt  <  2:391318/138‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
TCCGGCAGACAGATAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAGTGAT  >  NZ_CP009273/1026572‑1026839

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ 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
ACGTCACCCATTCCGGCAGACAGATAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAGTGATGTACGGTC  >  NZ_CP009273/1026561‑1026847
                                                                                                                                                |                                                                                                                                              
ACGTCACCCATTCCGGCAGACAGATAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCG                                                                                                                                          >  SRR3722139.398998/1‑151 (MQ=60)
       gggctcggagatgtgtataagagacagGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAA                                                                                                                                    <  SRR3722139.543991/123‑1 (MQ=60)
                   ACAGATAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGG                                                                                                                          >  SRR3722139.9269/1‑148 (MQ=60)
                   ACAGATAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGG                                                                                                                          >  SRR3722139.564999/1‑148 (MQ=60)
                             GTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATC                                                                                                             >  SRR3722139.135250/1‑151 (MQ=60)
                                  GAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCAC                                                                                                        >  SRR3722139.325857/1‑151 (MQ=60)
                                       agacaGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTAC                                                                                                   <  SRR3722139.61774/146‑1 (MQ=60)
                                               GATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCActg                                                                                           >  SRR3722139.63447/1‑148 (MQ=60)
                                                                  GTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGG                                                                        >  SRR3722139.183055/1‑151 (MQ=60)
                                                                       GGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGC                                                                   >  SRR3722139.158503/1‑151 (MQ=60)
                                                                                 GAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTC                                                          >  SRR3722139.203454/1‑150 (MQ=60)
                                                                                   GATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCctg                                                       >  SRR3722139.141179/1‑148 (MQ=60)
                                                                                   GATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACctgtctcttatacacatctccgagcccacgagaccgtact                                                           >  SRR3722139.33042/1‑106 (MQ=60)
                                                                                                 GCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCActgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgc                                         >  SRR3722139.458325/1‑98 (MQ=60)
                                                                                                     ATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGC                                     >  SRR3722139.142630/1‑151 (MQ=60)
                                                                                                          CCGTTTTATAGGTCAGGAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCActgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtcttctgcttgaaaaaaaaaaac                                 >  SRR3722139.262698/1‑72 (MQ=60)
                                                                                                                          GAACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGC                <  SRR3722139.566689/151‑1 (MQ=60)
                                                                                                                           ataagagacagATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGC                <  SRR3722139.559789/139‑1 (MQ=60)
                                                                                                                            ACATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGC                >  SRR3722139.395985/1‑149 (MQ=60)
                                                                                                                             ataagagacagACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAG              <  SRR3722139.131630/139‑1 (MQ=60)
                                                                                                                             CATCAAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTgtctcttatacacatct             >  SRR3722139.314851/1‑134 (MQ=60)
                                                                                                                                 AAACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCCGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAGTGATG         >  SRR3722139.247750/1‑151 (MQ=60)
                                                                                                                                   ACCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAGTGATGct       >  SRR3722139.329602/1‑149 (MQ=60)
                                                                                                                                    CCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAGTGATGTAC      >  SRR3722139.477386/1‑151 (MQ=60)
                                                                                                                                    CCATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGGATTACGCAGCACCTTATGGGTGCTAAGCAGTGATGTAC      >  SRR3722139.562973/1‑151 (MQ=60)
                                                                                                                                      ATACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCctgtctcttatacacatctccgagccccgagaccgtactagatctcgtatgccgtc    >  SRR3722139.383709/1‑95 (MQ=60)
                                                                                                                                        ACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAGTGATGTACGGTC  >  SRR3722139.125185/1‑151 (MQ=60)
                                                                                                                                        ACATACCCCCCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAGTGATGTACGGT   <  SRR3722139.64712/150‑1 (MQ=60)
                                                                                                                                                |                                                                                                                                              
ACGTCACCCATTCCGGCAGACAGATAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAGCAGTGATGTACGGTC  >  NZ_CP009273/1026561‑1026847

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