Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 1357 179.3 9500450 96.1% 9129932 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,983,890 T→C T88A (ACT→GCT)  araF ← L‑arabinose‑binding periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,983,8900TC98.2% 332.0 / ‑2.6 108T88A (ACT→GCT) araFL‑arabinose‑binding periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (53/53);  total (54/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

Only 100 of 118 total aligned reads displayed.

GTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGTACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAA  >  NC_000913/1983801‑1983979
                                                                                         |                                                                                         
gTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGcaca                                                                                        >  1:3205534/1‑93 (MQ=255)
gTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGcaca                                                                                        >  1:3205530/1‑93 (MQ=255)
gTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGcaca                                                                                        >  2:2298294/1‑93 (MQ=255)
gTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGcaca                                                                                        <  2:1402644/93‑1 (MQ=255)
 ttAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACaa                                                                                       >  2:4692767/1‑93 (MQ=255)
 ttAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACaa                                                                                       >  2:2650109/1‑93 (MQ=255)
  tAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGTACaaa                                                                                      <  2:4033835/93‑1 (MQ=255)
    aCAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATa                                                                                    <  2:3951241/93‑1 (MQ=255)
    aCAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGATTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATa                                                                                    >  2:2855989/1‑93 (MQ=255)
     cAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGCTTGGGGTCCGGAGCACAAATaa                                                                                   >  1:4425494/1‑93 (MQ=255)
       aaCTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACg                                                                                 >  2:3492267/1‑93 (MQ=255)
       aaCTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACg                                                                                 >  2:4377632/1‑93 (MQ=255)
       aaCTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACg                                                                                 >  1:4185626/1‑93 (MQ=255)
        aCTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGa                                                                                >  2:152064/1‑93 (MQ=255)
           ggTCATCCACGGCAATGACTTTCCTATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAAc                                                                             <  2:2591969/93‑1 (MQ=255)
           ggTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAAc                                                                             <  2:2141841/93‑1 (MQ=255)
             tCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCt                                                                           <  2:2978843/93‑1 (MQ=255)
             tCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCt                                                                           <  1:477784/93‑1 (MQ=255)
             tCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCt                                                                           <  2:1109441/93‑1 (MQ=255)
              cATCCGCGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCtt                                                                          <  1:2146545/93‑1 (MQ=255)
              cATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCtt                                                                          <  2:3486055/93‑1 (MQ=255)
                tCCAAGGCAGTGACGTTTATATCGTAGCCACGCGCGTACGCGACGATGGCAGAGCCGAGTTTGGGGGCCGGAGCACAAATAACGAAACCtttt                                                                        <  2:3838531/93‑1 (MQ=255)
                 ccACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTg                                                                       <  1:2280234/93‑1 (MQ=255)
                 ccACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTg                                                                       >  2:3643963/1‑93 (MQ=255)
                  cACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTgc                                                                      <  1:2316966/93‑1 (MQ=255)
                  cACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTgc                                                                      <  2:2720536/93‑1 (MQ=255)
                   aCGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTgcg                                                                     <  2:1196968/93‑1 (MQ=255)
                     ggCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGcc                                                                   <  2:1945172/93‑1 (MQ=255)
                      gCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGccc                                                                  >  1:3664619/1‑92 (MQ=255)
                       cAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCAc                                                                 <  2:4658685/93‑1 (MQ=255)
                         aTGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTg                                                               <  1:2115598/93‑1 (MQ=255)
                         aTGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTg                                                               >  2:3422356/1‑93 (MQ=255)
                          tGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTgg                                                              >  1:2806657/1‑93 (MQ=255)
                          tGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTgg                                                              >  1:504196/1‑93 (MQ=255)
                          tGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTgg                                                              >  1:7747/1‑93 (MQ=255)
                           gACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGc                                                             >  2:202569/1‑93 (MQ=255)
                           gACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGc                                                             <  1:1069391/93‑1 (MQ=255)
                            aCTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCa                                                            >  2:4728277/1‑93 (MQ=255)
                               ttCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGTAACCTTTTGCGCCACTGGCAGcc                                                         >  2:2583916/1‑93 (MQ=255)
                                  atatCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAgg                                                      >  1:1691377/1‑93 (MQ=255)
                                   aatCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGc                                                     <  2:3362525/92‑1 (MQ=255)
                                    atCGTAGGGACGCGCTTTCGAGTTGATGGCAGAGACGAGTTTGGGGTCCGGAGCACAAATAAAGCAACCTTTTGCGCCACTGGAAGCAAGGCt                                                    >  2:3067417/1‑93 (MQ=255)
                                    atCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCt                                                    >  1:3034109/1‑93 (MQ=255)
                                     tCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTg                                                   >  2:4134381/1‑93 (MQ=255)
                                     tCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTg                                                   <  1:4352849/93‑1 (MQ=255)
                                      cGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGt                                                  <  1:845114/93‑1 (MQ=255)
                                       gTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTc                                                 <  1:1580964/93‑1 (MQ=255)
                                          gCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGAt                                              <  2:4022986/93‑1 (MQ=255)
                                          gCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGAt                                              <  2:1155172/93‑1 (MQ=255)
                                            cACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATcg                                            >  2:1016163/1‑93 (MQ=255)
                                             ccgcgCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATcgc                                           <  1:1904828/92‑1 (MQ=255)
                                             aCGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGGCGATcgc                                           >  2:4307275/1‑93 (MQ=255)
                                              cgcgCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATcgcg                                          >  2:2100639/1‑93 (MQ=255)
                                                cgcTTTCGCGACGACGGCAGAGCCGAGGTTGGGGTCCGGCGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGtt                                        >  1:884595/1‑93 (MQ=255)
                                                    ttCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAAt                                    >  2:1597793/1‑93 (MQ=255)
                                                    ttCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAAt                                    <  2:3993571/93‑1 (MQ=255)
                                                      cgcgACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGt                                  <  2:2993955/93‑1 (MQ=255)
                                                       gcgACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGtt                                 <  2:79800/93‑1 (MQ=255)
                                                        cgacgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGttt                                <  1:4743991/93‑1 (MQ=255)
                                                          acgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGttttt                              <  1:4595111/93‑1 (MQ=255)
                                                          acgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGttttt                              <  2:2252333/93‑1 (MQ=255)
                                                          acgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGttttt                              <  2:2252469/93‑1 (MQ=255)
                                                           cgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGtttttt                             <  2:1145817/93‑1 (MQ=255)
                                                            gaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTc                            <  2:3312246/93‑1 (MQ=255)
                                                            gaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTc                            >  2:1805699/1‑93 (MQ=255)
                                                             aTGGCAGAGCCGAGTTTTGTGTCCGTAGCACAAATAACGAAACCTTTTGCGTCTCTTGTAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCg                           >  2:4337597/1‑93 (MQ=255)
                                                             aTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCg                           >  2:2011182/1‑93 (MQ=255)
                                                              tGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGc                          <  2:1613403/93‑1 (MQ=255)
                                                              tGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGc                          >  2:2354484/1‑93 (MQ=255)
                                                                gCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCa                        <  1:3931553/93‑1 (MQ=255)
                                                                  agagCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATc                      >  1:2468267/1‑93 (MQ=255)
                                                                  agagCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATc                      >  2:4703321/1‑93 (MQ=255)
                                                                   gagCCGAGTTTGGGGTCCGGCGCACAAATAACGAAGCCTTTTGCGCCTCGGGCAGCCAGGCTGTCGCTCGCGTTCAATGTTTTTTCGCCCTcc                     <  2:244098/93‑1 (MQ=255)
                                                                   gagCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATcc                     >  2:1971147/1‑93 (MQ=255)
                                                                     gCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcgg                   >  1:204496/1‑93 (MQ=255)
                                                                     gCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcgg                   >  1:411457/1‑93 (MQ=255)
                                                                       cGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggca                 <  2:1523397/93‑1 (MQ=255)
                                                                        gAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcgggac                >  1:443283/1‑93 (MQ=255)
                                                                        gAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcac                <  2:1656784/93‑1 (MQ=255)
                                                                        gAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcac                <  1:749748/93‑1 (MQ=255)
                                                                          gTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcacgg              >  1:2166167/1‑93 (MQ=255)
                                                                          gTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcacgg              >  2:188935/1‑93 (MQ=255)
                                                                          gTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcacgg              <  1:1449610/93‑1 (MQ=255)
                                                                           tttGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcacggc             <  2:856603/93‑1 (MQ=255)
                                                                           tttGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcacggc             >  1:769485/1‑93 (MQ=255)
                                                                             tGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCaa           <  2:104476/93‑1 (MQ=255)
                                                                             tGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCaa           <  1:4586716/93‑1 (MQ=255)
                                                                                ggTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCt        >  2:3262300/1‑93 (MQ=255)
                                                                                ggTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCt        >  1:2977106/1‑93 (MQ=255)
                                                                                ggTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCt        <  1:1143838/93‑1 (MQ=255)
                                                                                 gTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCtt       >  2:2877566/1‑93 (MQ=255)
                                                                                 gTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCtt       >  1:1434119/1‑93 (MQ=255)
                                                                                  tCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTta      >  2:4565518/1‑93 (MQ=255)
                                                                                   ccGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaa     <  2:376038/93‑1 (MQ=255)
                                                                                    cGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaat    <  1:1144038/93‑1 (MQ=255)
                                                                                    cGGAGCACAAATAACGAAACCTGTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaat    <  1:4556031/93‑1 (MQ=255)
                                                                                     cgAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaata   <  1:2712025/92‑1 (MQ=255)
                                                                                      gAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaataa  >  2:2158838/1‑93 (MQ=255)
                                                                                      gAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaataa  <  2:2073150/93‑1 (MQ=255)
                                                                                      gAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaataa  >  2:1663260/1‑93 (MQ=255)
                                                                                         |                                                                                         
GTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGTACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAA  >  NC_000913/1983801‑1983979

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCGTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGTACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCT  >  NC_000913/1983799‑1983982
                                                                                           |                                                                                            
GCGTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCAC                                                                                            <  SRR4302544.1115995/94‑1 (MQ=60)
 CGTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACA                                                                                           >  SRR4302544.3236103/1‑94 (MQ=60)
 CGTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACA                                                                                           >  SRR4302544.3236107/1‑94 (MQ=60)
      ACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGCTTGGGGTCCGGAGCACAAATAA                                                                                      >  SRR4302544.4470143/1‑94 (MQ=60)
        AAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACG                                                                                    >  SRR4302544.4227077/1‑94 (MQ=60)
               TCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTT                                                                             <  SRR4302544.2224847/94‑1 (MQ=60)
               TCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTT                                                                             <  SRR4302544.482461/94‑1 (MQ=60)
                CATCCGCGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTT                                                                            <  SRR4302544.2166998/94‑1 (MQ=60)
                 ATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTT                                                                           <  SRR4302544.4514974/94‑1 (MQ=60)
                   CCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGC                                                                         <  SRR4302544.2301719/94‑1 (MQ=60)
                    CACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCG                                                                        <  SRR4302544.2338897/94‑1 (MQ=60)
                       GGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCC                                                                     >  SRR4302544.3701163/1‑94 (MQ=60)
                           ATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGG                                                                 >  SRR4302544.7817/1‑94 (MQ=60)
                           ATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGG                                                                 >  SRR4302544.509152/1‑94 (MQ=60)
                           ATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGG                                                                 <  SRR4302544.2135921/94‑1 (MQ=60)
                           ATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGG                                                                 >  SRR4302544.2833719/1‑94 (MQ=60)
                             GACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCA                                                               <  SRR4302544.1080019/94‑1 (MQ=60)
                                   CATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGG                                                         >  SRR4302544.1708068/1‑94 (MQ=60)
                                     TATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCT                                                       >  SRR4302544.3063171/1‑94 (MQ=60)
                                      ATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGGCAGGCTG                                                      >  SRR4302544.1700667/1‑94 (MQ=60)
                                       TCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGT                                                     <  SRR4302544.4397062/94‑1 (MQ=60)
                                        CGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTC                                                    <  SRR4302544.853457/94‑1 (MQ=60)
                                         GTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCG                                                   <  SRR4302544.1596587/94‑1 (MQ=60)
                                               CCGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCG                                             <  SRR4302544.1923240/94‑1 (MQ=60)
                                                 GCGCTTTCGCGACGACGGCAGAGCCGAGGTTGGGGTCCGGCGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTT                                           >  SRR4302544.893358/1‑94 (MQ=60)
                                                          CGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTT                                  <  SRR4302544.4790741/94‑1 (MQ=60)
                                                            ACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTT                                <  SRR4302544.4640815/94‑1 (MQ=60)
                                                                 GGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTNCGCCACTGGCAGCCAGGCTGTNGNNNNCGNTCAATGTNTNTTCtnnn                           >  SRR4302544.1496168/1‑90 (MQ=60)
                                                                  GCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCAT                          <  SRR4302544.3970604/94‑1 (MQ=60)
                                                                   CAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATC                         >  SRR4302544.2491850/1‑94 (MQ=60)
                                                                     GAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCG                       <  SRR4302544.944634/94‑1 (MQ=60)
                                                                      AGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGG                      >  SRR4302544.415654/1‑94 (MQ=60)
                                                                      AGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGG                      >  SRR4302544.206580/1‑94 (MQ=60)
                                                                         CGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGGAC                   >  SRR4302544.447716/1‑94 (MQ=60)
                                                                          GAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACG                  <  SRR4302544.408019/94‑1 (MQ=60)
                                                                          GAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACG                  <  SRR4302544.757073/94‑1 (MQ=60)
                                                                           AGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGG                 >  SRR4302544.2186726/1‑94 (MQ=60)
                                                                            GTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGC                >  SRR4302544.2942444/1‑94 (MQ=60)
                                                                            GTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGC                >  SRR4302544.776976/1‑94 (MQ=60)
                                                                            GTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGC                <  SRR4302544.1464018/94‑1 (MQ=60)
                                                                               TGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAAT             <  SRR4302544.4632377/94‑1 (MQ=60)
                                                                               TGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAAT             >  SRR4302544.551099/1‑94 (MQ=60)
                                                                                 GGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCT           >  SRR4302544.3005662/1‑94 (MQ=60)
                                                                                  GGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTT          <  SRR4302544.1155212/94‑1 (MQ=60)
                                                                                  GGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTT          >  SRR4302544.1448388/1‑94 (MQ=60)
                                                                                      CGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATA      <  SRR4302544.1155414/94‑1 (MQ=60)
                                                                                      CGGAGCACAAATAACGAAACCTGTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATA      <  SRR4302544.4601496/94‑1 (MQ=60)
                                                                                          GCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCT  <  SRR4302544.4362013/94‑1 (MQ=60)
                                                                                          GCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCT  >  SRR4302544.1756338/1‑94 (MQ=60)
                                                                                           |                                                                                            
GCGTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGTACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCT  >  NC_000913/1983799‑1983982

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