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 4,367,563 T→C intergenic (‑24/+65) blc ← / ← ampC lipocalin Blc/BlaEC family class C beta‑lactamase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,367,5630TC96.5% 187.8 / ‑1.0 58intergenic (‑24/+65)blc/ampClipocalin Blc/BlaEC family class C beta‑lactamase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (23/33);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGA  >  NZ_CP009273/4367477‑4367646
                                                                                      |                                                                                   
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                  <  1:1069957/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                  <  2:539821/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                  <  2:490525/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                  <  2:470209/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                  <  1:1198621/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                  <  1:99601/90‑1 (MQ=255)
  aaCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAg                                                                                <  2:1061712/90‑1 (MQ=255)
  aaCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAg                                                                                <  1:247064/90‑1 (MQ=255)
  aaCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAg                                                                                <  1:568441/90‑1 (MQ=255)
  aaCTGCAGGCAACGACCAGAAATGCAGCTGTCGATGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAg                                                                                <  2:953092/90‑1 (MQ=255)
   aCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAgt                                                                               >  2:323103/1‑90 (MQ=255)
   aCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAgt                                                                               >  1:903606/1‑90 (MQ=255)
   aCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAgt                                                                               >  1:296297/1‑90 (MQ=255)
   aCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAgt                                                                               >  2:300036/1‑90 (MQ=255)
   aCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCTAAgtg                                                                              >  1:554963/1‑90 (MQ=255)
     tGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAgtgt                                                                             <  1:255950/90‑1 (MQ=255)
      gCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTa                                                                            >  2:397354/1‑90 (MQ=255)
      gCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTa                                                                            >  1:629447/1‑90 (MQ=255)
      gCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTa                                                                            >  1:950311/1‑90 (MQ=255)
        aGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGa                                                                          >  2:925292/1‑90 (MQ=255)
           cAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGa                                                                       >  2:257302/1‑90 (MQ=255)
           cAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGAGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGa                                                                       >  1:357112/1‑90 (MQ=255)
              cGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAg                                                                    <  1:718194/90‑1 (MQ=255)
                      aaTGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGaaa                                                            <  2:932443/90‑1 (MQ=255)
                                 cGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGaaaa                                                 >  2:671280/1‑90 (MQ=255)
                                      cGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTcc                                            >  1:762644/1‑90 (MQ=255)
                                      cGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTcc                                            >  1:355770/1‑90 (MQ=255)
                                      cGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGGGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTcc                                            >  1:225365/1‑90 (MQ=255)
                                         cAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGaa                                         <  1:552961/90‑1 (MQ=255)
                                             gagaGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAAtt                                     <  1:759927/90‑1 (MQ=255)
                                                  ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAc                                <  2:903606/90‑1 (MQ=255)
                                                    aGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAccc                              >  2:345180/1‑90 (MQ=255)
                                                       aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                           >  1:848866/1‑90 (MQ=255)
                                                       aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                           <  1:1025897/90‑1 (MQ=255)
                                                       aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                           <  1:397354/90‑1 (MQ=255)
                                                       aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                           <  1:1144396/90‑1 (MQ=255)
                                                       aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                           <  2:1078104/90‑1 (MQ=255)
                                                         gcgcATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATgg                         >  2:1114845/1‑90 (MQ=255)
                                                            cATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAt                      <  2:774708/90‑1 (MQ=255)
                                                            cATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAt                      <  1:968957/90‑1 (MQ=255)
                                                             aTAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAtt                     >  1:1072372/1‑90 (MQ=255)
                                                              tAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                    <  2:115589/90‑1 (MQ=255)
                                                              tAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                    <  2:115590/90‑1 (MQ=255)
                                                              tAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                    <  1:457613/90‑1 (MQ=255)
                                                              tAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                    <  1:185474/90‑1 (MQ=255)
                                                                  tGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACt                >  2:810792/1‑90 (MQ=255)
                                                                      tCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTag            >  1:908018/1‑90 (MQ=255)
                                                                      tCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTag            <  1:632998/90‑1 (MQ=255)
                                                                       ccTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTaga           <  1:1168189/90‑1 (MQ=255)
                                                                       ccTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTaga           <  2:440060/90‑1 (MQ=255)
                                                                       ccTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTaga           <  1:502162/90‑1 (MQ=255)
                                                                       ccTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTaga           <  1:94807/90‑1 (MQ=255)
                                                                         ttACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc         <  1:635465/90‑1 (MQ=255)
                                                                          tACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc         >  2:739978/1‑89 (MQ=255)
                                                                          tACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc         <  1:739978/89‑1 (MQ=255)
                                                                             tGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGc                                                         <  2:609924/38‑1 (MQ=255)
                                                                             tGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGc                                                         >  1:609924/1‑38 (MQ=255)
                                                                             tGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTa     >  1:862374/1‑90 (MQ=255)
                                                                                ttttttcCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGACGGAATTTTACTGTAGAGCGTTAAGa  >  2:1116243/1‑90 (MQ=255)
                                                                                      |                                                                                   
AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGA  >  NZ_CP009273/4367477‑4367646

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGG  >  NZ_CP009273/4367470‑4367659
                                                                                             |                                                                                                
GCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGT                                                                                            >  SRR3722089.916729/1‑100 (MQ=60)
GCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGT                                                                                            >  SRR3722089.300079/1‑100 (MQ=60)
GCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCTAAGTG                                                                                           >  SRR3722089.561407/1‑100 (MQ=60)
   TAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTA                                                                                         >  SRR3722089.637421/1‑100 (MQ=60)
   TAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTA                                                                                         >  SRR3722089.964274/1‑100 (MQ=60)
   TAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTA                                                                                         >  SRR3722089.1218845/1‑100 (MQ=60)
       AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATG                                                                                     <  SRR3722089.1086026/100‑1 (MQ=60)
       AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATG                                                                                     <  SRR3722089.100846/100‑1 (MQ=60)
       AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATG                                                                                     <  SRR3722089.1216847/100‑1 (MQ=60)
        GAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGAGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGA                                                                                    >  SRR3722089.361518/1‑100 (MQ=60)
         AACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGAC                                                                                   <  SRR3722089.575152/100‑1 (MQ=60)
         AACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGAC                                                                                   <  SRR3722089.250219/100‑1 (MQ=60)
            TGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGC                                                                                <  SRR3722089.259236/100‑1 (MQ=60)
                     CGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAG                                                                       <  SRR3722089.727532/100‑1 (MQ=60)
                                   GCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCC                                                         >  SRR3722089.360162/1‑100 (MQ=60)
                                   GCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCC                                                         >  SRR3722089.772859/1‑100 (MQ=60)
                                   GCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGGGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCC                                                         >  SRR3722089.228229/1‑100 (MQ=60)
                                                CAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACC                                            <  SRR3722089.559385/100‑1 (MQ=60)
                                                    GAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAT                                        >  SRR3722089.860935/1‑100 (MQ=60)
                                                    GAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAT                                        <  SRR3722089.770073/100‑1 (MQ=60)
                                                          GAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATT                                  >  SRR3722089.1088481/1‑100 (MQ=60)
                                                              AGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAC                              <  SRR3722089.402188/100‑1 (MQ=60)
                                                              AGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAC                              <  SRR3722089.1161811/100‑1 (MQ=60)
                                                              AGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAC                              <  SRR3722089.1041199/100‑1 (MQ=60)
                                                                   CATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAG                         >  SRR3722089.921218/1‑100 (MQ=60)
                                                                   CATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAG                         <  SRR3722089.983263/100‑1 (MQ=60)
                                                                     TAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAG                       <  SRR3722089.463091/100‑1 (MQ=60)
                                                                     TAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAG                       <  SRR3722089.187789/100‑1 (MQ=60)
                                                                          GTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGctgtctcttatacacatctgacgctgccgacgaccatcttag                  >  SRR3722089.617515/1‑58 (MQ=60)
                                                                          GTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTA                  >  SRR3722089.874714/1‑100 (MQ=60)
                                                                             TCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGA               <  SRR3722089.641048/100‑1 (MQ=60)
                                                                              CCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAA              <  SRR3722089.1185960/100‑1 (MQ=60)
                                                                              CCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAA              <  SRR3722089.508094/100‑1 (MQ=60)
                                                                              CCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAA              <  SRR3722089.95984/100‑1 (MQ=60)
                                                                                TTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATC            <  SRR3722089.643560/100‑1 (MQ=60)
                                                                                TTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATC            <  SRR3722089.749672/100‑1 (MQ=60)
                                                                                     GGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCtgtctcttatacaca       >  SRR3722089.371268/1‑85 (MQ=60)
                                                                                          TTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGG  >  SRR3722089.1146717/1‑100 (MQ=60)
                                                                                             |                                                                                                
GCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGG  >  NZ_CP009273/4367470‑4367659

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