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 F50 I1 R1 775 316.5 10245680 97.2% 9958800 138.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP025268 3,728,808 (C)10→8 pseudogene (1254‑1255/1457 nt) xylB ← xylulokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0252683,728,8000C.98.8% 389.0 / ‑2.7 83pseudogene (1263/1457 nt)xylBxylulokinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (49/33);  total (50/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
*NZ_CP0252683,728,8010C.95.1% 344.7 / 14.5 83pseudogene (1262/1457 nt)xylBxylulokinase
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base . (44/33);  total (50/33)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCC  >  NZ_CP025268/3728661‑3728932
                                                                                                                                           ||                                                                                                                                   
ggCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATc                                                                                                                                      <  1:1187416/140‑1 (MQ=255)
ggCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATc                                                                                                                                      >  1:257769/1‑140 (MQ=255)
   tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccc                                                                                                                                   <  1:2414153/140‑1 (MQ=255)
     aTAAGCGGCCTAACGCTGCGCATCTGGTAGAGGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccac                                                                                                                                 <  2:4389417/140‑3 (MQ=255)
     aTAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccc                                                                                                                                 <  2:975242/140‑1 (MQ=255)
      tAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccc                                                                                                                                >  2:975692/1‑140 (MQ=255)
      tAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccc                                                                                                                                >  2:1762623/1‑140 (MQ=255)
       aaGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccc                                                                                                                               >  1:1701202/1‑140 (MQ=255)
          cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccgta                                                                                                                            <  2:2553411/140‑4 (MQ=255)
          cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccgta                                                                                                                            <  2:2378003/140‑4 (MQ=255)
          cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccgta                                                                                                                            >  1:2553411/1‑137 (MQ=255)
          cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccgta                                                                                                                            >  2:1207932/1‑137 (MQ=255)
           ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccgtac                                                                                                                           >  2:2950362/1‑136 (MQ=255)
               tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTa                                                                                                                     >  2:5025937/1‑140 (MQ=255)
                    cTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGa                                                                                                                >  2:633482/1‑140 (MQ=255)
                    cTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGGAATCGa                                                                                                                >  2:193898/1‑140 (MQ=255)
                       cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAgct                                                                                                             >  1:3295992/1‑140 (MQ=255)
                       cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAgct                                                                                                             >  1:1746508/1‑140 (MQ=255)
                       cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGCACGGTACTCGAgct                                                                                                             >  2:4843947/1‑140 (MQ=255)
                         cATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAgctgc                                                                                                           >  2:1475808/1‑140 (MQ=255)
                            cTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGa                                                                                                        <  1:3449543/140‑1 (MQ=255)
                            cTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGa                                                                                                        >  1:4465510/1‑140 (MQ=255)
                               gTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTGCGGTTATCGAGCCGCTGACCg                                                                                                     >  1:4095483/1‑140 (MQ=255)
                               gTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAACCGAGCTGCTGACCg                                                                                                     >  2:4690000/1‑140 (MQ=255)
                                 aGATGCGAGTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTGGCCGCGATCTGCGCCAGCTTTGCTGCGCCCAGTGCGGGCCCCACAT‑‑GCCCCCCTGTACTGTTATCGAGCTGTTGACGGcg                                                                                                   >  2:3061135/1‑139 (MQ=255)
                                    tGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCCGACCGTAATCGAGCTGCTGAgcgcctat                                                                                                >  1:1945183/1‑136 (MQ=255)
                                    tGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGat                                                                                                >  2:2663343/1‑140 (MQ=255)
                                    tGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGat                                                                                                >  2:3758263/1‑140 (MQ=255)
                                    tGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACCGTAATCGAGCTGCTGACCGCCGaa                                                                                                >  2:2362773/1‑139 (MQ=255)
                                      cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCTTATGGTAATCGAGCCGCTGGCCGCTGatct                                                                                              >  2:1074331/1‑138 (MQ=255)
                                      cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCCTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGatat                                                                                              <  1:4052778/140‑1 (MQ=255)
                                        aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                            >  1:1032317/1‑140 (MQ=255)
                                         cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCTCCCATTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCCTCTGACCGCGGataaacg                                                                                           >  1:2529035/1‑136 (MQ=255)
                                            ttCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                        >  2:26658/1‑140 (MQ=255)
                                            ttCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                        >  2:97819/1‑140 (MQ=255)
                                                  cGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                  <  1:1476152/140‑1 (MQ=255)
                                                  cGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                  <  1:3742087/140‑1 (MQ=255)
                                                  cGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                  <  1:3776324/140‑1 (MQ=255)
                                                  cGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                  >  2:4286965/1‑140 (MQ=255)
                                                    gTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                                <  2:2740438/140‑1 (MQ=255)
                                                      aGTTGCGGCAACGATGCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCGGGCCCCACA‑‑AACCCGCCCGTTCGGCAAACGTGACGCTGCCCGCTGATATTCGCCAGCATCTAACg                                                                              >  2:3829022/1‑140 (MQ=255)
                                                      aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGGTCTGCGCCAACCTGTCTGCGCTCAGTGCTGGCCCCACAT‑‑CCCCCCCTGTAAGCGAATCGAGCTGCTGACCGCTGATATCCGGCAGCATATGACg                                                                              >  2:1837902/1‑140 (MQ=255)
                                                      aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCC‑‑TACGGTAAACCAGCTGCTGGCCGCTGATATCCGCCAGCCGCTGACg                                                                              >  1:4992005/1‑140 (MQ=255)
                                                      aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCGAGCATCTGGCg                                                                              >  2:2668899/1‑140 (MQ=255)
                                                      aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACg                                                                              >  2:1284248/1‑140 (MQ=255)
                                                      aGTTGCGGCAACAAGTCAATGAGCGATTTCTCTGGATGCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTATGTGAGTTGCTGGCCGGTGATATCGGGGAGCATCTGGCg                                                                              >  2:2668532/1‑140 (MQ=255)
                                                       gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                             <  2:1706754/140‑1 (MQ=255)
                                                        ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGcc                                                                            <  1:3171956/140‑1 (MQ=255)
                                                        ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ACCCGCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGcc                                                                            >  1:3645245/1‑140 (MQ=255)
                                                         tGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                           <  2:3394236/140‑1 (MQ=255)
                                                         tGCGGCAACAATTCAATGAGCGATTCCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                           <  2:1742985/140‑1 (MQ=255)
                                                            ggCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTa                                                                        <  1:1109029/140‑1 (MQ=255)
                                                            ggCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAAGCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTa                                                                        <  2:2519719/140‑1 (MQ=255)
                                                              caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACt                                                                      <  2:2962401/140‑1 (MQ=255)
                                                              caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACt                                                                      <  2:866439/140‑1 (MQ=255)
                                                               aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑‑GCCCCCCCCGTACGGTAATCGAGCTGCTGACCCCTGATATCCGCCAGCACCTGACGCCCGTACTc                                                                     >  2:178911/1‑140 (MQ=255)
                                                                 caaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCGCCCGTACGGGCAACGAGCTGCTGACCGCTGAGACCCGCCAGCATATGATTGCACTACTCat                                                                   >  2:711733/1‑139 (MQ=255)
                                                                 caaTTCAATGAGCGATTTCTCTGGATTCGACGCGATCTGCGCCAGCCTTTCTGCGCCCAGTGCTGGCCCCACAT‑‑GCCCCCCCGTACGGTAATCTAGCTGCTGACCGCTGATACCCGCCATCATCTGACACGAGTACTCAc                                                                   >  2:711761/1‑140 (MQ=255)
                                                                      cAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAcgc                                                              <  1:1344354/140‑1 (MQ=255)
                                                                      cAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAcgc                                                              <  1:1344562/140‑1 (MQ=255)
                                                                        aTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACgcgc                                                            >  1:3426446/1‑140 (MQ=255)
                                                                          gcgCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGTGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTCATCTAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccc                                                          <  2:3057584/138‑1 (MQ=255)
                                                                          gAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTTATCCACCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccc                                                          <  2:3057679/140‑1 (MQ=255)
                                                                          gAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccc                                                          <  1:2606210/140‑1 (MQ=255)
                                                                          gAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccc                                                          <  2:1847851/140‑1 (MQ=255)
                                                                           aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGGACGGTAATCGAGCTGCTGACCGCTGATATCCCCCAGCCTCTGACGCCCTTCCCCACCACGcgcccc                                                         >  2:3147191/1‑140 (MQ=255)
                                                                             cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgc                                                       >  1:1892119/1‑140 (MQ=255)
                                                                             cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCCAGCTGTTGGCCGCTTAAATCCGCCAGGCTCTAACGCGAGTACCCACTACGcgccccgc                                                       >  1:1892128/1‑140 (MQ=255)
                                                                                    tctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt                                                <  1:4286965/140‑1 (MQ=255)
                                                                                     ctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGCGCTTCTGACCGCTGATACCCGCCAGCATCTGCCGCCCGTACTCACTACGCGCCCCGCCCCCAATc                                               >  2:1200893/1‑140 (MQ=255)
                                                                                                   ctCTGCGCCAGCCTTGATGCGCCCAGTGCTGGCCCCCCAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTg                                 <  2:3162835/139‑1 (MQ=255)
                                                                                                   aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTg                                 <  1:2588516/140‑1 (MQ=255)
                                                                                                   aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTg                                 >  2:84524/1‑140 (MQ=255)
                                                                                                       gcgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                  <  1:2795090/87‑1 (MQ=255)
                                                                                                       gcgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                  >  2:2794625/1‑87 (MQ=255)
                                                                                                       gcgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                  >  2:2795090/1‑87 (MQ=255)
                                                                                                       gcgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                  <  1:2794625/87‑1 (MQ=255)
                                                                                                             gCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGGAACACTTTGCGGTTTAATa                       >  1:297364/1‑140 (MQ=255)
                                                                                                             gCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCCAGCTGCTGACCGCTGATAGCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTGACACTTCGCGTTTTAATa                       >  1:4284916/1‑140 (MQ=255)
                                                                                                                ttGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAACCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg                    <  2:2040913/140‑1 (MQ=255)
                                                                                                                        cccAGTGCTGGCCCCACATCCCCCCCCCTTACGGTAATCGGGCTGCTGACCGCTGAGATGCGCCAGCAGCTGACGGCAGTAGGCACTACGCGCCCCGGCCCCTATCAACGTAACACTTTTCAGTGTAATACCGCATGCaa             >  1:3908099/1‑139 (MQ=255)
                                                                                                                        cccAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAAACGAGCGGCTGACCGCTGATATCCGCCAACATCTGACGCCAGTACTCACTTCGCGCCCCTCCCCCCATCAACGTAACACTTTTGCGTTTAATACCGCAGGCAt             >  1:3908159/1‑140 (MQ=255)
                                                                                                                        cccAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGCTATCCGCCAGCATCTGACGCCAGTACTCCCTCCGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATg            >  1:464572/1‑140 (MQ=255)
                                                                                                                        cccAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCCCAGGCATg            >  1:464046/1‑140 (MQ=255)
                                                                                                                        cccAGTGCTGGCCCCACAT‑‑CCCCCCCAGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCacg            >  1:1826447/1‑138 (MQ=255)
                                                                                                                           aGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc         <  1:3147191/140‑1 (MQ=255)
                                                                                                                           aGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc         <  1:1475808/140‑1 (MQ=255)
                                                                                                                            gTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACg        <  1:2950362/140‑1 (MQ=255)
                                                                                                                            gTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACg        <  1:1762623/140‑1 (MQ=255)
                                                                                                                            gTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACg        <  2:1266999/140‑1 (MQ=255)
                                                                                                                            gTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACg        <  2:1207151/140‑1 (MQ=255)
                                                                                                                            gTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACg        <  1:2612958/140‑1 (MQ=255)
                                                                                                                             tGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCCGTACTCACTACGCGCCCCCCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGa       >  1:4450670/1‑140 (MQ=255)
                                                                                                                             tGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAACCGCGCTACTGACCGCTGATATCCGCCCGCATCTGACGCCAGGACTCACTAAGCCCCCCGCCTCACAACAACGTAAAACTTTGCGGTTTAATACCGCAGGCATACACGa       >  1:5039749/1‑140 (MQ=255)
                                                                                                                                  gCCCCACAT‑‑CCCCCCCCGTACGGTAATCGCGCTGCTGACGGCTGATATCCGCGAAGAGCTGACGCCAGGACTCACTGCGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGATTGATACCGCATGCGTGGACCACAAcc  >  2:2332560/1‑140 (MQ=255)
                                                                                                                                  gCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCCGCATCTGACGGCAGTACTCACTACGCGCCCCGCCCCCAATCAACCTAACACTTTGCGGGTTAATACCGCAGGCATGCACGGCATcc  >  2:2026891/1‑140 (MQ=255)
                                                                                                                                           ||                                                                                                                                   
GGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCC  >  NZ_CP025268/3728661‑3728932

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAC  >  NZ_CP025268/3728664‑3728958
                                                                                                                                                |                                                                                                                                                      
TGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGT                                                                                                                                                <  SRR14920169.2419527/151‑1 (MQ=60)
         GCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCT                                                                                                                                       >  SRR14920169.1750394/1‑151 (MQ=60)
         GCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCT                                                                                                                                       >  SRR14920169.3303734/1‑151 (MQ=60)
              ACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGA                                                                                                                                  >  SRR14920169.4476352/1‑151 (MQ=60)
                 CTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTGCGGTTATCGAGCCGCTGACCG                                                                                                                               >  SRR14920169.4105334/1‑151 (MQ=60)
                    CGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCATATCCCCCCCCCttcggtaatcgatcttttctccgttt                                                                                                                              >  SRR14920169.1609304/1‑125 (MQ=60)
                      CATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCC‑‑GACCGTAATCGAGCTGCTGAgcgcctat                                                                                                                          >  SRR14920169.1949600/1‑143 (MQ=60)
                         CTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATC                                                                                                                       <  SRR14920169.3457504/151‑1 (MQ=60)
                          TGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCC                                                                                                                      >  SRR14920169.1034652/1‑151 (MQ=60)
                           GGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCTCCCATTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCctctgaccgcggataaacg                                                                                                                     >  SRR14920169.2534654/1‑132 (MQ=60)
                                   CGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCCTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATC                                                                                                             <  SRR14920169.4062511/151‑1 (MQ=60)
                                        GTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCC‑‑TACGGTAAACCAGCTGCTGGCCGCTGATATCCGCCAGCCGCTGACG                                                                                                        >  SRR14920169.5003893/1‑151 (MQ=60)
                                          TCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ACCCGCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCC                                                                                                      >  SRR14920169.3653698/1‑151 (MQ=60)
                                               aggctccggtctcgtgggctcggagatgtgtataagagacagCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTAC                                                                                                 <  SRR14920169.2800857/109‑1 (MQ=60)
                                               aggctccggtctcgtgggctcggagatgtgtataagagacagCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTAC                                                                                                 <  SRR14920169.2801324/109‑1 (MQ=60)
                                               CGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTAC                                                                                                 <  SRR14920169.1479459/151‑1 (MQ=60)
                                               CGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTAC                                                                                                 <  SRR14920169.3785378/151‑1 (MQ=60)
                                               CGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTAC                                                                                                 <  SRR14920169.3750929/151‑1 (MQ=60)
                                                     TTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTA                                                                                           <  SRR14920169.3179478/151‑1 (MQ=60)
                                                         GGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCG                                                                                       <  SRR14920169.1111522/151‑1 (MQ=60)
                                                          GCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGC                                                                                      >  SRR14920169.3434383/1‑151 (MQ=60)
                                                            agaagacggcatacgagataggctccggtctcgtgggctcggagatgtgtataagagacagGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCC                                                                                    <  SRR14920169.1027722/90‑1 (MQ=60)
                                                               AATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGC                                                                                 >  SRR14920169.1896381/1‑151 (MQ=60)
                                                               AATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCCAGCTGTTGGCCGCTTAAATCCGCCAGGCTCTAACGCGAGTACCCACTACGCGCCCCGC                                                                                 >  SRR14920169.1896391/1‑151 (MQ=60)
                                                                   CAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCC                                                                             <  SRR14920169.1347542/151‑1 (MQ=60)
                                                                   CAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCC                                                                             <  SRR14920169.1347334/151‑1 (MQ=60)
                                                                       GAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATC                                                                         <  SRR14920169.2611980/151‑1 (MQ=60)
                                                                                 TCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACAC                                                               <  SRR14920169.4297385/151‑1 (MQ=60)
                                                                                               GATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGGAACACTTTGCGGTTTAATA                                                 >  SRR14920169.298080/1‑151 (MQ=60)
                                                                                               GATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCCAGCTGCTGACCGCTGATAGCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTGACACTTCGCGTTTTAATA                                                 >  SRR14920169.4295335/1‑151 (MQ=60)
                                                                                                ATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAC                                                <  SRR14920169.2594256/151‑1 (MQ=60)
                                                                                                          GCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCtacggtaaatgaccgtccgaccgccgatttccgcccgcatcttacgccaggcctcaccacccgccccgcccccaaacaccgtaaccccttccggtgtaatacccgcggcctg                                        >  SRR14920169.1191256/1‑39 (MQ=60)
                                                                                                          GCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCTTACGGTAATCGGGCTGCTGACCGCTGAGATGCGCCAGCAGCTGACGGCAGTAGGCACTACGCGCCCCGGCCCCTATCAACGTAACACTTTTCAGTGTAATACCGCAtgcaa                                       >  SRR14920169.3917452/1‑146 (MQ=60)
                                                                                                          GCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAAACGAGCGGCTGACCGCTGATATCCGCCAACATCTGACGCCAGTACTCACTTCGCGCCCCTCCCCCCATCAACGTAACACTTTTGCGTTTAATACCGCAGGCAT                                       >  SRR14920169.3917512/1‑151 (MQ=60)
                                                                                                          GCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGCTATCCGCCAGCATCTGACGCCAGTACTCCCTCCGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATG                                      >  SRR14920169.465669/1‑151 (MQ=60)
                                                                                                          GCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCCCAGGCATG                                      >  SRR14920169.465141/1‑151 (MQ=60)
                                                                                                          GCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCAGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCACG                                      >  SRR14920169.1830538/1‑151 (MQ=60)
                                                                                                               GCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCCGTACTCACTACGCGCCCCCCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGA                                 >  SRR14920169.4461484/1‑151 (MQ=60)
                                                                                                               GCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAACCGCGCTACTGACCGCTGATATCCGCCCGCATCTGACGCCAGGACTCACTAAGCCCCCCGCCTCACAACAACGTAAAACTTTGCGGTTTAATACCGCAGGCATACACGA                                 >  SRR14920169.5051701/1‑151 (MQ=60)
                                                                                                                        AGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGC                        <  SRR14920169.3154664/151‑1 (MQ=60)
                                                                                                                        AGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGC                        <  SRR14920169.1479114/151‑1 (MQ=60)
                                                                                                                         GTGCTGGCCCCACAT‑‑CCCCCCCCGTTCGGTAATCGAGCTGCTGACCGCTGGTATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAAAACTTTGCGGGTTAATACCGCAGGCATGCACGACATCCATGCC                       >  SRR14920169.1725718/1‑151 (MQ=60)
                                                                                                                         GTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCcccctgtctcttatacacatctccgagcccacgagaccggagcctatcgcgtatgccgtcttctgcttgtaaaaataaatctttgtgtgtgatat                       >  SRR14920169.289471/1‑56 (MQ=60)
                                                                                                                         GTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCctgtctcttatacacatctccgagcccacgagaccggagcctatctcgtatgccgtcttcttcttgcaaaaaatgataggtgtgggggttgg                       >  SRR14920169.289630/1‑59 (MQ=60)
                                                                                                                         GTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCC                       <  SRR14920169.2957261/151‑1 (MQ=60)
                                                                                                                         GTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCC                       <  SRR14920169.2618736/151‑1 (MQ=60)
                                                                                                                         GTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCC                       <  SRR14920169.1766568/151‑1 (MQ=60)
                                                                                                                         GTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCGGAAATCCGCCCGCATCTGACGCCAGTACTCACCACCCGCCCCGCCCCCAATCAACGTAACACTTTGCTGTTTAATACCGCCGGCCTGCACGACATCCATGCC                       >  SRR14920169.1010790/1‑151 (MQ=60)
                                                                                                                         GTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCCGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACctgtctcttatacacatctccgagcccacgagaccggagc                       >  SRR14920169.1069671/1‑111 (MQ=60)
                                                                                                                         GTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCGGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAACCAACATAACACTTTTCGGCTTAATACCGCAGGCCTGCACGACATCCATGCC                       >  SRR14920169.3460935/1‑151 (MQ=60)
                                                                                                                                CCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCC                >  SRR14920169.3860960/1‑151 (MQ=60)
                                                                                                                                    ACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACGCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCG            >  SRR14920169.1630385/1‑151 (MQ=60)
                                                                                                                                       catCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGCCGCCAGTACTCACTACGCGCCCCCCCCCCAATCAACGTCACACCTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGC           >  SRR14920169.3643367/4‑151 (MQ=60)
                                                                                                                                       catCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCCACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGC           >  SRR14920169.863290/4‑151 (MQ=60)
                                                                                                                                       catCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCCCGCCATCCATGCCATCTGCCAGCGC           >  SRR14920169.4538144/4‑151 (MQ=60)
                                                                                                                                            CCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAC      >  SRR14920169.4202662/1‑151 (MQ=60)
                                                                                                                                            CCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAC      <  SRR14920169.26760/151‑1 (MQ=60)
                                                                                                                                            CCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAC      <  SRR14920169.98085/151‑1 (MQ=60)
                                                                                                                                             CCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCtgtctcttatacacatctc     >  SRR14920169.4148102/1‑132 (MQ=60)
                                                                                                                                             CCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAGCAACGTAACACCTTGCGGTGTAATACCGCAGTCATGCACGGCATCCATGCtggctcttatacacatctc     >  SRR14920169.4147488/1‑132 (MQ=60)
                                                                                                                                               CCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCA   >  SRR14920169.378297/1‑151 (MQ=60)
                                                                                                                                                CCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAC  >  SRR14920169.1161737/1‑151 (MQ=60)
                                                                                                                                                |                                                                                                                                                      
TGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAC  >  NZ_CP025268/3728664‑3728958

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