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 F1 I1 R1 343 268.7 10873516 97.0% 10547310 150.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,031,017 Δ1 bp intergenic (+159/‑625) RPPX_17905 → / → RPPX_17915 23S rRNA methyltransferase/sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,031,0170A.91.7% 109.3 / 8.2 36intergenic (+159/‑625)RPPX_17905/RPPX_1791523S rRNA methyltransferase/sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base . (17/16);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGG  >  CP009974/4030871‑4031167
                                                                                                                                                  |                                                                                                                                                       
gtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGccc                                                                                                                                                    >  1:3106084/1‑151 (MQ=255)
gtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGccc                                                                                                                                                    >  1:3681641/1‑151 (MQ=16)
gtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGccc                                                                                                                                                    >  2:270408/1‑151 (MQ=16)
 ttctgcttctgcttctgcttctgcttctgcttctgcTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGcccc                                                                                                                                                   <  1:697406/151‑1 (MQ=255)
 ttctgcttctgcttctgcttctgcttctgcttctgcTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGcccc                                                                                                                                                   <  1:852286/151‑1 (MQ=255)
  tctgcttctgcttctgcttctgcttctgcttctgctTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCa                                                                                                                                                  <  2:5059353/151‑1 (MQ=16)
   ttgcttcggattctgcttctgcttctgcttctgcttTTGCTTCTCAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAg                                                                                                                                                <  1:2073309/150‑1 (MQ=11)
    tgcttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGCCGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAgg                                                                                                                                               >  2:2142805/1‑151 (MQ=255)
     gcttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGa                                                                                                                                              >  1:1763211/1‑151 (MQ=16)
     gctgctgcttctgcttctgcttctgcttctgcttTTGCTTCTCAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGACGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGa                                                                                                                                              <  2:3481242/151‑1 (MQ=11)
       ttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATg                                                                                                                                            >  2:1052902/1‑151 (MQ=16)
       ttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATg                                                                                                                                            <  2:3049066/151‑1 (MQ=16)
       ttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATg                                                                                                                                            <  1:5125066/151‑1 (MQ=16)
         ctgcttctgcttctgcttctgcttctgcttTTGCTTCTAGGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATggg                                                                                                                                          <  1:5295473/151‑1 (MQ=255)
         ctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATggg                                                                                                                                          >  1:323683/1‑151 (MQ=255)
         ctgcttctgcttcggcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATggg                                                                                                                                          <  1:3858899/151‑1 (MQ=255)
           gcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGc                                                                                                                                        <  2:1763211/151‑1 (MQ=16)
            cttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGcc                                                                                                                                       >  2:317539/1‑151 (MQ=255)
            cttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGcc                                                                                                                                       >  2:587619/1‑151 (MQ=12)
             ttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                                                                                                                                      <  2:714098/151‑1 (MQ=16)
             ttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                                                                                                                                      >  1:1211761/1‑151 (MQ=255)
             ttcggcttctgcttcggcttctgcttTTGCGTCAACGCGCGCGATGGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGCGCGCAGGGCTGGCGAAGGGGGGTGACGGGCTTTGATGCCCGTCACGCGCTGAGGACCCATGGATGGGGc                                                                                                                                        <  2:2858868/151‑1 (MQ=11)
              tctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtg                                                                                                                                     >  2:3273285/1‑150 (MQ=16)
              tctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtg                                                                                                                                     >  1:1182011/1‑150 (MQ=255)
                tgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGCGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgca                                                                                                                                   <  1:4093380/151‑4 (MQ=255)
                tgcgtctgcttctggttctgcttGTGCTTCTAGGGACGCGATAGGTGGGGCGACTCAGTTTGCTGCTTGAGGAGGCGGAGTGCAGGGCTTGTGAAGGGAGGTGACGTGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgca                                                                                                                                   <  2:4349027/151‑4 (MQ=2)
                  cttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagc                                                                                                                                 <  2:3074794/151‑6 (MQ=2)
                  cttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagc                                                                                                                                 >  1:3338203/1‑146 (MQ=255)
                   ttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagcg                                                                                                                                <  1:3939573/151‑7 (MQ=255)
                   ttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagcg                                                                                                                                >  2:346277/1‑145 (MQ=255)
                   ttctgctgcggcttcggcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagcg                                                                                                                                <  2:559177/151‑7 (MQ=1)
                      tgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCgg                                                                                                                             >  1:4216857/1‑151 (MQ=255)
                       gcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGt                                                                                                                            >  1:1951132/1‑151 (MQ=255)
                          tctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCt                                                                                                                         >  1:399776/1‑151 (MQ=255)
                                                                                                                                                 gAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAg   <  1:850631/151‑1 (MQ=11)
                                                                                                                                                 gAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAg   >  2:4130074/1‑151 (MQ=11)
                                                                                                                                                  aGGCCGGAGGATGGGGCGTGCGGGGCGGTCCTGCCGGGAGGAAGGCCGGAGTGAGGGGACGCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAgg  <  2:3162067/151‑1 (MQ=11)
                                                                                                                                                  |                                                                                                                                                       
GTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGG  >  CP009974/4030871‑4031167

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAG  >  CP009974/4030867‑4031166
                                                                                                                                                      |                                                                                                                                                     
aacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGG                                                                                                                                                       >  SRR10918230.3660964/39‑151 (MQ=0)
CGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGG                                                                                                                                                       <  SRR10918230.1246326/151‑1 (MQ=40)
 acgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGG                                                                                                                                                      >  SRR10918230.640075/38‑151 (MQ=0)
    GTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCC                                                                                                                                                  >  SRR10918230.3697619/1‑151 (MQ=22)
    GTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCC                                                                                                                                                  >  SRR10918230.3119618/1‑151 (MQ=22)
     TTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCC                                                                                                                                                 <  SRR10918230.855610/151‑1 (MQ=18)
     TTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCC                                                                                                                                                 <  SRR10918230.699962/151‑1 (MQ=18)
       ttgcttcggaTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTCAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAG                                                                                                                                               <  SRR10918230.2082066/141‑1 (MQ=16)
         GCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGA                                                                                                                                             >  SRR10918230.1770844/1‑151 (MQ=16)
           TTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATG                                                                                                                                           <  SRR10918230.5148601/151‑1 (MQ=16)
             CTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAGGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGG                                                                                                                                         <  SRR10918230.5320124/151‑1 (MQ=16)
             CTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGG                                                                                                                                         >  SRR10918230.324772/1‑151 (MQ=16)
             CTGCTTCTGCTTCGGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGG                                                                                                                                         <  SRR10918230.3875742/151‑1 (MQ=16)
                 TTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT                                                                                                                                     >  SRR10918230.1217039/1‑151 (MQ=16)
                  TCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTG                                                                                                                                    >  SRR10918230.1187106/1‑151 (MQ=11)
                    TGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGCGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTgca                                                                                                                                  <  SRR10918230.4111624/151‑4 (MQ=7)
                      CTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTgcagc                                                                                                                                >  SRR10918230.3352682/1‑146 (MQ=3)
                       TTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCG                                                                                                                               <  SRR10918230.3956859/151‑1 (MQ=1)
                                                                                                                     GAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTT                                  <  SRR10918230.3810202/151‑1 (MQ=0)
                                                                                                                       GGTGACGGGCATGGGTGCCCGTCAAGCGCGGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGT                                <  SRR10918230.474656/151‑1 (MQ=0)
                                                                                                                        GTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTC                               <  SRR10918230.840728/151‑1 (MQ=0)
                                                                                                                          GACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGC                             >  SRR10918230.3659067/1‑151 (MQ=0)
                                                                                                                              GGCTTGGTTGCCCGTCAAGCGCGGAGGCGCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTT                         <  SRR10918230.3333536/151‑1 (MQ=0)
                                                                                                                               GCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTG                        <  SRR10918230.745609/151‑1 (MQ=0)
                                                                                                                                 ATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGAC                      <  SRR10918230.4524896/151‑1 (MQ=0)
                                                                                                                                       GCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAG                <  SRR10918230.3801178/151‑1 (MQ=0)
                                                                                                                                         CCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTG              >  SRR10918230.3528011/1‑151 (MQ=0)
                                                                                                                                          CGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGA             >  SRR10918230.3729345/1‑151 (MQ=0)
                                                                                                                                           GTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGCCACAAAGTCAC            >  SRR10918230.2287252/1‑151 (MQ=0)
                                                                                                                                             CAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTC          <  SRR10918230.1940885/151‑1 (MQ=0)
                                                                                                                                             CAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTC          >  SRR10918230.3612786/1‑151 (MQ=0)
                                                                                                                                              AAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCG         <  SRR10918230.4262391/151‑1 (MQ=0)
                                                                                                                                              AAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCG         <  SRR10918230.794331/151‑1 (MQ=0)
                                                                                                                                                  GCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGT     <  SRR10918230.3001215/151‑1 (MQ=0)
                                                                                                                                                   CTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTA    >  SRR10918230.3570291/1‑151 (MQ=0)
                                                                                                                                                   CTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTA    >  SRR10918230.1157827/1‑151 (MQ=0)
                                                                                                                                                    TGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAA   <  SRR10918230.363219/151‑1 (MQ=0)
                                                                                                                                                     GAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAG  <  SRR10918230.853944/151‑1 (MQ=0)
                                                                                                                                                      |                                                                                                                                                     
CGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAG  >  CP009974/4030867‑4031166

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