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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,535,339 T→C F525S (TTC→TCC)  feoB → Fe(2+) transporter permease subunit FeoB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,535,3390TC100.0% 148.2 / NA 43F525S (TTC→TCC) feoBFe(2+) transporter permease subunit FeoB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/17);  total (26/17)

CTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGT  >  NZ_CP009273/3535263‑3535423
                                                                            |                                                                                    
cTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTc                                                                         >  1:1008835/1‑90 (MQ=255)
cTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTc                                                                         >  1:148708/1‑90 (MQ=255)
cTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTc                                                                         >  1:841235/1‑90 (MQ=255)
  gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCaa                                                                       >  2:319896/1‑90 (MQ=255)
  gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCaa                                                                       >  2:309754/1‑90 (MQ=255)
  gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCaa                                                                       >  2:1228060/1‑90 (MQ=255)
  gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGGTGTGCTTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCaa                                                                       >  2:868060/1‑90 (MQ=255)
   aTTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAAc                                                                      <  2:1108368/90‑1 (MQ=255)
           gACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTtctc                                                              >  2:1036619/1‑90 (MQ=255)
           gACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTtctc                                                              >  1:584328/1‑90 (MQ=255)
             ccTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGc                                                            <  1:1102361/90‑1 (MQ=255)
                ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                         >  2:362093/1‑90 (MQ=255)
                ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                         >  2:1049618/1‑90 (MQ=255)
                ggCAGCGTCGGAAAGGCTTCGTTCTGCGTGCTGGTAAACTGATGATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCtcg                                                         >  2:1102856/1‑88 (MQ=255)
                     cGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCggg                                                    >  1:1176519/1‑90 (MQ=255)
                      gTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGa                                                   >  2:536567/1‑90 (MQ=255)
                                cTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGa                                         <  2:808453/90‑1 (MQ=255)
                                 ttCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGAt                                        <  1:480965/90‑1 (MQ=255)
                                  tCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATa                                       <  1:416018/90‑1 (MQ=255)
                                    gTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAAc                                     >  1:264537/1‑90 (MQ=255)
                                    gTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAAc                                     >  1:654978/1‑90 (MQ=255)
                                         gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCaa                                >  1:346847/1‑90 (MQ=255)
                                              cTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACt                           >  2:1114862/1‑90 (MQ=255)
                                               tGGTAAGGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTc                          <  1:13192/90‑1 (MQ=255)
                                               tGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTc                          <  1:1114862/90‑1 (MQ=255)
                                               tGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTc                          <  1:799107/90‑1 (MQ=255)
                                               tGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTc                          <  1:75341/90‑1 (MQ=255)
                                               tGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTc                          <  1:527400/90‑1 (MQ=255)
                                                    aaGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGgcgc                     <  2:1176519/90‑1 (MQ=255)
                                                        gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAAc                               <  1:246941/76‑1 (MQ=255)
                                                        gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAAc                               >  2:246941/1‑76 (MQ=255)
                                                        gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc                 >  1:893058/1‑90 (MQ=255)
                                                          tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGGGCTGGcgt               >  2:748085/1‑90 (MQ=255)
                                                          tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt               >  1:1210354/1‑90 (MQ=255)
                                                           catcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtc              >  1:1232097/1‑90 (MQ=255)
                                                               atcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgtc          <  2:855142/90‑1 (MQ=255)
                                                                tcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCa         >  2:509124/1‑90 (MQ=255)
                                                                  gTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGc       >  2:294559/1‑90 (MQ=255)
                                                                   tCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGcc      <  1:567400/90‑1 (MQ=255)
                                                                   tCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGcc      <  1:300452/90‑1 (MQ=255)
                                                                     aGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCgg    <  1:1171312/90‑1 (MQ=255)
                                                                      gCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCggg   <  1:294559/90‑1 (MQ=255)
                                                                       cATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGt  >  2:361624/1‑90 (MQ=255)
                                                                            |                                                                                    
CTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGT  >  NZ_CP009273/3535263‑3535423

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ 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
ATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCT  >  NZ_CP009273/3535240‑3535437
                                                                                                   |                                                                                                  
ATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTC                                                                                                    >  SRR3722089.981065/1‑100 (MQ=60)
ATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTC                                                                                                    >  SRR3722089.1185541/1‑100 (MQ=60)
  cactaagatggtcgtcggcagcgtcagatgtgtataagagacaGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCC                                                                                                  <  SRR3722089.823464/57‑1 (MQ=60)
  CATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCC                                                                                                  <  SRR3722089.383344/100‑1 (MQ=60)
  CATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCC                                                                                                  <  SRR3722089.643795/100‑1 (MQ=60)
  CATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCC                                                                                                  <  SRR3722089.302126/100‑1 (MQ=60)
  CATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCC                                                                                                  <  SRR3722089.287435/100‑1 (MQ=60)
             CGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTC                                                                                       >  SRR3722089.1023860/1‑100 (MQ=60)
             CGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTC                                                                                       >  SRR3722089.150553/1‑100 (MQ=60)
             CGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTC                                                                                       >  SRR3722089.853161/1‑100 (MQ=60)
                        TGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTC                                                                            >  SRR3722089.591349/1‑100 (MQ=60)
                                  GACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGG                                                                  >  SRR3722089.1194449/1‑100 (MQ=60)
                                    CCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAA                                                                <  SRR3722089.1119026/100‑1 (MQ=60)
                                                 GAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAAC                                                   >  SRR3722089.663403/1‑100 (MQ=60)
                                                 GAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAAC                                                   >  SRR3722089.267923/1‑100 (MQ=60)
                                                      GCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAA                                              >  SRR3722089.351160/1‑100 (MQ=60)
                                                        TTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACG                                            <  SRR3722089.486677/100‑1 (MQ=60)
                                                         TCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGA                                           <  SRR3722089.421050/100‑1 (MQ=60)
                                                                 acaGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGC                                   <  SRR3722089.250094/97‑1 (MQ=60)
                                                                     CTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGC                               >  SRR3722089.905974/1‑100 (MQ=60)
                                                                      gatgtgtataagagacagCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCG                              <  SRR3722089.1246730/82‑1 (MQ=60)
                                                                      TGGTAAGGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCG                              <  SRR3722089.13352/100‑1 (MQ=60)
                                                                      TGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCG                              <  SRR3722089.1131745/100‑1 (MQ=60)
                                                                      TGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCG                              <  SRR3722089.533561/100‑1 (MQ=60)
                                                                      TGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCG                              <  SRR3722089.810149/100‑1 (MQ=60)
                                                                      TGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCG                              <  SRR3722089.76256/100‑1 (MQ=60)
                                                                       GGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGT                             >  SRR3722089.1228756/1‑100 (MQ=60)
                                                                        GTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTC                            >  SRR3722089.1250828/1‑100 (MQ=60)
                                                                                        CGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATC            >  SRR3722089.916664/1‑100 (MQ=60)
                                                                                          TCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCAC          <  SRR3722089.574085/100‑1 (MQ=60)
                                                                                          TCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCAC          <  SRR3722089.304271/100‑1 (MQ=60)
                                                                                            AGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCC        <  SRR3722089.1189137/100‑1 (MQ=60)
                                                                                             GCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCC       <  SRR3722089.298317/100‑1 (MQ=60)
                                                                                             GCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCC       >  SRR3722089.179029/1‑100 (MQ=60)
                                                                                                  TCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCT  <  SRR3722089.621208/100‑1 (MQ=60)
                                                                                                  TCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCT  >  SRR3722089.606445/1‑100 (MQ=60)
                                                                                                  TCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCT  >  SRR3722089.515351/1‑100 (MQ=60)
                                                                                                  TCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCT  >  SRR3722089.842059/1‑100 (MQ=60)
                                                                                                   |                                                                                                  
ATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCT  >  NZ_CP009273/3535240‑3535437

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