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,957,551 +T intergenic (‑120/‑741) ppiC ← / → rep peptidyl‑prolyl cis‑trans isomerase C/rep helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,957,5501.T98.5% 498.5 / ‑0.3 138intergenic (‑119/‑742)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (60/75);  total (63/75)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 139 total aligned reads displayed.

TAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTG‑‑CTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGT  >  NC_000913/3957470‑3957627
                                                                                     |                                                                               
tAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTcc                                                                     <  1:4646966/93‑1 (MQ=255)
   cATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTccctc                                                                  >  2:4558835/1‑90 (MQ=255)
    aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTccctca                                                                 <  1:3594481/93‑5 (MQ=255)
       aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAgg                                                               <  1:826261/93‑1 (MQ=255)
       aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAgg                                                               <  1:151896/93‑1 (MQ=255)
        tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAggg                                                              >  2:4325677/1‑93 (MQ=255)
        tCCGCAGATGCTAACATAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAggg                                                              >  2:2895915/1‑93 (MQ=255)
         ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAgggg                                                             <  2:5923120/93‑1 (MQ=255)
          cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGga                                                            <  2:3388717/93‑1 (MQ=255)
               aTGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGaaaa                                                       >  1:308986/1‑93 (MQ=255)
                tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                      >  2:4544786/1‑93 (MQ=255)
                tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                      >  1:2605310/1‑93 (MQ=255)
                 gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACt                                                     >  1:2938992/1‑93 (MQ=255)
                 gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CGTCGCAGGTCCCTCAAGGGGAGAAAACt                                                     >  2:1293094/1‑93 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                  <  1:1018440/93‑1 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                  >  1:1900277/1‑93 (MQ=255)
                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                  <  2:1925912/93‑1 (MQ=255)
                     aCGTAATTGTCAGTCAGGGAGTCGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                 <  1:3201179/93‑1 (MQ=255)
                       gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                               >  1:3249213/1‑93 (MQ=255)
                       gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                               <  2:1798259/93‑1 (MQ=255)
                       gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCTCTCAAGGGGAGAAAACTTAGGGc                                               >  1:3443540/1‑93 (MQ=255)
                        tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                              <  2:5958304/93‑1 (MQ=255)
                        tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                              >  2:1670348/1‑93 (MQ=255)
                        tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                              <  1:1056185/93‑1 (MQ=255)
                         aaTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                                             <  2:6071234/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                            <  2:4959993/93‑1 (MQ=255)
                          aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                            <  2:1815027/93‑1 (MQ=255)
                            tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                          <  2:3506523/93‑1 (MQ=255)
                              tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                        <  1:1129571/93‑1 (MQ=255)
                                agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                                      >  2:828683/1‑93 (MQ=255)
                                agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                                      >  2:1298989/1‑93 (MQ=255)
                                agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                                      >  2:1197914/1‑93 (MQ=255)
                                agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                                      <  1:521336/93‑1 (MQ=255)
                                agtcagGGAGTGGCTGACACATATATGGCAAAACCCCTGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                                      <  2:2739448/93‑1 (MQ=255)
                                  tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                    >  2:1566262/1‑93 (MQ=255)
                                  tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                    >  1:4757245/1‑93 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGGGAAAACTTAGGGCCCCCATggggtcc                                  >  2:1508563/1‑89 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                  >  1:4052741/1‑93 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                  <  2:2358848/93‑1 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                  <  1:619340/93‑1 (MQ=255)
                                    agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                  >  2:1508564/1‑93 (MQ=255)
                                     gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                 >  2:1727812/1‑93 (MQ=255)
                                     gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                 >  2:4555456/1‑93 (MQ=255)
                                      ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                                <  2:3183266/93‑1 (MQ=255)
                                      ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCCTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                                >  2:5165397/1‑93 (MQ=255)
                                       gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                               <  2:585246/93‑1 (MQ=255)
                                       gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                               <  1:5480111/93‑1 (MQ=255)
                                        aGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAg                              <  1:4886800/93‑1 (MQ=255)
                                        aGTGGCTGGTACATATAAGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAg                              <  1:1996586/93‑1 (MQ=255)
                                         gTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                             <  1:3299215/93‑1 (MQ=255)
                                          tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                            >  1:4814745/1‑93 (MQ=255)
                                            gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAt                          <  1:2671362/93‑1 (MQ=255)
                                             cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt                         <  1:907817/93‑1 (MQ=255)
                                              tGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                        >  1:1533937/1‑93 (MQ=255)
                                              tGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                        >  1:4211893/1‑93 (MQ=255)
                                              tGGTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                        <  1:5584220/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTTTATGGGCTACAAAGTCATTAt                       <  1:4971884/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                       <  1:5958066/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                       <  2:2162068/93‑1 (MQ=255)
                                               ggTACATATATGGCAAAAGCCCGAAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                       <  1:2561476/93‑1 (MQ=255)
                                                gTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                      >  1:255486/1‑93 (MQ=255)
                                                gTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                      >  2:1131429/1‑93 (MQ=255)
                                                gTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                      <  2:3238432/93‑1 (MQ=255)
                                                gTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                      >  1:3562583/1‑93 (MQ=255)
                                                gTACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGAGCCTCTATGGGCTACAAAGTCATTATc                      >  2:235248/1‑93 (MQ=255)
                                                 tACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                     <  1:4456508/93‑1 (MQ=255)
                                                 tACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                     >  2:4773745/1‑93 (MQ=255)
                                                 tACATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                     <  2:4320497/93‑1 (MQ=255)
                                                   cATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                   >  2:2286091/1‑93 (MQ=255)
                                                   cATATATGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                   >  1:5961297/1‑93 (MQ=255)
                                                    atatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                  >  1:5810927/1‑93 (MQ=255)
                                                    atatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                  >  1:3787/1‑93 (MQ=255)
                                                    atatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                  >  2:802404/1‑93 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                 <  1:686916/93‑1 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                 <  1:5811190/93‑1 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                 <  1:3881312/93‑1 (MQ=255)
                                                     tatatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                 >  2:260567/1‑93 (MQ=255)
                                                      atatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt                <  2:3942384/93‑1 (MQ=255)
                                                      atagGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt                <  2:158199/93‑1 (MQ=255)
                                                       tatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAACGTCATTATCGGGGATg               >  2:307469/1‑93 (MQ=255)
                                                       tatGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATg               <  2:5639335/93‑1 (MQ=255)
                                                         tGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt             <  1:5699089/93‑1 (MQ=255)
                                                         tGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt             <  2:651666/93‑1 (MQ=255)
                                                         tGGCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt             <  1:6090494/93‑1 (MQ=255)
                                                           gCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtta           <  1:870907/93‑1 (MQ=255)
                                                           gCAAAAGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtta           <  1:3366842/93‑1 (MQ=255)
                                                             aaaaGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt         >  1:2762356/1‑93 (MQ=255)
                                                              aaaGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattt        <  2:417131/93‑1 (MQ=255)
                                                               aaGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttt       <  2:3318719/93‑1 (MQ=255)
                                                               aaGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttt       >  1:3605180/1‑93 (MQ=255)
                                                                aGCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttta      >  2:5262722/1‑93 (MQ=255)
                                                                 gCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat     >  2:3086145/1‑93 (MQ=255)
                                                                 gCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat     >  2:2094205/1‑93 (MQ=255)
                                                                 gCCCGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat     <  2:1168370/93‑1 (MQ=255)
                                                                  cccGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt    >  1:469397/1‑93 (MQ=255)
                                                                  cccGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt    <  1:1553052/93‑1 (MQ=255)
                                                                   ccGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg   >  1:560076/1‑93 (MQ=255)
                                                                   ccGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg   <  2:1736995/93‑1 (MQ=255)
                                                                   ccGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg   <  2:3791210/93‑1 (MQ=255)
                                                                    cGGAAAATCCGGG‑‑CTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt  <  1:2959575/93‑1 (MQ=255)
                                                                                     |                                                                               
TAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTG‑‑CTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGT  >  NC_000913/3957470‑3957627

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTG‑‑CTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGAT  >  NC_000913/3957478‑3957642
                                                                           |                                                                                             
TCCGCAAATGCAAACGTAATGGCCAGCCAGGGAGTGGCTGGTAAATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGG                                                                            <  SRR4302532.1565089/94‑1 (MQ=34)
  CGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAG                                                                          <  SRR4302532.5422281/94‑1 (MQ=60)
     AGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCCAGGGGAGAAA                                                                       >  SRR4302532.175324/1‑94 (MQ=60)
      GATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAA                                                                      >  SRR4302532.319700/1‑94 (MQ=60)
       ATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAC                                                                     >  SRR4302532.2692445/1‑94 (MQ=60)
       ATGCTAACGTAATGGTCAGTCGGGGAGTGGCTGGTACATATAGGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAC                                                                     <  SRR4302532.3712800/94‑1 (MQ=57)
        TGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACT                                                                    >  SRR4302532.3036606/1‑94 (MQ=60)
           TAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAG                                                                 >  SRR4302532.1963294/1‑94 (MQ=60)
            AACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGG                                                                <  SRR4302532.1052998/94‑1 (MQ=60)
             ACGTAATTGTCAGTCAGGGAGTCGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGG                                                               <  SRR4302532.3307454/94‑1 (MQ=60)
              CGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGC                                                              >  SRR4302532.3357094/1‑94 (MQ=60)
              CGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCTCTCAAGGGGAGAAAACTTAGGGC                                                              >  SRR4302532.3558911/1‑94 (MQ=60)
                TAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCT                                                            <  SRR4302532.1091780/94‑1 (MQ=60)
                TAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCT                                                            <  SRR4302532.797056/94‑1 (MQ=60)
                    TGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAT                                                        <  SRR4302532.36448/94‑1 (MQ=60)
                    TGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAT                                                        <  SRR4302532.2961952/94‑1 (MQ=60)
                     GTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATG                                                       <  SRR4302532.3831612/94‑1 (MQ=60)
                      TCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGG                                                      <  SRR4302532.1167266/94‑1 (MQ=60)
                        AGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGC                                                    <  SRR4302532.539546/94‑1 (MQ=60)
                         GTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCT                                                   >  SRR4302532.4920009/1‑94 (MQ=60)
                         GTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCT                                                   <  SRR4302532.4626818/94‑1 (MQ=60)
                         GTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCT                                                   >  SRR4302532.3399002/1‑94 (MQ=60)
                           CAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAC                                                 >  SRR4302532.4189080/1‑94 (MQ=60)
                            AGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACA                                                <  SRR4302532.640927/94‑1 (MQ=60)
                               GAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAG                                             <  SRR4302532.5670115/94‑1 (MQ=60)
                                atgggaTGGTACATATATGGCAAAAGCCGGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGT                                            <  SRR4302532.754152/88‑1 (MQ=60)
                                AGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGT                                            <  SRR4302532.5054761/94‑1 (MQ=60)
                                AGTGGCTGGTACATATAAGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGT                                            <  SRR4302532.2063092/94‑1 (MQ=60)
                                 GTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTC                                           <  SRR4302532.3408990/94‑1 (MQ=60)
                                 GTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTC                                           >  SRR4302532.4980533/1‑94 (MQ=60)
                                    GCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATT                                        <  SRR4302532.2761018/94‑1 (MQ=60)
                                     CTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTA                                       >  SRR4302532.4353268/1‑94 (MQ=60)
                                     CTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTA                                       >  SRR4302532.4907254/1‑94 (MQ=60)
                                     CTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTA                                       <  SRR4302532.938826/94‑1 (MQ=60)
                                     CTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTA                                       >  SRR4302532.1585766/1‑94 (MQ=60)
                                      TGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAT                                      >  SRR4302532.6115746/1‑94 (MQ=60)
                                      TGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAT                                      <  SRR4302532.5778083/94‑1 (MQ=60)
                                       GGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTTTATGGGCTACAAAGTCATTATC                                     <  SRR4302532.5142330/94‑1 (MQ=60)
                                       GGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATC                                     >  SRR4302532.264363/1‑94 (MQ=60)
                                       GGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATC                                     <  SRR4302532.6164139/94‑1 (MQ=60)
                                       GGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATC                                     <  SRR4302532.2757391/94‑1 (MQ=60)
                                       GGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATC                                     >  SRR4302532.3681776/1‑94 (MQ=60)
                                       GGTACATATATGGCAAAAGCCCGAAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATC                                     <  SRR4302532.2646935/94‑1 (MQ=60)
                                         TACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGG                                   >  SRR4302532.5137121/1‑94 (MQ=60)
                                         TACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGG                                   <  SRR4302532.4608406/94‑1 (MQ=60)
                                          ACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGG                                  >  SRR4302532.6167455/1‑94 (MQ=60)
                                           CATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGG                                 >  SRR4302532.3944/1‑94 (MQ=60)
                                           CATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGG                                 >  SRR4302532.6012801/1‑94 (MQ=60)
                                           CATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGG                                 <  SRR4302532.4105087/94‑1 (MQ=60)
                                             TATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAT                               <  SRR4302532.4012325/94‑1 (MQ=60)
                                             TATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAT                               <  SRR4302532.6013074/94‑1 (MQ=60)
                                             TATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAT                               <  SRR4302532.710480/94‑1 (MQ=60)
                                               TATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGG                             <  SRR4302532.3264911/94‑1 (MQ=60)
                                                 TGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTT                           <  SRR4302532.5897878/94‑1 (MQ=60)
                                                 TGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTT                           <  SRR4302532.6300765/94‑1 (MQ=60)
                                                   GCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTAT                         <  SRR4302532.3479377/94‑1 (MQ=60)
                                                   GCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTAT                         <  SRR4302532.900730/94‑1 (MQ=60)
                                                    CAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATT                        >  SRR4302532.2855111/1‑94 (MQ=60)
                                                      AAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTT                      >  SRR4302532.3725681/1‑94 (MQ=60)
                                                         GCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATT                   >  SRR4302532.2583028/1‑94 (MQ=60)
                                                         GCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATT                   >  SRR4302532.485748/1‑94 (MQ=60)
                                                          CCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTG                  <  SRR4302532.1605565/94‑1 (MQ=60)
                                                          CCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTG                  >  SRR4302532.579600/1‑94 (MQ=60)
                                                            cggaaaatccgggctTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTC                <  SRR4302532.3057609/79‑1 (MQ=60)
                                                             ggaaaatccgggctTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCA               >  SRR4302532.1583036/15‑94 (MQ=60)
                                                                aaatccgggctTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATT            <  SRR4302532.5586381/83‑1 (MQ=60)
                                                                aaatccgggctTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATT            >  SRR4302532.5299622/12‑94 (MQ=60)
                                                                 aatccgggctTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTT           >  SRR4302532.5830244/11‑94 (MQ=60)
                                                                  atccgggctTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTA          >  SRR4302532.2379702/10‑94 (MQ=60)
                                                                  atccgggctTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTA          >  SRR4302532.4636875/10‑94 (MQ=60)
                                                                   tccgggctTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTAT         <  SRR4302532.6179839/86‑1 (MQ=60)
                                                                   tccgggctTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTAT         <  SRR4302532.5368813/86‑1 (MQ=60)
                                                                      gggctTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAG      <  SRR4302532.5038049/89‑1 (MQ=60)
                                                                         ctTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGA   >  SRR4302532.839555/3‑94 (MQ=60)
                                                                          TTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGAT  >  SRR4302532.1792266/1‑94 (MQ=60)
                                                                          TTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGAT  <  SRR4302532.2512962/94‑1 (MQ=60)
                                                                           |                                                                                             
TCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTG‑‑CTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGAT  >  NC_000913/3957478‑3957642

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