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 1548 369.5 20095359 96.9% 19472402 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,125,359 T→C K226K (AAA→AAG glcD ← glycolate oxidase subunit D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,125,3590TC99.6% 771.5 / NA 241K226K (AAA→AAGglcDglycolate oxidase subunit D
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (126/113);  total (127/114)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

GGCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCGTC  >  NC_000913/3125269‑3125448
                                                                                          |                                                                                         
ggCGATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGctt                                                                                         >  2:7846061/1‑93 (MQ=255)
ggCGATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGctt                                                                                         <  2:5046611/93‑1 (MQ=255)
 gCGATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGcttt                                                                                        <  2:2694742/93‑1 (MQ=255)
  cGATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTa                                                                                       >  2:2294120/1‑93 (MQ=255)
   gatgatGTCACCACCCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTAc                                                                                      <  1:909281/93‑1 (MQ=255)
   gatgatGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTAc                                                                                      <  1:9774745/93‑1 (MQ=255)
   gatgatGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTAc                                                                                      <  1:7543885/93‑1 (MQ=255)
    atgatgTCACTAACAGCAAGCCCTGCTTGTTCTACCGAGACAAACCTGACTAACAGAACCCGCGTGACGGTCGTCTTCGGCAGCAGCGTTAcc                                                                                     <  2:9901586/93‑1 (MQ=255)
    atgatgTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTAcc                                                                                     >  2:138709/1‑93 (MQ=255)
     tgatgTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCg                                                                                    >  1:5141435/1‑93 (MQ=255)
      gatgTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGt                                                                                   >  1:3135146/1‑93 (MQ=255)
      gatgGCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGt                                                                                   <  2:8641407/93‑1 (MQ=255)
       atgTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTc                                                                                  >  2:5542397/1‑93 (MQ=255)
       atgTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTc                                                                                  >  1:9830423/1‑93 (MQ=255)
        tgTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCa                                                                                 <  1:6558196/93‑1 (MQ=255)
         gTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCATGGGCGGCTTCGGCAGCAGCTTTACCGTCAc                                                                                <  1:5505811/93‑1 (MQ=255)
         gTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCAc                                                                                >  2:6809997/1‑93 (MQ=255)
         gTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCAc                                                                                >  1:4732634/1‑93 (MQ=255)
           cACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACtt                                                                              <  2:5173271/93‑1 (MQ=255)
            aCCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTc                                                                             >  1:1076896/1‑93 (MQ=255)
            aCCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTc                                                                             >  2:4608931/1‑93 (MQ=255)
            aCCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGATGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTc                                                                             >  2:9025942/1‑93 (MQ=255)
             ccAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCg                                                                            <  1:2855939/93‑1 (MQ=255)
               aaCCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTTGTTAACAGAACCCGCTCCACGGGCGGCTTCGGCAGGCGCTTTACCGTCACTTCggt                                                                          >  2:2548711/1‑93 (MQ=255)
                   gCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTc                                                                      >  2:21842/1‑93 (MQ=255)
                   gCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTc                                                                      <  1:852981/93‑1 (MQ=255)
                    cAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTcc                                                                     >  2:4637794/1‑92 (MQ=255)
                       gcccGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGc                                                                  <  1:2423331/91‑1 (MQ=255)
                       gcccGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGc                                                                  <  2:7623816/91‑1 (MQ=255)
                        cccGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGcc                                                                 <  1:6507229/92‑1 (MQ=255)
                        cccGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGcc                                                                 <  2:722858/92‑1 (MQ=255)
                         ccGGGTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCATCGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCg                                                                <  2:7365297/93‑1 (MQ=255)
                         ccGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCg                                                                >  2:1392752/1‑93 (MQ=255)
                          cGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCTCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGa                                                               <  2:2532309/93‑1 (MQ=255)
                          cGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGGAGCTTTACCGTCACTTCGGTGGGCACGCCgt                                                               >  1:8576875/1‑92 (MQ=255)
                          cGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGa                                                               >  1:4642725/1‑93 (MQ=255)
                          cGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGa                                                               >  1:874049/1‑93 (MQ=255)
                           cgCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAg                                                              <  1:8268086/92‑1 (MQ=255)
                              tttttCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCAt                                                           >  1:3235741/1‑93 (MQ=255)
                                tttCTACCGAGTCAAAGCTGGCTAACAGACCCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATAc                                                         <  1:2965516/93‑1 (MQ=255)
                                tttCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGATTTACCGTCACTTCGGTGGTCACGCCGCGCATAc                                                         >  1:7246028/1‑93 (MQ=255)
                                 ttCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATAcc                                                        <  1:8127390/93‑1 (MQ=255)
                                   cTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCtt                                                      >  2:1789358/1‑93 (MQ=255)
                                    tACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTc                                                     >  2:8796787/1‑93 (MQ=255)
                                      ccGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCg                                                   >  1:2760307/1‑93 (MQ=255)
                                         aGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATc                                                <  2:9308079/93‑1 (MQ=255)
                                          gTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATcc                                               <  1:2568451/93‑1 (MQ=255)
                                          gTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATcc                                               >  2:3212093/1‑93 (MQ=255)
                                          gTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATcc                                               <  2:8122094/93‑1 (MQ=255)
                                          gTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATcc                                               <  2:407744/93‑1 (MQ=255)
                                            cAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCgg                                             <  1:6032314/93‑1 (MQ=255)
                                            cAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCgg                                             >  2:8148106/1‑93 (MQ=255)
                                              aaGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTg                                           >  1:4833824/1‑93 (MQ=255)
                                              aaGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTg                                           <  2:2040651/93‑1 (MQ=255)
                                                gCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGaa                                         <  1:4177964/93‑1 (MQ=255)
                                                gCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGaa                                         <  1:6381249/93‑1 (MQ=255)
                                                gCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGaa                                         <  2:5768883/93‑1 (MQ=255)
                                                gCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGaa                                         <  2:6662731/93‑1 (MQ=255)
                                                acgtgCTAAGAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGaa                                         <  1:9935244/89‑1 (MQ=255)
                                                 cTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAAc                                        <  2:6185134/93‑1 (MQ=255)
                                                  tGGCTAACAGAACCCGCGCCACGAGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACa                                       <  2:243337/93‑1 (MQ=255)
                                                   ggCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAg                                      <  1:7101613/93‑1 (MQ=255)
                                                     cTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAgcg                                    >  1:5369652/1‑93 (MQ=255)
                                                      tAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAgcgc                                   >  1:2316735/1‑93 (MQ=255)
                                                        aCAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCca                                 >  2:3095473/1‑93 (MQ=255)
                                                        aCAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCca                                 >  1:182316/1‑93 (MQ=255)
                                                         cAGACCCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCcag                                <  1:5861643/93‑1 (MQ=255)
                                                          aGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCcagc                               <  1:337417/93‑1 (MQ=255)
                                                             aCCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAgg                            <  1:4354351/93‑1 (MQ=255)
                                                                cgcgCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCa                         >  2:7723821/1‑93 (MQ=255)
                                                                   gcCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCaaaa                      <  2:2809961/93‑1 (MQ=255)
                                                                   gcCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCaaaa                      >  1:9225761/1‑93 (MQ=255)
                                                                      ccGGGCGGCTGCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATCCCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCa                   <  2:1867239/92‑1 (MQ=255)
                                                                      aCGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCa                   >  1:6955537/1‑93 (MQ=255)
                                                                       cGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAg                  >  1:8109332/1‑93 (MQ=255)
                                                                        gggCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAgg                 <  1:1496064/93‑1 (MQ=255)
                                                                        gggCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAgg                 >  2:9161993/1‑93 (MQ=255)
                                                                         ggcggcTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAggg                >  1:9546796/1‑92 (MQ=255)
                                                                         ggcggcTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAggg                >  1:8308454/1‑92 (MQ=255)
                                                                         ggcggcTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGt                <  2:1245352/93‑1 (MQ=255)
                                                                         ggcggcTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGt                >  1:6536617/1‑93 (MQ=255)
                                                                         ggcggcTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGt                <  1:8982309/93‑1 (MQ=255)
                                                                         ggcggcGTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGt                <  1:3979741/93‑1 (MQ=255)
                                                                             gcTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGGAGGTCACAACCCGGAGAAt            >  2:6543857/1‑93 (MQ=255)
                                                                              cTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATc           <  1:7959308/93‑1 (MQ=255)
                                                                               ttCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATcc          <  2:8211597/93‑1 (MQ=255)
                                                                               ttCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATcc          <  1:2678922/93‑1 (MQ=255)
                                                                               ttCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATcc          >  1:2205787/1‑93 (MQ=255)
                                                                                 cGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAg        >  2:5918984/1‑93 (MQ=255)
                                                                                 cGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAg        <  2:4944730/93‑1 (MQ=255)
                                                                                 cGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAg        <  2:3926382/93‑1 (MQ=255)
                                                                                 cGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAg        <  2:321084/93‑1 (MQ=255)
                                                                                  ggCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAgc       >  2:10043754/1‑93 (MQ=255)
                                                                                  ggCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAgc       <  2:3697503/93‑1 (MQ=255)
                                                                                  ggCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAgc       >  1:572182/1‑93 (MQ=255)
                                                                                    cagcagCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAgcgc     >  1:6996741/1‑93 (MQ=255)
                                                                                     agcagcTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGcgcg    >  1:779976/1‑93 (MQ=255)
                                                                                     agcagcTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGcgcg    >  1:3698468/1‑93 (MQ=255)
                                                                                     agcagcTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGcgcg    >  1:7385639/1‑93 (MQ=255)
                                                                                       cagCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCGTc  <  2:5401511/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCGTC  >  NC_000913/3125269‑3125448

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTGGCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCGTCCG  >  NC_000913/3125267‑3125450
                                                                                            |                                                                                           
ttgcggatggtgtgacccacgGCAAGCCCGGCTTTTTCTACCGAGACAAAGCTGGCTAACAGAACCGGCGCCCCCGGCGGCTTCGGCAGCAGCT                                                                                            <  SRR4302540.3568387/73‑1 (MQ=58)
TTGGCGATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCT                                                                                            <  SRR4302540.3568386/94‑1 (MQ=60)
TTGGCGATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCT                                                                                            <  SRR4302540.3346374/94‑1 (MQ=60)
TTGGCGATGATGGCACCAACCGCAAGGCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCT                                                                                            <  SRR4302540.3568385/94‑1 (MQ=60)
 TGGCGATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTT                                                                                           >  SRR4302540.2407031/1‑94 (MQ=60)
   GCGATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTA                                                                                         <  SRR4302540.8665456/94‑1 (MQ=60)
    CGATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTAC                                                                                        <  SRR4302540.610460/94‑1 (MQ=60)
    CGATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTAC                                                                                        >  SRR4302540.2610219/1‑94 (MQ=60)
    CGATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTAC                                                                                        >  SRR4302540.8822922/1‑94 (MQ=60)
     GATGATGTCACCACCCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACC                                                                                       <  SRR4302540.918919/94‑1 (MQ=60)
     GATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACC                                                                                       <  SRR4302540.10000191/94‑1 (MQ=60)
     GATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACC                                                                                       >  SRR4302540.3356362/1‑94 (MQ=60)
     GATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACC                                                                                       <  SRR4302540.7628421/94‑1 (MQ=60)
     GATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACC                                                                                       <  SRR4302540.9889329/94‑1 (MQ=60)
      ATGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCG                                                                                      >  SRR4302540.5197123/1‑94 (MQ=60)
       TGATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGT                                                                                     >  SRR4302540.3169345/1‑94 (MQ=60)
        GATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTNCCGTC                                                                                    <  SRR4302540.8164709/94‑1 (MQ=60)
        GATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTC                                                                                    <  SRR4302540.515348/94‑1 (MQ=60)
        GATGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTC                                                                                    >  SRR4302540.9945398/1‑94 (MQ=60)
          TGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCAC                                                                                  >  SRR4302540.4783172/1‑94 (MQ=60)
          TGTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCAC                                                                                  <  SRR4302540.6630909/94‑1 (MQ=60)
           GTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCATGGGCGGCTTCGGCAGCAGCTTTACCGTCACT                                                                                 <  SRR4302540.5566445/94‑1 (MQ=60)
           GTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACT                                                                                 >  SRR4302540.9205885/1‑94 (MQ=60)
           GTCACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTAACCGTCACT                                                                                 >  SRR4302540.839859/1‑94 (MQ=60)
           GTCACCAACCGCAAGCCCGGCGCCCTCCACCGAGTCTAAGCGGGCTAACAGAACCCGCGCCAGGGGCGGGTTCGGCAGCAGCTTTACCGTCACT                                                                                 <  SRR4302540.3538296/94‑1 (MQ=36)
             CACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCCGCAGCTTTACCGTCACTTC                                                                               >  SRR4302540.8381618/1‑94 (MQ=60)
             CACCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTC                                                                               >  SRR4302540.1087902/1‑94 (MQ=60)
               CCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGG                                                                             <  SRR4302540.2886879/94‑1 (MQ=60)
               CCAACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGG                                                                             >  SRR4302540.5439068/1‑94 (MQ=60)
                 AACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGACCCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTG                                                                           <  SRR4302540.2263596/94‑1 (MQ=60)
                 AACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAAAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTG                                                                           <  SRR4302540.497044/94‑1 (MQ=60)
                  ACCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTCGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGG                                                                          >  SRR4302540.7461287/1‑94 (MQ=60)
                   CCGCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGT                                                                         <  SRR4302540.9118631/94‑1 (MQ=60)
                     GCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCA                                                                       <  SRR4302540.862037/94‑1 (MQ=60)
                     GCAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCA                                                                       >  SRR4302540.8212023/1‑94 (MQ=60)
                      CAAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCAC                                                                      <  SRR4302540.7304472/94‑1 (MQ=60)
                       AAGCCCGGCTTTTTCTACCGAGTCAAAGCTTGCTAACAGAACCCGCGCCACGTGCGGCTTTGGCAGCAGCTTTACCGTCACTTCGGTGGTCACG                                                                     >  SRR4302540.9076513/1‑94 (MQ=60)
                       AAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACG                                                                     >  SRR4302540.2926974/1‑94 (MQ=60)
                       AAGCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACG                                                                     >  SRR4302540.7235822/1‑94 (MQ=60)
                         GCCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCC                                                                   <  SRR4302540.2449439/94‑1 (MQ=60)
                          CCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCG                                                                  <  SRR4302540.6579544/94‑1 (MQ=60)
                          CCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCG                                                                  <  SRR4302540.6579545/94‑1 (MQ=60)
                           CCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGGAGCTTTACCGTCACTTCGGTGGGCACGCCGT                                                                 >  SRR4302540.8676559/1‑94 (MQ=60)
                           CCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGG                                                                 >  SRR4302540.6590537/1‑94 (MQ=60)
                           CCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGA                                                                 >  SRR4302540.883323/1‑94 (MQ=60)
                           CCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGA                                                                 >  SRR4302540.4692429/1‑94 (MQ=60)
                            CGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAG                                                                <  SRR4302540.8669442/94‑1 (MQ=60)
                            CGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAG                                                                <  SRR4302540.9728764/94‑1 (MQ=60)
                             GGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGC                                                               <  SRR4302540.8364155/94‑1 (MQ=60)
                             CGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGC                                                               <  SRR4302540.8364154/94‑1 (MQ=60)
                              GCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCA                                                              <  SRR4302540.3078522/94‑1 (MQ=60)
                               CTTTTTCTACCGAGTCAAAGCTGGCTAGCAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCCCTTCGGTGGTCACGCCGAGCAT                                                             <  SRR4302540.2161951/94‑1 (MQ=60)
                               CTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCAT                                                             >  SRR4302540.3270929/1‑94 (MQ=60)
                                TTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATA                                                            >  SRR4302540.4794335/1‑94 (MQ=60)
                                 TTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATAC                                                           <  SRR4302540.5993851/94‑1 (MQ=60)
                                 TTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGATTTACCGTCACTTCGGTGGTCACGCCGCGCATAC                                                           >  SRR4302540.7326871/1‑94 (MQ=60)
                                  TTTCTACCGAGTCAAAGCTGGCTAACAGACCCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACC                                                          <  SRR4302540.2998242/94‑1 (MQ=60)
                                   TTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCT                                                         >  SRR4302540.6611344/1‑94 (MQ=60)
                                   TTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCT                                                         <  SRR4302540.8221722/94‑1 (MQ=60)
                                    TCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTT                                                        >  SRR4302540.10048585/1‑94 (MQ=60)
                                    TCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTT                                                        <  SRR4302540.7175462/94‑1 (MQ=60)
                                     CTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTC                                                       >  SRR4302540.3244573/1‑94 (MQ=60)
                                       ACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCG                                                     >  SRR4302540.2790134/1‑94 (MQ=60)
                                       ACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCG                                                     >  SRR4302540.8107696/1‑94 (MQ=60)
                                         CGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGAT                                                   <  SRR4302540.4730927/94‑1 (MQ=60)
                                           AGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCC                                                 <  SRR4302540.1788679/94‑1 (MQ=60)
                                           AGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCC                                                 >  SRR4302540.1195208/1‑94 (MQ=60)
                                            GTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCG                                                <  SRR4302540.2596147/94‑1 (MQ=60)
                                            GTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCG                                                <  SRR4302540.3379246/94‑1 (MQ=60)
                                             TCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGG                                               <  SRR4302540.6559972/94‑1 (MQ=60)
                                             TCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTACGCTCCGG                                               >  SRR4302540.6097784/1‑94 (MQ=60)
                                              CAAAGCTGGCTAACAGAGCCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGT                                              <  SRR4302540.2784612/94‑1 (MQ=60)
                                              CAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGT                                              >  SRR4302540.3250892/1‑94 (MQ=60)
                                              CAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGT                                              <  SRR4302540.6099262/94‑1 (MQ=60)
                                               AAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTG                                             >  SRR4302540.9448874/1‑94 (MQ=60)
                                               AAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTG                                             >  SRR4302540.4885608/1‑94 (MQ=60)
                                                 AGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAA                                           <  SRR4302540.2230628/94‑1 (MQ=60)
                                                 AGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAA                                           <  SRR4302540.913007/94‑1 (MQ=60)
                                                  acgtGCTAAGAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAAC                                          <  SRR4302540.10051190/90‑1 (MQ=60)
                                                  GCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAAC                                          <  SRR4302540.4223323/94‑1 (MQ=60)
                                                  GCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAAC                                          >  SRR4302540.4257823/1‑94 (MQ=60)
                                                  GCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAAC                                          <  SRR4302540.6452292/94‑1 (MQ=60)
                                                   CTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACA                                         >  SRR4302540.5748307/1‑94 (MQ=60)
                                                    TGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAG                                        >  SRR4302540.6768115/1‑94 (MQ=60)
                                                     GGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGC                                       <  SRR4302540.7180628/94‑1 (MQ=60)
                                                      GCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCG                                      >  SRR4302540.5545894/1‑94 (MQ=60)
                                                      GCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCG                                      >  SRR4302540.5428667/1‑94 (MQ=60)
                                                      GCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCG                                      >  SRR4302540.9708636/1‑94 (MQ=60)
                                                       CTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGC                                     >  SRR4302540.2341677/1‑94 (MQ=60)
                                                         AACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCA                                   >  SRR4302540.184132/1‑94 (MQ=60)
                                                           CAGACCCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGC                                 <  SRR4302540.5926463/94‑1 (MQ=60)
                                                            AGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCA                                <  SRR4302540.340820/94‑1 (MQ=60)
                                                            AGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCA                                <  SRR4302540.4396996/94‑1 (MQ=60)
                                                            AGAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCA                                >  SRR4302540.9933761/1‑94 (MQ=60)
                                                             GAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAG                               <  SRR4302540.2988594/94‑1 (MQ=60)
                                                             GAACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGCTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAG                               <  SRR4302540.4071685/94‑1 (MQ=60)
                                                               ACCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGT                             <  SRR4302540.4401715/94‑1 (MQ=60)
                                                                CCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTC                            >  SRR4302540.907783/1‑94 (MQ=60)
                                                                CCCGCGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTC                            >  SRR4302540.8299113/1‑94 (MQ=60)
                                                                   GCTCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAA                         <  SRR4302540.708926/94‑1 (MQ=60)
                                                                    CGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAA                        >  SRR4302540.6961097/1‑94 (MQ=60)
                                                                    CGCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAA                        >  SRR4302540.9336375/1‑94 (MQ=60)
                                                                     GCCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAAC                       >  SRR4302540.794629/1‑94 (MQ=60)
                                                                      CCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACC                      >  SRR4302540.3875865/1‑94 (MQ=60)
                                                                      CCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACC                      <  SRR4302540.4581082/94‑1 (MQ=60)
                                                                      CCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACC                      <  SRR4302540.1834923/94‑1 (MQ=60)
                                                                      CCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACC                      >  SRR4302540.2496378/1‑94 (MQ=60)
                                                                      CCACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATACAGTGAACAGCGCCAGCACGTCAAAACC                      >  SRR4302540.562313/1‑94 (MQ=60)
                                                                       CACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCTGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGCCAAAACCA                     >  SRR4302540.7032926/1‑94 (MQ=60)
                                                                       CACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCA                     <  SRR4302540.4434967/94‑1 (MQ=60)
                                                                       CACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCA                     >  SRR4302540.5921136/1‑94 (MQ=60)
                                                                       CACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCA                     >  SRR4302540.7032924/1‑94 (MQ=60)
                                                                        ACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCGGCAGGTAAAAACCAG                    >  SRR4302540.8329335/1‑94 (MQ=60)
                                                                        ACGGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAG                    >  SRR4302540.8203454/1‑94 (MQ=60)
                                                                          GGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGT                  >  SRR4302540.6609164/1‑94 (MQ=60)
                                                                          GGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGT                  <  SRR4302540.1512529/94‑1 (MQ=60)
                                                                          GGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGG                  >  SRR4302540.9659740/1‑94 (MQ=60)
                                                                          GGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGG                  >  SRR4302540.8404916/1‑94 (MQ=60)
                                                                          GGGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGG                  >  SRR4302540.6609167/1‑94 (MQ=60)
                                                                           GGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTG                 <  SRR4302540.9089438/94‑1 (MQ=60)
                                                                           GGCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTG                 <  SRR4302540.5071110/94‑1 (MQ=60)
                                                                           GGCGGCGTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTG                 <  SRR4302540.4022610/94‑1 (MQ=60)
                                                                            GCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGA                >  SRR4302540.8267946/1‑94 (MQ=60)
                                                                            GCGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCAGACCTTGCGATCCGGTGAGCAGCGCCcgcgggtcaagacgcggaga                >  SRR4302540.180350/1‑74 (MQ=60)
                                                                             CGGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAA               >  SRR4302540.76191/1‑94 (MQ=60)
                                                                              GGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAAT              >  SRR4302540.7697148/1‑94 (MQ=60)
                                                                              GGCTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAAT              <  SRR4302540.5836436/94‑1 (MQ=60)
                                                                                CTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCC            >  SRR4302540.2375523/1‑94 (MQ=60)
                                                                                CTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCC            >  SRR4302540.2229618/1‑94 (MQ=60)
                                                                                CTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCC            <  SRR4302540.2806458/94‑1 (MQ=60)
                                                                                CTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCC            <  SRR4302540.8052076/94‑1 (MQ=60)
                                                                                CTTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCC            <  SRR4302540.8052075/94‑1 (MQ=60)
                                                                                CCTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCC            <  SRR4302540.9449808/94‑1 (MQ=60)
                                                                                 TTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCA           <  SRR4302540.2707906/94‑1 (MQ=60)
                                                                                 TTCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCA           <  SRR4302540.3749504/94‑1 (MQ=60)
                                                                                  TCGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAG          <  SRR4302540.2606137/94‑1 (MQ=60)
                                                                                   CGGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGC         >  SRR4302540.578393/1‑94 (MQ=60)
                                                                                    GGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCG        <  SRR4302540.5419252/94‑1 (MQ=60)
                                                                                    GGCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCG        >  SRR4302540.2802156/1‑94 (MQ=60)
                                                                                     GCAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGC       >  SRR4302540.7074397/1‑94 (MQ=60)
                                                                                      CAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCG      >  SRR4302540.7467927/1‑94 (MQ=60)
                                                                                      CAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCG      >  SRR4302540.788284/1‑94 (MQ=60)
                                                                                      CAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCG      >  SRR4302540.7570848/1‑94 (MQ=60)
                                                                                      CAGCAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCG      >  SRR4302540.3738841/1‑94 (MQ=60)
                                                                                         CAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCGTCC   >  SRR4302540.6017273/1‑94 (MQ=60)
                                                                                         CAGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCGTCC   <  SRR4302540.4798738/94‑1 (MQ=60)
                                                                                          AGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCGTCCG  >  SRR4302540.2730078/1‑94 (MQ=60)
                                                                                          AGCTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCGTCCG  >  SRR4302540.6270399/1‑94 (MQ=60)
                                                                                            |                                                                                           
TTGGCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGTCACTTCGGTGGTCACGCCGAGCATACCTTCCGATCCGGTGAACAGCGCCAGCAGGTCAAAACCAGGTGAATCCAGCGCGTCCG  >  NC_000913/3125267‑3125450

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