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 F5 I97 R1 275 44.0 1733530 96.5% 1672856 122.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,556,234 T→C E277G (GAA→GGA)  yjiN ← DUF445 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,556,2340TC95.9% 143.9 / ‑0.4 49E277G (GAA→GGA) yjiNDUF445 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (19/28);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

CACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTTCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATG  >  NZ_CP009273/4556103‑4556367
                                                                                                                                   |                                                                                                                                     
cacGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATc                                                                                                                                >  1:144917/1‑139 (MQ=255)
cacGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATc                                                                                                                                >  2:162393/1‑139 (MQ=255)
            cACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCa                                                                                                                        <  1:820272/135‑1 (MQ=255)
                 ttGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGa                                                                                                               >  2:178499/1‑139 (MQ=255)
                  tGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCa                                                                                                                        <  2:350000/129‑1 (MQ=255)
                  tGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCa                                                                                                                        >  1:350000/1‑129 (MQ=255)
                   gATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGAt                                                                                                              <  1:362845/138‑1 (MQ=255)
                    atatCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATc                                                                                                             <  2:324702/138‑1 (MQ=255)
                       tCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGAt                                                                                                              >  2:478677/1‑134 (MQ=255)
                       tCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGAt                                                                                                              <  1:478677/134‑1 (MQ=255)
                          aCTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGa                                                                                                      <  2:130152/139‑1 (MQ=255)
                                      tCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGAt                                                                                          >  1:203191/1‑139 (MQ=255)
                                        ccGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATc                                                                                         <  1:784197/138‑1 (MQ=255)
                                         cGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCaa                                                                                       <  1:162393/139‑1 (MQ=255)
                                                 cccccACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATgg                                                                               <  2:780963/139‑1 (MQ=255)
                                                 cccccACAATTCACGGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCGTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATgg                                                                               <  1:307515/139‑1 (MQ=255)
                                                          ttCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGgc                                                                                               >  2:653834/1‑114 (MQ=255)
                                                          ttCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGgc                                                                                               <  1:653834/114‑1 (MQ=255)
                                                            cACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGAt                                                                     <  1:120655/138‑1 (MQ=255)
                                                             aCTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCa                                                                                                                        <  2:491173/86‑1 (MQ=255)
                                                             aCTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCa                                                                                                                        >  1:491173/1‑86 (MQ=255)
                                                                 agaTAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATgcgc                                                                 <  1:719408/137‑1 (MQ=255)
                                                                 agaTAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCa                                                                <  1:83672/138‑1 (MQ=255)
                                                                 agaTAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCa                                                                <  2:203191/138‑1 (MQ=255)
                                                                    tAGCGATTTAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGAt                                                            <  1:670963/139‑1 (MQ=255)
                                                                    tAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGAt                                                            <  1:467301/139‑1 (MQ=255)
                                                                       cGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCCCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATcgc                                                         >  2:209131/1‑139 (MQ=255)
                                                                         aTTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGa                                                       >  1:504832/1‑139 (MQ=255)
                                                                         aTTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGAc                                                      >  2:805559/1‑139 (MQ=255)
                                                                          ttAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGAc                                                      <  2:504832/139‑1 (MQ=255)
                                                                           tAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCa                                                                                                                        >  1:716102/1‑72 (MQ=255)
                                                                           tAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCa                                                                                                                        <  2:716102/72‑1 (MQ=255)
                                                                            aaaaGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACt                                                     >  2:433113/1‑138 (MQ=255)
                                                                                    cATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGgcg                                                                                              >  1:259783/1‑89 (MQ=255)
                                                                                    cATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGgcg                                                                                              <  2:259783/89‑1 (MQ=255)
                                                                                    cATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATatcat                                             >  1:200942/1‑138 (MQ=255)
                                                                                     aTCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATatcatc                                            <  1:433113/138‑1 (MQ=255)
                                                                                          cttTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCaa                                          <  1:714965/135‑1 (MQ=255)
                                                                                                     cTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATc                                                                                                             >  2:474034/1‑57 (MQ=255)
                                                                                                     cTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATc                                                                                                             <  1:474034/57‑1 (MQ=255)
                                                                                                       ttttACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACcgcg                         >  1:492138/1‑139 (MQ=255)
                                                                                                              gCATCGGCTCGCGCTGCCATTCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAAc                 <  1:805559/139‑1 (MQ=255)
                                                                                                               cATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAAc                 >  1:503546/1‑139 (MQ=255)
                                                                                                                     cTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTc           <  1:714764/139‑1 (MQ=255)
                                                                                                                         cgcTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCcgcg       <  1:754697/139‑1 (MQ=255)
                                                                                                                          gcTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGa                                                       <  1:533083/90‑1 (MQ=255)
                                                                                                                          gcTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGa                                                       >  2:533083/1‑90 (MQ=255)
                                                                                                                          gcTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCgcgc      <  1:178499/139‑1 (MQ=255)
                                                                                                                              ccATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATg  >  2:786691/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                     
CACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTTCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATG  >  NZ_CP009273/4556103‑4556367

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ 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
GCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTTCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACC  >  NZ_CP009273/4556087‑4556378
                                                                                                                                                   |                                                                                                                                                
GCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCC                                                                                                                                                <  SRR3722120.108949/150‑1 (MQ=60)
GCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCC                                                                                                                                                <  SRR3722120.607056/150‑1 (MQ=60)
   ATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGAT                                                                                                                                            <  SRR3722120.470560/151‑1 (MQ=60)
     GCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCG                                                                                                                                          >  SRR3722120.146873/1‑151 (MQ=60)
                       ATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATC                                                                                                                        >  SRR3722120.355058/1‑151 (MQ=60)
                           TCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATC                                                                                                                        <  SRR3722120.831307/147‑1 (MQ=60)
                                  TGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGG                                                                                                              <  SRR3722120.368126/150‑1 (MQ=60)
                                  TGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGG                                                                                                              <  SRR3722120.485590/150‑1 (MQ=60)
                                   tttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGT                                                                                                             <  SRR3722120.480887/79‑1 (MQ=60)
                                           CTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATC                                                                                                    >  SRR3722120.206009/1‑151 (MQ=60)
                                                gtgtataagagacagTCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACG                                                                                               <  SRR3722120.662830/136‑1 (MQ=60)
                                                       CCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGC                                                                                         <  SRR3722120.794760/150‑1 (MQ=60)
                                                        CCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGG                                                                                       <  SRR3722120.164621/151‑1 (MQ=60)
                                                                CCCCCCACAATTCACGGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCGTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATG                                                                               <  SRR3722120.311975/151‑1 (MQ=60)
                                                                  CCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCctgtctcttatacacatctccgagcccacgagactaaggcga                                                                              >  SRR3722120.498228/1‑108 (MQ=60)
                                                                           TCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCG                                                                     <  SRR3722120.122249/150‑1 (MQ=60)
                                                                              CTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGAC                                                                 >  SRR3722120.512053/1‑151 (MQ=60)
                                                                                GAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtc                                                               >  SRR3722120.725861/1‑94 (MQ=60)
                                                                                GAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACT                                                                <  SRR3722120.84748/150‑1 (MQ=60)
                                                                                GAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGAC                                                                 <  SRR3722120.729197/149‑1 (MQ=60)
                                                                                   ATAGCGATTTAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATA                                                            <  SRR3722120.680154/151‑1 (MQ=60)
                                                                                   ATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATA                                                            <  SRR3722120.474077/151‑1 (MQ=60)
                                                                                        ccttagtctcgtgggctcggagatgtgtataagagacaGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATC                                                       <  SRR3722120.540650/113‑1 (MQ=60)
                                                                                         ATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGGAGGTGGCGCGATCAAACGCctgtctcttatacacatctccgagcccacgagactaaggc                                                      >  SRR3722120.263535/1‑111 (MQ=60)
                                                                                         ATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATC                                                       >  SRR3722120.203727/1‑150 (MQ=60)
                                                                                                GCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTTCCGGATCGTTTTTCAGCTTTTCGATCAGGGCGAAGGAGGTCGGATCAAACGTATGGCGGATCTGATGGGCACGATCGCGCCGAATATTATCAAGCAA                                                 >  SRR3722120.241811/1‑149 (MQ=60)
                                                                                                    CATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAAT                                            <  SRR3722120.439460/150‑1 (MQ=60)
                                                                                                         TCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTC                                          <  SRR3722120.724712/147‑1 (MQ=60)
                                                                                                            TTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGT                                   >  SRR3722120.499224/1‑151 (MQ=60)
                                                                                                                    GCTTTTTACGGCATCGGCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACC                           >  SRR3722120.510749/1‑151 (MQ=60)
                                                                                                                             GGCATCGGCTCGCGCTGCCATTCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGC                 <  SRR3722120.816422/151‑1 (MQ=60)
                                                                                                                                    GCTCGCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTT           <  SRR3722120.724508/151‑1 (MQ=60)
                                                                                                                                        GCGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCC       <  SRR3722120.764897/151‑1 (MQ=60)
                                                                                                                                         CGCTGCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCC      <  SRR3722120.180931/151‑1 (MQ=60)
                                                                                                                                             GCCATTCCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACC  >  SRR3722120.695344/1‑151 (MQ=60)
                                                                                                                                                   |                                                                                                                                                
GCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTTCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACC  >  NZ_CP009273/4556087‑4556378

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