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
A2 F5 I175 R1 328 64.7 2582366 96.2% 2484236 120.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 861,119 G→A R234R (CGC→CGT moeA ← molybdopterin molybdotransferase MoeA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273861,1190GA100.0% 163.5 / NA 43R234R (CGC→CGTmoeAmolybdopterin molybdotransferase MoeA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/26);  total (17/26)

GTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCGCGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCGGGCAACTG  >  NZ_CP009273/860983‑861253
                                                                                                                                        |                                                                                                                                      
gTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGgcacg                                                                                                                                      >  2:1168703/1‑139 (MQ=255)
gTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGgcacg                                                                                                                                      >  1:362994/1‑139 (MQ=255)
gTTTAATCGCCAGCTTCCAGAAGGCGACCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGgcacg                                                                                                                                      <  2:190551/139‑1 (MQ=255)
gTTTAATCGCCAGCGTCCAGAAGGCGATCTCCCCCAGCCCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGgcacg                                                                                                                                      <  2:933624/139‑1 (MQ=255)
 tttAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGgcacgc                                                                                                                                     <  2:1221665/139‑1 (MQ=255)
   tAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCa                                                                                                                                    <  1:872730/138‑1 (MQ=255)
   tAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCa                                                                                                                                    <  2:723487/138‑1 (MQ=255)
                       gCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGa                                                                                                               <  2:500503/139‑1 (MQ=255)
                       gCGATCGCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGa                                                                                                               <  2:736032/139‑1 (MQ=255)
                        cGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGAt                                                                                                              <  2:838441/139‑1 (MQ=255)
                        cGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGAt                                                                                                              <  2:1099563/139‑1 (MQ=255)
                               ccccAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATcc                                                                                                         <  2:441423/137‑1 (MQ=255)
                               ccccAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCt                                                                                                       <  1:416545/139‑1 (MQ=255)
                                 ccAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTaa                                                                                                     <  1:854692/139‑1 (MQ=255)
                                        ttCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTaa                                                                                                <  1:456904/137‑1 (MQ=255)
                                          cAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAAtt                                                                                              <  1:377745/137‑1 (MQ=255)
                                          cAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAAt                                                                                               <  1:236266/136‑1 (MQ=255)
                                              aaTCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCg                                                                                        <  1:1211628/139‑1 (MQ=255)
                                                tCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCt                                                                                                       >  2:890795/1‑122 (MQ=255)
                                                tCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCt                                                                                                       <  1:890795/122‑1 (MQ=255)
                                                  gTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATc                                                                                    >  1:1037375/1‑139 (MQ=255)
                                                  gTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATc                                                                                    >  2:885690/1‑139 (MQ=255)
                                                      tGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCa                                                                                 <  1:885690/138‑1 (MQ=255)
                                                       ggTGTAATCCGCCTCGCCCACTGAAACACCGCCGGAGCTGATCACCACATCCGCCTGGCTGTCGTCTTCAATAAATGCGGCACGCAGGGCATGCGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCaa                                                                                <  2:1136803/138‑1 (MQ=255)
                                                             aTCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTc                                                                          >  1:829673/1‑138 (MQ=255)
                                                             aTCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCt                                                                         <  1:87230/139‑1 (MQ=255)
                                                                  ccTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACAt                                                                    >  1:822476/1‑139 (MQ=255)
                                                                   cTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACAt                                                                    >  1:979337/1‑138 (MQ=255)
                                                                   cTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATc                                                                   <  1:316367/139‑1 (MQ=255)
                                                                                     gccgGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACgg                                                 >  2:173036/1‑139 (MQ=255)
                                                                                      ccgGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAg                                                                                                                                   >  2:736010/1‑56 (MQ=255)
                                                                                      ccgGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAg                                                                                                                                   <  1:736010/56‑1 (MQ=255)
                                                                                      ccgGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGt                                                >  1:271299/1‑139 (MQ=255)
                                                                                        ggAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTg                                              <  1:1168703/139‑1 (MQ=255)
                                                                                              gATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCg                                        >  1:541071/1‑139 (MQ=255)
                                                                                               aTCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCg                                        >  2:605232/1‑138 (MQ=255)
                                                                                                         ccGCCTGGCTGTCGGCTTCAATAAATGCGgcacgc                                                                                                                                     <  1:175389/35‑1 (MQ=255)
                                                                                                         ccGCCTGGCTGTCGGCTTCAATAAATGCGgcacgc                                                                                                                                     >  2:175389/1‑35 (MQ=255)
                                                                                                           gCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGt                           <  2:829673/139‑1 (MQ=255)
                                                                                                           gCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGt                           <  2:822476/139‑1 (MQ=255)
                                                                                                            ccTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTc                          >  2:120729/1‑139 (MQ=255)
                                                                                                              tGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGCGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGc                        >  2:88171/1‑139 (MQ=255)
                                                                                                                gCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGcc                       <  1:496907/138‑1 (MQ=255)
                                                                                                                gCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGcc                       <  2:189815/138‑1 (MQ=255)
                                                                                                                     cGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCgg                 >  1:996942/1‑139 (MQ=255)
                                                                                                                           cAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATccc                                                                                  <  2:343463/68‑1 (MQ=255)
                                                                                                                           cAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATccc                                                                                  >  1:343463/1‑68 (MQ=255)
                                                                                                                           cAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGAcc           >  1:573137/1‑139 (MQ=255)
                                                                                                                           cAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGAc            <  2:362994/138‑1 (MQ=255)
                                                                                                                                  tGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCGGGCAAc    >  2:1002877/1‑139 (MQ=255)
                                                                                                                                      gcgcgcAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCGGGCAACTg  >  2:218747/1‑137 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCGCGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCGGGCAACTG  >  NZ_CP009273/860983‑861253

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCGCGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCG  >  NZ_CP009273/860972‑861245
                                                                                                                                                   |                                                                                                                              
CGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGC                                                                                                                             >  SRR3722129.368239/1‑151 (MQ=60)
           atgtccgttgagtattttttttttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGG                                                                                                                  <  SRR3722129.177798/57‑1 (MQ=60)
             TTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGA                                                                                                                 <  SRR3722129.885189/150‑1 (MQ=60)
              ttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagCTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGAT                                                                                                                <  SRR3722129.746791/78‑1 (MQ=60)
                                         gagaCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTAC                                                                                    <  SRR3722129.903457/147‑1 (MQ=60)
                                         CCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTAC                                                                                    <  SRR3722129.422621/151‑1 (MQ=60)
                                           CCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCT                                                                                  <  SRR3722129.866890/151‑1 (MQ=60)
                                                  CTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCC                                                                           >  SRR3722129.1051962/1‑151 (MQ=60)
                                                  CTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCAT                                                                             <  SRR3722129.463503/149‑1 (MQ=60)
                                                    TCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCC                                                                           <  SRR3722129.383263/149‑1 (MQ=60)
                                                    TCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATC                                                                            <  SRR3722129.239632/148‑1 (MQ=60)
                                                        GAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTG                                                                     <  SRR3722129.1228518/151‑1 (MQ=60)
                                                             GTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCT                                                                 >  SRR3722129.841598/1‑150 (MQ=60)
                                                                TTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAAC                                                              <  SRR3722129.898300/150‑1 (MQ=60)
                                                                  GGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATC                                                           >  SRR3722129.834306/1‑151 (MQ=60)
                                                                   GTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATC                                                           >  SRR3722129.993190/1‑150 (MQ=60)
                                                                       AATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTG                                                      <  SRR3722129.88365/151‑1 (MQ=60)
                                                                             CCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGC                                                <  SRR3722129.320886/151‑1 (MQ=60)
                                                                                      CTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTT                                       >  SRR3722129.275147/1‑151 (MQ=60)
                                                                                              CCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGT                               >  SRR3722129.548989/1‑151 (MQ=60)
                                                                                                  CGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGAT                           <  SRR3722129.1184988/151‑1 (MQ=60)
                                                                                                                     CGCCTGGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGC        >  SRR3722129.1011001/1‑151 (MQ=60)
                                                                                                                          GGCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGAC    <  SRR3722129.504086/150‑1 (MQ=60)
                                                                                                                           GCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttct  >  SRR3722129.348420/1‑91 (MQ=60)
                                                                                                                           GCTGTCGGCTTCAATAAATGCGGCACGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCG  >  SRR3722129.581616/1‑151 (MQ=60)
                                                                                                                                                   |                                                                                                                              
CGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATAAATGCGGCGCGCAGGGCATGGGGATCGTCGCGGATAATCCCTAAGTTAATTACCTCGCATCCCAACTGTTCTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCG  >  NZ_CP009273/860972‑861245

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