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 1061 294.1 15065038 98.6% 14854127 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,077,514 Δ1 bp intergenic (‑165/+149) cmtB ← / ← tktA PTS system mannitol‑specific enzyme IIA component/transketolase 1 isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,077,5130C.100.0% 868.2 / NA 187intergenic (‑164/+150)cmtB/tktAPTS system mannitol‑specific enzyme IIA component/transketolase 1 isozyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (103/84);  total (103/84)

Only 100 of 193 total aligned reads displayed.

GCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCG  >  NC_000913/3077421‑3077600
                                                                                            |                                                                                       
gCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAc                                                                                         >  2:1027798/1‑93 (MQ=255)
   aTTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAcgca                                                                                      >  1:2802348/1‑90 (MQ=255)
   aTTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAcgca                                                                                      <  2:5366017/93‑4 (MQ=255)
      gAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCatat                                                                                  <  2:6263469/93‑1 (MQ=255)
       aGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATc                                                                                 <  1:4588532/93‑1 (MQ=255)
       aGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATc                                                                                 <  2:520339/93‑1 (MQ=255)
          aGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTg                                                                              <  2:405737/93‑1 (MQ=255)
          aGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTg                                                                              >  2:6601030/1‑93 (MQ=255)
          aGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTg                                                                              >  2:5883991/1‑93 (MQ=255)
           gggAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGa                                                                             <  2:7018932/93‑1 (MQ=255)
              aaaCCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGc                                                                          <  2:3153574/93‑1 (MQ=255)
              aaaCCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGc                                                                          <  1:3846136/93‑1 (MQ=255)
              aaaCCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGc                                                                          >  1:4767856/1‑93 (MQ=255)
              aaaCCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGc                                                                          <  1:2184619/93‑1 (MQ=255)
               aaCCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCa                                                                         <  2:5039122/93‑1 (MQ=255)
               aaCCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCa                                                                         >  2:1073698/1‑93 (MQ=255)
                 ccGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCata                                                                       <  1:6366685/93‑1 (MQ=255)
                  cGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCatat                                                                      >  1:4190115/1‑93 (MQ=255)
                  cGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCatat                                                                      <  1:5588525/93‑1 (MQ=255)
                  cGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCatat                                                                      <  2:4520584/93‑1 (MQ=255)
                   ggTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATAtt                                                                     <  2:757141/93‑1 (MQ=255)
                    gtgtTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTa                                                                    >  1:3870235/1‑93 (MQ=255)
                    gtgtTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTa                                                                    >  2:2960138/1‑93 (MQ=255)
                    gtgtTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTa                                                                    >  2:4732658/1‑93 (MQ=255)
                     tgtTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTAc                                                                   <  1:1106083/93‑1 (MQ=255)
                      gtTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACa                                                                  <  2:1173015/93‑1 (MQ=255)
                      gtTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACa                                                                  <  2:6189345/93‑1 (MQ=255)
                      gtTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACa                                                                  <  2:6385619/93‑1 (MQ=255)
                         ttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATaa                                                               <  2:3829601/93‑1 (MQ=255)
                            aaaaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTg                                                            >  1:1416879/1‑93 (MQ=255)
                             aaaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGa                                                           <  2:1432084/93‑1 (MQ=255)
                             aaaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGa                                                           >  1:4056445/1‑93 (MQ=255)
                              aaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGat                                                          <  1:1532606/93‑1 (MQ=255)
                              aaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGat                                                          >  1:1620019/1‑93 (MQ=255)
                                cAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGatat                                                        >  1:680258/1‑93 (MQ=255)
                                  gTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtattt                                                      >  1:1914372/1‑93 (MQ=255)
                                  gTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtattt                                                      >  2:4930028/1‑93 (MQ=255)
                                   tGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtattta                                                     <  1:7391789/93‑1 (MQ=255)
                                     aaTTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttatt                                                   <  2:2622556/93‑1 (MQ=255)
                                     aaTTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttatt                                                   <  1:2615987/93‑1 (MQ=255)
                                     aaTTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttatt                                                   <  1:701542/93‑1 (MQ=255)
                                       ttattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGc                                                 >  2:2260078/1‑93 (MQ=255)
                                       ttattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGc                                                 <  2:6684601/93‑1 (MQ=255)
                                       ttattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGc                                                 >  1:6231266/1‑93 (MQ=255)
                                       ttattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGc                                                 >  1:3012397/1‑93 (MQ=255)
                                        tattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCa                                                >  1:7194457/1‑93 (MQ=255)
                                        tattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCa                                                <  1:2858300/93‑1 (MQ=255)
                                        tattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCa                                                >  2:6622041/1‑93 (MQ=255)
                                        tattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCa                                                <  1:5426196/93‑1 (MQ=255)
                                          ttaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCata                                              >  1:2318657/1‑93 (MQ=255)
                                           taATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCatat                                             >  2:7387056/1‑93 (MQ=255)
                                           taATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCatat                                             <  2:2531590/93‑1 (MQ=255)
                                             aTTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaa                                           >  2:6766666/1‑93 (MQ=255)
                                              tttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaaa                                          >  2:6585227/1‑93 (MQ=255)
                                               ttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAat                                         <  2:4194210/93‑1 (MQ=255)
                                                tCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                                        <  1:5953696/93‑1 (MQ=255)
                                                tCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                                        >  2:6297684/1‑93 (MQ=255)
                                                  ttATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAtt                                      <  2:6088092/93‑1 (MQ=255)
                                                   tatataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAttt                                     <  2:5460348/93‑1 (MQ=255)
                                                    atataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtg                                    >  2:2696132/1‑93 (MQ=255)
                                                    atataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtg                                    <  1:4251418/93‑1 (MQ=255)
                                                     tataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtgt                                   >  2:1643555/1‑93 (MQ=255)
                                                     tataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtgt                                   >  2:7143728/1‑93 (MQ=255)
                                                       taACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGa                                 <  1:3385939/93‑1 (MQ=255)
                                                       taACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGa                                 >  1:5270556/1‑93 (MQ=255)
                                                        aaCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAt                                >  1:667020/1‑93 (MQ=255)
                                                        aaCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAt                                >  2:4977262/1‑93 (MQ=255)
                                                         aCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATc                               >  2:3615124/1‑93 (MQ=255)
                                                          cATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCt                              >  2:1124010/1‑93 (MQ=255)
                                                           aTTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTa                             >  1:5474906/1‑93 (MQ=255)
                                                           aTTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTa                             >  2:4185301/1‑93 (MQ=255)
                                                             tGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTaca                           >  1:4928148/1‑93 (MQ=255)
                                                              gAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTacaa                          >  2:3199206/1‑93 (MQ=255)
                                                               aGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTacaac                         <  1:6525553/93‑1 (MQ=255)
                                                                gTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTacaaca                        <  2:510757/93‑1 (MQ=255)
                                                                 ttATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAacac                       >  2:5600547/1‑93 (MQ=255)
                                                                 ttATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAacac                       >  2:3587293/1‑93 (MQ=255)
                                                                  tataGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACg                      >  1:3704135/1‑93 (MQ=255)
                                                                   ataGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGc                     >  1:3884250/1‑93 (MQ=255)
                                                                   ataGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGc                     <  2:7224981/93‑1 (MQ=255)
                                                                   ataGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGc                     >  2:3762422/1‑93 (MQ=255)
                                                                    taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGcc                    >  1:5546943/1‑93 (MQ=255)
                                                                    taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGcc                    >  1:7421074/1‑93 (MQ=255)
                                                                      gATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCtt                  >  2:1877781/1‑93 (MQ=255)
                                                                        tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTAt                >  2:2068776/1‑93 (MQ=255)
                                                                         ataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATc               >  1:5942478/1‑93 (MQ=255)
                                                                          taACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCt              >  2:2630296/1‑93 (MQ=255)
                                                                           aacaacAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTa             >  1:2711596/1‑93 (MQ=255)
                                                                            acaacaAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTAt            >  1:804956/1‑93 (MQ=255)
                                                                              aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTg          >  1:4170849/1‑93 (MQ=255)
                                                                               acaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGc         >  1:534941/1‑93 (MQ=255)
                                                                                 aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCtt       <  1:4904313/93‑1 (MQ=255)
                                                                                 aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCtt       <  1:3668491/93‑1 (MQ=255)
                                                                                  aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTg      >  1:2095399/1‑93 (MQ=255)
                                                                                    tGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTc    >  2:1585342/1‑93 (MQ=255)
                                                                                    tGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTc    <  2:3887766/93‑1 (MQ=255)
                                                                                     gTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTcc   >  2:3292635/1‑93 (MQ=255)
                                                                                      tCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCg  >  2:6634594/1‑93 (MQ=255)
                                                                                      tCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCg  >  2:3531191/1‑93 (MQ=255)
                                                                                      tCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCg  >  1:5908924/1‑93 (MQ=255)
                                                                                            |                                                                                       
GCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCG  >  NC_000913/3077421‑3077600

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGG  >  NC_000913/3077423‑3077607
                                                                                           |                                                                                             
AATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCA                                                                                            >  SRR4302542.2829331/1‑94 (MQ=60)
     AGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCG                                                                                       <  SRR4302542.4633409/94‑1 (MQ=60)
     AGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCG                                                                                       >  SRR4302542.976981/1‑94 (MQ=60)
     AGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCG                                                                                       >  SRR4302542.2604565/1‑94 (MQ=60)
      GCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGT                                                                                      <  SRR4302542.6725167/94‑1 (MQ=60)
       CAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTG                                                                                     <  SRR4302542.4105500/94‑1 (MQ=60)
         GGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGAT                                                                                   <  SRR4302542.1561704/94‑1 (MQ=60)
           GAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGC                                                                                 >  SRR4302542.4814331/1‑94 (MQ=60)
            AAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCA                                                                                <  SRR4302542.3884144/94‑1 (MQ=60)
            AAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCA                                                                                <  SRR4302542.2203440/94‑1 (MQ=60)
               CCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATAT                                                                             <  SRR4302542.6429385/94‑1 (MQ=60)
               CCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATAT                                                                             >  SRR4302542.4231434/1‑94 (MQ=60)
                CGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATT                                                                            <  SRR4302542.5642827/94‑1 (MQ=60)
                 GGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTA                                                                           >  SRR4302542.3908471/1‑94 (MQ=60)
                   TGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACA                                                                         <  SRR4302542.1115541/94‑1 (MQ=60)
                   TGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACA                                                                         <  SRR4302542.3885247/94‑1 (MQ=60)
                    GTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACAT                                                                        <  SRR4302542.4233116/94‑1 (MQ=60)
                     TTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATA                                                                       <  SRR4302542.3295849/94‑1 (MQ=60)
                         GAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTG                                                                   >  SRR4302542.1429146/1‑94 (MQ=60)
                          AAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGA                                                                  >  SRR4302542.4096449/1‑94 (MQ=60)
                           AAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAT                                                                 >  SRR4302542.1634014/1‑94 (MQ=60)
                            AACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATA                                                                <  SRR4302542.1545842/94‑1 (MQ=60)
                            AACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATA                                                                >  SRR4302542.1438647/1‑94 (MQ=60)
                            AACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATA                                                                <  SRR4302542.7063068/94‑1 (MQ=60)
                             ACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATAT                                                               >  SRR4302542.686132/1‑94 (MQ=60)
                              CAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATT                                                              <  SRR4302542.5862734/94‑1 (MQ=60)
                               AGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTT                                                             <  SRR4302542.4672301/94‑1 (MQ=60)
                               AGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTT                                                             >  SRR4302542.1930917/1‑94 (MQ=60)
                               AGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTT                                                             >  SRR4302542.6348818/1‑94 (MQ=60)
                                 TGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT                                                           <  SRR4302542.7464142/94‑1 (MQ=60)
                                 TGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT                                                           >  SRR4302542.5368740/1‑94 (MQ=60)
                                 TGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT                                                           <  SRR4302542.390461/94‑1 (MQ=60)
                                   AATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTG                                                         <  SRR4302542.707618/94‑1 (MQ=60)
                                   AATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTG                                                         <  SRR4302542.2641260/94‑1 (MQ=60)
                                    ATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGC                                                        >  SRR4302542.6292619/1‑94 (MQ=60)
                                    ATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGC                                                        >  SRR4302542.3041192/1‑94 (MQ=60)
                                     TTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCA                                                       >  SRR4302542.7265062/1‑94 (MQ=60)
                                      TATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCAT                                                      <  SRR4302542.5478924/94‑1 (MQ=60)
                                      TATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCAT                                                      <  SRR4302542.2885764/94‑1 (MQ=60)
                                       ATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATA                                                     >  SRR4302542.2338542/1‑94 (MQ=60)
                                       ATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATA                                                     >  SRR4302542.2873869/1‑94 (MQ=60)
                                         TAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATA                                                   >  SRR4302542.5538604/1‑94 (MQ=60)
                                          AATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAA                                                  >  SRR4302542.7485570/1‑94 (MQ=60)
                                              TCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAT                                              <  SRR4302542.6011356/94‑1 (MQ=60)
                                              TCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAT                                              <  SRR4302542.6011355/94‑1 (MQ=60)
                                              TCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAT                                              >  SRR4302542.5996893/1‑94 (MQ=60)
                                              TCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAT                                              >  SRR4302542.7187803/1‑94 (MQ=60)
                                                 TATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTG                                           <  SRR4302542.2969598/94‑1 (MQ=60)
                                                  ATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGT                                          >  SRR4302542.1850001/1‑94 (MQ=60)
                                                  ATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGT                                          <  SRR4302542.4293298/94‑1 (MQ=60)
                                                   TATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTG                                         >  SRR4302542.7581007/1‑94 (MQ=60)
                                                   TATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTG                                         <  SRR4302542.3699582/94‑1 (MQ=60)
                                                    ATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGA                                        >  SRR4302542.5321773/1‑94 (MQ=60)
                                                    ATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGA                                        <  SRR4302542.4593241/94‑1 (MQ=60)
                                                    ATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGA                                        <  SRR4302542.5603085/94‑1 (MQ=60)
                                                     TAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAT                                       <  SRR4302542.4925668/94‑1 (MQ=60)
                                                     TAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAT                                       >  SRR4302542.672783/1‑94 (MQ=60)
                                                     TAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAT                                       <  SRR4302542.6738147/94‑1 (MQ=60)
                                                     TAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAT                                       <  SRR4302542.3417727/94‑1 (MQ=60)
                                                      AACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATC                                      >  SRR4302542.4970801/1‑94 (MQ=60)
                                                        CATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTA                                    >  SRR4302542.5528095/1‑94 (MQ=60)
                                                         ATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTAC                                   <  SRR4302542.5164002/94‑1 (MQ=60)
                                                          TTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACA                                  >  SRR4302542.4976053/1‑94 (MQ=60)
                                                             AGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACA                               <  SRR4302542.6589856/94‑1 (MQ=60)
                                                             AGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACA                               >  SRR4302542.3927503/1‑94 (MQ=60)
                                                               TTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACG                             <  SRR4302542.4221815/94‑1 (MQ=60)
                                                               TTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACG                             >  SRR4302542.3740787/1‑94 (MQ=60)
                                                                TATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGC                            >  SRR4302542.2315942/1‑94 (MQ=60)
                                                                TATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGC                            >  SRR4302542.3922609/1‑94 (MQ=60)
                                                                 ATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCC                           >  SRR4302542.5600848/1‑94 (MQ=60)
                                                                 ATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCC                           >  SRR4302542.266913/1‑94 (MQ=60)
                                                                 ATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCC                           >  SRR4302542.7493685/1‑94 (MQ=60)
                                                                      TATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATC                      >  SRR4302542.6000046/1‑94 (MQ=60)
                                                                      TATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATC                      >  SRR4302542.3841139/1‑94 (MQ=60)
                                                                        TAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTA                    >  SRR4302542.2818121/1‑94 (MQ=60)
                                                                        TAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTA                    >  SRR4302542.2737749/1‑94 (MQ=60)
                                                                         AACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTAT                   <  SRR4302542.6472140/94‑1 (MQ=60)
                                                                         AACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTAT                   >  SRR4302542.811878/1‑94 (MQ=60)
                                                                          ACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATT                  >  SRR4302542.2190265/1‑94 (MQ=60)
                                                                           CAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTG                 >  SRR4302542.4211984/1‑94 (MQ=60)
                                                                            AACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGC                >  SRR4302542.5494884/1‑94 (MQ=60)
                                                                            AACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGC                >  SRR4302542.539563/1‑94 (MQ=60)
                                                                              CAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTT              <  SRR4302542.203137/94‑1 (MQ=60)
                                                                               AAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTG             >  SRR4302542.2113477/1‑94 (MQ=60)
                                                                               AAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTG             <  SRR4302542.4952017/94‑1 (MQ=60)
                                                                               AAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTG             <  SRR4302542.3703597/94‑1 (MQ=60)
                                                                                  TGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCC          <  SRR4302542.6185389/94‑1 (MQ=60)
                                                                                   GTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCG         >  SRR4302542.5966191/1‑94 (MQ=60)
                                                                                      AATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAA      >  SRR4302542.3199529/1‑94 (MQ=60)
                                                                                         taCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAAC    <  SRR4302542.2486739/92‑1 (MQ=60)
                                                                                           CGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGG  <  SRR4302542.4787068/94‑1 (MQ=60)
                                                                                           CGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGG  >  SRR4302542.6363285/1‑94 (MQ=60)
                                                                                           |                                                                                             
AATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGG  >  NC_000913/3077423‑3077607

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