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 3,273,721 A→G E49G (GAA→GGA)  agaV → PTS N‑acetylgalactosamine transporter subunit IIB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,273,7210AG98.2% 172.6 / ‑4.1 54E49G (GAA→GGA) agaVPTS N‑acetylgalactosamine transporter subunit IIB
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (22/31);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

TTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAG  >  NZ_CP009273/3273634‑3273806
                                                                                       |                                                                                     
ttCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGgaa                                                                                     <  2:1529671/90‑1 (MQ=255)
ttCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGgaa                                                                                     <  2:1556973/90‑1 (MQ=255)
 tCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAAt                                                                                    <  1:101666/90‑1 (MQ=255)
 tCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAAt                                                                                    <  2:600810/90‑1 (MQ=255)
  cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                   >  1:1164595/1‑90 (MQ=255)
  cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                   <  1:997412/90‑1 (MQ=255)
  cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                   <  1:713011/90‑1 (MQ=255)
  cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                   <  1:471854/90‑1 (MQ=255)
  cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                   >  2:791014/1‑90 (MQ=255)
  cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                   <  1:1466363/90‑1 (MQ=255)
      gggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTAc                                                                               <  2:1152989/90‑1 (MQ=255)
      gggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTAc                                                                               <  1:210938/90‑1 (MQ=255)
      gggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTAATGGGAATGGTAc                                                                               <  1:1494419/90‑1 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              >  2:131683/1‑90 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              <  2:218281/90‑1 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              >  2:676112/1‑90 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              >  1:355165/1‑90 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              >  1:409716/1‑90 (MQ=255)
            gATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAg                                                                         <  2:601635/90‑1 (MQ=255)
                 gCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAggg                                                                    >  2:738555/1‑90 (MQ=255)
                    ggggCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATc                                                                 <  1:243973/90‑1 (MQ=255)
                    ggggCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATc                                                                 <  1:1145616/90‑1 (MQ=255)
                     gggCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCg                                                                >  1:195907/1‑90 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  1:1360068/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  2:808210/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  1:330174/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  1:1229268/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  2:355165/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  2:334641/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  1:1507169/90‑1 (MQ=255)
                           aTCTGGGGCTGGTAGCCAACGATGAGGTTGCCGAAAGTCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  2:754954/90‑1 (MQ=255)
                            tCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACg                                                         <  2:809/90‑1 (MQ=255)
                                  gCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTg                                                   >  2:583750/1‑90 (MQ=255)
                                    tGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              <  2:1029251/61‑1 (MQ=255)
                                    tGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              >  1:1029251/1‑61 (MQ=255)
                                     ggTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAc                                                >  2:961996/1‑90 (MQ=255)
                                     ggTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAc                                                >  2:926468/1‑90 (MQ=255)
                                     ggTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAc                                                >  1:947801/1‑90 (MQ=255)
                                      gTAGCCAACGATGAGGTTGTCGAAGATCCGGTACAACAAAACCTGATGGGAAAGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                               <  1:528335/90‑1 (MQ=255)
                                      gTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                               <  2:197689/90‑1 (MQ=255)
                                      gTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                               <  2:245303/90‑1 (MQ=255)
                                      gTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                               <  1:647974/90‑1 (MQ=255)
                                       tAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGc                                              <  1:375732/90‑1 (MQ=255)
                                         gCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTg                                            >  1:1226235/1‑90 (MQ=255)
                                                       tGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGa                              >  1:698318/1‑90 (MQ=255)
                                                         ccGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGaca                            <  2:332063/90‑1 (MQ=255)
                                                                  cGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTTATc                   >  1:76965/1‑90 (MQ=255)
                                                                        aacaaAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGcgccg             >  1:387651/1‑90 (MQ=255)
                                                                         acaaAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGcgccgc            >  2:1563049/1‑90 (MQ=255)
                                                                              aCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATc       >  1:1552970/1‑90 (MQ=255)
                                                                              aCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATc       <  2:1244452/90‑1 (MQ=255)
                                                                              aCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATc       >  1:1013607/1‑90 (MQ=255)
                                                                              aCCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATc       >  2:10671/1‑90 (MQ=255)
                                                                                cTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGa     <  1:738555/90‑1 (MQ=255)
                                                                                 tGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGAc    >  2:1146515/1‑90 (MQ=255)
                                                                                   aTGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAg  >  1:780117/1‑90 (MQ=255)
                                                                                       |                                                                                     
TTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAG  >  NZ_CP009273/3273634‑3273806

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGGTCGGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGT  >  NZ_CP009273/3273625‑3273820
                                                                                                |                                                                                                   
AGGTCGGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAAT                                                                                                  <  SRR3722078.1256269/100‑1 (MQ=60)
 GGTCGGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATG                                                                                                 >  SRR3722078.1182676/1‑100 (MQ=60)
     cgtcggcagcgtcagatgtgtataagagacagTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTAC                                                                                             <  SRR3722078.447445/68‑1 (MQ=60)
     GGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTAC                                                                                             <  SRR3722078.707984/100‑1 (MQ=60)
      GCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACT                                                                                            >  SRR3722078.359492/1‑100 (MQ=60)
      GCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACT                                                                                            >  SRR3722078.414670/1‑100 (MQ=60)
          TCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCA                                                                                        <  SRR3722078.102871/100‑1 (MQ=60)
           CAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAG                                                                                       <  SRR3722078.1011942/100‑1 (MQ=60)
           CAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAG                                                                                       <  SRR3722078.1490288/100‑1 (MQ=60)
           CAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAG                                                                                       <  SRR3722078.477453/100‑1 (MQ=60)
           CAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAG                                                                                       <  SRR3722078.721754/100‑1 (MQ=60)
               GGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGG                                                                                   <  SRR3722078.213491/100‑1 (MQ=60)
               GGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTAATGGGAATGGTACTGGCAGAAGG                                                                                   <  SRR3722078.1518784/100‑1 (MQ=60)
                    GGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCG                                                                              >  SRR3722078.198242/1‑100 (MQ=60)
                             GGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTT                                                                     <  SRR3722078.1163257/100‑1 (MQ=60)
                             GGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTT                                                                     <  SRR3722078.246961/100‑1 (MQ=60)
                                   AATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGctgtctcttatacacatct                                                               >  SRR3722078.1044529/1‑81 (MQ=60)
                                    ATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAC                                                              >  SRR3722078.961077/1‑100 (MQ=60)
                                    ATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAC                                                              <  SRR3722078.334213/100‑1 (MQ=60)
                                    ATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAC                                                              <  SRR3722078.1531743/100‑1 (MQ=60)
                                    ATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAC                                                              <  SRR3722078.1381921/100‑1 (MQ=60)
                                    ATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAC                                                              <  SRR3722078.1248633/100‑1 (MQ=60)
                                        GGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTG                                                          >  SRR3722078.1245549/1‑100 (MQ=60)
                                               GTAGCCAACGATGAGGTTGTCGAAGATCCGGTACAACAAAACCTGATGGGAAAGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAG                                                   <  SRR3722078.534578/100‑1 (MQ=60)
                                               GTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAG                                                   <  SRR3722078.655419/100‑1 (MQ=60)
                                                TAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGT                                                  <  SRR3722078.380322/100‑1 (MQ=60)
                                                      ACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGA                                            >  SRR3722078.706704/1‑100 (MQ=60)
                                                                 GCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTTATC                                 >  SRR3722078.77873/1‑100 (MQ=60)
                                                                       GATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCG                           >  SRR3722078.392374/1‑100 (MQ=60)
                                                                             GTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATC                     >  SRR3722078.1028541/1‑100 (MQ=60)
                                                                             GTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATC                     >  SRR3722078.1578256/1‑100 (MQ=60)
                                                                                  ACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAG                >  SRR3722078.790246/1‑100 (MQ=60)
                                                                                        CCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATC          >  SRR3722078.746341/1‑100 (MQ=60)
                                                                                         CTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCC         <  SRR3722078.747878/100‑1 (MQ=60)
                                                                                                GAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGT  >  SRR3722078.1035760/1‑100 (MQ=60)
                                                                                                |                                                                                                   
AGGTCGGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGT  >  NZ_CP009273/3273625‑3273820

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