Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F2 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,593,447 T→C S52P (TCC→CCC)  yjjQ → response regulator transcription factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,593,4470TC98.4% 202.2 / ‑3.9 63S52P (TCC→CCC) yjjQresponse regulator transcription factor
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (43/19);  total (44/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

AGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCAC  >  NZ_CP009273/4593364‑4593528
                                                                                   |                                                                                 
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  1:1079057/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  1:643497/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  1:655733/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  1:666355/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  1:275076/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  1:810233/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  2:1213870/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  2:442260/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  2:560889/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  1:1545950/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  2:564399/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  2:765651/1‑90 (MQ=255)
 ggTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAtt                                                                            <  1:1159442/90‑1 (MQ=255)
 ggTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAtt                                                                            >  2:1244757/1‑90 (MQ=255)
     aTGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTc                                                                        <  1:1197595/90‑1 (MQ=255)
                   cccTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGc                                                          >  2:426764/1‑90 (MQ=255)
                      tGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCgg                                                       >  2:293110/1‑90 (MQ=255)
                        aaTATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTg                                                     <  2:671248/90‑1 (MQ=255)
                        aaTATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTg                                                     <  2:206693/90‑1 (MQ=255)
                        aaTATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTg                                                     >  2:1348318/1‑90 (MQ=255)
                          tatGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGaa                                                   >  2:854467/1‑90 (MQ=255)
                            tGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAg                                                 <  1:564399/90‑1 (MQ=255)
                                   atatCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGAt                                          >  2:1162248/1‑90 (MQ=255)
                                       ccAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCAc                                      >  1:223923/1‑90 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                   >  2:1134598/1‑90 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAg                                                 <  2:259551/76‑1 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAg                                                 >  1:259551/1‑76 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                   >  2:1347688/1‑90 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                   >  1:356327/1‑90 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                   >  1:672134/1‑90 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                   >  2:1216736/1‑90 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                   >  2:295797/1‑90 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                   >  1:1440049/1‑90 (MQ=255)
                                            gcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCt                                 <  1:560294/90‑1 (MQ=255)
                                             ccTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTg                                >  1:1082880/1‑90 (MQ=255)
                                              ctctGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGt                               >  1:1562800/1‑90 (MQ=255)
                                              ctctGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGt                               >  1:497071/1‑90 (MQ=255)
                                               tctGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGtt                              >  1:372819/1‑90 (MQ=255)
                                                 tGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTtg                            <  2:864397/90‑1 (MQ=255)
                                                    ggaggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGa                         <  2:655733/90‑1 (MQ=255)
                                                    ggaggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGa                         <  2:1529287/90‑1 (MQ=255)
                                                    ggaggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGa                         >  2:212213/1‑90 (MQ=255)
                                                    ggaggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGa                         <  1:732306/90‑1 (MQ=255)
                                                      atgaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAc                       <  1:1577285/88‑1 (MQ=255)
                                                      aggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAc                       <  2:643497/90‑1 (MQ=255)
                                                      aggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAc                       <  2:1562800/90‑1 (MQ=255)
                                                      aggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAc                       <  2:497071/90‑1 (MQ=255)
                                                      aggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAc                       <  1:182146/90‑1 (MQ=255)
                                                       ggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACa                      >  2:384113/1‑90 (MQ=255)
                                                         accttaccttATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACatt                    >  1:1517144/1‑90 (MQ=255)
                                                          ccttaccttaTTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACatta                   >  1:1124596/1‑90 (MQ=255)
                                                            ttaccttaTTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACattatt                 >  1:623655/1‑90 (MQ=255)
                                                                   attaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACattattat               >  1:218444/1‑85 (MQ=255)
                                                                   attaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACattattat               <  2:218444/85‑1 (MQ=255)
                                                                   attaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCttt          >  2:1133761/1‑90 (MQ=255)
                                                                   attaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCttt          >  2:678741/1‑90 (MQ=255)
                                                                    ttaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTa         >  2:1390108/1‑90 (MQ=255)
                                                                    ttaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTa         >  1:649961/1‑90 (MQ=255)
                                                                      aCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAt       <  1:1267966/90‑1 (MQ=255)
                                                                         aTTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAtctc    >  1:482083/1‑90 (MQ=255)
                                                                         aTTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAtctc    >  1:52885/1‑90 (MQ=255)
                                                                         aTTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAtctc    >  1:1198884/1‑90 (MQ=255)
                                                                           tACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTcac  <  1:924496/90‑1 (MQ=255)
                                                                                   |                                                                                 
AGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCAC  >  NZ_CP009273/4593364‑4593528

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ 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
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCA  >  NZ_CP009273/4593354‑4593544
                                                                                             |                                                                                                 
GGGTTAAAGGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAT                                                                                             >  SRR3722078.1571134/1‑100 (MQ=60)
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAT                                                                                             >  SRR3722078.278430/1‑100 (MQ=60)
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAT                                                                                             >  SRR3722078.820949/1‑100 (MQ=60)
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAT                                                                                             >  SRR3722078.1095382/1‑100 (MQ=60)
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAT                                                                                             >  SRR3722078.650889/1‑100 (MQ=60)
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAT                                                                                             >  SRR3722078.673984/1‑100 (MQ=60)
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAT                                                                                             >  SRR3722078.663243/1‑100 (MQ=60)
   agagacagGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAG                                                                                          <  SRR3722078.991345/92‑1 (MQ=60)
     AAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTC                                                                                        <  SRR3722078.571596/100‑1 (MQ=60)
     AAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTC                                                                                        <  SRR3722078.21380/100‑1 (MQ=60)
     AAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTC                                                                                        <  SRR3722078.324108/100‑1 (MQ=60)
           GGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCT                                                                                  <  SRR3722078.1177397/100‑1 (MQ=60)
               ATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATT                                                                              <  SRR3722078.1216380/100‑1 (MQ=60)
                                      TGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCA                                                       <  SRR3722078.571018/100‑1 (MQ=60)
                                       GAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCAC                                                      >  SRR3722078.226645/1‑100 (MQ=60)
                                          ATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACtgt                                                   >  SRR3722078.262721/1‑97 (MQ=60)
                                          ATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTT                                                   >  SRR3722078.679823/1‑100 (MQ=60)
                                          ATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTT                                                   >  SRR3722078.360673/1‑100 (MQ=60)
                                          ATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTT                                                   >  SRR3722078.1463488/1‑100 (MQ=60)
                                             ATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTG                                                >  SRR3722078.1099263/1‑100 (MQ=60)
                                              TATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGT                                               >  SRR3722078.1588240/1‑100 (MQ=60)
                                              TATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGT                                               >  SRR3722078.502954/1‑100 (MQ=60)
                                               ATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTT                                              >  SRR3722078.377379/1‑100 (MQ=60)
                                                      GCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAC                                       <  SRR3722078.566884/100‑1 (MQ=60)
                                                         TCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATT                                    >  SRR3722078.1541887/1‑100 (MQ=60)
                                                          CTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTA                                   >  SRR3722078.1141793/1‑100 (MQ=60)
                                                            GCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATT                                 >  SRR3722078.630839/1‑100 (MQ=60)
                                                              GGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTAT                               <  SRR3722078.741477/100‑1 (MQ=60)
                                                                ATGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTC                             <  SRR3722078.1602962/100‑1 (MQ=60)
                                                                AGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTC                             <  SRR3722078.184307/100‑1 (MQ=60)
                                                                   ACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTT                          >  SRR3722078.221085/1‑100 (MQ=60)
                                                                    CCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTA                         >  SRR3722078.657422/1‑100 (MQ=60)
                                                                         CCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTC                    >  SRR3722078.53507/1‑100 (MQ=60)
                                                                         CCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTC                    >  SRR3722078.487796/1‑100 (MQ=60)
                                                                         CCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTC                    >  SRR3722078.1217697/1‑100 (MQ=60)
                                                                                ACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATAT             <  SRR3722078.1288081/100‑1 (MQ=60)
                                                                                 CAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATC            >  SRR3722078.247333/1‑100 (MQ=60)
                                                                                 CAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATC            >  SRR3722078.751482/1‑100 (MQ=60)
                                                                                     TACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGA        <  SRR3722078.937165/100‑1 (MQ=60)
                                                                                           GTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCA  >  SRR3722078.818231/1‑100 (MQ=60)
                                                                                             |                                                                                                 
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCA  >  NZ_CP009273/4593354‑4593544

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: