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 F41 I0 R1 212 10.4 395026 92.3% 364608 148.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 9,058 C→A 70.5% A274E (GCG→GAG)  talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092739,0580CA70.5% 53.5 / 26.7 37A274E (GCG→GAG) talBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base C (3/8);  new base A (12/14);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTC  >  CP009273/8911‑9197
                                                                                                                                                   |                                                                                                                                           
cAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGgagc                                                                                                                                           >  1:157288/1‑150 (MQ=255)
 aaGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGagcg                                                                                                                                          <  1:114722/150‑1 (MQ=255)
  aGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGagcg                                                                                                                                          >  1:54782/1‑149 (MQ=255)
    cTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTAt                                                                                                                                       >  2:26170/1‑150 (MQ=255)
       ccGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATc                                                                                                                                      <  1:160757/148‑1 (MQ=255)
            acgttggCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGaa                                                                                                                                <  2:76814/144‑2 (MQ=255)
             cgttggCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATc                                                                                                                              <  2:176142/146‑1 (MQ=255)
               tCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCg                                                                                                                            <  1:5890/150‑1 (MQ=255)
                 ggCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCCTCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGCAGGGGCTATCGAACGTAAACTGTCCTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAg                                                                                                                          >  2:159742/1‑150 (MQ=255)
                        ttCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCtgt                                                                                                                   >  1:92404/1‑150 (MQ=255)
                          cTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCtgtg                                                                                                                  <  1:139394/149‑1 (MQ=255)
                                            gACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATcc                                                                                                >  1:10300/1‑149 (MQ=255)
                                            gACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATcc                                                                                                >  2:173408/1‑149 (MQ=255)
                                             aCCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCaa                                                                                              <  2:122902/150‑1 (MQ=255)
                                                         tCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGaca                                                                                  <  1:14118/150‑3 (MQ=255)
                                                              ccGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATccga                                                                                              >  1:110610/1‑131 (MQ=255)
                                                                        tgAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTa                                                                   >  1:73313/1‑150 (MQ=255)
                                                                           aaGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATcc                                                                <  1:26170/150‑1 (MQ=255)
                                                                                   tccgaaagCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGtt                                                         <  2:112853/143‑1 (MQ=255)
                                                                                    cGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTg                                                       <  1:147552/150‑1 (MQ=255)
                                                                                       aaagCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGcga                                                    >  1:99271/3‑148 (MQ=255)
                                                                                         agCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGAtc                                                  <  1:134147/150‑2 (MQ=255)
                                                                                            gAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGAc                                               >  2:123044/1‑150 (MQ=255)
                                                                                                 ggCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGa                                          <  1:72400/150‑1 (MQ=255)
                                                                                                  gCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGaa                                         <  1:125336/150‑1 (MQ=255)
                                                                                                      tCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGaaaaa                                      <  1:145340/149‑1 (MQ=255)
                                                                                                                    gtatatacCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATc                       >  1:10501/7‑150 (MQ=255)
                                                                                                                     tatatacCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCg                      <  1:138557/145‑1 (MQ=255)
                                                                                                                       tatacCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCgac                    <  1:190215/147‑3 (MQ=255)
                                                                                                                         tacCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACg                   >  2:152619/2‑149 (MQ=255)
                                                                                                                            cGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATct               <  1:159742/150‑1 (MQ=255)
                                                                                                                             ggtgAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActg              >  1:71862/4‑150 (MQ=255)
                                                                                                                             ggtgAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActg              <  2:183261/147‑1 (MQ=255)
                                                                                                                               tgAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActgct            <  2:92404/149‑1 (MQ=255)
                                                                                                                               cGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATctgct            >  2:160600/1‑150 (MQ=255)
                                                                                                                                gAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActgctg           >  1:123210/1‑150 (MQ=255)
                                                                                                                                 aaGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAACTGCTGt          <  2:82091/150‑1 (MQ=255)
                                                                                                                                         gcgcGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAACTGCTGTaagagaat  <  1:63609/150‑7 (MQ=255)
                                                                                                                                             gTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActgct            >  2:64725/1‑136 (MQ=255)
                                                                                                                                              tCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActgctg           <  1:64725/136‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
CAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTC  >  CP009273/8911‑9197

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ 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
TCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTC  >  CP009273/8917‑9197
                                                                                                                                             |                                                                                                                                           
TCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATC                                                                                                                                      <  SRR12117902.163690/149‑1 (MQ=60)
        ATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCG                                                                                                                            <  SRR12117902.5974/151‑1 (MQ=60)
                                      gaccgcctgaccatcgccccggcgctgctgaaagagctgtccgaaagcgaaggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCG                                                                                               >  SRR12117902.10479/79‑150 (MQ=33)
                                                  atcgccccggcgctgctgaaagagctgtccgaaagcgaaggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACA                                                                                  <  SRR12117902.14350/85‑1 (MQ=33)
                                                    cgccccggcgctgctgaaagagctgtccgaaagcgaaggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAA                                                                                <  SRR12117902.27670/87‑1 (MQ=32)
                                                     gccccggcgctgctgaaagagctgtccgaaagcgaaggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAAC                                                                               <  SRR12117902.176547/88‑1 (MQ=32)
                                                       cccggcgctgctgaaagagctgtccgaaagcgaaggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTG                                                                             <  SRR12117902.81740/90‑1 (MQ=35)
                                                       cccggcgctgctgaaagagctgtccgaaagcgaaggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTG                                                                             <  SRR12117902.139681/90‑1 (MQ=35)
                                                        ccggcgctgctgaaagagctgtccgaaagcgaaggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGAT                                                                                             >  SRR12117902.112607/61‑134 (MQ=32)
                                                                  tgaaagagctgtccgaaagcgaaggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTAT                                                                  >  SRR12117902.74652/51‑151 (MQ=54)
                                                                    aaagagctgtccgaaagcgaaggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCC                                                                <  SRR12117902.26636/103‑1 (MQ=57)
                                                                             GCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTG                                                       <  SRR12117902.150250/151‑1 (MQ=60)
                                                                                 aaagcgaaggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGAT                                                   >  SRR12117902.101092/36‑151 (MQ=60)
                                                                                  aagcgaaggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCgatc                                                  <  SRR12117902.136594/117‑5 (MQ=60)
                                                                                        aggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACC                                              >  SRR12117902.30191/29‑149 (MQ=60)
                                                                                         ggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGC                                                                                          <  SRR12117902.187461/77‑1 (MQ=33)
                                                                                          GGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGA                                          <  SRR12117902.73720/151‑1 (MQ=60)
                                                                                           GGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAA                                         <  SRR12117902.127624/151‑1 (MQ=60)
                                                                                               gttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAAC                                     <  SRR12117902.43878/130‑1 (MQ=60)
                                                                                               ATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAA                                      <  SRR12117902.147998/150‑1 (MQ=60)
                                                                                                     cgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAA                                >  SRR12117902.88019/16‑150 (MQ=60)
                                                                                                          gctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATG                          >  SRR12117902.93627/11‑151 (MQ=60)
                                                                                                              gtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCG                      <  SRR12117902.141090/145‑1 (MQ=60)
                                                                                                              gtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCG                      >  SRR12117902.10682/7‑151 (MQ=60)
                                                                                                                atatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGAC                    <  SRR12117902.193673/147‑1 (MQ=60)
                                                                                                                     CCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCT               <  SRR12117902.162664/151‑1 (MQ=60)
                                                                                                                       GGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAACTGC             >  SRR12117902.73171/1‑151 (MQ=60)
                                                                                                                          GAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAACTGCTGT          >  SRR12117902.125456/1‑151 (MQ=60)
                                                                                                                                  AGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAACTGCTGTAAgagaat  <  SRR12117902.64755/151‑7 (MQ=60)
                                                                                                                                       GTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAACTGCTG           <  SRR12117902.65891/137‑1 (MQ=60)
                                                                                                                                             |                                                                                                                                           
TCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTC  >  CP009273/8917‑9197

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