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
A1 F1 I1 R1 66 46.7 5280542 98.1% 5180211 47.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_017635 2,407,182 Δ1 bp coding (15/1011 nt) mglC ← galactose/methyl galactoside ABC transporter permease MglC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176352,407,1820A.100.0% 431.4 / NA 88coding (15/1011 nt)mglCgalactose/methyl galactoside ABC transporter permease MglC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (18/70);  total (18/70)

ACGTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTATTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCA  >  NC_017635/2407138‑2407222
                                            |                                        
acgTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTtttta                                       <  2:205906/48‑2 (MQ=255)
  gTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAAc                                    >  2:1213169/1‑48 (MQ=39)
  gTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAAc                                    <  1:2296630/48‑1 (MQ=39)
  gTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAAc                                    <  2:1446980/48‑1 (MQ=39)
  gTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAAc                                    <  1:2631785/48‑1 (MQ=39)
  gTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAAc                                    >  1:550544/1‑48 (MQ=39)
    aaaTACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACg                                   >  1:771832/1‑47 (MQ=39)
     aaTACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCa                                 <  1:446717/48‑1 (MQ=39)
     aaTACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCa                                 <  2:870152/48‑1 (MQ=39)
     aaTACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCa                                 <  1:200031/48‑1 (MQ=39)
      aTACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCAc                                <  1:824062/48‑1 (MQ=39)
      aTACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCAc                                >  1:1748487/1‑48 (MQ=39)
      aTACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCAc                                >  1:2011746/1‑48 (MQ=39)
       tACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACt                               >  2:938750/1‑48 (MQ=39)
         ccGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTc                              <  1:602867/47‑1 (MQ=39)
         ccGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTc                              <  1:1185946/47‑1 (MQ=39)
         ccGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTc                              <  1:114765/47‑1 (MQ=39)
         ccGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCa                             <  2:2047106/48‑1 (MQ=39)
           gCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATg                           <  2:1012281/48‑1 (MQ=39)
           gCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATg                           <  2:146856/48‑1 (MQ=39)
           gCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATg                           <  1:1545483/48‑1 (MQ=39)
            cccTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATgg                          <  2:1434743/48‑1 (MQ=39)
                  ttCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcc                    >  1:1604017/1‑48 (MQ=255)
                  ttCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcc                    >  2:597807/1‑48 (MQ=255)
                  ttCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcc                    >  2:2332590/1‑48 (MQ=255)
                  ttCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcc                    >  2:370550/1‑48 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAg                      <  2:310674/45‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAg                      <  2:59419/45‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcc                    <  1:298617/47‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcc                    <  1:1883203/47‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcc                    <  2:1915900/47‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcc                    <  2:1270238/47‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcc                    <  1:863720/47‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcc                    <  2:1073895/47‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcc                    <  1:985334/47‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcc                    <  1:1292895/47‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGc                     <  2:363027/46‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGc                     <  2:650395/46‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGc                     <  2:160028/46‑1 (MQ=255)
                   tCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGc                     <  2:360736/46‑1 (MQ=255)
                    cAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGccc                   <  2:1612089/47‑1 (MQ=255)
                    cAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGccc                   <  2:1583116/47‑1 (MQ=255)
                    cAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGccc                   <  1:1003425/47‑1 (MQ=255)
                    cAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGccc                   <  1:640844/47‑1 (MQ=255)
                    cAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGccc                   <  1:967477/47‑1 (MQ=255)
                    cAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGccc                   <  1:854487/47‑1 (MQ=255)
                    cAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGccc                   <  2:2284870/47‑1 (MQ=255)
                    cAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGccc                   <  1:2550808/47‑1 (MQ=255)
                    cAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGccc                   <  1:572368/47‑1 (MQ=255)
                     aGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGcccc                  <  2:1844033/47‑1 (MQ=255)
                     aGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCt                 <  1:657792/48‑1 (MQ=255)
                      ggTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTa                <  2:2264793/48‑1 (MQ=255)
                      ggTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTa                <  1:1333441/48‑1 (MQ=255)
                      ggTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTa                <  1:1478226/48‑1 (MQ=255)
                      ggTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTa                <  1:209682/48‑1 (MQ=255)
                      ggTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTa                <  2:2243373/48‑1 (MQ=255)
                      ggTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTa                <  2:1003680/48‑1 (MQ=255)
                       gtaagtaaGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTaa               <  2:1423051/48‑1 (MQ=255)
                       gtaagtaaGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTaa               <  1:2038358/48‑1 (MQ=255)
                       gtaagtaaGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTaa               <  2:2440336/48‑1 (MQ=255)
                        taagtaagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAAt              <  2:1032610/48‑1 (MQ=255)
                        taagtaagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAAt              <  2:1908099/48‑1 (MQ=255)
                        taagtaagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAAt              <  1:1333918/48‑1 (MQ=255)
                        taagtaagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAAt              <  1:1344586/48‑1 (MQ=255)
                        taagtaagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAAt              <  1:1989998/48‑1 (MQ=255)
                        taagtaagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAAt              <  1:37113/48‑1 (MQ=255)
                        taagtaagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAAt              <  1:680615/48‑1 (MQ=255)
                        taagtaagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAAt              <  1:684956/48‑1 (MQ=255)
                          agtaagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATc             <  1:2214563/47‑1 (MQ=255)
                            taagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCtt           >  2:2310385/1‑47 (MQ=255)
                             aagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTa          >  1:1028674/1‑47 (MQ=255)
                             aagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTa          >  2:1772342/1‑47 (MQ=255)
                             aagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTa          >  2:1746057/1‑47 (MQ=255)
                             aagAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTa          >  2:1111012/1‑47 (MQ=255)
                              agAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTaa         <  2:1101127/47‑1 (MQ=255)
                              agAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTaa         <  2:233996/47‑1 (MQ=255)
                              agAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTaa         <  2:1178188/47‑1 (MQ=255)
                              agAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTaa         <  2:583905/47‑1 (MQ=255)
                               gAAAACTTTTCTTTTTAACGCACTCATGGGAATCCCCTAATCTTAAAg       >  1:1078430/1‑48 (MQ=255)
                               gAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAg       <  2:1761859/48‑1 (MQ=255)
                               gAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAg       <  2:916225/48‑1 (MQ=255)
                               gAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAg       <  2:6476/48‑1 (MQ=255)
                                aaaaCTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAgg      >  2:1349806/1‑48 (MQ=255)
                                aaaaCTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAgg      <  2:921137/48‑1 (MQ=255)
                                   aCTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGc   <  1:1809953/48‑1 (MQ=255)
                                   aCTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGc   <  1:2529212/48‑1 (MQ=255)
                                   aCTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGc   <  2:2118704/48‑1 (MQ=255)
                                   aCTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGc   <  1:997284/48‑1 (MQ=255)
                                    cttttcttttTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGc   >  2:2110367/1‑47 (MQ=255)
                                     ttttCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCa  >  1:814180/1‑47 (MQ=255)
                                            |                                        
ACGTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTATTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCA  >  NC_017635/2407138‑2407222

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGAACGACGTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTATTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCAAAGACGC  >  NC_017635/2407132‑2407229
                                                  |                                               
AGAACGACGTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTT                                                 <  M01420:84:000000000‑ABVEH:1:2114:22017:20383/51‑1 (MQ=60)
 GAACGACGTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTT                                                <  M01420:84:000000000‑ABVEH:1:2109:17096:13220/51‑1 (MQ=60)
  AACGACGTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTT                                               >  M01420:84:000000000‑ABVEH:1:2112:11266:13955/1‑51 (MQ=60)
      ACGTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACG                                          >  M01420:84:000000000‑ABVEH:1:1105:19171:20681/1‑51 (MQ=60)
      ACGTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACG                                          >  M01420:84:000000000‑ABVEH:1:2112:22031:4414/1‑51 (MQ=60)
       CGTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGC                                         <  M01420:84:000000000‑ABVEH:1:2109:16803:14329/51‑1 (MQ=60)
       CGTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGC                                         <  M01420:84:000000000‑ABVEH:1:2106:18790:18104/51‑1 (MQ=60)
        GTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGC                                         >  M01420:84:000000000‑ABVEH:1:1107:2913:18097/1‑50 (MQ=60)
          AAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACT                                      >  M01420:84:000000000‑ABVEH:1:2104:10560:7374/1‑51 (MQ=60)
          AAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACT                                      >  M01420:84:000000000‑ABVEH:1:2101:7381:23829/1‑51 (MQ=60)
          AAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACT                                      <  M01420:84:000000000‑ABVEH:1:2112:27747:14485/51‑1 (MQ=60)
          AAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACT                                      <  M01420:84:000000000‑ABVEH:1:1102:18847:20529/51‑1 (MQ=60)
          AAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACT                                      <  M01420:84:000000000‑ABVEH:1:1104:6341:23394/51‑1 (MQ=60)
           AATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTC                                     <  M01420:84:000000000‑ABVEH:1:2114:18008:6226/51‑1 (MQ=60)
           AATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTC                                     <  M01420:84:000000000‑ABVEH:1:1108:13600:4715/51‑1 (MQ=60)
           AATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTC                                     <  M01420:84:000000000‑ABVEH:1:2114:20329:15632/51‑1 (MQ=60)
              ACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCAT                                   <  M01420:84:000000000‑ABVEH:1:1111:20990:6730/50‑1 (MQ=60)
              ACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCAT                                   <  M01420:84:000000000‑ABVEH:1:1101:14366:28602/50‑1 (MQ=60)
              ACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCAT                                   <  M01420:84:000000000‑ABVEH:1:1106:9578:8012/50‑1 (MQ=60)
                CGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGG                                <  M01420:84:000000000‑ABVEH:1:1114:26054:6914/51‑1 (MQ=60)
                      CTTTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCC                          >  M01420:84:000000000‑ABVEH:1:1114:25601:18223/1‑51 (MQ=60)
                        TTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCC                         <  M01420:84:000000000‑ABVEH:1:2103:13865:3657/50‑1 (MQ=60)
                        TTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCC                         <  M01420:84:000000000‑ABVEH:1:1109:16580:14027/50‑1 (MQ=60)
                        TTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCC                         <  M01420:84:000000000‑ABVEH:1:1108:17981:13262/50‑1 (MQ=60)
                        TTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCC                         <  M01420:84:000000000‑ABVEH:1:1112:11991:3321/50‑1 (MQ=60)
                        TTCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCC                         <  M01420:84:000000000‑ABVEH:1:1103:24869:16986/50‑1 (MQ=60)
                         TCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCT                        <  M01420:84:000000000‑ABVEH:1:2108:18967:21609/50‑1 (MQ=60)
                         TCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCT                        <  M01420:84:000000000‑ABVEH:1:2113:24361:5998/50‑1 (MQ=60)
                         TCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCT                        <  M01420:84:000000000‑ABVEH:1:2112:11002:20966/50‑1 (MQ=60)
                         TCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCT                        <  M01420:84:000000000‑ABVEH:1:1109:8102:17424/50‑1 (MQ=60)
                         TCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCT                        <  M01420:84:000000000‑ABVEH:1:1109:5720:10595/50‑1 (MQ=60)
                         TCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCT                        <  M01420:84:000000000‑ABVEH:1:1108:26174:11485/50‑1 (MQ=60)
                         TCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCT                        <  M01420:84:000000000‑ABVEH:1:1105:9295:25731/50‑1 (MQ=60)
                         TCAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCT                        <  M01420:84:000000000‑ABVEH:1:1106:5730:15481/50‑1 (MQ=60)
                          CAGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAA                      <  M01420:84:000000000‑ABVEH:1:1106:11359:18770/51‑1 (MQ=60)
                           AGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAAT                     >  M01420:84:000000000‑ABVEH:1:2114:23403:6250/1‑51 (MQ=60)
                           AGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAAT                     <  M01420:84:000000000‑ABVEH:1:1113:18521:17157/51‑1 (MQ=60)
                           AGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAAT                     <  M01420:84:000000000‑ABVEH:1:1112:8108:12453/51‑1 (MQ=60)
                           AGGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAAT                     <  M01420:84:000000000‑ABVEH:1:1102:19098:22572/51‑1 (MQ=60)
                            GGTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATC                    <  M01420:84:000000000‑ABVEH:1:2104:19378:12848/51‑1 (MQ=60)
                             GTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCT                   <  M01420:84:000000000‑ABVEH:1:1112:21030:14542/51‑1 (MQ=60)
                             GTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCT                   <  M01420:84:000000000‑ABVEH:1:1101:17062:11311/51‑1 (MQ=60)
                             GTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCT                   <  M01420:84:000000000‑ABVEH:1:2103:15507:27142/51‑1 (MQ=60)
                             GTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCT                   <  M01420:84:000000000‑ABVEH:1:1106:11002:23520/51‑1 (MQ=60)
                             GTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCT                   <  M01420:84:000000000‑ABVEH:1:1106:21856:24540/51‑1 (MQ=60)
                             GTAAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCT                   <  M01420:84:000000000‑ABVEH:1:1112:6716:12541/51‑1 (MQ=60)
                               AAGTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTT                  <  M01420:84:000000000‑ABVEH:1:2105:18638:25520/50‑1 (MQ=60)
                                 GTAAGAAAACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAA                >  M01420:84:000000000‑ABVEH:1:1109:8777:22460/1‑50 (MQ=60)
                                   AAGAAAACTTTTCTTTTTAACGCACTCATGGGAATCCCCTAATCTTAAAGG             >  M01420:84:000000000‑ABVEH:1:1110:20989:9104/1‑51 (MQ=57)
                                        AACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCAA        <  M01420:84:000000000‑ABVEH:1:1109:5526:16264/51‑1 (MQ=60)
                                        AACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCAA        <  M01420:84:000000000‑ABVEH:1:2108:27285:17434/51‑1 (MQ=60)
                                        AACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCAA        <  M01420:84:000000000‑ABVEH:1:2102:8451:13120/51‑1 (MQ=60)
                                         ACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCAA        >  M01420:84:000000000‑ABVEH:1:1107:12964:28188/1‑50 (MQ=60)
                                         ACTTTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCAA        >  M01420:84:000000000‑ABVEH:1:2112:27843:8612/1‑50 (MQ=60)
                                            TTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCAAAG      <  M01420:84:000000000‑ABVEH:1:1106:8855:21591/49‑1 (MQ=60)
                                            TTTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCAAAG      <  M01420:84:000000000‑ABVEH:1:2108:23934:22599/49‑1 (MQ=60)
                                             TTCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCAAAGAC    <  M01420:84:000000000‑ABVEH:1:1104:20857:25478/50‑1 (MQ=60)
                                              TCTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCAAAGACGC  >  M01420:84:000000000‑ABVEH:1:2106:15708:20333/1‑51 (MQ=60)
                                               CTTTTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCAAAGACGC  >  M01420:84:000000000‑ABVEH:1:1104:23529:9279/1‑50 (MQ=60)
                                                  |                                               
AGAACGACGTAAATACCGCCCTCTTTCAGGTAAGTAAGAAAACTTTTCTTATTTAACGCACTCATGGGAAGCCCCTAATCTTAAAGGTGCAAAGACGC  >  NC_017635/2407132‑2407229

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