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,728,024 A→C S315S (TCA→TCC xylH → xylose ABC transporter permease XylH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,728,0240AC95.8% 226.9 / 0.9 71S315S (TCA→TCCxylHxylose ABC transporter permease XylH
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (30/38);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCG  >  NZ_CP009273/3727940‑3728109
                                                                                    |                                                                                     
cACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGCTTAATCCTTAGTACTCCACTTGGCGCTGGTTCCCCTTc                                                                                  <  1:1191787/90‑1 (MQ=255)
cACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTc                                                                                  >  1:509238/1‑90 (MQ=255)
cACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTc                                                                                  >  1:165274/1‑90 (MQ=255)
         tGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGaaa                                                                         >  2:1237983/1‑90 (MQ=255)
          gCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAat                                                                        >  1:164287/1‑90 (MQ=255)
          gCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAat                                                                        <  1:1019678/90‑1 (MQ=255)
          gCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAat                                                                        >  2:1476128/1‑90 (MQ=255)
          gCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAat                                                                        >  2:1573186/1‑90 (MQ=255)
          gCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGc                                                                               >  2:1558503/1‑83 (MQ=255)
          gCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGc                                                                               <  1:1558503/83‑1 (MQ=255)
               gtTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGc                                                                   >  1:270358/1‑90 (MQ=255)
                    cgATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAAc                                                              <  1:128095/90‑1 (MQ=255)
                         aaCGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         <  1:364939/90‑1 (MQ=255)
                         aaCGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         <  2:973686/90‑1 (MQ=255)
                         aaCGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         <  1:881402/90‑1 (MQ=255)
                         aaCGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         <  2:1395737/90‑1 (MQ=255)
                         aaCGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         <  2:1437312/90‑1 (MQ=255)
                         aaCGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         <  1:599817/90‑1 (MQ=255)
                         aaCGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         <  2:1569570/90‑1 (MQ=255)
                         aaCGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         <  1:1288273/90‑1 (MQ=255)
                         aaCGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         >  2:250858/1‑90 (MQ=255)
                         aaCGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         <  2:285005/90‑1 (MQ=255)
                           cGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGc                                                       <  2:164179/90‑1 (MQ=255)
                           cGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGc                                                       <  2:697704/90‑1 (MQ=255)
                           cGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGc                                                       <  1:1517804/90‑1 (MQ=255)
                           cGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGc                                                       <  2:270358/90‑1 (MQ=255)
                              aTTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAAt                                                    >  1:1191667/1‑90 (MQ=255)
                              aTTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAAt                                                    >  1:282158/1‑90 (MQ=255)
                                tAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTg                                                  >  2:787450/1‑90 (MQ=255)
                                  aTGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTgca                                                <  2:47166/90‑1 (MQ=255)
                                      tAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGGAATATCGCCGAACTGGACGCAATTGCAGCAt                                            <  2:458648/90‑1 (MQ=255)
                                       aGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         <  1:1287207/76‑1 (MQ=255)
                                       aGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         >  2:1287207/1‑76 (MQ=255)
                                       aGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATg                                           >  1:711588/1‑90 (MQ=255)
                                           aTCGCCGGATTAATCCTTAGTTCTCGACTTGTCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTg                                       >  2:1193346/1‑90 (MQ=255)
                                            tCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGa                                      >  2:788203/1‑90 (MQ=255)
                                            tCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGa                                      >  1:1194396/1‑90 (MQ=255)
                                                cGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTgg                                  <  2:153254/90‑1 (MQ=255)
                                                 ggATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTggc                                 <  2:954858/90‑1 (MQ=255)
                                                 ggATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTggc                                 <  2:722900/90‑1 (MQ=255)
                                                 ggATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTggc                                 <  2:509726/90‑1 (MQ=255)
                                                 ggATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTggc                                 <  1:1321652/90‑1 (MQ=255)
                                                  gATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTggcg                                >  2:768200/1‑90 (MQ=255)
                                                      aaTCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCAc                            >  1:1136282/1‑90 (MQ=255)
                                                      aaTCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCAc                            <  1:1237983/90‑1 (MQ=255)
                                                      aaTCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCAc                            <  2:223451/90‑1 (MQ=255)
                                                      aaTCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCAc                            <  1:1573186/90‑1 (MQ=255)
                                                       aTCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCAcc                           <  2:509238/90‑1 (MQ=255)
                                                       aTCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCAcc                           <  2:188801/90‑1 (MQ=255)
                                                           ttAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCCCCAGcc                       <  1:880688/90‑1 (MQ=255)
                                                           ttAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGcc                       <  1:934894/90‑1 (MQ=255)
                                                           ttAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGcc                       <  1:634537/90‑1 (MQ=255)
                                                            tAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCct                      >  2:1532871/1‑90 (MQ=255)
                                                            tAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCct                      >  2:1299914/1‑90 (MQ=255)
                                                            tAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCct                      >  2:757744/1‑90 (MQ=255)
                                                                tctcGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCctggct                  <  1:1548244/90‑1 (MQ=255)
                                                                tctcGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCctggct                  >  1:672168/1‑90 (MQ=255)
                                                                tctcGACTTGGCGCTGGCTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCctggct                  >  2:802328/1‑90 (MQ=255)
                                                                 ctcGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAAtt                                                   <  2:149381/56‑1 (MQ=255)
                                                                 ctcGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAAtt                                                   >  1:149381/1‑56 (MQ=255)
                                                                 ctcGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCctggctg                 >  2:181153/1‑90 (MQ=255)
                                                                 ctcGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCctggctg                 >  1:613592/1‑90 (MQ=255)
                                                                   cGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCtggctggc               <  2:377751/90‑1 (MQ=255)
                                                                   cGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCtggctggc               >  1:1302581/1‑90 (MQ=255)
                                                                    gACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGTTTGGCGGCACCAGCCTGGCTGGCg              <  2:314959/90‑1 (MQ=255)
                                                                       ttGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGgtg           <  1:585848/90‑1 (MQ=255)
                                                                        tGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGgtgt          >  1:458670/1‑90 (MQ=255)
                                                                        tGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGgtgt          >  1:485196/1‑90 (MQ=255)
                                                                        tGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGgtgt          >  2:890591/1‑90 (MQ=255)
                                                                        tGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGgggt          >  1:1445019/1‑90 (MQ=255)
                                                                                gTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCg  <  1:431911/90‑1 (MQ=255)
                                                                                    |                                                                                     
CACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCG  >  NZ_CP009273/3727940‑3728109

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ 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
ACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGC  >  NZ_CP009273/3727930‑3728119
                                                                                              |                                                                                               
ACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTC                                                                                            >  SRR3722078.515265/1‑100 (MQ=60)
ACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTC                                                                                            <  SRR3722078.852919/100‑1 (MQ=60)
ACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTC                                                                                            <  SRR3722078.358162/100‑1 (MQ=60)
ACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTC                                                                                            >  SRR3722078.167200/1‑100 (MQ=60)
          CACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGCTTAATCCTTAGTACTCCACTTGGCGCTGGTTCCCCTTCTGCGGGAAAT                                                                                  <  SRR3722078.1210436/100‑1 (MQ=60)
          CACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAAT                                                                                  >  SRR3722078.166200/1‑100 (MQ=60)
             CAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATC                                                                               <  SRR3722078.1583872/100‑1 (MQ=60)
               AACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGC                                                                             >  SRR3722078.273660/1‑100 (MQ=60)
                    GCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAAC                                                                        <  SRR3722078.1034740/100‑1 (MQ=60)
                              CGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAAT                                                              >  SRR3722078.1210313/1‑100 (MQ=60)
                              CGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAAT                                                              <  SRR3722078.129573/100‑1 (MQ=60)
                              CGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAAT                                                              >  SRR3722078.285599/1‑100 (MQ=60)
                                   acaGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAG                                                         <  SRR3722078.1307698/97‑1 (MQ=60)
                                   AACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAG                                                         <  SRR3722078.606759/100‑1 (MQ=60)
                                   AACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAG                                                         <  SRR3722078.893251/100‑1 (MQ=60)
                                   AACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAG                                                         <  SRR3722078.369397/100‑1 (MQ=60)
                                   AACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAG                                                         <  SRR3722078.1308780/100‑1 (MQ=60)
                                     CGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCA                                                       <  SRR3722078.1542560/100‑1 (MQ=60)
                                       GATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATG                                                     >  SRR3722078.720291/1‑100 (MQ=60)
                                            ATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGA                                                >  SRR3722078.1213114/1‑100 (MQ=60)
                                                      TCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCAC                                      >  SRR3722078.1153707/1‑100 (MQ=60)
                                                           GGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCC                                 <  SRR3722078.1342773/100‑1 (MQ=60)
                                                            GATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGtctcttatacacatctgacgctgccgacgatggcggtcgtgtagatctcggtggtcgc                                >  SRR3722078.376044/1‑42 (MQ=60)
                                                                AATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCT                            <  SRR3722078.1257527/100‑1 (MQ=60)
                                                                AATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCT                            <  SRR3722078.1598789/100‑1 (MQ=60)
                                                                AATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCT                            >  SRR3722078.679857/1‑100 (MQ=60)
                                                                 ATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGctgtctcttatacacatctgacgc                           >  SRR3722078.151115/1‑76 (MQ=60)
                                                                 ATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTG                           >  SRR3722078.620670/1‑100 (MQ=60)
                                                                   CCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGC                         >  SRR3722078.1323352/1‑100 (MQ=60)
                                                                     TTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCCCCAGCCTGGCTGGCGG                       <  SRR3722078.892531/100‑1 (MQ=60)
                                                                     TTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGG                       <  SRR3722078.641829/100‑1 (MQ=60)
                                                                     TTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGG                       <  SRR3722078.947831/100‑1 (MQ=60)
                                                                        GTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGT                    >  SRR3722078.490959/1‑100 (MQ=60)
                                                                        GTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGT                    >  SRR3722078.464148/1‑100 (MQ=60)
                                                                        GTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGGGT                    >  SRR3722078.1468569/1‑100 (MQ=60)
                                                                          TCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGG                  <  SRR3722078.1573467/100‑1 (MQ=60)
                                                                                 TTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGT           <  SRR3722078.592676/100‑1 (MQ=60)
                                                                                  TGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGGAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTT          >  SRR3722078.1093736/1‑100 (MQ=60)
                                                                                         GGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGGGTGGGAAGCGTTGCCGGAG   >  SRR3722078.743064/1‑100 (MQ=60)
                                                                                          GTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGC  <  SRR3722078.437108/100‑1 (MQ=60)
                                                                                          GTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGC  >  SRR3722078.374525/1‑100 (MQ=60)
                                                                                              |                                                                                               
ACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGC  >  NZ_CP009273/3727930‑3728119

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