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
A0 F0 I1 R1 19 24.3 471414 95.8% 451614 257.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,543 C→A P1100Q (CCG→CAG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5430CA91.7% 34.2 / NA 12P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (6/5);  total (6/6)
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

TCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCAAAGGTATGCCAATCGCAACGCCGGTGTTCGACGGTGC  >  NC_000913/4184318‑4184813
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                              
tCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                                                                                                                                                               <  2:128161/275‑1 (MQ=255)
       gaaACTCGAAGCGAAACGCCGCAAAATAACCCAGGGCGACGATCTGCCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTa                                                                                                                                                                                                                        <  2:81307/275‑1 (MQ=255)
                                                           ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagcag                                                                                                                                                                      >  1:197377/1‑273 (MQ=255)
                                                                              aGGTGTCTCTTGCGGGTAAACGCCGGATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCTGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcg                                                                                                                                                  <  2:53785/274‑1 (MQ=255)
                                                                                 tataTCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa                                                                                                                                              <  2:192579/275‑1 (MQ=255)
                                                                                           ggTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgc                                                                                                                                    >  2:31658/1‑275 (MQ=255)
                                                                                                aatgccttactcAGCCTGGTGCCAAGAGGGGGGGTCGTGTTGGTAACAGGGGTGTAATTTCTAAGAGCAACACGAGCGAAGATAAGCCTTACGATGAAAACGGCGTGCCGGTGGACATCGCACCGAACCAGCGGAGGGGACCGTCCGGTATCACCAGCGGTCGGATCCTCGGGACCCACCTGGGTAGGGCTGCGAAAGGTATCGGCGACGAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAGCTGGGCCGATTCATCGGgcgtgcg                                                                                                                               <  2:205220/264‑1 (MQ=255)
                                                                                                                                                                     aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGAATCGGGCGCTGACGTTCGGCAGACAGTTGCACCGAGTCCCTTCCGCGAGGGAGGAGCTATg                                                             >  2:202747/1‑272 (MQ=255)
                                                                                                                                                                                         tACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGAGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCGGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGAGTGCGTACGAGCCGGGCGCTGACGTTCGGCACACAGTTGACCTGCGTACCTTCAGCGATGAAGAAGTTATGCGTCTGTCTGCACACCTgcgaa                                       >  1:210751/1‑272 (MQ=255)
                                                                                                                                                                                          aCGATGAAAACGTTACGCCGGTAGACATCGTACTGAACCCGATGGGTGTACCGTTTCTTAGGAAAGTCGGTCAGATCCTGGAAACCCACCTGGGTAGGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCTCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCaaa                                     <  1:112573/275‑1 (MQ=255)
                                                                                                                                                                                                                 gACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCAAAGGTATGCCAATCGCAACGCCGgt              >  2:187468/1‑275 (MQ=255)
                                                                                                                                                                                                                             aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCTCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACATGAGTCCCTCAAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCCAAGGTATGCCAATCGCACCGCCGGTGTTCGACGGTGc  >  1:117630/1‑275 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                              
TCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCAAAGGTATGCCAATCGCAACGCCGGTGTTCGACGGTGC  >  NC_000913/4184318‑4184813

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 11 ≤ 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
ACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCAAAGGTATGCCAATCGCAACGCCGGTGTTCGACGGTGCGAAAGAAGCAGAAAT  >  NC_000913/4184280‑4184828
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                             
ACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCCACCCGATCGAAGGTATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAACGTCTCGTATcaccgtctcttatcca                                                                                                                                                                                                                                                          >  M01186:109:000000000‑A8RG1:1:1108:27142:12570/1‑285 (MQ=60)
         atagagagaagaagaacgaacgcacgaggtcggcaacaaacccgaatccaaacggcacaaaatctccgcgggtGACGATCTGGCACCGGGCCTGCTTAAGATTGTTAAGGTATGTCTGTCGGTCAATCGCCGTATCCAGCCTTGTGACAAGGTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCGAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACCTCGGTCAGATCCTCGAAAC                                                                                                                                                                                                                                                  <  M01186:109:000000000‑A8RG1:1:2101:17478:14266/227‑1 (MQ=60)
           tacaggaagagtaaaagaaaaacatgtgcgaaaatacacgcgaagagagatgccgtaaagcctcgccgggcgtagatgcggcgCCGGGCGTGGTGAAGATGGAGAAGGTATATCTGGCGGTTAAACGCCTGATCCAGCCTGGTGACAATATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGAGGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACC                                                                                                                                                                                                                                                 <  M01186:109:000000000‑A8RG1:1:2106:17475:1751/216‑1 (MQ=60)
                  tgatgaacggatacacgagtctcgaccgaaacaccagacgaaacgcgccaaaataacccctggaggcgagctgggccctggcaggctgacgatggttaaagtattaCTGGCAGTTAAACGCCGTATCCATCCTGGTGACATGATGGCATGTCGTCACGGTAACAAGGGGGTAATTTCTACGATCCACCCGATCGAAGATACGCCTTACGATGAAAACGGTACTCCGGTAGACATCGTACTGACCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAATCCCACCTGGGT                                                                                                                                                                                                                                        <  M01186:109:000000000‑A8RG1:1:2109:10157:3661/195‑1 (MQ=60)
                         ATGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAACTCCCCCAGGGCGACGATCTGGCACCGTGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAAGTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGTGCGTACCG                                                                                                                                                                                                                                                                                  <  M01186:109:000000000‑A8RG1:1:2103:9434:4488/252‑1 (MQ=60)
                                                               CCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCTCGGTAAAAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCAGTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACAAGCTGGGCTTACCGTCGAGTATGAACAGCGGTGAGAACCGCGAAACCCACCTGGGTATGGCTGagacaggagtcggcggaaagatcaacggcaagctgcaa                                                                                                                                                                                           >  M01186:109:000000000‑A8RG1:1:2109:17509:7495/1‑263 (MQ=60)
                                                                                      ATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGCAGTCGCGAAA                                                                                                                                                                      >  M01186:109:000000000‑A8RG1:1:2109:27080:6296/1‑299 (MQ=60)
                                                                                      ATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCGAAGATCAACCCGATCGAAGATATGCCTTACGATGAAACCGGTAAGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCCGATCCTCGAACCCCACCTGGGTACGGTTGCGGAAGGTAATGGCGACAAGAACAACGCAATGGTGAATACGCAGCAAGAAGTCGCGAAACA                                                                                                                                                                    >  M01186:109:000000000‑A8RG1:1:2112:19652:15970/1‑301 (MQ=60)
                                                                                          GGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCAAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAAAAAGGGAGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATCAAAACGGTACGCCGGTCGACAACGTACTGCAACAGCTGGGCGTACCGTCTCGTATGAACAACGGTCAGCTCCTCGAAAACCACCTGGGTCTGACGGCAACAGGTATCGGCGACAAGATCCACGCCctactgcaccaacaacaagaagacgacaaaatgac                                                                                                                                                                 >  M01186:109:000000000‑A8RG1:1:1101:5931:20850/1‑265 (MQ=60)
                                                                                                                                                                              CGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGATAAGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAGGAAAATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTAATCCAGCGTGAGTACGATCTGGGCGCTGACGTTCGTCAGACACTTGACCTGAGTACCTTCAGCGATGCAGAAGTTATA                                                                            >  M01186:109:000000000‑A8RG1:1:2109:17106:5726/1‑301 (MQ=60)
                                                                                                                                                                                                      GATCAACCCGATAGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGCCACGATACACGCCAGGCTGAAACAGCAGCACGAAGTCGCGCAAATGCGCGACTTCATCCAGCGTGCGGACGAACTGGCCGCTGACGTTCGTCAGAAAGTTGACCTGAGTCCCTTCAGCGActacgaagcaatgcgtctagctgaacacatggccaca                                                    >  M01186:109:000000000‑A8RG1:1:1112:13958:5302/1‑264 (MQ=60)
                                                                                                                                                                                                                 TAGAAGATATGCCTTACGATGAAAACGTTACGCCGGTAGACATCGTACTGAACCCGATGGGTGTACCGTTTCTTAGGAAAGTCGGTCAGATCCTGGAAACCCACCTGGGTAGGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCTCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCAAAGGTATGCCAAT                                         <  M01186:109:000000000‑A8RG1:1:1119:12614:21950/301‑1 (MQ=60)
                                                                                                                                                                                                                    AAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGAGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCGGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGAGTGCGTACGAGCCGGGCGCTGACGTTCGGCACACAGTTGACCTGCGTACCTTCAGCGATGAAGAAGTTATGCGTCTGTCTGCACACCTGCGAAAAGGTAGGCCAATCG                                       >  M01186:109:000000000‑A8RG1:1:2111:12220:12692/1‑300 (MQ=60)
                                                                                                                                                                                                                                              ccgccttgAGACATCGTACTGCACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCTAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCAACAAGAGCAACGCCATGCTGAAGCAGCAGCAAGAAGTCGCGAAGCTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCAAAGGTATGCCAATCGCAACGCCGGTGTTCGACGGTGCGAAAG            <  M01186:109:000000000‑A8RG1:1:1104:28914:10397/293‑1 (MQ=60)
                                                                                                                                                                                                                                                        ACATCGTACTGACCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCCCCTGGTTATGGCTGCGAAAGGTATCGGCGACAAGACCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCAAAGGTATGCCAATCGCAACGCCGGTGTTCGACGGTGCGAAAGAAGCAGAAA   <  M01186:109:000000000‑A8RG1:1:1103:12447:17750/300‑1 (MQ=60)
                                                                                                                                                                                                                                                        ACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCTCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACATGAGTCCCTCAAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCCAAGGTATGCCAATCGCACCGCCGGTGTTCGACGGTGCGAcagacacagaaaa  >  M01186:109:000000000‑A8RG1:1:2101:11424:11196/1‑288 (MQ=60)
                                                                                                                                                                                                                                                             GTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAACCCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCAAAGG                                                  <  M01186:109:000000000‑A8RG1:1:2102:11916:15554/248‑1 (MQ=60)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                             
ACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCAAAGGTATGCCAATCGCAACGCCGGTGTTCGACGGTGCGAAAGAAGCAGAAAT  >  NC_000913/4184280‑4184828

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑