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 F11 I1 R1 15 43.1 1622188 96.5% 1565411 127.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,820,792 T→C D407G (GAC→GGC)  alaS ← alanyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,820,7920TC86.0% 127.9 / 10.2 50D407G (GAC→GGC) alaSalanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base C (21/22);  total (25/25)
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

ATGCACTGTCAACACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTT  >  NC_000913/2820646‑2820924
                                                                                                                                                  |                                                                                                                                    
aTGCACTGTCAACACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCa                                                                                                                                    >  4:184663/1‑149 (MQ=255)
 tGCACTGTCAACACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCaa                                                                                                                                   >  5:66303/1‑149 (MQ=255)
           acacGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGcc                                                                                                                         <  5:106232/149‑1 (MQ=255)
             acGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCa                                                                                                                        >  3:43972/1‑148 (MQ=255)
               ggATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAg                                                                                                                     <  6:129760/149‑1 (MQ=255)
                    aTTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTc                                                                                                                >  5:51285/1‑149 (MQ=255)
                     ttGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCg                                                                                                               <  6:66303/149‑1 (MQ=255)
                        cGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTa                                                                                                             >  7:82792/1‑148 (MQ=255)
                         gTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAAc                                                                                                                                  <  6:36717/126‑1 (MQ=255)
                         gTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAAc                                                                                                                                  >  5:36717/1‑126 (MQ=255)
                                  ggCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGaaa                                                                                                   <  7:86492/148‑1 (MQ=255)
                                    cACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCa                                                                                                <  4:43972/149‑1 (MQ=255)
                                        aaaGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGtt                                                                                             >  6:24803/1‑148 (MQ=255)
                                        aaaGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGtt                                                                                             <  5:24803/148‑1 (MQ=255)
                                          aGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCa                                                                                          <  3:30917/149‑1 (MQ=255)
                                            ccGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAg                                                                                                                     <  2:210693/120‑1 (MQ=255)
                                            ccGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAg                                                                                                                     <  2:207491/120‑1 (MQ=255)
                                            ccGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAg                                                                                                                     >  1:210693/1‑120 (MQ=255)
                                            ccGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAg                                                                                                                     >  1:207491/1‑120 (MQ=255)
                                                    ttCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGc                                                                                 >  2:113977/1‑148 (MQ=255)
                                                      cgcgcgcgcgACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGc                                                                                 <  1:113977/146‑1 (MQ=255)
                                                       gcgcgcgcgACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTa                                                                                                             >  1:14036/1‑117 (MQ=255)
                                                       gcgcgcgcgACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTa                                                                                                             <  2:14036/117‑1 (MQ=255)
                                                          cgcgcgACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCAc                                                                          <  2:212572/149‑1 (MQ=255)
                                                                     gctgctCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGCCGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGtttt                                                               >  3:138473/1‑149 (MQ=255)
                                                                     gctgctCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGtttt                                                               <  8:82792/149‑1 (MQ=255)
                                                                        gctCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGc                                                             >  2:62344/1‑148 (MQ=255)
                                                                            ttCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAg                                                          <  1:62344/147‑1 (MQ=255)
                                                                             tCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGctc                                                       >  2:163400/1‑149 (MQ=255)
                                                                                aTTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCtt                                                     >  8:42227/1‑148 (MQ=255)
                                                                                       cTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCATATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGc                                             <  3:184663/149‑1 (MQ=255)
                                                                                        ttCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGtttt                                                               >  6:76933/1‑130 (MQ=255)
                                                                                         tCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCaa                                           >  5:146191/1‑149 (MQ=255)
                                                                                            aaaCCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGc                                        <  1:163400/149‑1 (MQ=255)
                                                                                                   cTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCa                                                                                          >  6:35872/1‑92 (MQ=255)
                                                                                                   cTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCa                                                                                          <  5:35872/92‑1 (MQ=255)
                                                                                                        tCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTc                            <  6:146191/149‑1 (MQ=255)
                                                                                                        tCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTc                            >  4:6739/1‑149 (MQ=255)
                                                                                                            cTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCaga                        <  3:6739/149‑1 (MQ=255)
                                                                                                                       cgcTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCaaa             <  7:42227/149‑1 (MQ=255)
                                                                                                                        gcTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCaaa             >  4:120374/1‑148 (MQ=255)
                                                                                                                        gcTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCaaa             <  3:120374/148‑1 (MQ=255)
                                                                                                                        gcTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAAc            <  6:51285/149‑1 (MQ=255)
                                                                                                                          tCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGAACAGACGCAAAGAAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTg          >  7:54281/1‑149 (MQ=255)
                                                                                                                                   aCATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTcttctt  >  5:42912/1‑148 (MQ=255)
                                                                                                                                    cATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTcttctt  >  6:93927/1‑147 (MQ=255)
                                                                                                                                     aTCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGa                                    <  7:22192/112‑1 (MQ=255)
                                                                                                                                     aTCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGa                                    >  8:22192/1‑112 (MQ=255)
                                                                                                                                        aGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCa                          <  2:64807/119‑1 (MQ=255)
                                                                                                                                        aGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCa                          >  1:64807/1‑119 (MQ=255)
                                                                                                                                                  |                                                                                                                                    
ATGCACTGTCAACACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTT  >  NC_000913/2820646‑2820924

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATGCACTGTCAACACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTC  >  NC_000913/2820646‑2820925
                                                                                                                                                  |                                                                                                                                     
ATGCACTGTCAACACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAAC                                                                                                                                   >  NB501016:228:HFK2FAFXY:3:11511:20470:5500/1‑151 (MQ=60)
          AACACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCA                                                                                                                         <  NB501016:228:HFK2FAFXY:3:11612:11087:6696/151‑1 (MQ=60)
            CACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCAT                                                                                                                        >  NB501016:228:HFK2FAFXY:2:11206:22718:12382/1‑150 (MQ=60)
                   CATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCG                                                                                                                >  NB501016:228:HFK2FAFXY:3:11506:10058:5864/1‑151 (MQ=60)
                       GCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTAC                                                                                                             >  NB501016:228:HFK2FAFXY:4:11608:14932:14691/1‑150 (MQ=60)
                        CGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACC                                                                                                                                  >  NB501016:228:HFK2FAFXY:3:11501:10334:8143/1‑128 (MQ=60)
                                 CGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAG                                                                                                   <  NB501016:228:HFK2FAFXY:4:11610:11370:5780/150‑1 (MQ=60)
                                       CAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTT                                                                                             <  NB501016:228:HFK2FAFXY:3:11409:16149:8234/150‑1 (MQ=60)
                                         AAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCAC                                                                                          <  NB501016:228:HFK2FAFXY:2:11201:10169:8727/151‑1 (MQ=60)
                                           GCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCACCTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGG                                                                                                                     >  NB501016:228:HFK2FAFXY:1:21307:26750:1412/1‑122 (MQ=60)
                                           GCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGG                                                                                                                     >  NB501016:228:HFK2FAFXY:1:21306:1164:1437/1‑122 (MQ=60)
                                                     TCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCG                                                                                 <  NB501016:228:HFK2FAFXY:1:11312:7197:12378/148‑1 (MQ=60)
                                                      CGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTAC                                                                                                             >  NB501016:228:HFK2FAFXY:1:11105:2291:6380/1‑119 (MQ=60)
                                                                    cgtgtgcttcttccgatctCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTG                                                               <  NB501016:228:HFK2FAFXY:3:11602:22300:15198/132‑1 (MQ=60)
                                                                    CGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGCCGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTG                                                               >  NB501016:228:HFK2FAFXY:2:21207:15074:4624/1‑151 (MQ=60)
                                                                           CTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGC                                                          <  NB501016:228:HFK2FAFXY:1:11208:15411:6369/149‑1 (MQ=60)
                                                                                      GCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCATATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCA                                             <  NB501016:228:HFK2FAFXY:2:21312:21177:9010/151‑1 (MQ=60)
                                                                                        TTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAAC                                           >  NB501016:228:HFK2FAFXY:3:21502:22604:5631/1‑151 (MQ=60)
                                                                                           AAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCC                                        <  NB501016:228:HFK2FAFXY:1:21204:24093:2805/151‑1 (MQ=60)
                                                                                                  GCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCAC                                                                                          <  NB501016:228:HFK2FAFXY:3:11501:6041:2465/94‑1 (MQ=60)
                                                                                                           ACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAG                        <  NB501016:228:HFK2FAFXY:2:11103:3338:13009/151‑1 (MQ=60)
                                                                                                                      GCGCTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAAC             <  NB501016:228:HFK2FAFXY:4:11504:18530:8318/151‑1 (MQ=60)
                                                                                                                       CGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAAC             <  NB501016:228:HFK2FAFXY:2:21111:22619:12829/150‑1 (MQ=60)
                                                                                                                         CTCACGACAAACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGAACAGACGCAAAGAAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGC          >  NB501016:228:HFK2FAFXY:4:11509:10003:6659/1‑151 (MQ=60)
                                                                                                                                  AACATCAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTC  >  NB501016:228:HFK2FAFXY:3:11503:7803:9810/1‑150 (MQ=60)
                                                                                                                                    CATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGAC                                    <  NB501016:228:HFK2FAFXY:4:11409:7047:2397/114‑1 (MQ=60)
                                                                                                                                       CAGCCGTCAGGCCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAG                          >  NB501016:228:HFK2FAFXY:1:11209:13312:2085/1‑121 (MQ=60)
                                                                                                                                                  |                                                                                                                                     
ATGCACTGTCAACACGGATCATTGCGTTGTAATCGGCACCAAAGCCGCTGGCTTCGCGCGCGCGACGACGCTGCTCTTCCATTGCAGCTTCAAAACCAGCTTCGTCAACTTTGATGTTGCGCTCACGACAAACATCAGCCGTCAGGTCAACCGGGAAGCCATAGGTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTC  >  NC_000913/2820646‑2820925

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