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 F27 I0 R1 60 79.9 2053078 82.8% 1699948 215.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,305 C→T 100% R395C (CGT→TGT)  icd → isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,3050CT99.7% 196.1 / ‑7.1 63R395C (CGT→TGT) icdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/26);  minor base G (2/0);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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.

GGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCC  >  NC_000913/1196070‑1196534
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     
ggTTGGCGGGATCGGGATCGCCCCTGGTGCAAACAGCGGTGACGAAGGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTAGTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAgtgt                                                                                                                                                                                                                                     <  1:741615/238‑1 (MQ=255)
        gtatcggtatcgCCCCTGGTGCAAACATCGGTTACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgatgga                                                                                                                                                                                                                             <  1:658456/238‑1 (MQ=255)
        gtatcggtatcgCCCATGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCGCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATTGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgatgga                                                                                                                                                                                                                             <  1:74356/238‑1 (MQ=255)
                atcgCCCCTGGTGCAAACCTCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTa                                                                                                                                                                                                                     <  2:231905/238‑1 (MQ=255)
                                aCATCGGTGAAGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCa                                                                                                                                                                                                    <  1:436931/239‑1 (MQ=255)
                                     gggtacgtaagcgcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGtt                                                                                                                                                                                               <  1:237769/229‑1 (MQ=255)
                                     ggTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAACCAACGCGAAAAACGTAACCTATGACTTCGAGGGTCTGGTGGATGGCGCTAAACTGCTGAAAATTTCAGGGtt                                                                                                                                                                                               >  2:936018/1‑239 (MQ=255)
                                               gcgcCCTGTTTGAAGTCCCCCACGGTTCTTCGCTGAAATAGGCCGGTCAGGACAAGGTAAATCCTGGTCCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgcg                                                                                                                                                                                     <  2:296763/239‑1 (MQ=25)
                                                     tGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGgcg                                                                                                                                                                                                                        >  1:751926/1‑198 (MQ=255)
                                                     tGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGgcg                                                                                                                                                                                                                        <  2:751926/198‑1 (MQ=255)
                                                               cacccacGGTACGGCGCCGAAAGATGCCGGTCAAGACAAATGAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGAGGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                                                                                                                                      <  1:108768/238‑1 (MQ=255)
                                                                             cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTg                                                                                                                                                       >  1:262455/1‑239 (MQ=255)
                                                                                  aaaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACTTATGCCTTCTAGTGTCGGATGGATGGCGATAAACTGCTGAAATGTTCAAAGTTAGGTGACGCGATACTCGAAAAAATTGAAACACGGAGTGTGTTaaa                                                                                                                                                  >  2:1025647/1‑239 (MQ=21)
                                                                                  aaaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGt                                                                                                                                                          <  1:572405/231‑1 (MQ=255)
                                                                                  aaaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGt                                                                                                                                                          >  2:572405/1‑231 (MQ=255)
                                                                                  aaaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGAATGGCAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGGCCTAAACTGCGGAAATGTTCACAGTTTGGTGCCGGGACCATCGAAAACATGTAATTCCGAAGTTTTTTaaa                                                                                                                                                  >  2:65020/1‑239 (MQ=255)
                                                                                              cAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACgg                                                                                                                                       <  1:152839/238‑1 (MQ=255)
                                                                                                      aGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTCGTTTGTTAATTTTATTAACgggggggaaac                                                                                                                              >  2:558116/1‑231 (MQ=255)
                                                                                                      aGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTTTTAACGGGAGCGTAAc                                                                                                                              >  1:256343/1‑239 (MQ=255)
                                                                                                           aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACGTGGGTTGGACAGAAGCGGCTGACTTAATTGTTAAAGGAATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAACTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAGTTTATTAACGGGAGCGTAACGCTc                                                                                                                          >  1:591062/1‑238 (MQ=255)
                                                                                                           aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGaa                                                                                                                                                                                                           >  1:997537/1‑157 (MQ=255)
                                                                                                           aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTc                                                                                                                          >  1:638169/1‑238 (MQ=255)
                                                                                                           aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTc                                                                                                                          >  1:77704/1‑238 (MQ=255)
                                                                                                                ggCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGGATGGAAGGCCCAATCAAAGCGAAAACCGTAACCTAGGACTTCGAGGGTCTGATGGATGGAGCTAAACTGCTGAAATGTTCAAAGTTTGGTGACGCGATCCTCGAAAACATGTAATCCCGCAGTTTGTTAAATTTATTAACGGAAGGGTAACGCTCCCcttg                                                                                                                    >  1:407319/1‑239 (MQ=255)
                                                                                                                 gCTCTATTATTCTCTCCGCTGAGCTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGCGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAAAGTTTGGTGACGCGATCATCGAAAACATGTAATTCCGTAGTTTGTTAAATTTATTAACGGGCGCGTAACGCTACCgatgt                                                                                                                   >  1:703117/1‑239 (MQ=25)
                                                                                                                           tttctcCGCGGAGATTATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAg                                                                                                         <  2:77704/237‑1 (MQ=255)
                                                                                                                            ctctcGGCTGAGATTATGGTGCGGACAATGGGTTGGCCAGAAGGGGCTGACTTAATTGTTAAAGGTATGGAAGGGGCAATCAACGCGAAAACCGTAACCTATGGCTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAgg                                                                                                        <  2:936473/239‑1 (MQ=255)
                                                                                                                               tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGCCTTAATTGTTACAGGTATGGAAGGCGCAATCAACGCGGAACCCGTAACCTATGACTTCAAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACCCGACACTCCAAAACATGTAATGACGGAGTTTGTTCAATTATTTAACGGGAGATTAACGCTCCCGTAGTtttttttttcgatg                                                                                                     >  2:992626/1‑233 (MQ=255)
                                                                                                                               tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGATGCGATCATCGAAAACATGCAATGCCGTAGTTTGTTAAATTTATTAACGGAAGCGTAACGCTCCCGTTGTTTTTTGTTAgggt                                                                                                      >  1:530950/1‑236 (MQ=255)
                                                                                                                               tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgct                                                                                                      >  1:391205/1‑238 (MQ=255)
                                                                                                                                          gatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAGCCGTAACCTATGGCTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTAt                                                                                          <  2:228068/239‑1 (MQ=255)
                                                                                                                                           atgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACgg                                                                                              <  1:99300/234‑1 (MQ=255)
                                                                                                                                           atgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACgg                                                                                              >  2:99300/1‑234 (MQ=255)
                                                                                                                                           atgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAAAATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGGGTAAAGCTCCCGTTGTTTTTTGTTAGGTTGTTAAGGGTTaac                                                                                         >  2:580851/1‑237 (MQ=255)
                                                                                                                                             gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGCTGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGAAATGCCGTAGTTTGTTAAATTTATTAATGGGAGCGTAACGCTCCCGTTGTTTTTGTTTGGCCGGGTAACGGTAATCaa                                                                                       >  1:179992/1‑239 (MQ=21)
                                                                                                                                                           gTTGGAGCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGCAAACCGTAACCTATGACTTCGAGTGTCTGATGGATAGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGCACACATGCAATGCCGTCGTTTGTTACATTTATTACCGGGGGCGTCAAGCTCCCGCTGTTTTTTGTTAGGCTGCTAAAGGGTTTAAAAAATTTTCCaaata                                                                         >  2:191849/1‑237 (MQ=21)
                                                                                                                                                              ggACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAg                                                                      >  2:177908/1‑239 (MQ=255)
                                                                                                                                                              ggACCGAAGCGGCTGACTTAATTGGTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaaaaaa                                                                       <  1:558116/238‑1 (MQ=255)
                                                                                                                                                                     aGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGGGAGTCTGATGGATGGCGCTAAACTGCTGAACTGTTCAAAGGTTGGTGACGCGATCATCGAAACCATGTATTGCCGTAGTTTGTTAAATTTCTTATAGGGTGCGTAACGGCCCCGTTGTTTTTTGTTATGCTGCTAACGGTTATACAAATTTAATAAAAAAAAGtttttaa                                                               >  2:888263/1‑234 (MQ=21)
                                                                                                                                                                     aGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCaa                                                               >  1:220668/1‑239 (MQ=255)
                                                                                                                                                                     aGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGGGACGCGATCATCGAAAACCTGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCaa                                                               >  2:205011/1‑239 (MQ=255)
                                                                                                                                                                         gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAc                                                            >  1:98758/1‑238 (MQ=255)
                                                                                                                                                                         gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAc                                                            >  1:225887/1‑238 (MQ=255)
                                                                                                                                                                              cTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTTTCAAAATTTTTTCAAAAAAAGTTATAAAAACCCCACg                                                      >  2:568347/1‑239 (MQ=255)
                                                                                                                                                                               tttaTTTTTAAAGGTTTGGGAGGGCGAATAAACGCGAAAACCGTATCCTATGGCTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCCGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCgg                                                     <  2:220668/236‑1 (MQ=37)
                                                                                                                                                                               ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCgg                                                     <  1:580851/239‑1 (MQ=255)
                                                                                                                                                                                                ggAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAATCCCCTCGGTAGtt                                                <  2:524324/227‑1 (MQ=255)
                                                                                                                                                                                                ggAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGtt                                                >  1:524324/1‑227 (MQ=255)
                                                                                                                                                                                                    ggCGCAAACAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGGTAAACTGCTGAAATGTTCAGAGTTTGGTTACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGc                                 <  2:638169/238‑1 (MQ=255)
                                                                                                                                                                                                      cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCgg                              >  2:984833/1‑239 (MQ=255)
                                                                                                                                                                                                       gcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTAATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCgg                              >  1:108554/1‑238 (MQ=255)
                                                                                                                                                                                                        cAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGGCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGt                             <  1:568347/238‑1 (MQ=255)
                                                                                                                                                                                                          aTCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAg                          >  1:41993/1‑239 (MQ=255)
                                                                                                                                                                                                            cAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAggt                        <  1:177908/239‑1 (MQ=255)
                                                                                                                                                                                                             aaCGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAggtgg                      <  1:1025647/239‑1 (MQ=255)
                                                                                                                                                                                                               cgcgAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAATTTTTATCAAAAAAAGTTATCAAAACCCCTCgg                                                     >  2:313324/1‑207 (MQ=255)
                                                                                                                                                                                                               cgcgAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCgg                                                     <  1:313324/207‑1 (MQ=255)
                                                                                                                                                                                                               cgcgAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAggtggt                     >  2:645379/1‑239 (MQ=255)
                                                                                                                                                                                                                         gTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGAAGTTttacta           >  1:241107/1‑239 (MQ=255)
                                                                                                                                                                                                                              cTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGt      <  1:469069/239‑1 (MQ=255)
                                                                                                                                                                                                                                  taCTTCGAGTGTCGGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAATAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTcc  <  1:191849/238‑1 (MQ=255)
                                                                                                                                                                                                                                  gACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTTGGGTAGGCTGGCCGGTCAGGCGGTAGTTCTACTACTAGtctc   >  1:175816/1‑238 (MQ=255)
                                                                                                                                                                                                                                       cGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGt                             <  2:89388/207‑1 (MQ=255)
                                                                                                                                                                                                                                       cGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGt                             >  1:89388/1‑207 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     
GGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCC  >  NC_000913/1196070‑1196534

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGC  >  NC_000913/1196069‑1196553
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
AGGTTGGCGGGATCGGGATCGCCCCTGGTGCAAACAGCGGTGACGAAGGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTAGTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGG                                                                                                                                                                                                                                             <  M01186:129:000000000‑A5M8V:1:2107:9778:6154/250‑1 (MQ=60)
        GGTATCGGTATCGCCCCTGGTGCAAACATCGGTTACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAAC                                                                                                                                                                                                                                     <  M01186:129:000000000‑A5M8V:1:2104:14541:21330/250‑1 (MQ=60)
        GGTATCGGTATCGCCCATGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCGCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATTGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAAC                                                                                                                                                                                                                                     <  M01186:129:000000000‑A5M8V:1:1102:11978:26415/250‑1 (MQ=60)
                                AACATCGGTGAAGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGA                                                                                                                                                                                                            <  M01186:129:000000000‑A5M8V:1:1112:14276:14618/251‑1 (MQ=60)
                                     tgggtacgtaaGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGA                                                                                                                                                                                                       <  M01186:129:000000000‑A5M8V:1:1107:14734:9572/240‑1 (MQ=60)
                                           CGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGctggcccttttacacatttccgagcccacga                                                                                                                                                                                                 >  M01186:129:000000000‑A5M8V:1:2107:27719:12952/1‑220 (MQ=60)
                                                               CCACCCACGGTACGGCGCCGAAAGATGCCGGTCAAGACAAATGAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGAGGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAG                                                                                                                                                                              <  M01186:129:000000000‑A5M8V:1:1103:25148:23313/250‑1 (MQ=60)
                                                                   CCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGT                                                                                                                                                                         >  M01186:129:000000000‑A5M8V:1:1107:23329:24337/1‑251 (MQ=60)
                                                                          TCTGCCCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTT                                                                                                                                                                  <  M01186:129:000000000‑A5M8V:1:2102:13560:11296/251‑1 (MQ=60)
                                                                                            GGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTTTTAACGGGAGCGTAACG                                                                                                                                                >  M01186:129:000000000‑A5M8V:1:1107:19788:20111/1‑251 (MQ=60)
                                                                                              TCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGC                                                                                                                                               <  M01186:129:000000000‑A5M8V:1:1105:8808:3686/250‑1 (MQ=60)
                                                                                                 GGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACGTGGGTTGGACAGAAGCGGCTGACTTAATTGTTAAAGGAATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAACTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAGTTTATTAACGGGAGCGTAACGCTCC                                                                                                                                            >  M01186:129:000000000‑A5M8V:1:2102:14605:22720/1‑250 (MQ=60)
                                                                                                 GGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCC                                                                                                                                            >  M01186:129:000000000‑A5M8V:1:2104:18033:9393/1‑250 (MQ=60)
                                                                                                 GGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCC                                                                                                                                            >  M01186:129:000000000‑A5M8V:1:1103:20970:4646/1‑250 (MQ=60)
                                                                                                 GGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAA                                                                                                                                                                                                                             >  M01186:129:000000000‑A5M8V:1:2114:16618:8580/1‑169 (MQ=60)
                                                                                                      AAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGGATGGAAGGCCCAATCAAAGCGAAAACCGTAACCTAGGACTTCGAGGGTCTGATGGATGGAGCTAAACTGCTGAAATGTTCAAAGTTTGGTGACGCGATCCTCGAAAACATGTAATCCCGCAGTTTGTTAAATTTATTAACGGAAGGGTAACGCTCCCCTTGT                                                                                                                                      >  M01186:129:000000000‑A5M8V:1:1111:16751:18931/1‑251 (MQ=60)
                                                                                                       AGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGCTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGCGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAAAGTTTGGTGACGCGATCATCGAAAACATGTAATTCCGTAGTTTGTTAAATTTATTAACGGGCGCGTAACGCTACCGATGTT                                                                                                                                     >  M01186:129:000000000‑A5M8V:1:2106:12045:3254/1‑251 (MQ=60)
                                                                                                                     CTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGATGCGATCATCGAAAACATGCAATGCCGTAGTTTGTTAAATTTATTAACGGAAGCGTAACGCTCCCGTTGTTTTTTGTTAGGgtc                                                                                                                        >  M01186:129:000000000‑A5M8V:1:2101:13262:6392/1‑247 (MQ=60)
                                                                                                                     CTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTG                                                                                                                        >  M01186:129:000000000‑A5M8V:1:1111:20543:10001/1‑250 (MQ=60)
                                                                                                                                   GCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGCTGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGAAATGCCGTAGTTTGTTAAATTTATTAATGGGAGCGTAACGCTCCCGTTGTTTTTgtttggccgggtaacggtaatcaaa                                                                                                         >  M01186:129:000000000‑A5M8V:1:1105:8395:19750/1‑226 (MQ=60)
                                                                                                                                      TAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATT                                                                                                      <  M01186:129:000000000‑A5M8V:1:1103:12970:17455/251‑1 (MQ=60)
                                                                                                                                                           GGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAA                                                                                 >  M01186:129:000000000‑A5M8V:1:1106:15971:21081/1‑251 (MQ=60)
                                                                                                                                                              TGGACCGAAGCGGCTGACTTAATTGGTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAC                                                                               <  M01186:129:000000000‑A5M8V:1:2101:16152:23871/250‑1 (MQ=60)
                                                                                                                                                               GGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACC                                                                              >  M01186:129:000000000‑A5M8V:1:1106:23123:24703/1‑250 (MQ=60)
                                                                                                                                                               GGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACC                                                                              >  M01186:129:000000000‑A5M8V:1:1103:23946:17152/1‑250 (MQ=60)
                                                                                                                                                                               ATTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGG                                                             <  M01186:129:000000000‑A5M8V:1:2102:19905:16174/251‑1 (MQ=60)
                                                                                                                                                                                      GTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGC                                                        >  M01186:129:000000000‑A5M8V:1:1114:13743:24636/1‑249 (MQ=60)
                                                                                                                                                                                             GTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTAATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGT                                                >  M01186:129:000000000‑A5M8V:1:1103:22989:23161/1‑250 (MQ=60)
                                                                                                                                                                                                TGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGG                                            >  M01186:129:000000000‑A5M8V:1:1102:22139:5821/1‑251 (MQ=60)
                                                                                                                                                                                                     GGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGG                                                             <  M01186:129:000000000‑A5M8V:1:1109:23716:9418/229‑1 (MQ=60)
                                                                                                                                                                                                        GCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGGCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGT                                     <  M01186:129:000000000‑A5M8V:1:2102:14535:8792/250‑1 (MQ=60)
                                                                                                                                                                                                            TCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTAC                                <  M01186:129:000000000‑A5M8V:1:1105:16248:18599/251‑1 (MQ=60)
                                                                                                                                                                                                             CAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTA                              <  M01186:129:000000000‑A5M8V:1:2114:18070:27133/251‑1 (MQ=60)
                                                                                                                                                                                                               ACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGAAGTTTTACTAA                             >  M01186:129:000000000‑A5M8V:1:1107:22230:11540/1‑251 (MQ=60)
                                                                                                                                                                                                                        CCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTTGGGTAGGCTGGCCGGTCAGGCGGTAGTTCTACTACTAGTCTCC                     >  M01186:129:000000000‑A5M8V:1:1105:7446:17495/1‑250 (MQ=60)
                                                                                                                                                                                                                             ACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTctgtctcttatacccatttccg               >  M01186:129:000000000‑A5M8V:1:1103:16247:12134/1‑229 (MQ=60)
                                                                                                                                                                                                                              CCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAG              <  M01186:129:000000000‑A5M8V:1:1113:11004:12467/251‑1 (MQ=60)
                                                                                                                                                                                                                               CTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGA             >  M01186:129:000000000‑A5M8V:1:1103:27946:21476/1‑251 (MQ=60)
                                                                                                                                                                                                                                  TTACTTCGAGTGTCGGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAATAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATAT          <  M01186:129:000000000‑A5M8V:1:1106:16542:2944/251‑1 (MQ=60)
                                                                                                                                                                                                                                     CTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGACGTAGTTTGTTACATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAATTTATCAAAACCCCTCGTTAGTTTTGGGGTAGGCTGGCGGGGCAGGGGGTAGTTCTACTACTAGTCTCCCACCTAGAAATT         >  M01186:129:000000000‑A5M8V:1:1111:22806:20340/1‑249 (MQ=60)
                                                                                                                                                                                                                                           GTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGC  <  M01186:129:000000000‑A5M8V:1:2113:11065:22103/250‑1 (MQ=60)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
AGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGC  >  NC_000913/1196069‑1196553

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