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
A3 F36 I1 R1 22 37.6 1496218 97.1% 1452827 127.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,768,131 A→T pseudogene (911/945 nt) rhsJ → Rhs‑family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,768,1310AT100.0% 144.1 / NA 41pseudogene (911/945 nt)rhsJRhs‑family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (20/21);  total (20/21)

GGAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCAAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTAT  >  NC_000913/3767986‑3768265
                                                                                                                                                 |                                                                                                                                      
ggAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAgg                                                                                                                                     >  7:113287/1‑149 (MQ=255)
                               tGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAg                                                                                                                                      >  4:160306/1‑117 (MQ=255)
                               tGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAg                                                                                                                                      <  3:160306/117‑1 (MQ=255)
                                 cTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGTCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGt                                                                                                    <  7:36052/149‑1 (MQ=255)
                                       gCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGa                                                                                              <  8:113287/149‑1 (MQ=255)
                                        cACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGAt                                                                                             <  5:166619/149‑1 (MQ=255)
                                                aaCCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGtt                                                                                                   <  8:24222/135‑1 (MQ=255)
                                                       aTTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGt                                                                                                                  <  6:99840/113‑1 (MQ=255)
                                                       aTTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGt                                                                                                                  >  5:99840/1‑113 (MQ=255)
                                                        ttCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGa                                                                                                                            >  5:2867/1‑102 (MQ=255)
                                                        ttCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGa                                                                                                                            <  6:2867/102‑1 (MQ=255)
                                                        ttCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCtt                                                                             >  7:100472/1‑149 (MQ=255)
                                                                                aacaacAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAg                                                                                                          >  8:163341/1‑96 (MQ=255)
                                                                                aacaacAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAg                                                                                                          <  7:163341/96‑1 (MQ=255)
                                                                                 acaacaAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCtttt                                                    <  8:100472/149‑1 (MQ=255)
                                                                                  caacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTa                                                   >  1:48229/1‑149 (MQ=255)
                                                                                   aacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGt                                                                                           >  1:26971/1‑108 (MQ=255)
                                                                                   aacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGt                                                                                           <  2:26971/108‑1 (MQ=255)
                                                                                   aacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGaaa                                                                                     <  8:83478/114‑1 (MQ=255)
                                                                                   aacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGaaa                                                                                     >  7:83478/1‑114 (MQ=255)
                                                                                     caaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTtctc                                                        >  4:40735/1‑141 (MQ=255)
                                                                                     caaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTtctc                                                        <  3:40735/141‑1 (MQ=255)
                                                                                                       ggATAACATATGGGCCATCTGGGCACCCAATTTCTAATATGCTAGGTACGCCAGAAAAACCAAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTTTCAATTGCATGGTGCt                               >  8:131567/1‑148 (MQ=255)
                                                                                                        gATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTg                                                                                                  <  7:109696/80‑1 (MQ=255)
                                                                                                        gATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTg                                                                                                  >  8:109696/1‑80 (MQ=255)
                                                                                                        gATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGa                             >  6:22589/1‑149 (MQ=255)
                                                                                                            aCATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGAtgt                          <  5:22589/148‑1 (MQ=255)
                                                                                                               tatGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATg                                                                                            >  5:59501/1‑79 (MQ=255)
                                                                                                               tatGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATg                                                                                            <  6:59501/79‑1 (MQ=255)
                                                                                                                atGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGaaaa                     >  5:61396/1‑149 (MQ=255)
                                                                                                                   ggCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaa                  <  2:48229/149‑1 (MQ=255)
                                                                                                                   ggCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTAttttt                                                                 <  1:158097/101‑1 (MQ=255)
                                                                                                                   ggCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTAttttt                                                                 >  2:158097/1‑101 (MQ=255)
                                                                                                                            ggCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGt                                              >  5:62395/1‑112 (MQ=255)
                                                                                                                            ggCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGt                                              <  6:62395/112‑1 (MQ=255)
                                                                                                                                 cccATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGaaa                      <  7:32/131‑1 (MQ=255)
                                                                                                                                 cccATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGaaa                      >  8:32/1‑131 (MQ=255)
                                                                                                                                 cccATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaaa                >  5:85818/1‑137 (MQ=255)
                                                                                                                                 cccATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaaa                <  6:85818/137‑1 (MQ=255)
                                                                                                                                  ccATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTa   >  5:122023/1‑149 (MQ=255)
                                                                                                                                   cATTTCTAATATGCTAGGTACTCCAGATAAAACTATTTATATAAGATATGTTGTAGATGTATGATAGCATGCTTTTTACTATTTTTTTATTTCTATTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTat  <  7:131567/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                      
GGAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCAAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTAT  >  NC_000913/3767986‑3768265

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
GCGAGGAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCAAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATG  >  NC_000913/3767982‑3768266
                                                                                                                                                     |                                                                                                                                       
GCGAGGAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTA                                                                                                                                        <  NB501016:228:HFK2FAFXY:3:11501:22208:11658/151‑1 (MQ=60)
   AGGAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGT                                                                                                                                     >  NB501016:228:HFK2FAFXY:4:21411:4586:17876/1‑151 (MQ=60)
                      AAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCT                                                                                                                                         <  NB501016:228:HFK2FAFXY:3:11510:17109:15936/128‑1 (MQ=60)
                              ggcatactagataagaattcgtgactggattcagacgtgtgatcttccgatctTAACAACAAAAGACCTTTGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGA                                                                                                          <  NB501016:228:HFK2FAFXY:4:21610:8387:2678/98‑1 (MQ=60)
                                  GTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGG                                                                                                                                      <  NB501016:228:HFK2FAFXY:2:21308:8789:14215/119‑1 (MQ=60)
                                    GCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGTCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTT                                                                                                    <  NB501016:228:HFK2FAFXY:4:11503:16131:2195/151‑1 (MQ=60)
                                           GCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATG                                                                                             <  NB501016:228:HFK2FAFXY:3:21602:7174:14424/151‑1 (MQ=60)
                                                   GAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGagatcggaaggcac                                                                                     >  NB501016:228:HFK2FAFXY:4:11410:6045:8666/1‑137 (MQ=60)
                                                          AATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTG                                                                                                                  >  NB501016:228:HFK2FAFXY:3:21401:21123:6729/1‑115 (MQ=60)
                                                           ATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCCCCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTT                                                                             >  NB501016:228:HFK2FAFXY:4:21406:8403:20059/1‑151 (MQ=60)
                                                           ATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAA                                                                                                                            >  NB501016:228:HFK2FAFXY:3:11402:16010:2021/1‑104 (MQ=60)
                                                                                     ACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTAT                                                   >  NB501016:228:HFK2FAFXY:1:11204:9221:19076/1‑151 (MQ=60)
                                                                                      CAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAG                                                                                     >  NB501016:228:HFK2FAFXY:4:11611:16304:20385/1‑116 (MQ=60)
                                                                                      CAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTA                                                                                           >  NB501016:228:HFK2FAFXY:1:11109:12604:17297/1‑110 (MQ=60)
                                                                                        ACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCT                                                        <  NB501016:228:HFK2FAFXY:2:11206:9569:12021/143‑1 (MQ=60)
                                                                                                           GGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGT                                                                                                  <  NB501016:228:HFK2FAFXY:4:21410:3208:11710/82‑1 (MQ=60)
                                                                                                               AACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTG                          <  NB501016:228:HFK2FAFXY:3:11409:14780:9839/150‑1 (MQ=60)
                                                                                                                  ATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGT                                                                                            >  NB501016:228:HFK2FAFXY:3:11510:25412:17471/1‑81 (MQ=60)
                                                                                                                   TATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAAT                     >  NB501016:228:HFK2FAFXY:3:11511:4679:11361/1‑151 (MQ=60)
                                                                                                                      GGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTT                                                                 <  NB501016:228:HFK2FAFXY:1:21206:10729:4810/103‑1 (MQ=60)
                                                                                                                               GGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTC                                              >  NB501016:228:HFK2FAFXY:3:11511:7184:18730/1‑114 (MQ=60)
                                                                                                                                    CCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAAC                >  NB501016:228:HFK2FAFXY:3:11608:10885:4275/1‑139 (MQ=60)
                                                                                                                                    CCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAA                      <  NB501016:228:HFK2FAFXY:4:11401:14042:1346/133‑1 (MQ=60)
                                                                                                                                     CCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTAT   >  NB501016:228:HFK2FAFXY:3:21409:17855:14731/1‑151 (MQ=60)
                                                                                                                                      CCATTTCTAATATGCTAGGTACTCCAGATAAAACTATTTATATAAGATATGTTGTAGATGTATGATAGCATGCTTTTTACTATTTTTTTATTTCTATTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATG  <  NB501016:228:HFK2FAFXY:4:21507:19803:18116/151‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                       
GCGAGGAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCAAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATG  >  NC_000913/3767982‑3768266

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