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,068 T→C 69.6% T277T (ACT→ACC talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092739,0680TC69.6% 47.0 / 24.2 36T277T (ACT→ACCtalBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base T (4/7);  new base C (11/14);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

ACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGT  >  CP009273/8923‑9215
                                                                                                                                                 |                                                                                                                                                   
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                                <  2:183261/147‑1 (MQ=255)
                                                                                                                 ggtgAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActg                                >  1:71862/4‑150 (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)
                                                                                                                                               aCTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGt  >  1:191113/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGT  >  CP009273/8923‑9215

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ 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
ATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTC  >  CP009273/8925‑9216
                                                                                                                                               |                                                                                                                                                    
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.139681/90‑1 (MQ=35)
                                               cccggcgctgctgaaagagctgtccgaaagcgaaggcgccgttgagcgtaagctggtatatACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTG                                                                                                <  SRR12117902.81740/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)
                                                                                                                                             ACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTC  >  SRR12117902.194589/1‑151 (MQ=60)
                                                                                                                                               |                                                                                                                                                    
ATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTC  >  CP009273/8925‑9216

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: