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 1548 369.5 20095359 96.9% 19472402 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,887,232 A→G G297G (GGT→GGC)  cysI ← sulfite reductase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,887,2320AG94.9% 823.3 / 22.0 294G297G (GGT→GGC) cysIsulfite reductase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (11/4);  new base G (151/126);  total (163/131)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

Only 100 of 308 total aligned reads displayed.

CGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGC  >  NC_000913/2887142‑2887321
                                                                                          |                                                                                         
   gcgACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCa                                                                                      >  2:6036006/1‑93 (MQ=255)
   gcgACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCa                                                                                      >  1:3544187/1‑93 (MQ=255)
     gACCCTCCATTTCCGCTTGAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGt                                                                                    <  1:7578163/93‑1 (MQ=255)
      aCGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTc                                                                                   >  2:7606599/1‑93 (MQ=255)
          tCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGc                                                                               <  1:10003970/93‑1 (MQ=255)
          tCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGc                                                                               >  2:2189228/1‑93 (MQ=255)
          tCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGc                                                                               <  1:4589687/93‑1 (MQ=255)
            cACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTg                                                                             >  1:2955759/1‑93 (MQ=255)
             aCTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGa                                                                            <  2:1315898/93‑1 (MQ=255)
             aCTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGa                                                                            >  1:7856335/1‑93 (MQ=255)
               ttCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGt                                                                          <  2:1574020/93‑1 (MQ=255)
                  cGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGt                                                                       >  2:8022686/1‑93 (MQ=255)
                  cGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTgg                                                                       >  2:3830721/1‑92 (MQ=255)
                    cTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCa                                                                     <  1:7001481/93‑1 (MQ=255)
                     tttAAACGTCTCACCCCCCGCGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCac                                                                    <  2:9654038/93‑1 (MQ=255)
                     tttAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCac                                                                    <  2:3093263/93‑1 (MQ=255)
                      ttAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCacg                                                                   >  2:4596226/1‑93 (MQ=255)
                       tAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCacga                                                                  <  1:6907279/93‑1 (MQ=255)
                        aaaCGTCTCACCCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCacgac                                                                 <  1:5802197/93‑1 (MQ=255)
                        aaaCGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCacgac                                                                 <  1:6316440/93‑1 (MQ=255)
                        aaaCGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCacgac                                                                 <  1:1483683/93‑1 (MQ=255)
                         aaCGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCacgacg                                                                >  1:9660286/1‑93 (MQ=255)
                           cGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGc                                                              <  2:349391/93‑1 (MQ=255)
                           cGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGc                                                              <  1:7325016/93‑1 (MQ=255)
                            gTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCt                                                             >  2:9264187/1‑93 (MQ=255)
                            gTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGCTTCGGATGACCCAGTCACGCTGAGTTGTCACGACGGCt                                                             >  2:6945351/1‑93 (MQ=255)
                             tctcAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCtt                                                            <  2:5479309/93‑1 (MQ=255)
                             tctcAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCtt                                                            >  2:1806701/1‑93 (MQ=255)
                             tctcAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCtc                                                            >  1:7215389/1‑92 (MQ=255)
                               tcAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCg                                                          >  1:3926350/1‑93 (MQ=255)
                                 aaCCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGc                                                        <  1:8868621/93‑1 (MQ=255)
                                  aCCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGcc                                                       <  2:5257991/93‑1 (MQ=255)
                                   cccccACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCa                                                      >  1:1414183/1‑93 (MQ=255)
                                     cccACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCAcc                                                    >  1:847971/1‑93 (MQ=255)
                                     cccACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCAcc                                                    <  1:6151867/93‑1 (MQ=255)
                                     cccACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCAcc                                                    >  2:4911919/1‑93 (MQ=255)
                                      ccACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCg                                                   >  1:4499186/1‑93 (MQ=255)
                                       cACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGc                                                  >  1:3655397/1‑93 (MQ=255)
                                        aCGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGcc                                                 >  2:2326730/1‑93 (MQ=255)
                                          gcgcTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAg                                               <  2:1233810/93‑1 (MQ=255)
                                            gcTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCg                                             <  2:8167863/93‑1 (MQ=255)
                                            gcTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCg                                             >  2:7240614/1‑93 (MQ=255)
                                             cTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGt                                            >  2:5334835/1‑93 (MQ=255)
                                              tCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTa                                           <  1:8272585/93‑1 (MQ=255)
                                              tCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTa                                           >  2:7454340/1‑93 (MQ=255)
                                              tCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTa                                           <  2:9324102/93‑1 (MQ=255)
                                               ccAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTAt                                          >  2:4220756/1‑93 (MQ=255)
                                               ccAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTAt                                          <  2:6209187/93‑1 (MQ=255)
                                               ccAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTAt                                          >  2:3401017/1‑93 (MQ=255)
                                                cAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATg                                         >  2:1947212/1‑93 (MQ=255)
                                                cAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATg                                         <  2:306043/93‑1 (MQ=255)
                                                 aGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGc                                        <  2:8804841/93‑1 (MQ=255)
                                                 aGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGc                                        <  2:8758221/93‑1 (MQ=255)
                                                  gCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCt                                       >  1:7509749/1‑93 (MQ=255)
                                                   cGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTc                                      <  1:1967818/93‑1 (MQ=255)
                                                    gtgtATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTcc                                     <  1:7095510/93‑1 (MQ=255)
                                                     tgtATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTccc                                    >  2:1978038/1‑92 (MQ=255)
                                                     tgtATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCa                                    <  1:3420920/93‑1 (MQ=255)
                                                      gtATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAg                                   >  2:6865230/1‑93 (MQ=255)
                                                      gtATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAg                                   >  2:2453611/1‑93 (MQ=255)
                                                      gtATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAg                                   <  2:773407/93‑1 (MQ=255)
                                                      gtATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAg                                   >  1:9250301/1‑93 (MQ=255)
                                                       tATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAgg                                  >  1:2312030/1‑92 (MQ=255)
                                                       tATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGc                                  >  1:4871644/1‑93 (MQ=255)
                                                       tATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACAGCCAGCGTATGCTCCCGc                                  >  2:2771367/1‑93 (MQ=255)
                                                        aTTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCg                                 <  1:12448/93‑1 (MQ=255)
                                                         tttGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCgg                                >  2:6598268/1‑93 (MQ=255)
                                                          ttGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGc                               >  1:6581161/1‑93 (MQ=255)
                                                             gTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGa                            >  1:8653951/1‑93 (MQ=255)
                                                              ttttGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGAt                           <  2:67551/93‑1 (MQ=255)
                                                              ttttGGCATTTTTACGATCGGGTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGAt                           <  2:3179598/93‑1 (MQ=255)
                                                               tttGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCAGGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATa                          <  2:2893973/93‑1 (MQ=255)
                                                               tttGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATa                          <  2:1054627/93‑1 (MQ=255)
                                                               ttgggCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATa                          <  2:8561781/90‑1 (MQ=255)
                                                                ttGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAg                         >  2:4708510/1‑93 (MQ=255)
                                                                ttGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAg                         >  2:6345508/1‑93 (MQ=255)
                                                                   gCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCa                      >  2:7071201/1‑93 (MQ=255)
                                                                    cATTTTTACGATCGGTTCGGTTGCCGCAGTCTCGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCaa                     <  1:2138768/93‑1 (MQ=255)
                                                                        tttACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAactc                 >  2:1652304/1‑93 (MQ=255)
                                                                        tttACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAactc                 <  1:574063/93‑1 (MQ=255)
                                                                         ttACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAactca                <  1:2860522/93‑1 (MQ=255)
                                                                         ttACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAactca                >  2:7919115/1‑93 (MQ=255)
                                                                           aCGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAactcact              <  2:8928370/93‑1 (MQ=255)
                                                                             gATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTcg            >  2:1723625/1‑93 (MQ=255)
                                                                              aTCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTcgc           >  1:3237781/1‑93 (MQ=255)
                                                                              aTCGGTGCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTcgc           <  1:9029851/93‑1 (MQ=255)
                                                                              aTCGGTGCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTcgc           <  1:5184613/93‑1 (MQ=255)
                                                                                cggttcggttGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCgg         <  1:1814298/93‑1 (MQ=255)
                                                                                cggttcggttCCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCgg         <  2:5689199/93‑1 (MQ=255)
                                                                                 ggttcggttGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGCCCACCGCCAGCGTATGCTCCAGCTGCAGATAGCCAAACTCACTCGCGGt        <  2:4516359/93‑1 (MQ=255)
                                                                                  gttcggttGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTg       >  1:9159903/1‑93 (MQ=255)
                                                                                  gttcggttGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTg       <  1:3293940/93‑1 (MQ=255)
                                                                                  gtgcggttGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTg       <  2:5218580/90‑1 (MQ=255)
                                                                                   ttcggttGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGc      <  1:5264944/93‑1 (MQ=255)
                                                                                   ttcggttGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGc      <  1:7336217/93‑1 (MQ=255)
                                                                                   ttcggttGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGc      >  1:3961232/1‑93 (MQ=255)
                                                                                   ttcggttGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGc      >  2:9555439/1‑93 (MQ=255)
                                                                                     cggttACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCgg    <  1:55410/93‑1 (MQ=255)
                                                                                      ggttGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCggg   >  1:9425374/1‑93 (MQ=255)
                                                                                       gttGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGc  <  2:3661678/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGC  >  NC_000913/2887142‑2887321

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAG  >  NC_000913/2887139‑2887325
                                                                                             |                                                                                             
CCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTG                                                                                               <  SRR4302540.1344368/94‑1 (MQ=60)
 CCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGC                                                                                              <  SRR4302540.7633181/94‑1 (MQ=60)
  CCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCC                                                                                             >  SRR4302540.1104414/1‑94 (MQ=60)
    GCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCC                                                                                           <  SRR4302540.9382200/94‑1 (MQ=60)
    GCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCC                                                                                           <  SRR4302540.139188/94‑1 (MQ=60)
    GCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCC                                                                                           <  SRR4302540.8410886/94‑1 (MQ=60)
     CGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCA                                                                                          <  SRR4302540.5816182/94‑1 (MQ=60)
     CGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCA                                                                                          >  SRR4302540.3582902/1‑94 (MQ=60)
      GCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAG                                                                                         >  SRR4302540.2116099/1‑94 (MQ=60)
       CGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGT                                                                                        >  SRR4302540.159873/1‑94 (MQ=60)
        GACCCTCCATTTCCGCTTGAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTC                                                                                       <  SRR4302540.7663204/94‑1 (MQ=60)
           GCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACG                                                                                    <  SRR4302540.2963034/94‑1 (MQ=60)
           GCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACG                                                                                    >  SRR4302540.2144166/1‑94 (MQ=60)
           GCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACG                                                                                    <  SRR4302540.8487061/94‑1 (MQ=60)
             TCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCT                                                                                  >  SRR4302540.7210816/1‑94 (MQ=60)
             TCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCT                                                                                  <  SRR4302540.4639096/94‑1 (MQ=60)
             TCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCT                                                                                  >  SRR4302540.392995/1‑94 (MQ=60)
             TCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCT                                                                                  <  SRR4302540.10120496/94‑1 (MQ=60)
             CCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCT                                                                                  <  SRR4302540.2232986/94‑1 (MQ=60)
              CCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTG                                                                                 >  SRR4302540.2988315/1‑94 (MQ=60)
               CACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCCGTCACGCTGA                                                                                >  SRR4302540.3979457/1‑94 (MQ=60)
               CACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTag                                                                                >  SRR4302540.1307362/1‑92 (MQ=60)
               CACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGA                                                                                >  SRR4302540.7948133/1‑94 (MQ=60)
                CCTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAG                                                                               <  SRR4302540.1036118/94‑1 (MQ=60)
                 CTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGT                                                                              >  SRR4302540.4664861/1‑94 (MQ=60)
                  TTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTT                                                                             >  SRR4302540.7383299/1‑94 (MQ=60)
                  TTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTT                                                                             <  SRR4302540.7466291/94‑1 (MQ=60)
                     CGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTC                                                                          <  SRR4302540.2665543/94‑1 (MQ=60)
                     CGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTC                                                                          <  SRR4302540.4395076/94‑1 (MQ=60)
                     CGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTC                                                                          >  SRR4302540.906381/1‑94 (MQ=60)
                       CTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCAC                                                                        >  SRR4302540.3542296/1‑94 (MQ=60)
                       CTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCAC                                                                        >  SRR4302540.3128880/1‑94 (MQ=60)
                       CTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCAC                                                                        <  SRR4302540.7079174/94‑1 (MQ=60)
                       CTTTAAACGTCGCCACCCCCACGCGCTCCAGCGTGTATTTGGTTTTTGCATTTTTACGATCGGTTCGGTTTCCCCAGTCACGCTGAGTTGTCAC                                                                        <  SRR4302540.9129569/94‑1 (MQ=60)
                        TTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACG                                                                       <  SRR4302540.6103014/94‑1 (MQ=60)
                          TAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGAC                                                                     <  SRR4302540.6984098/94‑1 (MQ=60)
                          TAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGAC                                                                     >  SRR4302540.2184865/1‑94 (MQ=60)
                          TAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGAC                                                                     <  SRR4302540.1848442/94‑1 (MQ=60)
                          TAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGAC                                                                     >  SRR4302540.8218675/1‑94 (MQ=60)
                           AAACGTCTCACCCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACG                                                                    <  SRR4302540.5866305/94‑1 (MQ=60)
                           AAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACG                                                                    >  SRR4302540.9773998/1‑94 (MQ=60)
                           AAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACG                                                                    <  SRR4302540.6620674/94‑1 (MQ=60)
                           AAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACG                                                                    <  SRR4302540.6386778/94‑1 (MQ=60)
                           AAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACG                                                                    <  SRR4302540.1499971/94‑1 (MQ=60)
                            AACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGG                                                                   <  SRR4302540.3678974/94‑1 (MQ=60)
                             ACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGC                                                                  <  SRR4302540.5191241/94‑1 (MQ=60)
                              CGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCT                                                                 <  SRR4302540.7406809/94‑1 (MQ=60)
                               GTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTC                                                                >  SRR4302540.7295819/1‑94 (MQ=60)
                                 CTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCG                                                              >  SRR4302540.3968824/1‑94 (MQ=60)
                                  TCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGG                                                             >  SRR4302540.9686199/1‑94 (MQ=60)
                                   CACCCCCCCCGCGCCCCAGCGGGTATTTGGTTGTTGCATTTGTTCGATCCGTTCGGTTACCCCAGTCCCGCGGAGTTGTCACGACGGCTTCGGC                                                            <  SRR4302540.7684868/94‑1 (MQ=36)
                                   CAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGTC                                                            >  SRR4302540.7408965/1‑94 (MQ=60)
                                   CAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGC                                                            >  SRR4302540.7408964/1‑94 (MQ=60)
                                   CAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCACCAGTCACGCTGAGTTGTCACGACGGCTTCGGC                                                            >  SRR4302540.9804111/1‑94 (MQ=60)
                                    AACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCC                                                           <  SRR4302540.8971246/94‑1 (MQ=60)
                                     ACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCA                                                          >  SRR4302540.1429441/1‑94 (MQ=60)
                                       CCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACC                                                        >  SRR4302540.856965/1‑94 (MQ=60)
                                        CCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCG                                                       >  SRR4302540.4547962/1‑94 (MQ=60)
                                        CCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCG                                                       >  SRR4302540.4463253/1‑94 (MQ=60)
                                        CCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCG                                                       <  SRR4302540.6220390/94‑1 (MQ=60)
                                         CCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGC                                                      >  SRR4302540.841032/1‑94 (MQ=60)
                                         CCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGC                                                      >  SRR4302540.3695430/1‑94 (MQ=60)
                                         CCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGC                                                      >  SRR4302540.9505161/1‑94 (MQ=60)
                                          CACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGATGGCTTCGGCCACCGCC                                                     >  SRR4302540.7076350/1‑94 (MQ=60)
                                          CACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCC                                                     >  SRR4302540.7715797/1‑94 (MQ=60)
                                            CGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAG                                                   <  SRR4302540.3008723/94‑1 (MQ=60)
                                             CCGCTCCAGCGTTTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGC                                                  <  SRR4302540.8622181/94‑1 (MQ=60)
                                               GCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATTGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGT                                                <  SRR4302540.6903567/94‑1 (MQ=60)
                                                 TCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTAT                                              <  SRR4302540.8368692/94‑1 (MQ=60)
                                                 TCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTAT                                              <  SRR4302540.3999984/94‑1 (MQ=60)
                                                 TCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTAT                                              <  SRR4302540.6357676/94‑1 (MQ=60)
                                                  CCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATG                                             <  SRR4302540.5732126/94‑1 (MQ=60)
                                                   CAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGC                                            <  SRR4302540.8595300/94‑1 (MQ=60)
                                                    AGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCT                                           >  SRR4302540.7593689/1‑94 (MQ=60)
                                                      CGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCC                                         <  SRR4302540.6554805/94‑1 (MQ=60)
                                                      CGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCC                                         >  SRR4302540.2903933/1‑94 (MQ=60)
                                                      CGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCC                                         <  SRR4302540.1989191/94‑1 (MQ=60)
                                                      CGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCC                                         >  SRR4302540.8095755/1‑94 (MQ=60)
                                                       GTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCA                                        <  SRR4302540.7174471/94‑1 (MQ=60)
                                                       GTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCA                                        <  SRR4302540.7642224/94‑1 (MQ=60)
                                                       GTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCA                                        >  SRR4302540.9832430/1‑94 (MQ=60)
                                                       GTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCA                                        >  SRR4302540.5260031/1‑94 (MQ=60)
                                                       GTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGACACCGACAGCGTATGCTCCC                                        >  SRR4302540.484313/1‑94 (MQ=60)
                                                        TGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAG                                       >  SRR4302540.2684607/1‑94 (MQ=60)
                                                        TGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAG                                       <  SRR4302540.3458125/94‑1 (MQ=60)
                                                        TGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAG                                       >  SRR4302540.9361036/1‑94 (MQ=60)
                                                         GTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGG                                      >  SRR4302540.2336915/1‑94 (MQ=60)
                                                         GTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGC                                      >  SRR4302540.4923910/1‑94 (MQ=60)
                                                          TATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCG                                     <  SRR4302540.3485970/94‑1 (MQ=60)
                                                          TATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCG                                     >  SRR4302540.7849867/1‑94 (MQ=60)
                                                          TATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCG                                     >  SRR4302540.1484185/1‑94 (MQ=60)
                                                           ATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGG                                    <  SRR4302540.12583/94‑1 (MQ=60)
                                                           ATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGG                                    >  SRR4302540.1857480/1‑94 (MQ=60)
                                                           ATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTAACCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACAGCCAGAGTATGCTCCAGCGG                                    >  SRR4302540.6896228/1‑94 (MQ=60)
                                                            TTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGC                                   >  SRR4302540.6654061/1‑94 (MQ=60)
                                                            TTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGC                                   >  SRR4302540.3274894/1‑94 (MQ=60)
                                                              TGGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAG                                 >  SRR4302540.3705996/1‑94 (MQ=60)
                                                               GGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGA                                >  SRR4302540.8754480/1‑94 (MQ=60)
                                                               GGTTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCAACCGCCAGCGTATGCTCCACCGGCAGA                                >  SRR4302540.8754481/1‑94 (MQ=60)
                                                                GTTGTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGAT                               <  SRR4302540.9655805/94‑1 (MQ=60)
                                                                 TTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATA                              <  SRR4302540.7082682/94‑1 (MQ=60)
                                                                 GTTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCGCGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATA                              <  SRR4302540.1618591/94‑1 (MQ=60)
                                                                  TTTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAG                             <  SRR4302540.4601306/94‑1 (MQ=60)
                                                                   TTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGC                            <  SRR4302540.9236443/94‑1 (MQ=60)
                                                                   TTGGCATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGC                            >  SRR4302540.6343838/1‑94 (MQ=60)
                                                                       CATTTTTACGATCGGTTCGGTTGCCGCAGTCTCGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAA                        <  SRR4302540.2161750/94‑1 (MQ=60)
                                                                       CATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAA                        >  SRR4302540.5940822/1‑94 (MQ=60)
                                                                        ATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAAC                       <  SRR4302540.5144497/94‑1 (MQ=60)
                                                                        ATTTTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAAC                       >  SRR4302540.5473927/1‑94 (MQ=60)
                                                                        ATTTTTACGATCGGTTCGGGTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAAC                       <  SRR4302540.3348469/94‑1 (MQ=60)
                                                                           TTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGTCAGATAGCCAAACTCA                    >  SRR4302540.9028481/1‑94 (MQ=60)
                                                                           TTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCA                    <  SRR4302540.580290/94‑1 (MQ=60)
                                                                           TTTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCA                    >  SRR4302540.8795025/1‑94 (MQ=60)
                                                                            TTACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC                   <  SRR4302540.2891530/94‑1 (MQ=60)
                                                                              ACGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTC                 >  SRR4302540.7273617/1‑94 (MQ=60)
                                                                               CGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCG                <  SRR4302540.9278870/94‑1 (MQ=60)
                                                                               CGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCG                <  SRR4302540.1485405/94‑1 (MQ=60)
                                                                               CGATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCG                >  SRR4302540.7365337/1‑94 (MQ=60)
                                                                                GATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGC               >  SRR4302540.3273000/1‑94 (MQ=60)
                                                                                 ATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCCGCGGCAGATAGCCACACTCACTCGCG              >  SRR4302540.929841/1‑94 (MQ=60)
                                                                                 ATCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCG              <  SRR4302540.3544428/94‑1 (MQ=60)
                                                                                 ATCGGTGCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCG              <  SRR4302540.9138212/94‑1 (MQ=60)
                                                                                 ATCGGTGCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCG              <  SRR4302540.5241051/94‑1 (MQ=60)
                                                                                  acaGTTCGGTTCCCCCAGTCACGATGAGTTGACACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGG             <  SRR4302540.5783463/91‑1 (MQ=60)
                                                                                  TCGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGG             <  SRR4302540.3663996/94‑1 (MQ=60)
                                                                                   CGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGT            <  SRR4302540.922647/94‑1 (MQ=60)
                                                                                   CGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGT            <  SRR4302540.7030436/94‑1 (MQ=60)
                                                                                   CGGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGT            <  SRR4302540.1834107/94‑1 (MQ=60)
                                                                                   CGGTTCGGTTGCCCCAGTCACGCTGAGTTCTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGT            <  SRR4302540.600896/94‑1 (MQ=60)
                                                                                    GGTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTG           >  SRR4302540.9269844/1‑94 (MQ=60)
                                                                                     GTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTcg          >  SRR4302540.7615532/1‑92 (MQ=60)
                                                                                     GTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGC          >  SRR4302540.4003955/1‑94 (MQ=60)
                                                                                     GTTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGC          <  SRR4302540.3329984/94‑1 (MQ=60)
                                                                                     GTTCGGTTGCCCCAGTCACGCTGAGTGGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGC          <  SRR4302540.4616901/94‑1 (MQ=60)
                                                                                      TTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCG         <  SRR4302540.5322398/94‑1 (MQ=60)
                                                                                      TTCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCG         <  SRR4302540.7418093/94‑1 (MQ=60)
                                                                                      TTCGGTTGCCCCACTCACGCTCAGTTGTCACGACGGTTTCGGCCCCCGCCAGCGTATGCTCGAGCGGCAGATAGCCAAACTCACTCGCGGTGCG         <  SRR4302540.5322397/94‑1 (MQ=60)
                                                                                      TTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCG         <  SRR4302540.3307041/94‑1 (MQ=60)
                                                                                       TCGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGG        >  SRR4302540.3387178/1‑94 (MQ=60)
                                                                                        CGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGG       >  SRR4302540.7158910/1‑94 (MQ=60)
                                                                                        CGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGG       <  SRR4302540.536128/94‑1 (MQ=60)
                                                                                        CGGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGG       >  SRR4302540.9537371/1‑94 (MQ=60)
                                                                                        CGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGG       <  SRR4302540.55966/94‑1 (MQ=60)
                                                                                         GGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGG      >  SRR4302540.7421052/1‑94 (MQ=60)
                                                                                         GGTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGC      >  SRR4302540.5846160/1‑94 (MQ=60)
                                                                                          GTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCG     >  SRR4302540.2674014/1‑94 (MQ=60)
                                                                                          GTTGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCG     >  SRR4302540.7706330/1‑94 (MQ=60)
                                                                                            TGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTA   <  SRR4302540.1194854/94‑1 (MQ=60)
                                                                                            TGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTA   >  SRR4302540.9308940/1‑94 (MQ=60)
                                                                                            TGCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTA   >  SRR4302540.3988391/1‑94 (MQ=60)
                                                                                             GCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTCG  >  SRR4302540.4852733/1‑94 (MQ=60)
                                                                                             GCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAG  >  SRR4302540.6122381/1‑94 (MQ=60)
                                                                                             GCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGCATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAG  >  SRR4302540.873875/1‑94 (MQ=60)
                                                                                             |                                                                                             
CCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAG  >  NC_000913/2887139‑2887325

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 2 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑