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 521 238.2 12299412 96.0% 11807435 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,554,113 Δ1 bp coding (372/396 nt) yhgKJ ← RNA phosphage cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,554,1130C.99.1% 989.7 / ‑0.6 217coding (372/396 nt)yhgKJRNA phosphage cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (104/111);  minor base G (1/1);  total (105/112)
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.95e-01

Only 100 of 237 total aligned reads displayed.

TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGC  >  NC_000913/3554026‑3554199
                                                                                         |                                                                                       
tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgc                                                                                   <  2:5106125/93‑1 (MQ=255)
 ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAACTCCGCAGGCGGCGCCGACGGGTTATCGGtgccgcc                                                                                  <  1:6116065/93‑2 (MQ=255)
 ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgcc                                                                                  >  1:5745745/1‑92 (MQ=255)
 ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgcc                                                                                  <  1:3271546/93‑2 (MQ=255)
  gTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccg                                                                                 >  2:5960364/1‑91 (MQ=255)
      gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCa                                                                              <  2:368950/93‑1 (MQ=255)
      gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCa                                                                              >  1:3595083/1‑93 (MQ=255)
       cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCAc                                                                             <  2:552618/93‑1 (MQ=255)
       cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCAc                                                                             <  2:140838/93‑1 (MQ=255)
        tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACt                                                                            <  1:3614906/93‑1 (MQ=255)
        tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACt                                                                            <  2:2803082/93‑1 (MQ=255)
        tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCAGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACt                                                                            >  1:4649155/1‑93 (MQ=255)
         gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGGGCCGACGGGTTATCGGTGCCGCCGCTCACtt                                                                           >  2:2456030/1‑93 (MQ=255)
         gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCTGCCGCTCACtt                                                                           >  1:2262459/1‑93 (MQ=255)
         gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACtt                                                                           <  2:1704437/93‑1 (MQ=255)
         gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGACGACGGGTTGTCGGTGCCGCCGCTCACtt                                                                           >  1:2146772/1‑93 (MQ=255)
           tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGGCGCCGCTCACTTCa                                                                         >  2:72051/1‑93 (MQ=255)
            gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCaa                                                                        >  1:4330893/1‑93 (MQ=255)
             aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAac                                                                       <  2:2518328/93‑1 (MQ=255)
              aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAaca                                                                      >  1:4044213/1‑93 (MQ=255)
              aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAaca                                                                      >  1:5152306/1‑93 (MQ=255)
              aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAaca                                                                      <  1:5792925/93‑1 (MQ=255)
              aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAaca                                                                      >  1:1093900/1‑93 (MQ=255)
                tCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcg                                                                    >  1:2621689/1‑93 (MQ=255)
                 ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcgc                                                                   >  1:1369888/1‑93 (MQ=255)
                   tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGa                                                                 >  2:2341466/1‑93 (MQ=255)
                    aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGaa                                                                <  2:2846502/93‑1 (MQ=255)
                    aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGaa                                                                <  2:908857/93‑1 (MQ=255)
                    aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGaa                                                                >  1:3229691/1‑93 (MQ=255)
                     ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAg                                                               <  1:2528060/93‑1 (MQ=255)
                      tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAgg                                                              <  1:3819659/93‑1 (MQ=255)
                      tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCCAAgg                                                              >  2:4206211/1‑93 (MQ=255)
                       ttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGt                                                             <  1:203210/93‑1 (MQ=255)
                       ttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGt                                                             <  2:3843900/93‑1 (MQ=255)
                        tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTc                                                            >  1:2823921/1‑93 (MQ=255)
                        tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGAGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAAAACGCGAAGGTc                                                            >  2:1128360/1‑93 (MQ=255)
                          gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCa                                                          >  1:4117577/1‑93 (MQ=255)
                          gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCa                                                          >  1:5446225/1‑93 (MQ=255)
                          gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCa                                                          <  1:4616587/93‑1 (MQ=255)
                           ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                         <  1:5036629/93‑1 (MQ=255)
                             agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                                                       >  2:1935657/1‑93 (MQ=255)
                             agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                                                       <  2:398096/93‑1 (MQ=255)
                             agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                                                       >  1:975540/1‑93 (MQ=255)
                              gcagcGGCTCCAGCACCCGGCGGATAAAACCCGCAGGCGGCGGCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                      <  2:4904803/93‑1 (MQ=255)
                               cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                                                     >  2:2525851/1‑93 (MQ=255)
                               cagcGGCCCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                                                     <  1:1214809/93‑1 (MQ=255)
                                 ggcgCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACGTCAACACGCGAAGGTCCATCGGCaa                                                   <  2:4243033/90‑1 (MQ=255)
                                    gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGCAGGTCCATCGGCAAAcc                                                >  1:3325581/1‑93 (MQ=255)
                                    gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCCAAcc                                                >  2:2364643/1‑93 (MQ=255)
                                     cTCCAGCACCCGGCGGATAAAATCTGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                                               >  1:3767623/1‑93 (MQ=255)
                                     cTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                                               <  1:4729534/93‑1 (MQ=255)
                                      tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                                              <  1:5925167/93‑1 (MQ=255)
                                       ccAGCCCCCGGCGGAGAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                                             <  2:619616/93‑1 (MQ=255)
                                       ccAGCACCCGGCGGATAAAATCCTCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                                             <  1:1774345/93‑1 (MQ=255)
                                       ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                                             <  2:3942790/93‑1 (MQ=255)
                                        cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGCCGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                                            <  1:5668838/93‑1 (MQ=255)
                                        cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                                            <  2:5788909/93‑1 (MQ=255)
                                         agaccccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                           <  1:367236/89‑1 (MQ=255)
                                         aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                           >  2:5374872/1‑93 (MQ=255)
                                           cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                                         <  2:1081765/93‑1 (MQ=255)
                                            ccccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                        <  2:3572140/92‑1 (MQ=255)
                                            aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                        <  1:6129880/93‑1 (MQ=255)
                                              ccGGCGGATAAAATCCGCAGTCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCggg                                      <  1:2060123/93‑1 (MQ=255)
                                              ccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCggg                                      >  2:1023053/1‑93 (MQ=255)
                                               cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                                     <  2:5291463/93‑1 (MQ=255)
                                               cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                                     <  2:4634965/93‑1 (MQ=255)
                                                ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca                                    <  1:5822580/93‑1 (MQ=255)
                                                 gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag                                   >  1:5541446/1‑93 (MQ=255)
                                                  cggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagc                                  >  1:5646944/1‑93 (MQ=255)
                                                   caaTAAAATCCGCAGGCGGCGCCGAAGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                 <  1:5201813/91‑1 (MQ=255)
                                                    gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                >  2:2149066/1‑93 (MQ=255)
                                                    gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                <  2:839045/93‑1 (MQ=255)
                                                    gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                >  1:180724/1‑93 (MQ=255)
                                                    gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                >  1:5382920/1‑93 (MQ=255)
                                                        aaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt                            >  1:2376141/1‑93 (MQ=255)
                                                        aaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt                            <  2:3132004/93‑1 (MQ=255)
                                                             ccgcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca                        <  2:5683959/89‑1 (MQ=255)
                                                               gcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagc                      >  1:4540780/3‑93 (MQ=255)
                                                               gcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagc                      <  1:3128080/91‑1 (MQ=255)
                                                                 aGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAc                    <  2:422001/93‑1 (MQ=255)
                                                                 aGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAc                    >  1:4209814/1‑93 (MQ=255)
                                                                 aGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAc                    >  1:2737028/1‑93 (MQ=255)
                                                                  ggcggcGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAcc                   <  1:289894/93‑1 (MQ=255)
                                                                  ggcggcGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAcc                   <  1:796386/93‑1 (MQ=255)
                                                                    cggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAg                 <  1:2462803/93‑1 (MQ=255)
                                                                    cggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAg                 <  2:461439/93‑1 (MQ=255)
                                                                       cgcCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGt              <  1:1934202/93‑1 (MQ=255)
                                                                         ccGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGGCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAc            >  2:3672220/1‑93 (MQ=255)
                                                                          cGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACa           <  2:4322286/93‑1 (MQ=255)
                                                                           gACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACaa          >  2:4789332/1‑93 (MQ=255)
                                                                           gACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACaa          <  1:2694482/93‑1 (MQ=255)
                                                                           gACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACaa          >  2:1310621/1‑93 (MQ=255)
                                                                            aCGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc         <  2:3912960/93‑1 (MQ=255)
                                                                             cGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACt        >  2:4467410/1‑93 (MQ=255)
                                                                              gggTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACtt       >  1:4944580/1‑93 (MQ=255)
                                                                              gggTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACtt       >  1:2721045/1‑93 (MQ=255)
                                                                                gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  1:722257/93‑1 (MQ=255)
                                                                                  tATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCgg   <  2:3010776/93‑1 (MQ=255)
                                                                                  tATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATAGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCgg   <  1:519353/93‑1 (MQ=255)
                                                                                   aTCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGgc  <  2:4191921/93‑1 (MQ=255)
                                                                                         |                                                                                       
TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGC  >  NC_000913/3554026‑3554199

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGT‑‑CGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGCG  >  NC_000913/3554025‑3554200
                                                                                            |                                                                                        
GTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGC                                                                                    >  SRR4302532.1167454/1‑94 (MQ=60)
 TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCC                                                                                       >  SRR4302532.5945758/1‑94 (MQ=60)
 TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCC                                                                                       <  SRR4302532.4341421/94‑1 (MQ=60)
  GGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAACTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCG                                                                                      <  SRR4302532.6327140/94‑1 (MQ=60)
  GGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCG                                                                                      <  SRR4302532.3380335/94‑1 (MQ=60)
      TGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCA                                                                               >  SRR4302532.3715307/1‑94 (MQ=60)
       GCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCAC                                                                              <  SRR4302532.4186234/94‑1 (MQ=60)
       GCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCAC                                                                              >  SRR4302532.283092/1‑94 (MQ=60)
        CTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACT                                                                             <  SRR4302532.3724081/94‑1 (MQ=60)
        CTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCAGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACT                                                                             >  SRR4302532.4807665/1‑94 (MQ=60)
         TGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTT                                                                            <  SRR4302532.3735806/94‑1 (MQ=60)
         TGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTT                                                                            >  SRR4302532.6000525/1‑94 (MQ=60)
         TGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGACGACGGGTTGTCGGT‑‑GCCGCCGCTCACTT                                                                            >  SRR4302532.2218508/1‑94 (MQ=60)
          GATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTC                                                                           <  SRR4302532.286394/94‑1 (MQ=60)
           ATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCA                                                                          >  SRR4302532.2536094/1‑94 (MQ=60)
           ATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATACGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCA                                                                          >  SRR4302532.6348394/1‑94 (MQ=60)
            TGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAA                                                                         >  SRR4302532.4477207/1‑94 (MQ=60)
            TGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAA                                                                         >  SRR4302532.5491375/1‑94 (MQ=60)
             GAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAAC                                                                        >  SRR4302532.3786782/1‑94 (MQ=60)
              AATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACA                                                                       >  SRR4302532.4180387/1‑94 (MQ=60)
              AATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACA                                                                       >  SRR4302532.5329556/1‑94 (MQ=60)
              AATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACA                                                                       >  SRR4302532.1130481/1‑94 (MQ=60)
               ATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACAC                                                                      >  SRR4302532.380866/1‑94 (MQ=60)
               ATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACAC                                                                      <  SRR4302532.5994311/94‑1 (MQ=60)
                TTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACG                                                                     >  SRR4302532.4081427/1‑94 (MQ=60)
                TTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACG                                                                     >  SRR4302532.2709467/1‑94 (MQ=60)
                 TCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGC                                                                    >  SRR4302532.1415778/1‑94 (MQ=60)
                   CTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGA                                                                  >  SRR4302532.1747960/1‑94 (MQ=60)
                   CTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGA                                                                  >  SRR4302532.5236094/1‑94 (MQ=60)
                    TATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAA                                                                 >  SRR4302532.3381084/1‑94 (MQ=60)
                    TATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAA                                                                 >  SRR4302532.3336897/1‑94 (MQ=60)
                    TATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAA                                                                 >  SRR4302532.4764109/1‑94 (MQ=60)
                     ATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAG                                                                >  SRR4302532.5601616/1‑94 (MQ=60)
                     ATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAG                                                                <  SRR4302532.1306136/94‑1 (MQ=60)
                     ATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAG                                                                <  SRR4302532.1612057/94‑1 (MQ=60)
                      TTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGG                                                               <  SRR4302532.3240674/94‑1 (MQ=60)
                      TTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGG                                                               <  SRR4302532.2612416/94‑1 (MQ=60)
                       TTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGT                                                              <  SRR4302532.3948611/94‑1 (MQ=60)
                       TTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGT                                                              <  SRR4302532.2483741/94‑1 (MQ=60)
                        TTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTC                                                             <  SRR4302532.210134/94‑1 (MQ=60)
                        TTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTC                                                             >  SRR4302532.2918402/1‑94 (MQ=60)
                        TTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTC                                                             <  SRR4302532.577012/94‑1 (MQ=60)
                          CGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCA                                                           >  SRR4302532.5634994/1‑94 (MQ=60)
                          CGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCA                                                           >  SRR4302532.4255777/1‑94 (MQ=60)
                          CGCCACCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCA                                                           <  SRR4302532.2597435/94‑1 (MQ=60)
                           GCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAT                                                          <  SRR4302532.4774000/94‑1 (MQ=60)
                           GCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAT                                                          <  SRR4302532.225898/94‑1 (MQ=60)
                            CCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATC                                                         <  SRR4302532.5209373/94‑1 (MQ=60)
                            CCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATC                                                         <  SRR4302532.5808196/94‑1 (MQ=60)
                            CCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACCCGAAGGTCCATC                                                         >  SRR4302532.3919333/1‑94 (MQ=60)
                             CAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCG                                                        >  SRR4302532.1008842/1‑94 (MQ=60)
                             CAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCG                                                        <  SRR4302532.4909568/94‑1 (MQ=60)
                             CAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCG                                                        <  SRR4302532.6303167/94‑1 (MQ=60)
                              AGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGG                                                       >  SRR4302532.3212424/1‑94 (MQ=60)
                                CAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCA                                                     <  SRR4302532.2126483/94‑1 (MQ=60)
                                CAGCGGCCCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCA                                                     <  SRR4302532.1255469/94‑1 (MQ=60)
                                 AGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAA                                                    >  SRR4302532.5310616/1‑94 (MQ=60)
                                 AGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGGCCATCGGCAA                                                    >  SRR4302532.5981412/1‑94 (MQ=60)
                                  GCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAA                                                   >  SRR4302532.6327782/1‑94 (MQ=60)
                                    GGCTCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACC                                                 <  SRR4302532.2718328/94‑1 (MQ=60)
                                    GGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGCAGGTCCATCGGCAAACC                                                 >  SRR4302532.3436480/1‑94 (MQ=60)
                                     GCTCCAGCACCCGGCGGATAAAATCTGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCA                                                >  SRR4302532.3894418/1‑94 (MQ=60)
                                      CTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCAC                                               <  SRR4302532.4890966/94‑1 (MQ=60)
                                      CTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCAC                                               <  SRR4302532.5894063/94‑1 (MQ=60)
                                       cccacCACCCGGCGGATAAAATCCGCAGGCGGGGCCGACGGGTTAGCGGT‑‑GCCGCCGCTAACTTCAACACGCGAAGGTCCATCGGCAAACCACA                                              <  SRR4302532.329794/89‑1 (MQ=60)
                                       TCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACA                                              <  SRR4302532.3664470/94‑1 (MQ=60)
                                       TCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACA                                              >  SRR4302532.2770200/1‑94 (MQ=60)
                                       TCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACA                                              >  SRR4302532.2667519/1‑94 (MQ=60)
                                       TCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACA                                              <  SRR4302532.6130334/94‑1 (MQ=60)
                                       TCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACA                                              >  SRR4302532.4761851/1‑94 (MQ=60)
                                       TCCAGCACCCGGCGGATAAAAGCCGCAGGCGGCGCCGACGGATTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACA                                              <  SRR4302532.679527/94‑1 (MQ=60)
                                        CCAGCACCCGGCGGATAAAATCCTCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAG                                             <  SRR4302532.1833771/94‑1 (MQ=60)
                                        CCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAG                                             >  SRR4302532.1866899/1‑94 (MQ=60)
                                        CCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAG                                             >  SRR4302532.5919235/1‑94 (MQ=60)
                                         CAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGCCGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGC                                            <  SRR4302532.5866500/94‑1 (MQ=60)
                                         CAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGC                                            <  SRR4302532.784513/94‑1 (MQ=60)
                                         CAGCACCCGGCGAATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGC                                            <  SRR4302532.3753405/94‑1 (MQ=60)
                                          aggcCCCGGCGGATAAAATCGGCAGCCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCG                                           <  SRR4302532.2059126/90‑1 (MQ=60)
                                          agacCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCG                                           <  SRR4302532.379997/90‑1 (MQ=60)
                                          AGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCG                                           <  SRR4302532.5350210/94‑1 (MQ=60)
                                            CACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCG                                         >  SRR4302532.6063291/1‑94 (MQ=60)
                                            CACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCG                                         >  SRR4302532.5491309/1‑94 (MQ=60)
                                             CCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGG                                        <  SRR4302532.5853777/94‑1 (MQ=60)
                                             ACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGG                                        <  SRR4302532.6341392/94‑1 (MQ=60)
                                              CCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGG                                       >  SRR4302532.909408/1‑94 (MQ=60)
                                               CCGGCGGATAAAATCCGCAGTCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGC                                      <  SRR4302532.2128745/94‑1 (MQ=60)
                                                 GGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAG                                    <  SRR4302532.6024780/94‑1 (MQ=60)
                                                 GGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAG                                    <  SRR4302532.4598361/94‑1 (MQ=60)
                                                 GGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAG                                    >  SRR4302532.5733771/1‑94 (MQ=60)
                                                 GGAGGAGAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑CCCGCCGCTCACGGCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAG                                    <  SRR4302532.3966696/94‑1 (MQ=60)
                                                  GCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGC                                   >  SRR4302532.5843361/1‑94 (MQ=60)
                                                   CGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCC                                  >  SRR4302532.4825452/1‑94 (MQ=60)
                                                    caATAAAATCCGCAGGCGGCGCCGAAGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAC                                 <  SRR4302532.5381178/92‑1 (MQ=60)
                                                    GGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAC                                 <  SRR4302532.5381177/94‑1 (MQ=60)
                                                    GGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAC                                 >  SRR4302532.186805/1‑94 (MQ=60)
                                                    GGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAC                                 >  SRR4302532.5569367/1‑94 (MQ=60)
                                                     GATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACC                                <  SRR4302532.520318/94‑1 (MQ=60)
                                                     GATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACC                                <  SRR4302532.6088966/94‑1 (MQ=60)
                                                      ATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCG                               <  SRR4302532.3390266/94‑1 (MQ=60)
                                                        AAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTT                             >  SRR4302532.2361573/1‑94 (MQ=60)
                                                        AAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTT                             >  SRR4302532.2455340/1‑94 (MQ=60)
                                                         AAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT                            <  SRR4302532.1980937/94‑1 (MQ=60)
                                                           agtccgcAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTG                           <  SRR4302532.2712122/87‑1 (MQ=60)
                                                               cgcAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGC                       >  SRR4302532.2046186/4‑94 (MQ=60)
                                                               cgcAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGC                       >  SRR4302532.4695782/4‑94 (MQ=60)
                                                                gcAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCA                      <  SRR4302532.3231920/92‑1 (MQ=60)
                                                                CAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAC                     >  SRR4302532.4351133/1‑94 (MQ=60)
                                                                CAGGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAC                     >  SRR4302532.2829016/1‑94 (MQ=60)
                                                                   GGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCA                   <  SRR4302532.299884/94‑1 (MQ=60)
                                                                   GGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCA                   <  SRR4302532.823620/94‑1 (MQ=60)
                                                                   GGCGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCA                   <  SRR4302532.3367374/94‑1 (MQ=60)
                                                                     CGGCGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGC                 <  SRR4302532.2544741/94‑1 (MQ=60)
                                                                        CGCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTA              <  SRR4302532.1998561/94‑1 (MQ=60)
                                                                         GCCGACGGGTTATCGGT‑‑GCCTCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAC             <  SRR4302532.2256815/94‑1 (MQ=60)
                                                                         GCCGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCCGGCAGCACCGTTTGCAGAACCAGCGTCC             >  SRR4302532.2885167/1‑94 (MQ=60)
                                                                           CGACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAA           >  SRR4302532.4785175/1‑94 (MQ=60)
                                                                            GGCGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAC          <  SRR4302532.2784986/94‑1 (MQ=60)
                                                                            GACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAC          <  SRR4302532.2784982/94‑1 (MQ=60)
                                                                             ACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACT         <  SRR4302532.4426581/94‑1 (MQ=60)
                                                                             ACGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACT         >  SRR4302532.4313674/1‑94 (MQ=60)
                                                                              CGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTT        >  SRR4302532.2812558/1‑94 (MQ=60)
                                                                              CGGGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTT        >  SRR4302532.5114291/1‑94 (MQ=60)
                                                                                GGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCC      <  SRR4302532.2133651/94‑1 (MQ=60)
                                                                                GGTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCC      <  SRR4302532.5661784/94‑1 (MQ=60)
                                                                                 GTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCG     <  SRR4302532.2963875/94‑1 (MQ=60)
                                                                                 GTTATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCG     <  SRR4302532.746897/94‑1 (MQ=60)
                                                                                   TATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGC   >  SRR4302532.3036012/1‑94 (MQ=60)
                                                                                   TATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATAGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGC   <  SRR4302532.537487/94‑1 (MQ=60)
                                                                                    ATCGGT‑‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGCG  >  SRR4302532.3646904/1‑94 (MQ=60)
                                                                                            |                                                                                        
GTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGT‑‑CGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGCG  >  NC_000913/3554025‑3554200

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