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 3,934,461 A→G P186P (CCA→CCG rbsB → D‑ribose periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,934,4610AG99.6% 767.9 / ‑4.6 237P186P (CCA→CCGrbsBD‑ribose periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (121/115);  minor base T (0/1);  total (121/116)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 270 total aligned reads displayed.

CATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCAGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGTTCAGGCTGTATTCGCGCAGAATGA  >  NC_000913/3934371‑3934550
                                                                                          |                                                                                         
cATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCcagccagccggc                                                                                         <  1:2613394/93‑1 (MQ=255)
 aTTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCcggca                                                                                        <  1:2567611/93‑1 (MQ=255)
 aTTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCcggca                                                                                        <  1:275895/93‑1 (MQ=255)
 aTTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCcggca                                                                                        <  2:1100749/93‑1 (MQ=255)
  ttGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCcggcag                                                                                       >  2:874370/1‑93 (MQ=255)
      tGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAttt                                                                                   >  1:4704072/1‑93 (MQ=255)
       ggTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAtttt                                                                                  >  2:4419695/1‑93 (MQ=255)
       ggTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAtttt                                                                                  <  2:5496384/93‑1 (MQ=255)
         tACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGa                                                                                <  2:5396015/93‑1 (MQ=255)
         tACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGa                                                                                <  1:5779198/93‑1 (MQ=255)
          aCATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGAt                                                                               >  2:1761000/1‑93 (MQ=255)
          aCATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGAt                                                                               <  1:1066009/93‑1 (MQ=255)
           cATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATc                                                                              >  2:6041711/1‑93 (MQ=255)
             tCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGc                                                                            >  2:934514/1‑93 (MQ=255)
             tCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGc                                                                            >  2:394809/1‑93 (MQ=255)
             tCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGc                                                                            >  2:5984065/1‑93 (MQ=255)
                gCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCAtt                                                                         >  1:905329/1‑93 (MQ=255)
                 cAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTCTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTa                                                                        >  1:354408/1‑93 (MQ=255)
                  aGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTaa                                                                       <  1:2197073/93‑1 (MQ=255)
                  aGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTaa                                                                       >  1:4049223/1‑93 (MQ=255)
                    cccGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAg                                                                     <  2:5718892/93‑1 (MQ=255)
                    cccGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAg                                                                     >  1:4447769/1‑93 (MQ=255)
                      cGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGt                                                                   <  1:768834/93‑1 (MQ=255)
                      cGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGt                                                                   <  1:3047824/93‑1 (MQ=255)
                       gTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGtt                                                                  >  1:298795/1‑93 (MQ=255)
                       gTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGtt                                                                  <  1:5463884/93‑1 (MQ=255)
                       gTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGtt                                                                  >  1:2899926/1‑93 (MQ=255)
                        tGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGttt                                                                 <  2:2621301/93‑1 (MQ=255)
                        tGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGttt                                                                 >  2:4884453/1‑93 (MQ=255)
                           aCGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGaa                                                              <  2:2743035/93‑1 (MQ=255)
                           aCGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGaa                                                              >  1:3853232/1‑93 (MQ=255)
                             gTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACg                                                            <  1:240532/93‑1 (MQ=255)
                             gTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACg                                                            >  1:4309016/1‑93 (MQ=255)
                             gTGGCGAAGGCTTCAAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACg                                                            >  2:3848630/1‑93 (MQ=255)
                              tGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAATTTGATCGCATTAAAGGTTTGAACGt                                                           >  1:3070318/1‑93 (MQ=255)
                               ggCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTa                                                          >  1:2907197/1‑93 (MQ=255)
                               ggCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTa                                                          <  1:2030590/93‑1 (MQ=255)
                                gCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTaa                                                         <  2:2377543/93‑1 (MQ=255)
                                 cGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAAt                                                        <  1:5133875/93‑1 (MQ=255)
                                  gAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATg                                                       >  2:1447606/1‑93 (MQ=255)
                                   aaGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGc                                                      >  1:5977171/1‑93 (MQ=255)
                                   aaGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGc                                                      <  1:1294398/93‑1 (MQ=255)
                                    aGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCa                                                     <  2:3715227/93‑1 (MQ=255)
                                     ggCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAg                                                    >  2:479216/1‑93 (MQ=255)
                                     ggCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAg                                                    <  1:5598689/93‑1 (MQ=255)
                                      gCTTCCAGCAGGCCTTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGa                                                   <  2:2008397/93‑1 (MQ=255)
                                      gCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGa                                                   <  1:5083137/93‑1 (MQ=255)
                                        ttCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAAc                                                 <  2:860129/93‑1 (MQ=255)
                                         tCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAAcc                                                <  2:3656663/93‑1 (MQ=255)
                                         tCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAAcc                                                <  1:2789108/93‑1 (MQ=255)
                                          ccAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGCTTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAAACt                                               >  1:3606651/1‑93 (MQ=255)
                                          ccAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGCACCt                                               >  2:4620363/1‑93 (MQ=255)
                                          ccAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCt                                               >  1:3363279/1‑93 (MQ=255)
                                          ccAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCt                                               <  2:1253094/93‑1 (MQ=255)
                                           cagcagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTg                                              >  1:5555603/1‑93 (MQ=255)
                                           cagcagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTg                                              >  1:625470/1‑93 (MQ=255)
                                           cagcagGCCGTTGCTGCTCACAAGTTTAACGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTg                                              <  2:6187191/93‑1 (MQ=255)
                                            agcagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGt                                             <  2:3182474/93‑1 (MQ=255)
                                               agGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGa                                          >  1:1928461/1‑93 (MQ=255)
                                               agGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGa                                          >  2:3004152/1‑93 (MQ=255)
                                               agGCCGTTGCTGCTCACAAGTCTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGa                                          <  2:354408/93‑1 (MQ=255)
                                                ggCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGAc                                         <  1:1549214/93‑1 (MQ=255)
                                                ggCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGAc                                         <  2:3133544/93‑1 (MQ=255)
                                                 gctgTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGATCCTGTTGAcc                                        <  2:3244514/90‑1 (MQ=255)
                                                   cGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGATTGTACGTAATGCAGAACCAGTTGAACGc                                      >  2:4308655/1‑93 (MQ=255)
                                                        ctgctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATc                                 <  2:1336919/93‑1 (MQ=255)
                                                        ctgctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATc                                 <  1:119206/93‑1 (MQ=255)
                                                        cggctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATc                                 <  2:21797/91‑1 (MQ=255)
                                                         tgctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATcc                                <  2:6078593/93‑1 (MQ=255)
                                                         tgctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATcc                                <  2:503188/93‑1 (MQ=255)
                                                         tgctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATcc                                >  1:2934464/1‑93 (MQ=255)
                                                         tgctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATcc                                <  2:1645032/93‑1 (MQ=255)
                                                          gctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCg                               >  2:1638129/1‑93 (MQ=255)
                                                           ctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCgg                              >  1:493810/1‑93 (MQ=255)
                                                           ctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCgg                              >  1:5849675/1‑93 (MQ=255)
                                                            tCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGa                             <  2:5558233/93‑1 (MQ=255)
                                                             cacaAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGAt                            >  2:5269230/1‑93 (MQ=255)
                                                              acaAGTTTAATGTTCTTGCCAGCCAGCAGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATg                           <  2:4514717/93‑1 (MQ=255)
                                                                  gTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCa                       <  1:3227318/93‑1 (MQ=255)
                                                                  gTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCa                       <  1:1601242/93‑1 (MQ=255)
                                                                   tttAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAg                      >  1:1483069/1‑93 (MQ=255)
                                                                    ttAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAgg                     >  1:4113135/1‑93 (MQ=255)
                                                                    ttAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAgg                     >  2:47537/1‑93 (MQ=255)
                                                                     tAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCGTTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGc                    >  2:2543262/1‑93 (MQ=255)
                                                                      aaTGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCt                   <  2:1495007/93‑1 (MQ=255)
                                                                      aaTGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCt                   >  2:1586273/1‑93 (MQ=255)
                                                                       aTGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTg                  <  2:5650831/93‑1 (MQ=255)
                                                                        tGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGt                 >  1:666824/1‑93 (MQ=255)
                                                                        tGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGt                 >  2:1862796/1‑93 (MQ=255)
                                                                             ttGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTcg            >  1:3451422/1‑93 (MQ=255)
                                                                              tGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTcgc           >  2:4185380/1‑93 (MQ=255)
                                                                               gccagccagccgGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTcgcg          >  2:4023686/1‑93 (MQ=255)
                                                                                  agccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAg       >  2:5781422/1‑93 (MQ=255)
                                                                                  agccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAg       >  2:2255074/1‑93 (MQ=255)
                                                                                   gccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGa      >  2:2835893/1‑93 (MQ=255)
                                                                                   gccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGa      >  1:6078646/1‑93 (MQ=255)
                                                                                    ccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGaa     >  2:727058/1‑93 (MQ=255)
                                                                                      agccggcAGATTTTGATCGCATTAAAGGCTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatg   >  1:5936097/1‑93 (MQ=255)
                                                                                       gccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatga  <  2:4116362/93‑1 (MQ=255)
                                                                                       gccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatga  <  1:1354045/93‑1 (MQ=255)
                                                                                          |                                                                                         
CATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCAGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGTTCAGGCTGTATTCGCGCAGAATGA  >  NC_000913/3934371‑3934550

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
AGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCAGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGTTCAGGCTGTATTCGCGCAGAATGATGAA  >  NC_000913/3934368‑3934554
                                                                                             |                                                                                             
AGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCG                                                                                               <  SRR4302533.5179469/94‑1 (MQ=60)
 GGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGG                                                                                              <  SRR4302533.6242151/94‑1 (MQ=60)
   CATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCA                                                                                            <  SRR4302533.2629435/94‑1 (MQ=60)
    ATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAG                                                                                           <  SRR4302533.2583413/94‑1 (MQ=60)
    ATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAG                                                                                           <  SRR4302533.277574/94‑1 (MQ=60)
        CTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTT                                                                                       >  SRR4302533.4732586/1‑94 (MQ=60)
        CTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTT                                                                                       >  SRR4302533.2616685/1‑94 (MQ=60)
          GGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTG                                                                                     <  SRR4302533.3615091/94‑1 (MQ=60)
            TACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGAT                                                                                   <  SRR4302533.5814103/94‑1 (MQ=60)
            TACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGAT                                                                                   <  SRR4302533.3509480/94‑1 (MQ=60)
            TACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGAT                                                                                   >  SRR4302533.3105817/1‑94 (MQ=60)
             TCATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATC                                                                                  <  SRR4302533.2292368/94‑1 (MQ=60)
             ACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATC                                                                                  >  SRR4302533.1975779/1‑94 (MQ=60)
             ACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATC                                                                                  <  SRR4302533.1072271/94‑1 (MQ=60)
                TCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCA                                                                               <  SRR4302533.3427373/94‑1 (MQ=60)
                  CGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATT                                                                             >  SRR4302533.910635/1‑94 (MQ=60)
                   GCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTCTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTA                                                                            >  SRR4302533.356568/1‑94 (MQ=60)
                    CAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAA                                                                           >  SRR4302533.4073885/1‑94 (MQ=60)
                     AGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAA                                                                          >  SRR4302533.3089073/1‑94 (MQ=60)
                     AGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAA                                                                          >  SRR4302533.1597632/1‑94 (MQ=60)
                     AGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAA                                                                          <  SRR4302533.2210609/94‑1 (MQ=60)
                      GCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAG                                                                         >  SRR4302533.4474667/1‑94 (MQ=60)
                      GCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAG                                                                         <  SRR4302533.613494/94‑1 (MQ=60)
                       CCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGG                                                                        <  SRR4302533.154008/94‑1 (MQ=60)
                        CCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGT                                                                       >  SRR4302533.2155982/1‑94 (MQ=60)
                        CCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGT                                                                       >  SRR4302533.2922555/1‑94 (MQ=60)
                         CGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTT                                                                      >  SRR4302533.300622/1‑94 (MQ=60)
                         CGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTT                                                                      >  SRR4302533.2917745/1‑94 (MQ=60)
                         CGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTT                                                                      <  SRR4302533.3066483/94‑1 (MQ=60)
                         CGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTT                                                                      <  SRR4302533.773365/94‑1 (MQ=60)
                          GTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTT                                                                     <  SRR4302533.2138876/94‑1 (MQ=60)
                          GTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTT                                                                     >  SRR4302533.2449327/1‑94 (MQ=60)
                          GTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTT                                                                     <  SRR4302533.5496810/94‑1 (MQ=60)
                           TGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTG                                                                    <  SRR4302533.5840970/94‑1 (MQ=60)
                             AACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAA                                                                  >  SRR4302533.2986554/1‑94 (MQ=60)
                             AACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAA                                                                  >  SRR4302533.3876770/1‑94 (MQ=60)
                             AACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAA                                                                  <  SRR4302533.1316478/94‑1 (MQ=60)
                               CGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACG                                                                >  SRR4302533.4334953/1‑94 (MQ=60)
                                GTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGT                                                               >  SRR4302533.3089169/1‑94 (MQ=60)
                                GTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGT                                                               <  SRR4302533.242025/94‑1 (MQ=60)
                                GTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGT                                                               >  SRR4302533.5536685/1‑94 (MQ=60)
                                GTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGT                                                               <  SRR4302533.197225/94‑1 (MQ=60)
                                GTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAATTTGATCGCATTAAAGGTTTGAACGT                                                               >  SRR4302533.3089170/1‑94 (MQ=60)
                                 TGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTA                                                              >  SRR4302533.2925046/1‑94 (MQ=60)
                                  GGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAA                                                             <  SRR4302533.2043123/94‑1 (MQ=60)
                                    CGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATG                                                           <  SRR4302533.5164889/94‑1 (MQ=60)
                                    CGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATG                                                           <  SRR4302533.6090861/94‑1 (MQ=60)
                                     GAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGC                                                          <  SRR4302533.4414592/94‑1 (MQ=60)
                                     GAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGC                                                          >  SRR4302533.6013375/1‑94 (MQ=60)
                                      AAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCA                                                         >  SRR4302533.744993/1‑94 (MQ=60)
                                      AAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCA                                                         <  SRR4302533.1302191/94‑1 (MQ=60)
                                        GGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGA                                                       <  SRR4302533.5632514/94‑1 (MQ=60)
                                         GCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAA                                                      >  SRR4302533.4658891/1‑94 (MQ=60)
                                         GCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAA                                                      <  SRR4302533.5113858/94‑1 (MQ=60)
                                          gttccaccagtccggTGCTGCTCACAAGTATAATGTTCTTGCCCGCGAGCCGGCAGATTTTGATCACATTAAAGGTTTGAACGTAATGCAGAAC                                                     <  SRR4302533.2007923/79‑1 (MQ=60)
                                          CTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAAC                                                     <  SRR4302533.3004019/94‑1 (MQ=60)
                                          CTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAAC                                                     >  SRR4302533.283268/1‑94 (MQ=60)
                                            TCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGCTTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAAACT                                                   >  SRR4302533.3628855/1‑94 (MQ=60)
                                            TCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAGCCT                                                   >  SRR4302533.5171090/1‑94 (MQ=60)
                                            TCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCT                                                   <  SRR4302533.2806254/94‑1 (MQ=60)
                                            TCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCT                                                   >  SRR4302533.3384060/1‑94 (MQ=60)
                                             CCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTG                                                  >  SRR4302533.5589167/1‑94 (MQ=60)
                                             CCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTG                                                  <  SRR4302533.4128143/94‑1 (MQ=60)
                                             CCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTG                                                  >  SRR4302533.629142/1‑94 (MQ=60)
                                              CAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGT                                                 >  SRR4302533.5967181/1‑94 (MQ=60)
                                               AGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTT                                                >  SRR4302533.5684063/1‑94 (MQ=60)
                                                 CAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGA                                              >  SRR4302533.1940228/1‑94 (MQ=60)
                                                 CAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGA                                              >  SRR4302533.5991754/1‑94 (MQ=60)
                                                  AGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGAC                                             <  SRR4302533.3373633/94‑1 (MQ=60)
                                                   GGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACC                                            >  SRR4302533.3583474/1‑94 (MQ=60)
                                                   GGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACC                                            >  SRR4302533.1478734/1‑94 (MQ=60)
                                                   GGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACC                                            <  SRR4302533.1558552/94‑1 (MQ=60)
                                                      CGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCT                                         >  SRR4302533.5517101/1‑94 (MQ=60)
                                                      CGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCT                                         >  SRR4302533.5930607/1‑94 (MQ=60)
                                                       GTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTC                                        >  SRR4302533.1311428/1‑94 (MQ=60)
                                                       GTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTC                                        >  SRR4302533.230513/1‑94 (MQ=60)
                                                        TTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCA                                       >  SRR4302533.3805818/1‑94 (MQ=60)
                                                        TTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCA                                       <  SRR4302533.1579666/94‑1 (MQ=60)
                                                        TTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCA                                       >  SRR4302533.476347/1‑94 (MQ=60)
                                                           CTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCC                                    <  SRR4302533.2676078/94‑1 (MQ=60)
                                                           CTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCC                                    >  SRR4302533.2929614/1‑94 (MQ=60)
                                                           CTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCC                                    >  SRR4302533.2952406/1‑94 (MQ=60)
                                                           CTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCC                                    >  SRR4302533.3014502/1‑94 (MQ=60)
                                                           CTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCC                                    <  SRR4302533.119962/94‑1 (MQ=60)
                                                            TGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCG                                   <  SRR4302533.5679977/94‑1 (MQ=60)
                                                            TGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCG                                   >  SRR4302533.1255430/1‑94 (MQ=60)
                                                            TGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCG                                   >  SRR4302533.3444764/1‑94 (MQ=60)
                                                             GCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGG                                  <  SRR4302533.3350504/94‑1 (MQ=60)
                                                             GCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGG                                  >  SRR4302533.496756/1‑94 (MQ=60)
                                                             GCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGG                                  >  SRR4302533.5885042/1‑94 (MQ=60)
                                                              CTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGA                                 >  SRR4302533.63576/1‑94 (MQ=60)
                                                               TCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGAT                                <  SRR4302533.3830470/94‑1 (MQ=60)
                                                               TCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGAT                                <  SRR4302533.560218/94‑1 (MQ=60)
                                                                CACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATG                               <  SRR4302533.3062107/94‑1 (MQ=60)
                                                                CACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATG                               <  SRR4302533.2596462/94‑1 (MQ=60)
                                                                  CAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCT                             >  SRR4302533.2921827/1‑94 (MQ=60)
                                                                  CAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCT                             <  SRR4302533.743229/94‑1 (MQ=60)
                                                                  CAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCT                             >  SRR4302533.5377689/1‑94 (MQ=60)
                                                                    AGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCA                           >  SRR4302533.409594/1‑94 (MQ=60)
                                                                     GTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAG                          <  SRR4302533.3247215/94‑1 (MQ=60)
                                                                     GTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAG                          >  SRR4302533.1492023/1‑94 (MQ=60)
                                                                     GTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAG                          <  SRR4302533.1610934/94‑1 (MQ=60)
                                                                      TTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGG                         >  SRR4302533.4138131/1‑94 (MQ=60)
                                                                        TAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCT                       <  SRR4302533.782407/94‑1 (MQ=60)
                                                                        TAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCT                       >  SRR4302533.2801312/1‑94 (MQ=60)
                                                                        TAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCT                       <  SRR4302533.1699361/94‑1 (MQ=60)
                                                                         AATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTG                      <  SRR4302533.1192472/94‑1 (MQ=60)
                                                                          ATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGT                     >  SRR4302533.670763/1‑94 (MQ=60)
                                                                          ATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGT                     <  SRR4302533.1243110/94‑1 (MQ=60)
                                                                            GTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTAT                   <  SRR4302533.6031807/94‑1 (MQ=60)
                                                                             TTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATT                  <  SRR4302533.784883/94‑1 (MQ=60)
                                                                             TTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATT                  <  SRR4302533.1499852/94‑1 (MQ=60)
                                                                               CTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCG                >  SRR4302533.3472732/1‑94 (MQ=60)
                                                                               CTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCG                <  SRR4302533.3975469/94‑1 (MQ=60)
                                                                                TTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGC               >  SRR4302533.4231956/1‑94 (MQ=60)
                                                                                   CCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCA            >  SRR4302533.5883605/1‑94 (MQ=60)
                                                                                   CCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCA            >  SRR4302533.4270727/1‑94 (MQ=60)
                                                                                     AGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGA          >  SRR4302533.4868169/1‑94 (MQ=60)
                                                                                     AGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGA          >  SRR4302533.6115409/1‑94 (MQ=60)
                                                                                      GCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAA         <  SRR4302533.6192466/94‑1 (MQ=60)
                                                                                       CCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAAT        >  SRR4302533.515727/1‑94 (MQ=60)
                                                                                       CCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAAT        >  SRR4302533.2713353/1‑94 (MQ=60)
                                                                                       CCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAAT        <  SRR4302533.2842071/94‑1 (MQ=60)
                                                                                       CCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAAT        <  SRR4302533.4996797/94‑1 (MQ=60)
                                                                                       CCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAAT        >  SRR4302533.5072211/1‑94 (MQ=60)
                                                                                        CAGCCGGCAGATTTTGATCGCATTAAAGGCTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAATG       >  SRR4302533.5971989/1‑94 (MQ=60)
                                                                                         AGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAATGA      >  SRR4302533.4356247/1‑94 (MQ=60)
                                                                                          GCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGACTGAT     >  SRR4302533.2084656/1‑94 (MQ=60)
                                                                                          GCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAATGAT     <  SRR4302533.1362184/94‑1 (MQ=60)
                                                                                          GCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAATGAT     >  SRR4302533.3019763/1‑94 (MQ=60)
                                                                                           CCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAATGATG    >  SRR4302533.3911488/1‑94 (MQ=60)
                                                                                            CGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAATGATGA   <  SRR4302533.6008160/94‑1 (MQ=60)
                                                                                            CGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAATGATGA   >  SRR4302533.5673801/1‑94 (MQ=60)
                                                                                            CGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAATGATGA   >  SRR4302533.755825/1‑94 (MQ=60)
                                                                                            CGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAATGATGA   >  SRR4302533.5542887/1‑94 (MQ=60)
                                                                                             GGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAATGATGAA  <  SRR4302533.4258434/94‑1 (MQ=60)
                                                                                             GGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAATGATGAA  >  SRR4302533.1500799/1‑94 (MQ=60)
                                                                                             GGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAATGATGAA  <  SRR4302533.2164931/94‑1 (MQ=60)
                                                                                             |                                                                                             
AGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCAGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGTTCAGGCTGTATTCGCGCAGAATGATGAA  >  NC_000913/3934368‑3934554

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