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 I2 R1 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,857,508 T→C E79G (GAA→GGA)  yidK ← putative cotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,857,5080TC99.6% 696.0 / ‑4.7 223E79G (GAA→GGA) yidKputative cotransporter
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (121/101);  total (122/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-01

Only 100 of 230 total aligned reads displayed.

AGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAA  >  NC_000913/3857418‑3857597
                                                                                          |                                                                                         
aGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTccc                                                                                         >  2:3478484/1‑93 (MQ=255)
aGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTccc                                                                                         >  2:1689018/1‑93 (MQ=255)
 gAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTcccc                                                                                        <  2:1661293/93‑1 (MQ=255)
  aaaaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCa                                                                                       >  2:4476225/1‑93 (MQ=255)
    aaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGc                                                                                     >  1:1004715/1‑93 (MQ=255)
    aaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGc                                                                                     >  2:2685130/1‑93 (MQ=255)
     aTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGcc                                                                                    >  1:3098505/1‑93 (MQ=255)
      tCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGccc                                                                                   >  1:2151122/1‑93 (MQ=255)
        ggggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCAt                                                                                 >  1:658241/1‑93 (MQ=255)
        ggggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCAt                                                                                 <  1:3373522/93‑1 (MQ=255)
           gATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCAc                                                                              <  2:3689560/93‑1 (MQ=255)
             tggtggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCg                                                                            >  1:3087573/1‑93 (MQ=255)
               gtggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGAc                                                                          >  1:2059568/1‑93 (MQ=255)
                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  1:3514045/93‑1 (MQ=255)
                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  2:1379366/93‑1 (MQ=255)
                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  2:4002050/93‑1 (MQ=255)
                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  1:1423851/93‑1 (MQ=255)
                 ggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACAt                                                                        <  1:1832419/93‑1 (MQ=255)
                 ggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACAt                                                                        >  2:2044387/1‑93 (MQ=255)
                    aaTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATGACCGACATGcc                                                                     >  1:4213119/1‑93 (MQ=255)
                    aaTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGcc                                                                     >  2:2179568/1‑93 (MQ=255)
                     aTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCg                                                                    >  2:2422075/1‑93 (MQ=255)
                      tGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGc                                                                   >  1:866405/1‑93 (MQ=255)
                       gCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCt                                                                  >  2:1275378/1‑93 (MQ=255)
                         cgcgCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCttt                                                                <  2:3227078/93‑1 (MQ=255)
                           cgcTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTg                                                              <  1:3371548/93‑1 (MQ=255)
                             cTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTa                                                            <  2:2164979/93‑1 (MQ=255)
                             cTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTa                                                            <  2:1623203/93‑1 (MQ=255)
                              ttCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAg                                                           >  1:3045542/1‑93 (MQ=255)
                                  gATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCt                                                       >  2:915167/1‑93 (MQ=255)
                                    tAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGc                                                     >  2:1586534/1‑93 (MQ=255)
                                      gcgcGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGccc                                                   <  1:2272712/93‑1 (MQ=255)
                                        gcgGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGTCTGCCCgg                                                 >  1:1220799/1‑93 (MQ=255)
                                           gTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaaa                                              <  1:4478960/93‑1 (MQ=255)
                                            tAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAg                                             <  1:1000592/93‑1 (MQ=255)
                                             aaaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACGTGCCGCTTTTGTAGGACTGCCCGGAAAgg                                            >  2:1533727/1‑93 (MQ=255)
                                             aaaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAgg                                            >  2:2281458/1‑93 (MQ=255)
                                             aaaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAgg                                            >  2:4169260/1‑93 (MQ=255)
                                             aaaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAgg                                            >  2:4210806/1‑93 (MQ=255)
                                              aaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGc                                           >  1:4064979/1‑93 (MQ=255)
                                              aaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGc                                           <  1:1662580/93‑1 (MQ=255)
                                               aaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGcc                                          >  2:486229/1‑93 (MQ=255)
                                               aaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGcc                                          <  1:3143322/93‑1 (MQ=255)
                                                aaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCg                                         >  1:4568093/1‑93 (MQ=255)
                                                aaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCg                                         >  1:2561434/1‑93 (MQ=255)
                                                  gATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGAc                                       >  1:604887/1‑93 (MQ=255)
                                                  aaTTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGAc                                       <  2:3089266/92‑1 (MQ=255)
                                                   aTTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGAcc                                      <  2:3747564/93‑1 (MQ=255)
                                                    ttAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCCGCCCATCACCGACATGCAGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCa                                     >  1:2100422/1‑93 (MQ=255)
                                                    ttAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCa                                     <  1:1938035/93‑1 (MQ=255)
                                                     tAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAg                                    <  2:2859294/93‑1 (MQ=255)
                                                      aGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgt                                   >  1:3773266/1‑93 (MQ=255)
                                                       gcgcGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgtt                                  >  1:4670129/1‑93 (MQ=255)
                                                         gcgAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgttgt                                >  1:2038400/1‑93 (MQ=255)
                                                          cgAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgttgtt                               <  2:3171027/93‑1 (MQ=255)
                                                           gAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTc                              <  1:2173105/93‑1 (MQ=255)
                                                            aGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTcc                             <  1:3607000/93‑1 (MQ=255)
                                                            aGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTcc                             >  2:2486960/1‑93 (MQ=255)
                                                             ggAAGATCAGCGTAACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCg                            >  2:384043/1‑93 (MQ=255)
                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           >  1:2456081/1‑93 (MQ=255)
                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           >  1:4528115/1‑93 (MQ=255)
                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           <  1:4299708/93‑1 (MQ=255)
                                                               aaGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTg                          <  1:1118334/93‑1 (MQ=255)
                                                                 gATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGa                        >  2:140980/1‑93 (MQ=255)
                                                                  aTCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaa                       >  2:1466219/1‑93 (MQ=255)
                                                                   tCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaaa                      <  1:2879871/93‑1 (MQ=255)
                                                                   tCAGCATCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaaa                      <  2:4312215/93‑1 (MQ=255)
                                                                     aGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAgg                    >  2:948744/1‑93 (MQ=255)
                                                                     aGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAgg                    <  2:1857318/93‑1 (MQ=255)
                                                                     aGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAgg                    >  2:1007729/1‑93 (MQ=255)
                                                                      gCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggt                   >  2:2607610/1‑93 (MQ=255)
                                                                      gCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggt                   >  1:1506154/1‑93 (MQ=255)
                                                                      gCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggt                   >  1:3905845/1‑93 (MQ=255)
                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 >  1:1333322/1‑93 (MQ=255)
                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 >  2:74922/1‑93 (MQ=255)
                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 <  2:4351405/93‑1 (MQ=255)
                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 >  1:3688849/1‑93 (MQ=255)
                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 >  2:3505257/1‑93 (MQ=255)
                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTACGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 <  1:2516157/93‑1 (MQ=255)
                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCCACATGCCGCTTTTGTAGGACTGCCCGGGAACGCCGACCAGGTGGTCCGTGGAAAggttg                 >  1:3496112/1‑93 (MQ=255)
                                                                        gTCACCGCTCAAGTCGCTCCCCAGCCCATCCGCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 <  1:1746345/93‑1 (MQ=255)
                                                                         tCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttgg                >  1:3984409/1‑93 (MQ=255)
                                                                          cACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGCAAggctggt               >  2:3436928/1‑93 (MQ=255)
                                                                            ccGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTa             >  2:2373833/1‑93 (MQ=255)
                                                                            ccGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTa             <  1:1550071/93‑1 (MQ=255)
                                                                             cGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAg            >  2:627778/1‑93 (MQ=255)
                                                                              gCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGc           >  1:2994908/1‑93 (MQ=255)
                                                                               cTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCa          <  1:258707/93‑1 (MQ=255)
                                                                                tGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAt         >  2:289627/1‑93 (MQ=255)
                                                                                tGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAt         <  1:4306469/93‑1 (MQ=255)
                                                                                tGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAt         <  2:4565694/93‑1 (MQ=255)
                                                                                 gAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAtt        <  2:85419/93‑1 (MQ=255)
                                                                                    gTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     >  1:1353385/1‑93 (MQ=255)
                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGTAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     <  2:2649901/93‑1 (MQ=255)
                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     >  1:1867707/1‑93 (MQ=255)
                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     <  1:797927/93‑1 (MQ=255)
                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     >  2:3270175/1‑93 (MQ=255)
                                                                                      cACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGa   <  1:3845702/93‑1 (MQ=255)
                                                                                      cACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGa   <  2:3387807/93‑1 (MQ=255)
                                                                                       aCTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCAGTGGAAAGGTTGGTTAGCATTAACGaa  <  2:1130823/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAA  >  NC_000913/3857418‑3857597

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAAGCG  >  NC_000913/3857415‑3857600
                                                                                             |                                                                                             
TCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTC                                                                                               <  SRR4302535.3242068/94‑1 (MQ=60)
 CCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCC                                                                                              >  SRR4302535.2495023/1‑94 (MQ=60)
 CCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCC                                                                                              >  SRR4302535.2366794/1‑94 (MQ=60)
 CCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCC                                                                                              >  SRR4302535.1719013/1‑94 (MQ=60)
    GAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCA                                                                                          <  SRR4302535.4671164/94‑1 (MQ=60)
     AAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAG                                                                                         <  SRR4302535.4333924/94‑1 (MQ=60)
     AAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAG                                                                                         <  SRR4302535.81298/94‑1 (MQ=60)
      AAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGC                                                                                        >  SRR4302535.1343210/1‑94 (MQ=60)
      AAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGC                                                                                        >  SRR4302535.1011852/1‑94 (MQ=60)
       AATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCC                                                                                       >  SRR4302535.3120077/1‑94 (MQ=60)
        ATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCC                                                                                      >  SRR4302535.2166587/1‑94 (MQ=60)
          CGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCAT                                                                                    >  SRR4302535.662837/1‑94 (MQ=60)
           GGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATC                                                                                   >  SRR4302535.3200768/1‑94 (MQ=60)
           GGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATC                                                                                   <  SRR4302535.3396985/94‑1 (MQ=60)
              GATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCCCCGCTGAAGTCACTCCCCAGCCCATCACC                                                                                <  SRR4302535.598409/94‑1 (MQ=60)
              GATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACC                                                                                >  SRR4302535.4638223/1‑94 (MQ=60)
               ATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCG                                                                               >  SRR4302535.3109078/1‑94 (MQ=60)
               ATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCG                                                                               >  SRR4302535.4379717/1‑94 (MQ=60)
                 GGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGAC                                                                             >  SRR4302535.2074407/1‑94 (MQ=60)
                   TGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACAT                                                                           <  SRR4302535.3064647/94‑1 (MQ=60)
                   TGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACAT                                                                           <  SRR4302535.1434019/94‑1 (MQ=60)
                   TGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACAT                                                                           <  SRR4302535.3538464/94‑1 (MQ=60)
                    GGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATG                                                                          >  SRR4302535.141741/1‑94 (MQ=60)
                    GGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATG                                                                          <  SRR4302535.1845642/94‑1 (MQ=60)
                     GCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGC                                                                         >  SRR4302535.1529668/1‑94 (MQ=60)
                      CAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATGACCGACATGCC                                                                        >  SRR4302535.4241953/1‑94 (MQ=60)
                      CAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCC                                                                        >  SRR4302535.3138392/1‑94 (MQ=60)
                       AATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCG                                                                       >  SRR4302535.1627569/1‑94 (MQ=60)
                        ATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGC                                                                      >  SRR4302535.872532/1‑94 (MQ=60)
                             GCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTG                                                                 <  SRR4302535.672321/94‑1 (MQ=60)
                             GCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTG                                                                 <  SRR4302535.2582997/94‑1 (MQ=60)
                              CGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGT                                                                <  SRR4302535.3395005/94‑1 (MQ=60)
                              CGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGA                                                                >  SRR4302535.2069463/1‑94 (MQ=60)
                               GCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTA                                                               >  SRR4302535.4548968/1‑94 (MQ=60)
                               GCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTA                                                               >  SRR4302535.148231/1‑94 (MQ=60)
                                CTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAG                                                              >  SRR4302535.3066876/1‑94 (MQ=60)
                                CTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAG                                                              >  SRR4302535.2032788/1‑94 (MQ=60)
                                    AGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCT                                                          <  SRR4302535.1752810/94‑1 (MQ=60)
                                    AGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCT                                                          >  SRR4302535.3511663/1‑94 (MQ=60)
                                    AGATAGCGCGCTAAAAAGATTAGCGCTAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCTCCGACATGCCGCTTTTGTAGGCCT                                                          <  SRR4302535.650558/94‑1 (MQ=60)
                                     GATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTG                                                         <  SRR4302535.689451/94‑1 (MQ=60)
                                      ATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGC                                                        >  SRR4302535.1870694/1‑94 (MQ=60)
                                         GCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCG                                                     <  SRR4302535.408754/94‑1 (MQ=60)
                                         GCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCG                                                     <  SRR4302535.2288925/94‑1 (MQ=60)
                                          CGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGTCTGCCCGG                                                    >  SRR4302535.1229548/1‑94 (MQ=60)
                                           GCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGA                                                   >  SRR4302535.1488388/1‑94 (MQ=60)
                                              GTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAG                                                <  SRR4302535.3853820/94‑1 (MQ=60)
                                              GTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAG                                                <  SRR4302535.4509711/94‑1 (MQ=60)
                                              GTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAG                                                <  SRR4302535.991615/94‑1 (MQ=60)
                                               TAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGG                                               <  SRR4302535.1007703/94‑1 (MQ=60)
                                               TAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGG                                               >  SRR4302535.3101227/1‑94 (MQ=60)
                                                AAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGC                                              >  SRR4302535.4092795/1‑94 (MQ=60)
                                                 AAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCC                                             <  SRR4302535.1674592/94‑1 (MQ=60)
                                                 AAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCC                                             <  SRR4302535.2284330/94‑1 (MQ=60)
                                                  AAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCG                                            <  SRR4302535.3165189/94‑1 (MQ=60)
                                                  AAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCG                                            <  SRR4302535.2169167/94‑1 (MQ=60)
                                                  AAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCG                                            >  SRR4302535.2579876/1‑94 (MQ=60)
                                                  AAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCG                                            >  SRR4302535.4599249/1‑94 (MQ=60)
                                                  AAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAACCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCG                                            <  SRR4302535.3649487/94‑1 (MQ=60)
                                                    AGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGAC                                          >  SRR4302535.609117/1‑94 (MQ=60)
                                                    AGATTAGCGCGAGGAAGATCAGCGTCACCCCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGACTGCCCGGACAAGCCGAC                                          >  SRR4302535.2040539/1‑94 (MQ=60)
                                                     GATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACC                                         >  SRR4302535.2514067/1‑94 (MQ=60)
                                                      ATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCCGCCCATCACCGACATGCAGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCA                                        >  SRR4302535.2115569/1‑94 (MQ=60)
                                                      ATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGGCCGGAAAGGCCGACCA                                        >  SRR4302535.2115568/1‑94 (MQ=60)
                                                       TTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAG                                       <  SRR4302535.1951988/94‑1 (MQ=60)
                                                        TAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGT                                      >  SRR4302535.3799292/1‑94 (MQ=60)
                                                        TAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGT                                      <  SRR4302535.3914655/94‑1 (MQ=60)
                                                         AGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTT                                     >  SRR4302535.4701676/1‑94 (MQ=60)
                                                           CGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGT                                   >  SRR4302535.2053079/1‑94 (MQ=60)
                                                              GAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCC                                <  SRR4302535.3938801/94‑1 (MQ=60)
                                                              GAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCC                                >  SRR4302535.3216624/1‑94 (MQ=60)
                                                              GAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCC                                <  SRR4302535.1009919/94‑1 (MQ=60)
                                                              GAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCC                                <  SRR4302535.2188650/94‑1 (MQ=60)
                                                               AGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCG                               >  SRR4302535.1012944/1‑94 (MQ=60)
                                                               AGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCG                               <  SRR4302535.3632020/94‑1 (MQ=60)
                                                                GGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGT                              >  SRR4302535.2473789/1‑94 (MQ=60)
                                                                GGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGT                              >  SRR4302535.4559125/1‑94 (MQ=60)
                                                                 GAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTG                             <  SRR4302535.4329131/94‑1 (MQ=60)
                                                                  AAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGG                            <  SRR4302535.1126312/94‑1 (MQ=60)
                                                                    GATCAGCGTCACCGCTGAAGTCCCTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAA                          <  SRR4302535.3396580/94‑1 (MQ=60)
                                                                     ATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAA                         >  SRR4302535.3068742/1‑94 (MQ=60)
                                                                     ATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAA                         >  SRR4302535.3527343/1‑94 (MQ=60)
                                                                      TCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAG                        <  SRR4302535.2900260/94‑1 (MQ=60)
                                                                      ACAGCGTCACCGCTGAAGTCCCTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAG                        <  SRR4302535.3727016/94‑1 (MQ=60)
                                                                       CAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGG                       >  SRR4302535.12344/1‑94 (MQ=60)
                                                                        AGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGCAAGGT                      >  SRR4302535.1751308/1‑94 (MQ=60)
                                                                        AGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGT                      >  SRR4302535.3932605/1‑94 (MQ=60)
                                                                        AGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGT                      >  SRR4302535.1516983/1‑94 (MQ=60)
                                                                        AGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGT                      >  SRR4302535.212869/1‑94 (MQ=60)
                                                                         GCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTT                     <  SRR4302535.4587843/94‑1 (MQ=60)
                                                                         GCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTT                     <  SRR4302535.1929825/94‑1 (MQ=60)
                                                                         GCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTT                     <  SRR4302535.4046084/94‑1 (MQ=60)
                                                                         GCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTT                     <  SRR4302535.283851/94‑1 (MQ=60)
                                                                          CGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTG                    >  SRR4302535.3714377/1‑94 (MQ=60)
                                                                          CGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTG                    >  SRR4302535.1342831/1‑94 (MQ=60)
                                                                          CGTCACCGCTGAAGTCACTCCCCAGCCCATCACCCACATGCCGCTTTTGTAGGACTGCCCGGGAACGCCGACCAGGTGGTCCGTGGAAAGGTTG                    >  SRR4302535.3520409/1‑94 (MQ=60)
                                                                           GTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGG                   >  SRR4302535.4011633/1‑94 (MQ=60)
                                                                           GTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTACGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGG                   <  SRR4302535.2534293/94‑1 (MQ=60)
                                                                           GTCACCGCTCAAGTCGCTCCCCAGCCCATCCGCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGG                   <  SRR4302535.1758979/94‑1 (MQ=60)
                                                                            TCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGT                  >  SRR4302535.4236791/1‑94 (MQ=60)
                                                                            TCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGT                  >  SRR4302535.1778977/1‑94 (MQ=60)
                                                                             CACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTT                 >  SRR4302535.3110073/1‑94 (MQ=60)
                                                                              ACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTA                <  SRR4302535.4291330/94‑1 (MQ=60)
                                                                               CCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAG               >  SRR4302535.2557062/1‑94 (MQ=60)
                                                                               CCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAG               <  SRR4302535.1561237/94‑1 (MQ=60)
                                                                                CGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGC              >  SRR4302535.3015904/1‑94 (MQ=60)
                                                                                 GCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCA             >  SRR4302535.3561618/1‑94 (MQ=60)
                                                                                  CTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAT            >  SRR4302535.2374048/1‑94 (MQ=60)
                                                                                  CTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAT            <  SRR4302535.260497/94‑1 (MQ=60)
                                                                                   TGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATT           <  SRR4302535.4335933/94‑1 (MQ=60)
                                                                                    GAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTA          <  SRR4302535.2939852/94‑1 (MQ=60)
                                                                                     AAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAA         <  SRR4302535.2134809/94‑1 (MQ=60)
                                                                                     AAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAA         <  SRR4302535.2399968/94‑1 (MQ=60)
                                                                                     AAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAA         >  SRR4302535.2653374/1‑94 (MQ=60)
                                                                                     AAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAA         >  SRR4302535.1303967/1‑94 (MQ=60)
                                                                                      AGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAC        >  SRR4302535.1363069/1‑94 (MQ=60)
                                                                                      AGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAC        >  SRR4302535.3815336/1‑94 (MQ=60)
                                                                                      AGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAC        >  SRR4302535.1881169/1‑94 (MQ=60)
                                                                                      AGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAACGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAC        >  SRR4302535.4406288/1‑94 (MQ=60)
                                                                                       GTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACG       >  SRR4302535.3616461/1‑94 (MQ=60)
                                                                                       GTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACG       <  SRR4302535.803527/94‑1 (MQ=60)
                                                                                         CACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAA     <  SRR4302535.3872181/94‑1 (MQ=60)
                                                                                          ACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAAG    <  SRR4302535.1202406/94‑1 (MQ=60)
                                                                                          ACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAAG    >  SRR4302535.3106532/1‑94 (MQ=60)
                                                                                           CTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAAGC   <  SRR4302535.824826/94‑1 (MQ=60)
                                                                                           CTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAAGC   >  SRR4302535.4587367/1‑94 (MQ=60)
                                                                                           CTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAAGC   >  SRR4302535.4572447/1‑94 (MQ=60)
                                                                                            TCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAAGCG  >  SRR4302535.546954/1‑94 (MQ=60)
                                                                                            TCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAAGCG  <  SRR4302535.1298772/94‑1 (MQ=60)
                                                                                            TCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAAGCG  <  SRR4302535.2323358/94‑1 (MQ=60)
                                                                                            TCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAAGCG  <  SRR4302535.3595984/94‑1 (MQ=60)
                                                                                             TCCCAAGCCCATCAACGACATGCCGCTTTAGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCATGGAAAGGTTGGTTAGCATTAACGAAGCG  <  SRR4302535.4712515/94‑1 (MQ=60)
                                                                                             |                                                                                             
TCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAAGCG  >  NC_000913/3857415‑3857600

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