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
A4 F4 I1 R1 17 37.4 1246952 98.2% 1224506 142.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,964,929 C→T E234E (GAG→GAA)  rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,964,9290CT100.0% 112.2 / NA 34E234E (GAG→GAA) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/21);  total (13/21)

ACACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGG  >  NC_000913/3964783‑3965072
                                                                                                                                                  |                                                                                                                                               
acacGCTGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGGGTTTGGGGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGttc                                                                                                                                               <  4:104597/149‑1 (MQ=255)
   cGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGttcttc                                                                                                                                            <  2:76978/149‑1 (MQ=255)
     aTGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTcttctt                                                                                                                                          <  1:46791/149‑1 (MQ=255)
     aTGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTcttctt                                                                                                                                          <  1:8924/149‑1 (MQ=255)
                  gcCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGcc                                                                                                                             >  6:94707/1‑149 (MQ=255)
                     aGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGttcttc                                                                                                                                            >  4:20112/1‑131 (MQ=255)
                     aGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGttcttc                                                                                                                                            <  3:20112/131‑1 (MQ=255)
                            cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGt                                                                                                                   >  1:58890/1‑149 (MQ=255)
                                        cacaACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCAc                                                                                                        <  5:113244/148‑1 (MQ=255)
                                        cacaACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACt                                                                                                       >  4:77277/1‑149 (MQ=255)
                                             cGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAat                                                                                                  <  1:32837/149‑1 (MQ=255)
                                                tgtTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAatata                                                                                               <  4:72251/149‑1 (MQ=255)
                                                      gtgtTGGCGTAAATAATCGCTCTGTCTGTCCTCTCTGCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGGTAGAAGGGTATATAATTTCTTCTTTAATACGGTGGCACGTTTTCTGTTCCGGTTCCACTTCAATATATTCgg                                                                                          <  7:137720/148‑1 (MQ=255)
                                                       tgtTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCa                                                                                        >  7:35882/1‑149 (MQ=255)
                                                       tgtTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCa                                                                                        <  6:43018/149‑1 (MQ=255)
                                                       tgtTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCa                                                                                        <  1:50673/149‑1 (MQ=255)
                                                           ggCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGGTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTg                                                                                     <  8:92738/148‑1 (MQ=255)
                                                            gCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGtt                                                                                   >  5:67831/1‑149 (MQ=255)
                                                              gAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAAt                                                                                                                                      >  3:70693/1‑96 (MQ=255)
                                                              gAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAAt                                                                                                                                      <  4:70693/96‑1 (MQ=255)
                                                               aaaaTAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCAt                                                                                <  2:58890/149‑1 (MQ=255)
                                                                              tctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGGATTTTTTCTTCGTTAGAAGTGTAGAAAAGTTCTTCTTTAATACGGTTGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCa                                                                  <  6:34652/148‑1 (MQ=255)
                                                                              tctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCa                                                                  <  5:94707/148‑1 (MQ=255)
                                                                              tctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCa                                                                  <  3:77277/148‑1 (MQ=255)
                                                                                   ccACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCAc                                                            <  3:77156/149‑1 (MQ=255)
                                                                                              cGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTCTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTAc                                                 >  7:49197/1‑149 (MQ=255)
                                                                                                  atgCGTTTGCTGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGATAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGa                                               >  1:47638/3‑147 (MQ=255)
                                                                                                       tttGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTg                                        <  8:35882/149‑1 (MQ=255)
                                                                                                            agcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCg                                     >  5:155307/1‑147 (MQ=255)
                                                                                                            agcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGga                                   <  6:155307/149‑1 (MQ=255)
                                                                                                                 aCGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTCTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACa                              <  8:49197/149‑1 (MQ=255)
                                                                                                                   gCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGc                            <  6:67831/149‑1 (MQ=255)
                                                                                                                     aTTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCTGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCAGGTACGAAAGCGTGGCGGAGAACAGCAt                          >  7:32384/1‑149 (MQ=255)
                                                                                                                                gTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCt               >  6:36646/1‑149 (MQ=255)
                                                                                                                                 ttAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACa                              <  3:127960/133‑1 (MQ=255)
                                                                                                                                 ttAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACa                              >  4:127960/1‑133 (MQ=255)
                                                                                                                                          tAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAgg     >  3:33779/1‑149 (MQ=255)
                                                                                                                                             aaaaGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggcgg  <  5:36646/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ACACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGG  >  NC_000913/3964783‑3965072

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGC  >  NC_000913/3964787‑3965073
                                                                                                                                              |                                                                                                                                                
GATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTT                                                                                                                                          <  NB501016:228:HFK2FAFXY:1:11104:23508:8225/151‑1 (MQ=60)
GATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTT                                                                                                                                          <  NB501016:228:HFK2FAFXY:1:11207:5496:7886/151‑1 (MQ=60)
                CAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCT                                                                                                                                            <  NB501016:228:HFK2FAFXY:2:11111:24218:12627/133‑1 (MQ=60)
                       CCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTT                                                                                                                   >  NB501016:228:HFK2FAFXY:1:11212:5735:5289/1‑151 (MQ=60)
                                   TCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACT                                                                                                        <  NB501016:228:HFK2FAFXY:3:21503:6571:8630/150‑1 (MQ=60)
                                        ACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATA                                                                                                  <  NB501016:228:HFK2FAFXY:1:11201:20113:15664/151‑1 (MQ=60)
                                                 GGTGTTGGCGTAAATAATCGCTCTGTCTGTCCTCTCTGCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGGTAGAAGGGTATATAATTTCTTCTTTAATACGGTGGCACGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGC                                                                                          <  NB501016:228:HFK2FAFXY:4:21611:20866:19515/150‑1 (MQ=60)
                                                  GTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCAT                                                                                        <  NB501016:228:HFK2FAFXY:1:11208:24843:18084/151‑1 (MQ=60)
                                                  GTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCATCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCAT                                                                                        >  NB501016:228:HFK2FAFXY:4:11506:24877:6647/1‑151 (MQ=60)
                                                       GGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTC                                                                                   >  NB501016:228:HFK2FAFXY:3:11606:23351:14208/1‑151 (MQ=60)
                                                         CGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATA                                                                                                                                      >  NB501016:228:HFK2FAFXY:2:11312:14661:4271/1‑98 (MQ=60)
                                                                         GTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAG                                                                  <  NB501016:228:HFK2FAFXY:2:21103:16033:9844/150‑1 (MQ=60)
                                                                         GTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAG                                                                  <  NB501016:228:HFK2FAFXY:3:21406:8586:13840/150‑1 (MQ=60)
                                                                              GCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACG                                                            <  NB501016:228:HFK2FAFXY:2:21103:17126:8822/151‑1 (MQ=60)
                                                                                CACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTC                                                                                                                                                >  NB501016:228:HFK2FAFXY:1:21303:14665:1081/1‑65 (MQ=60)
                                                                                         TCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTCTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACG                                                 >  NB501016:228:HFK2FAFXY:4:11601:22990:8781/1‑151 (MQ=60)
                                                                                             aatGCGTTTGCTGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGATAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAA                                               >  NB501016:228:HFK2FAFXY:1:11207:22739:14381/4‑149 (MQ=60)
                                                                                                       CAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGG                                     >  NB501016:228:HFK2FAFXY:3:21610:14345:13063/1‑149 (MQ=60)
                                                                                                                CATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCTGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCAGGTACGAAAGCGTGGCGGAGAACAGCATG                          >  NB501016:228:HFK2FAFXY:4:11504:1727:8424/1‑151 (MQ=60)
                                                                                                                        CTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTATTCTGTTCCGGTTCCACTTCAATATATTCGGCATAGTTAATCAGCTCGAACGCCAGTTCACGTACCAGGTACGAAAGAGTGGCGGAGAACAGCATGTTGAGGCC                  >  NB501016:228:HFK2FAFXY:1:21105:24380:18542/1‑151 (MQ=60)
                                                                                                                            GTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAG                              <  NB501016:228:HFK2FAFXY:2:21306:8780:1094/135‑1 (MQ=60)
                                                                                                                                     GTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGC     >  NB501016:228:HFK2FAFXY:2:11206:12268:13803/1‑151 (MQ=60)
                                                                                                                                        GAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGC  <  NB501016:228:HFK2FAFXY:3:11505:12385:8208/151‑1 (MQ=60)
                                                                                                                                              |                                                                                                                                                
GATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGC  >  NC_000913/3964787‑3965073

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: ‑