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 1624 250.6 12620048 98.0% 12367647 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,460,914 T→C Y457C (TAC→TGC)  treC ← trehalase 6‑P hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,460,9140TC99.7% 924.7 / ‑5.0 283Y457C (TAC→TGC) treCtrehalase 6‑P hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (155/127);  minor base A (1/0);  total (156/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 287 total aligned reads displayed.

ATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCCAGCCGT  >  NC_000913/4460825‑4461002
                                                                                         |                                                                                        
  aCAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAac                                                                                     >  2:6015058/1‑93 (MQ=255)
   cAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAaca                                                                                    >  2:4807839/1‑93 (MQ=255)
    aGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAacac                                                                                   >  2:4988619/1‑93 (MQ=255)
     gggCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACAcc                                                                                  >  1:2950533/1‑93 (MQ=255)
     gggCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACAcc                                                                                  <  2:3751052/93‑1 (MQ=255)
     gggCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACAcc                                                                                  <  1:2476915/93‑1 (MQ=255)
     gggCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACAcc                                                                                  <  2:39721/93‑1 (MQ=255)
      ggCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCg                                                                                 <  1:5373402/93‑1 (MQ=255)
      ggCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCg                                                                                 <  1:3193440/93‑1 (MQ=255)
       gCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGa                                                                                >  2:1515335/1‑93 (MQ=255)
       gCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGa                                                                                >  2:5529065/1‑93 (MQ=255)
       gCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGa                                                                                <  1:2714137/93‑1 (MQ=255)
        cTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAg                                                                               <  1:4345196/93‑1 (MQ=255)
         tGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAgg                                                                              <  2:3676240/93‑1 (MQ=255)
          gTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGa                                                                             >  1:1834515/1‑93 (MQ=255)
          gTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGa                                                                             <  1:4840938/93‑1 (MQ=255)
          gTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGa                                                                             >  1:4032184/1‑93 (MQ=255)
           tttGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGaa                                                                            >  1:946013/1‑93 (MQ=255)
           tttGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGGAGGTGCAAAACACCGAGGaa                                                                            >  1:2098725/1‑93 (MQ=255)
                 agcagATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCACAACACCGAGGAATCATCg                                                                      >  2:3702465/1‑93 (MQ=255)
                 agcagATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCg                                                                      >  2:4435334/1‑93 (MQ=255)
                     gATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCa                                                                  >  2:5291917/1‑93 (MQ=255)
                     gATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCa                                                                  >  2:1325528/1‑93 (MQ=255)
                       tCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAgc                                                                <  1:1397236/93‑1 (MQ=255)
                         cTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAgcgc                                                              >  2:4754988/1‑93 (MQ=255)
                         cTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAgcgc                                                              >  1:1543900/1‑93 (MQ=255)
                          tGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGcgcg                                                             >  2:3362482/1‑93 (MQ=255)
                           ggTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCgg                                                            >  2:1021498/1‑93 (MQ=255)
                           ggTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCgg                                                            >  2:6172717/1‑93 (MQ=255)
                            gTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGc                                                           <  1:6087733/93‑1 (MQ=255)
                            gTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGc                                                           <  1:4194619/93‑1 (MQ=255)
                              aaTTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCtt                                                         <  1:1152032/93‑1 (MQ=255)
                              aaTTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCtt                                                         <  1:425405/93‑1 (MQ=255)
                              aaTTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCtt                                                         >  1:4162005/1‑93 (MQ=255)
                               aTTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTc                                                        >  2:2394044/1‑93 (MQ=255)
                               aTTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTc                                                        <  1:2852191/93‑1 (MQ=255)
                                 tGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACTGAGGAATCATCGGCCAGCGCGGCTTCTa                                                      <  1:3527911/93‑1 (MQ=255)
                                    cccATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGt                                                   >  1:221800/1‑93 (MQ=255)
                                     ccATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGtt                                                  >  2:3768301/1‑93 (MQ=255)
                                      cATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTg                                                 >  1:4243533/1‑93 (MQ=255)
                                      cATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTg                                                 >  2:825312/1‑93 (MQ=255)
                                       aTGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGa                                                <  1:1597787/93‑1 (MQ=255)
                                       aTGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGa                                                >  2:5207747/1‑93 (MQ=255)
                                       aTGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGa                                                >  1:1999197/1‑93 (MQ=255)
                                       aTGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGa                                                >  2:231689/1‑93 (MQ=255)
                                         gTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGAtt                                              <  1:5462728/93‑1 (MQ=255)
                                         gTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGAtt                                              >  2:3408603/1‑93 (MQ=255)
                                          tCAGGATGGCTTCCTTCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGAttt                                             <  1:1757924/93‑1 (MQ=255)
                                          tCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGAttt                                             <  2:5903307/93‑1 (MQ=255)
                                          tCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGATGAttt                                             >  1:4515427/1‑93 (MQ=255)
                                           cAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATttg                                            <  1:4898775/93‑1 (MQ=255)
                                             ggATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATttgtt                                          <  2:1282792/93‑1 (MQ=255)
                                             ggATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCAGCCAGCGCGGCTTCTACGTTGATttgtt                                          >  1:101018/1‑93 (MQ=255)
                                                 ggCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATa                                      <  2:1864724/93‑1 (MQ=255)
                                                  gCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAg                                     <  1:5250961/93‑1 (MQ=255)
                                                  gCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAg                                     >  1:578219/1‑93 (MQ=255)
                                                  gCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAg                                     >  1:636248/1‑93 (MQ=255)
                                                   cTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGt                                    <  1:3636535/93‑1 (MQ=255)
                                                   cTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGt                                    <  1:3784976/93‑1 (MQ=255)
                                                   cTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGt                                    <  2:3037944/93‑1 (MQ=255)
                                                   cTTCCTGCTTACGCAGTGCGATTAACTTTTGGGAGGTGCAAAACACCCAGGAATCATCGGCCAGCGCGGCTTCTACGTTGCTTTGTTGATAGt                                    >  1:477787/1‑93 (MQ=255)
                                                    ttCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGtt                                   >  2:5870098/1‑93 (MQ=255)
                                                    tgccTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGtt                                   <  2:3909628/91‑1 (MQ=255)
                                                     tCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTa                                  >  2:2726063/1‑93 (MQ=255)
                                                     tCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTa                                  >  1:2184692/1‑93 (MQ=255)
                                                     tCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTa                                  <  1:5718560/93‑1 (MQ=255)
                                                     tCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTa                                  >  2:3853017/1‑93 (MQ=255)
                                                     tCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTa                                  <  1:3278366/93‑1 (MQ=255)
                                                     tCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTa                                  >  1:2902277/1‑93 (MQ=255)
                                                      ccTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTAt                                 <  2:4070595/93‑1 (MQ=255)
                                                      ccTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTAt                                 >  2:3787310/1‑93 (MQ=255)
                                                       cTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTTGTGATAGTTATc                                >  1:5623926/1‑93 (MQ=255)
                                                       cTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATc                                >  1:1477661/1‑93 (MQ=255)
                                                       cTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATc                                >  1:5200071/1‑93 (MQ=255)
                                                       cTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATc                                >  1:5118564/1‑93 (MQ=255)
                                                       cTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATc                                >  2:1759677/1‑93 (MQ=255)
                                                         gCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGc                              <  2:4681820/93‑1 (MQ=255)
                                                          cTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGcc                             <  2:5673944/93‑1 (MQ=255)
                                                            tACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCa                           >  2:3600752/1‑93 (MQ=255)
                                                             aCGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAg                          <  2:2925888/93‑1 (MQ=255)
                                                             aCGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAg                          <  2:748201/93‑1 (MQ=255)
                                                             aCGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAg                          <  2:455048/93‑1 (MQ=255)
                                                             aCGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAg                          >  2:2821751/1‑93 (MQ=255)
                                                              cGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAgg                         >  2:3252574/1‑93 (MQ=255)
                                                               gCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGc                        >  1:2245177/1‑93 (MQ=255)
                                                               gCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGc                        >  2:6263732/1‑93 (MQ=255)
                                                               gCAATGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGc                        >  1:772177/1‑93 (MQ=255)
                                                                cAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGcc                       >  1:622074/1‑93 (MQ=255)
                                                                 aGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCa                      >  2:4648801/1‑93 (MQ=255)
                                                                  gTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCaa                     >  1:3943662/1‑93 (MQ=255)
                                                                   tGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAAt                    >  1:3993488/1‑93 (MQ=255)
                                                                      gATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCa                 <  1:1185556/93‑1 (MQ=255)
                                                                      gATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAAGCCAATCCa                 >  2:3553560/1‑93 (MQ=255)
                                                                       aTTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCAc                >  2:2582904/1‑93 (MQ=255)
                                                                        ttAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACg               >  2:156295/1‑93 (MQ=255)
                                                                        ttAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACg               >  2:24391/1‑93 (MQ=255)
                                                                         tAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACgg              >  2:2869978/1‑93 (MQ=255)
                                                                          aaCTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGt             <  1:3078043/93‑1 (MQ=255)
                                                                           aCTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGtt            <  2:3192721/93‑1 (MQ=255)
                                                                            cTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTc           <  2:1853261/93‑1 (MQ=255)
                                                                                         |                                                                                        
ATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCCAGCCGT  >  NC_000913/4460825‑4461002

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ 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
CACAATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCCAGCCGTAAACC  >  NC_000913/4460821‑4461007
                                                                                             |                                                                                             
CACAATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGC                                                                                               >  SRR4302533.4888039/1‑94 (MQ=60)
 ACAATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCA                                                                                              <  SRR4302533.4929332/94‑1 (MQ=60)
 ACAATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCA                                                                                              <  SRR4302533.3733483/94‑1 (MQ=60)
 ACAATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCA                                                                                              <  SRR4302533.5472178/94‑1 (MQ=60)
   AATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAA                                                                                            >  SRR4302533.4304106/1‑94 (MQ=60)
    ATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAA                                                                                           >  SRR4302533.2309003/1‑94 (MQ=60)
      ACAGCGCTGTTTGGCAGCAGATCCTGGTAATTGCCCAATGTCAGGATGGCTTCCTGATTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACA                                                                                         <  SRR4302533.3150939/94‑1 (MQ=60)
        AGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACC                                                                                       >  SRR4302533.2968574/1‑94 (MQ=60)
        AGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACC                                                                                       >  SRR4302533.4686742/1‑94 (MQ=60)
         GGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCG                                                                                      >  SRR4302533.2205327/1‑94 (MQ=60)
         GGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCG                                                                                      <  SRR4302533.2492191/94‑1 (MQ=60)
         GGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCG                                                                                      <  SRR4302533.5580661/94‑1 (MQ=60)
          GGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGA                                                                                     <  SRR4302533.3213084/94‑1 (MQ=60)
          GGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGA                                                                                     <  SRR4302533.5405689/94‑1 (MQ=60)
           GCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAG                                                                                    <  SRR4302533.2730803/94‑1 (MQ=60)
            CTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGG                                                                                   <  SRR4302533.4371382/94‑1 (MQ=60)
            CTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGG                                                                                   <  SRR4302533.3124584/94‑1 (MQ=60)
            CTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGG                                                                                   >  SRR4302533.1036095/1‑94 (MQ=60)
             TGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGA                                                                                  >  SRR4302533.1845711/1‑94 (MQ=60)
             TGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGA                                                                                  >  SRR4302533.6044242/1‑94 (MQ=60)
             TGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGA                                                                                  <  SRR4302533.641103/94‑1 (MQ=60)
             TGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGA                                                                                  >  SRR4302533.4056730/1‑94 (MQ=60)
              GTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAA                                                                                 <  SRR4302533.4870280/94‑1 (MQ=60)
              GTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAA                                                                                 >  SRR4302533.172259/1‑94 (MQ=60)
              GTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAA                                                                                 >  SRR4302533.951575/1‑94 (MQ=60)
              GTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAA                                                                                 <  SRR4302533.5457727/94‑1 (MQ=60)
              GTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGGAGGTGCAAAACACCGAGGAA                                                                                 >  SRR4302533.2111673/1‑94 (MQ=60)
               TTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAAT                                                                                <  SRR4302533.1426249/94‑1 (MQ=60)
                  GGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCAT                                                                             >  SRR4302533.5313296/1‑94 (MQ=60)
                  GGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCAT                                                                             <  SRR4302533.651894/94‑1 (MQ=60)
                  GGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCAT                                                                             >  SRR4302533.2878240/1‑94 (MQ=60)
                   GCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATC                                                                            <  SRR4302533.43111/94‑1 (MQ=60)
                    CAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATAG                                                                           >  SRR4302533.401720/1‑94 (MQ=60)
                      GCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGC                                                                         >  SRR4302533.5749987/1‑94 (MQ=60)
                        AGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCA                                                                       >  SRR4302533.2556895/1‑94 (MQ=60)
                        AGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCA                                                                       >  SRR4302533.2430234/1‑94 (MQ=60)
                         GATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAG                                                                      <  SRR4302533.4034826/94‑1 (MQ=60)
                           TCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCG                                                                    >  SRR4302533.4859415/1‑94 (MQ=60)
                           TCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCG                                                                    <  SRR4302533.1405657/94‑1 (MQ=60)
                            CCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGC                                                                   >  SRR4302533.4982393/1‑94 (MQ=60)
                            CCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGC                                                                   >  SRR4302533.1553200/1‑94 (MQ=60)
                              TGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGG                                                                 >  SRR4302533.842240/1‑94 (MQ=60)
                              TGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGG                                                                 <  SRR4302533.4906087/94‑1 (MQ=60)
                              TGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGG                                                                 <  SRR4302533.1405965/94‑1 (MQ=60)
                               GGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGC                                                                >  SRR4302533.3230459/1‑94 (MQ=60)
                                GTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCT                                                               <  SRR4302533.4220009/94‑1 (MQ=60)
                                GTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCT                                                               >  SRR4302533.1180437/1‑94 (MQ=60)
                                GTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCT                                                               <  SRR4302533.6124533/94‑1 (MQ=60)
                                 TAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTT                                                              >  SRR4302533.4187225/1‑94 (MQ=60)
                                 TAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTT                                                              <  SRR4302533.1262348/94‑1 (MQ=60)
                                 TAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTT                                                              >  SRR4302533.2743105/1‑94 (MQ=60)
                                  AATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTC                                                             <  SRR4302533.427943/94‑1 (MQ=60)
                                  AATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTC                                                             <  SRR4302533.158913/94‑1 (MQ=60)
                                  AATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTC                                                             >  SRR4302533.3870703/1‑94 (MQ=60)
                                  AATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTC                                                             >  SRR4302533.6214570/1‑94 (MQ=60)
                                  AATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTC                                                             <  SRR4302533.1158931/94‑1 (MQ=60)
                                   ATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCT                                                            <  SRR4302533.2869736/94‑1 (MQ=60)
                                   ATTGCCCCATGTAAGTATTGCTTCCTGCTGTCGTAGGGCTATTATCTGTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCT                                                            <  SRR4302533.4822859/94‑1 (MQ=36)
                                     TGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACTGAGGAATCATCGGCCAGCGCGGCTTCTAC                                                          <  SRR4302533.3549656/94‑1 (MQ=60)
                                     TGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTAC                                                          <  SRR4302533.4855603/94‑1 (MQ=60)
                                      GCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACG                                                         >  SRR4302533.966851/1‑94 (MQ=60)
                                      GCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACG                                                         >  SRR4302533.5764082/1‑94 (MQ=60)
                                      GCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACG                                                         <  SRR4302533.2929774/94‑1 (MQ=60)
                                      GCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACG                                                         <  SRR4302533.2256014/94‑1 (MQ=60)
                                       CCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGT                                                        >  SRR4302533.3992073/1‑94 (MQ=60)
                                       CCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGT                                                        >  SRR4302533.223188/1‑94 (MQ=60)
                                       CCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGT                                                        >  SRR4302533.2174652/1‑94 (MQ=60)
                                         CCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTG                                                      >  SRR4302533.4269159/1‑94 (MQ=60)
                                          CATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGA                                                     >  SRR4302533.2011500/1‑94 (MQ=60)
                                           ATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGAT                                                    <  SRR4302533.1607461/94‑1 (MQ=60)
                                           ATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGAT                                                    <  SRR4302533.5484362/94‑1 (MQ=60)
                                            TGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATT                                                   >  SRR4302533.4045261/1‑94 (MQ=60)
                                             GTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTT                                                  <  SRR4302533.5495639/94‑1 (MQ=60)
                                             GTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTT                                                  <  SRR4302533.4586293/94‑1 (MQ=60)
                                             GTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGATGATTT                                                  >  SRR4302533.4542768/1‑94 (MQ=60)
                                              TCAGGATGGCTTCCTTCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTG                                                 <  SRR4302533.1768638/94‑1 (MQ=60)
                                              TCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTG                                                 >  SRR4302533.2175734/1‑94 (MQ=60)
                                               CAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGT                                                <  SRR4302533.4928456/94‑1 (MQ=60)
                                               CAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGT                                                <  SRR4302533.3225849/94‑1 (MQ=60)
                                               CAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGT                                                <  SRR4302533.4509250/94‑1 (MQ=60)
                                                AGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCAGCCAGCGCGGCTTCTACGTTGATTTGTT                                               >  SRR4302533.101646/1‑94 (MQ=60)
                                                    TGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATA                                           <  SRR4302533.299162/94‑1 (MQ=60)
                                                    TGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGGTGATTTGTTGATA                                           >  SRR4302533.4222421/1‑94 (MQ=60)
                                                     GGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAG                                          >  SRR4302533.639995/1‑94 (MQ=60)
                                                     GGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAG                                          >  SRR4302533.581629/1‑94 (MQ=60)
                                                      GCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGT                                         >  SRR4302533.1782883/1‑94 (MQ=60)
                                                      GCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGT                                         <  SRR4302533.5282514/94‑1 (MQ=60)
                                                      GCTTCCTGCTTACGCAGTGCGATTAACTTTTGGGAGGTGCAAAACACCCAGGAATCATCGGCCAGCGCGGCTTCTACGTTGCTTTGTTGATAGT                                         >  SRR4302533.480640/1‑94 (MQ=60)
                                                       CTTGCTGCTGACGCAGTGAGATTAACTTTTGGTAGGTGCAAAACGCCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTT                                        <  SRR4302533.1173566/94‑1 (MQ=60)
                                                       CTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTT                                        <  SRR4302533.3658896/94‑1 (MQ=60)
                                                       CTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTT                                        <  SRR4302533.3808150/94‑1 (MQ=60)
                                                       CTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTT                                        <  SRR4302533.3552378/94‑1 (MQ=60)
                                                       CTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTT                                        >  SRR4302533.1554791/1‑94 (MQ=60)
                                                        TTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTA                                       >  SRR4302533.2920105/1‑94 (MQ=60)
                                                        TTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTA                                       <  SRR4302533.5238092/94‑1 (MQ=60)
                                                        TTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTA                                       <  SRR4302533.4273882/94‑1 (MQ=60)
                                                        TTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTA                                       >  SRR4302533.1535605/1‑94 (MQ=60)
                                                        TTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTA                                       >  SRR4302533.2198159/1‑94 (MQ=60)
                                                        TTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTA                                       <  SRR4302533.2933688/94‑1 (MQ=60)
                                                         TCCTGCTTACGCAGTGCGATTATCTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTAT                                      <  SRR4302533.1646002/94‑1 (MQ=60)
                                                         TCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTAT                                      <  SRR4302533.5753091/94‑1 (MQ=60)
                                                         TCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTAT                                      <  SRR4302533.2096172/94‑1 (MQ=60)
                                                         TCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTAT                                      <  SRR4302533.4207322/94‑1 (MQ=60)
                                                         TCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTAT                                      <  SRR4302533.3298609/94‑1 (MQ=60)
                                                          CCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTTGTGATAGTTATC                                     >  SRR4302533.5657919/1‑94 (MQ=60)
                                                          CCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATC                                     >  SRR4302533.5231420/1‑94 (MQ=60)
                                                          CCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATC                                     >  SRR4302533.1486590/1‑94 (MQ=60)
                                                          CCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATC                                     >  SRR4302533.5149490/1‑94 (MQ=60)
                                                             GCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATGGTTATCGCC                                  >  SRR4302533.4141606/1‑94 (MQ=60)
                                                              CTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCC                                 >  SRR4302533.2548412/1‑94 (MQ=60)
                                                              CTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCC                                 >  SRR4302533.2871150/1‑94 (MQ=60)
                                                              CTTACGCACTTCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCC                                 <  SRR4302533.2200981/94‑1 (MQ=60)
                                                                TACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAG                               <  SRR4302533.415700/94‑1 (MQ=60)
                                                                  CGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGC                             >  SRR4302533.6197644/1‑94 (MQ=60)
                                                                  CGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGC                             >  SRR4302533.2259051/1‑94 (MQ=60)
                                                                  CGCAATGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGC                             >  SRR4302533.776720/1‑94 (MQ=60)
                                                                   GCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCC                            <  SRR4302533.4921230/94‑1 (MQ=60)
                                                                   GCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCC                            >  SRR4302533.625729/1‑94 (MQ=60)
                                                                    CAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCA                           >  SRR4302533.6300933/1‑94 (MQ=60)
                                                                    CAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCA                           <  SRR4302533.2240396/94‑1 (MQ=60)
                                                                     AGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAA                          >  SRR4302533.3800046/1‑94 (MQ=60)
                                                                     AGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAA                          >  SRR4302533.4951060/1‑94 (MQ=60)
                                                                     AGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAA                          >  SRR4302533.3967718/1‑94 (MQ=60)
                                                                      GTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAAT                         >  SRR4302533.4389492/1‑94 (MQ=60)
                                                                      GTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAAT                         >  SRR4302533.4017801/1‑94 (MQ=60)
                                                                         CGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGCTGATTTGTTGATAGTTATCGCCCAGGCCAATCCA                      >  SRR4302533.3562344/1‑94 (MQ=60)
                                                                          GATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCAC                     <  SRR4302533.5721933/94‑1 (MQ=60)
                                                                          GATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCAC                     <  SRR4302533.3095926/94‑1 (MQ=60)
                                                                          GATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCAC                     <  SRR4302533.1192672/94‑1 (MQ=60)
                                                                          GATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCAC                     <  SRR4302533.4222089/94‑1 (MQ=60)
                                                                          GATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCAC                     <  SRR4302533.5579574/94‑1 (MQ=60)
                                                                           ATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACG                    >  SRR4302533.2106410/1‑94 (MQ=60)
                                                                           ATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACG                    <  SRR4302533.6051508/94‑1 (MQ=60)
                                                                            TTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGG                   >  SRR4302533.4192389/1‑94 (MQ=60)
                                                                            TTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGG                   <  SRR4302533.5810862/94‑1 (MQ=60)
                                                                             TAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGT                  <  SRR4302533.5297191/94‑1 (MQ=60)
                                                                             TAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGT                  >  SRR4302533.3305756/1‑94 (MQ=60)
                                                                              AACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTT                 <  SRR4302533.3096932/94‑1 (MQ=60)
                                                                                CTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCG               <  SRR4302533.1497305/94‑1 (MQ=60)
                                                                                CTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCG               >  SRR4302533.3821387/1‑94 (MQ=60)
                                                                                CTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCG               >  SRR4302533.1909007/1‑94 (MQ=60)
                                                                                CTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCG               >  SRR4302533.1272327/1‑94 (MQ=60)
                                                                                CTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCG               <  SRR4302533.517572/94‑1 (MQ=60)
                                                                                  TTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCC             <  SRR4302533.1978257/94‑1 (MQ=60)
                                                                                  TTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCC             <  SRR4302533.3334738/94‑1 (MQ=60)
                                                                                    TGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCCAG           <  SRR4302533.5198781/94‑1 (MQ=60)
                                                                                        AGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCCAGCCGT       >  SRR4302533.802353/1‑94 (MQ=60)
                                                                                            GCAAACCACCGAGGAATCATCGCCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCCAGCCGTAAAC   <  SRR4302533.2983537/94‑1 (MQ=60)
                                                                                            GCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCCAGCCGTAAAC   >  SRR4302533.2941055/1‑94 (MQ=60)
                                                                                             CAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCCAGCCGTAAACC  >  SRR4302533.6345496/1‑94 (MQ=60)
                                                                                             |                                                                                             
CACAATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCCAGCCGTAAACC  >  NC_000913/4460821‑4461007

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