Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A2 F17 I0 R1 43 55.2 1410130 81.3% 1146435 215.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,560,455:1 +G 100% intergenic (‑2/+1) glpR ← / ← glpR pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,4551.G96.4% 74.1 / ‑4.5 27intergenic (‑2/+1)glpR/glpRpseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (21/5);  total (22/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCC  >  NC_000913/3560229‑3560661
                                                                                                                                                                                                                                    |                                                                                                                                                                                                              
aGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGGGCCGCACCGCCATGAGGGcg                                                                                                                                                                                                                >  1:423568/1‑227 (MQ=255)
    aCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATAGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGGGCCGCACGCCCATGCTGGcgcag                                                                                                                                                                                                             >  1:485732/1‑223 (MQ=255)
       ttGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgccggat                                                                                                                                                                                                          >  1:672036/1‑220 (MQ=255)
                ttgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACGGCCATGATGGCGCAGGAGCGGGtctgg                                                                                                                                                                                                 >  1:98828/1‑223 (MQ=255)
                 tgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACGGCCATGATGGCGCAGGATCAGGTTTTGCt                                                                                                                                                                                               >  1:540484/1‑228 (MQ=255)
                              cAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCCTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                                                                                                  >  1:556367/1‑228 (MQ=255)
                                           ttGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGCAGGCAGGGCCCCCCCGCCATAGTGGGCCGGGATAGGGTTTTGCCCCGCCGCCCCCTTGGGGGCGCGGGGa                                                                                                                                                                     >  2:222770/1‑228 (MQ=255)
                                                  gtgcgtgcGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt                                                                                                                                                              >  1:199225/1‑228 (MQ=255)
                                                   tgcgtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTg                                                                                                                                                             >  1:691380/1‑228 (MQ=255)
                                                   tgcgtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCAACGGCGTGTTAACCGAAGTGGAAAGCAGCGCCGCACCGCCAGGATGGCGCAGAATAAGGTTTGGCTCCGCCAGCCCATGTGGGCCGCGGGGAATAGtacg                                                                                                                                                             >  2:201846/1‑225 (MQ=255)
                                                                            gtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCGGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTTACGGGGAAAtccctgcc                                                                                                                                     >  1:354442/1‑220 (MQ=255)
                                                                                          aaaCGGCGTCCGGCCTTGGGGGATTTGCCCCGCCACGTGGCGGGCGATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                                                                                       <  1:222770/227‑1 (MQ=255)
                                                                                                  tCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATGATGGCGCGGGATCAGGTTTTGTTCCGCCAGCTCATTGAGGTCGCGGCGAAAAGTCTGCGGGCTGACGGGGCAAGGGCCTACCCGCTCTTCGGGCTTGACATaaa                                                                                                              >  2:423339/1‑227 (MQ=255)
                                                                                                                    gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGACACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGGCGCAGGATCCGGTTTTGATCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGGGGGCTGCCGGGGAAATTCTCTACCAGCTCTCCGGACCTGCCATAACCCCTGTGTTTAACCAAtt                                                                                            >  2:270163/1‑228 (MQ=255)
                                                                                                                     cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTACGGTCGCGGCGATTAGTCTGCGGCCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCGGTTTAACCAGTTc                                                                                           >  1:339024/1‑228 (MQ=255)
                                                                                                                          ccACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGGGAATAGTCTGCGGCCTGTCGGGGAAATGCTCTACCACCTCACCGGTCCTGACCGAAACCTGCTGTTTAACCCGTTcaagaa                                                                                      >  2:115913/1‑223 (MQ=255)
                                                                                                                               tggCGGGCGATGCGCTGTTTTTCTCGGGTCTGGGTGCCCGTGCGATCGTGCCACGGTGTGTTAGCCGAACTGGAGGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                                                                                 <  1:115913/226‑1 (MQ=255)
                                                                                                                                               ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAGGCTCTACCGGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCCGTTCGATAATACCGCTGGGACCTTGtgttt                                                                 >  2:442282/1‑228 (MQ=255)
                                                                                                                                               ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGCGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGTTATTACCGTTTTGACGTTGtgttt                                                                 >  1:461646/1‑228 (MQ=255)
                                                                                                                                                   tgcttgGGTCTGGGTGGCCTGGCGATCGGGCCCCGGGGTGTTAACGGAACTGGAAGGCAGCGCCGAACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgtt                                                              <  1:14363/225‑1 (MQ=255)
                                                                                                                                                               ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGTTGACGGAGAACTGCCCTACCAGCTCTTCGGTACTGCCATAACCCCGCTGTTTAACCAGTTCGCTAATACCCTTGTTACGTATTGTTTTTTTCCTTTTTAAATc                                                 >  2:182105/1‑228 (MQ=255)
                                                                                                                                                               ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAACTGCTCTACCAGCTCTTCGGTACTTACATAACCCTGCTGTTTAACCAGTTCGACAATACCGTTGTGACGTTGTGTTTGTTTCATTTTTaaaac                                                 >  1:638548/1‑226 (MQ=255)
                                                                                                                                                                     cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTTCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGGCATAACCCTGCTGTTTAACCAGTTCGGTCATACCGTTGTGACGTTGGTTTTGTTTCATTTATAAATCCCTGGa                                           >  1:654995/1‑228 (MQ=255)
                                                                                                                                                                                         tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTCCTGCAATAACCCTGCCGTTTAACCCGTTCGATATACCCCTTTTTACGTTGTGTTTGTTTCCTTTATATATCCCTGGAATTTttttcgttttccggca                       >  2:34552/1‑223 (MQ=255)
                                                                                                                                                                                         tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTTTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCgcgcc                       >  1:444955/1‑227 (MQ=255)
                                                                                                                                                                                                 gAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCCTTGAGCg                >  2:374195/1‑227 (MQ=255)
                                                                                                                                                                                                            cAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAGTTATTTTCGTTTTCGCGCATTGAGCGCATCAACaaaa     >  1:583555/1‑227 (MQ=255)
                                                                                                                                                                                                               cgccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATCGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc  <  1:442282/227‑1 (MQ=255)
                                                                                                                                                                                                               cgccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                                                                                         >  1:403243/1‑108 (MQ=255)
                                                                                                                                                                                                                gccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                                                                                        <  2:403243/108‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                              
AGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCC  >  NC_000913/3560229‑3560661

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGC  >  NC_000913/3560216‑3560696
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
CCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGG                                                                                                                                                                                                                                         >  M01186:129:000000000‑A5M8V:1:1105:15160:10660/1‑251 (MQ=60)
  ATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGGGCCGCACCGCCATGAGGGCGCGGGATCAGGTT                                                                                                                                                                                                                                       >  M01186:129:000000000‑A5M8V:1:2103:22336:16605/1‑251 (MQ=60)
      ACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATAGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGGGCCGCACGCCCATGCTGGCGCAGGATCAGGTTTTG                                                                                                                                                                                                                                    >  M01186:129:000000000‑A5M8V:1:2105:6090:15191/1‑250 (MQ=60)
         TGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTGCTC                                                                                                                                                                                                                                 >  M01186:129:000000000‑A5M8V:1:2112:22449:17920/1‑250 (MQ=60)
                  CGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACGGCCATGATGGCgcaggagcgggtctggctcttatacaca                                                                                                                                                                                                                        >  M01186:129:000000000‑A5M8V:1:1103:22447:26710/1‑222 (MQ=60)
                   GTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACGGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAT                                                                                                                                                                                                                      >  M01186:129:000000000‑A5M8V:1:2107:18873:10806/1‑251 (MQ=60)
                          tgtgtggcgcaaagcgcaaattgtggttgattgacaaggcggggccgtcccgcttccgggtgggtgcgaattccgaaaaccgccgtggggcttTTGGGGTTTTGCTCCGCCATTTTGCGGGCTATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGT                                                                                                                                                                                                                <  M01186:129:000000000‑A5M8V:1:1104:19257:27963/157‑1 (MQ=60)
                                GTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCCTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCG                                                                                                                                                                                                         >  M01186:129:000000000‑A5M8V:1:2107:15981:22070/1‑251 (MQ=60)
                                      gtgataaagggagcgggctgggggagtgcgtcggctcgcgctgcgggtgggggcgcgatacgtATAAACAGCGTCGATCTATTGGGAATTGGCGCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCTTTGCGATCGTGCCACGGCGTGTTAACGGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAG                                                                                                                                                                                                    <  M01186:129:000000000‑A5M8V:1:2103:11957:2169/187‑1 (MQ=60)
                                       tgcgaaatgggctggcgagaggctggggcgctgttccctttgcgggtgggtgccctatcgcgtagaacagcgccggcctttgggggttttgcgccgcccgttcgcgggcgagtcgccttgTTTCTCCGGTCTGGGGGCCCTTGCGATCGGCCCACGGGTTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTC                                                                                                                                                                                                  <  M01186:129:000000000‑A5M8V:1:1114:3560:13528/131‑1 (MQ=60)
                                                    GTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGG                                                                                                                                                                                     >  M01186:129:000000000‑A5M8V:1:1107:17111:9804/1‑251 (MQ=60)
                                                     TGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGA                                                                                                                                                                                    >  M01186:129:000000000‑A5M8V:1:2113:17862:25279/1‑251 (MQ=60)
                                                                  tgtgcaggtggtgtgtgcgtgtgacagtggttgtcgtttgtttccaagaaaaagccggcatacggagtcctcaggggtctcgggggctcggagaggtgtataagagccaGGCCTTGCGATCGTGCCAGGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACC                                                                                                                                                                       <  M01186:129:000000000‑A5M8V:1:1104:16611:28267/142‑1 (MQ=60)
                                                                              CTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCGGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTTACGGGGAAATccctgccccttttacacat                                                                                                                                                            >  M01186:129:000000000‑A5M8V:1:2101:15764:13457/1‑231 (MQ=60)
                                                                               tccgggcgggggcccataacggaaaacagtgtgcggcccgttggggtttggccCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCGTCGGTCGGGGGGGCCTTGCGATCGGGCCACGGGGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAC                                                                                                                                                          <  M01186:129:000000000‑A5M8V:1:2101:21005:21744/198‑1 (MQ=60)
                                                                                           gcggatttgggtaaacggcgtccggccttgGGGGATTTGCCCCGCCACGTGGCGGGCGATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCC                                                                                                                                               <  M01186:129:000000000‑A5M8V:1:1108:12969:7705/220‑1 (MQ=60)
                                                                                                 ATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAGGTCCCGGCGAATAGTCTGCGGGCTGACGCGGAAATGCTCTcccaccccttggggcctcacctaaccctgctgtt                                                                                                                                        >  M01186:129:000000000‑A5M8V:1:1109:26202:11144/1‑217 (MQ=60)
                                                                                                  TCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGGGAAATCCTCTACCtgtctcttttacccatctccgagcccacgaga                                                                                                                                       >  M01186:129:000000000‑A5M8V:1:1102:11374:23511/1‑219 (MQ=60)
                                                                                                      aaaacagggtcgtccttttggagtttgGCTCCGCCCCTTTGCGGGCGATGCGCTTTTTTTCTTCGGTCTGGGTGGCCTTGCGACCGTGCCAGGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACC                                                                                                                                   <  M01186:129:000000000‑A5M8V:1:2106:24773:6754/224‑1 (MQ=60)
                                                                                                                 aatcttttgggtattggctccgcccattggcgggcgatgggctttttttcctggggcgGGGTGGCCTGGCGATGGTGCCACGGGTTGTTAACCGAACGGGAAGGCAGCCCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA                                                                                                                         <  M01186:129:000000000‑A5M8V:1:2109:18094:10848/192‑1 (MQ=60)
                                                                                                                      ttgggggatctgctcagccctgttgcgggaggtgagcactttttctcaggtttggtgggtctggcaattgtgcgcgGGCGTGTTAACCTAACTGGAAGCCCGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCG                                                                                                                    <  M01186:129:000000000‑A5M8V:1:2110:21056:27899/174‑1 (MQ=60)
                                                                                                                       TTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTACGGTCGCGGCGATTAGTCTGCGGCCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCGGTTTAACCAGTTCGATtactttctt                                                                                                                  >  M01186:129:000000000‑A5M8V:1:1114:25171:22451/1‑242 (MQ=60)
                                                                                                                             ttttgccccgacctttgggggcggagtcgcttttttccttggggCTGGGTGGCCTTGCGATCGGGCCCCGGGGTTGTAACCGAACTGTAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACG                                                                                                            <  M01186:129:000000000‑A5M8V:1:1104:10072:10955/207‑1 (MQ=60)
                                                                                                                                TGCCCCGCCTCTTGGCGGGCGATGCGCTGTTTTTCTCGGGTCTGGGTGCCCGTGCGATCGTGCCACGGTGTGTTAGCCGAACTGGAGGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTG                                                                                                         <  M01186:129:000000000‑A5M8V:1:1104:4313:15012/251‑1 (MQ=60)
                                                                                                                                                 GGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGCGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGTTATTACCGTTTTGACGTTGTGTTTGTTTCATTTATA                                                                                        >  M01186:129:000000000‑A5M8V:1:2104:27177:20604/1‑251 (MQ=60)
                                                                                                                                                    GAGGCGCTTTTTTGCTTGGGTCTGGGTGGCCTGGCGATCGGGCCCCGGGGTGTTAACGGAACTGGAAGGCAGCGCCGAACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAA                                                                                      <  M01186:129:000000000‑A5M8V:1:1101:21765:13045/250‑1 (MQ=60)
                                                                                                                                                                 TCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAACTGCTCTACCAGCTCTTCGGTACTTACATAACCCTGCTGTTTAACCAGTTCGACAATACCGTTGTGACGTTGTGTTTGTTTCATTTTTAAAACCCTGGCCTTATT                                                                        >  M01186:129:000000000‑A5M8V:1:2110:17341:26251/1‑251 (MQ=60)
                                                                                                                                                                       GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTTCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGGCATAACCCTGCTGTTTAACCAGTTCGGTCATACCGTTGTGACGTTGGTTTTGTTTCATTTATAAATCCCTGGAGTTATTTTCGTT                                                                  >  M01186:129:000000000‑A5M8V:1:2111:17695:19741/1‑251 (MQ=60)
                                                                                                                                                                                           GTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTTTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCcttgacgcattca                                              >  M01186:129:000000000‑A5M8V:1:2104:9655:9993/1‑238 (MQ=60)
                                                                                                                                                                                                           ACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGATCTAACAGATCTTCGGTACTGACATAACCCCGCTGTTTAACCAGTTCGATAATACCGTTAGGAAGTTGGGTTTGGTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGAGCCTTGGGCGAATAAACCAAAGCCAAGGCTA                              >  M01186:129:000000000‑A5M8V:1:1103:26641:9636/1‑251 (MQ=60)
                                                                                                                                                                                                              GAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAGTTATTTTCGTTTTCGCGCATTGAGCGCATCAACAAAAGCCATCGCTAAA                            >  M01186:129:000000000‑A5M8V:1:2108:7830:21404/1‑250 (MQ=60)
                                                                                                                                                                                                                cgTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATCGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACC                          <  M01186:129:000000000‑A5M8V:1:2104:24041:8116/248‑1 (MQ=60)
                                                                                                                                                                                                                 CTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACC                                                                                                                                                >  M01186:129:000000000‑A5M8V:1:2102:10355:25379/1‑131 (MQ=60)
                                                                                                                                                                                                                                     GCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGgat    >  M01186:129:000000000‑A5M8V:1:2105:5814:16970/1‑248 (MQ=60)
                                                                                                                                                                                                                                       CATGATGGCGCAGGGTCAGGTTTTACGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGC  <  M01186:129:000000000‑A5M8V:1:1103:16799:4762/251‑1 (MQ=60)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
CCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGC  >  NC_000913/3560216‑3560696

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