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 2,037,676 (C)5→6 coding (177/1005 nt) b1971 → putative reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,037,6711.C93.5% 226.2 / 9.6 78coding (172/1005 nt)b1971putative reductase
Reads supporting (aligned to +/- strand):  ref base . (0/5);  new base C (42/30);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 1.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

ACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCG‑‑CCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTT  >  NC_000913/2037583‑2037758
                                                                                          |                                                                                        
aCTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGcccg‑‑cccc                                                                                     >  2:4468539/1‑93 (MQ=255)
  cgggCATCAGTGCACCTGCCCTTGCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGcccg‑‑cccccc                                                                                   <  2:1539808/92‑2 (MQ=255)
   gggCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGcccg‑‑ccccccg                                                                                  >  2:1443492/1‑91 (MQ=255)
        tCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGa                                                                              >  2:2728668/1‑93 (MQ=255)
        tCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGa                                                                              >  1:647321/1‑93 (MQ=255)
          aGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGaaa                                                                            >  2:4186843/1‑93 (MQ=255)
          aGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGaaa                                                                            >  2:3725009/1‑93 (MQ=255)
          aGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGaaa                                                                            <  2:542106/93‑1 (MQ=255)
          aGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGaaa                                                                            >  1:973113/1‑93 (MQ=255)
               aaCTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGTCCCCCGCCGGAAAAgcgc                                                                       <  2:4513210/93‑1 (MQ=255)
               aaCTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAgcgc                                                                       >  1:4071945/1‑93 (MQ=255)
                  tGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTgg                                                                    >  1:4384507/1‑93 (MQ=255)
                  tGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCG‑‑CCCCCGCCGGAAAAGCGCTGGa                                                                   <  2:1094696/93‑1 (MQ=255)
                    cACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCG‑‑CCCCCGCCGGAAAAGCGCTGGAGt                                                                 <  2:4532211/93‑1 (MQ=255)
                      ctttctttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGtt                                                                >  2:87669/1‑93 (MQ=255)
                       tttctttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTc                                                               <  1:34172/93‑1 (MQ=255)
                        ttctttGCCTCACGCTGCGCATGCCGATCTGTGTAGCTGGTTTAAAGGAAACGATCTCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCa                                                              <  2:3667867/93‑1 (MQ=255)
                        ttctttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCa                                                              <  1:592560/93‑1 (MQ=255)
                        ttctttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCa                                                              <  2:3667868/93‑1 (MQ=255)
                          ctttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGc                                                            >  1:4589845/1‑93 (MQ=255)
                            ttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCaa                                                          <  2:835653/93‑1 (MQ=255)
                            ttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCaa                                                          <  1:3612862/93‑1 (MQ=255)
                              gCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGc                                                        >  2:4395140/1‑93 (MQ=255)
                               ccTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGcc                                                       <  2:4589845/93‑1 (MQ=255)
                                cTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCt                                                      <  2:3278718/93‑1 (MQ=255)
                                 tCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTg                                                     >  2:851322/1‑93 (MQ=255)
                                  cACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTgc                                                    <  1:3861626/93‑1 (MQ=255)
                                  cACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTgc                                                    <  1:1475537/93‑1 (MQ=255)
                                   aCGCTGCGGATGCCGATCTGGTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTgcc                                                   <  2:1180172/93‑1 (MQ=255)
                                     gCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTgccgc                                                 >  1:4524610/1‑93 (MQ=255)
                                     gCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTgccgc                                                 >  1:1756676/1‑93 (MQ=255)
                                       tGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCt                                               >  2:2301025/1‑93 (MQ=255)
                                        gcgcATGCCGATCTGCTTCGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTg                                              <  2:2967308/93‑1 (MQ=255)
                                        gcgcATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTg                                              >  1:4220571/1‑93 (MQ=255)
                                          gcATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGc                                            >  1:3797226/1‑93 (MQ=255)
                                            aTGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCaa                                          <  1:1181531/93‑1 (MQ=255)
                                             tGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCaca                                         >  2:2840132/1‑91 (MQ=255)
                                                cGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATa                                      >  1:4053232/1‑93 (MQ=255)
                                                 gATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATaa                                     >  2:183103/1‑93 (MQ=255)
                                                 gATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATaa                                     >  1:3818842/1‑93 (MQ=255)
                                                 gATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATaa                                     >  2:684888/1‑93 (MQ=255)
                                                 gATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATaa                                     >  2:684889/1‑93 (MQ=255)
                                                 gATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATaa                                     >  1:4299643/1‑93 (MQ=255)
                                                  aTCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAAc                                    >  1:3195382/1‑93 (MQ=255)
                                                   tCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAAcc                                   <  2:2930920/93‑1 (MQ=255)
                                                    cTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCACGCCTGCCGCCTGGCAAAATAACCt                                  >  1:551561/1‑93 (MQ=255)
                                                    cTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCt                                  >  1:1950820/1‑93 (MQ=255)
                                                    cTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCCGGCAAAATAACCt                                  >  1:564596/1‑93 (MQ=255)
                                                     tGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTg                                 >  2:3379092/1‑93 (MQ=255)
                                                      gCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCG‑‑CCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGcc                               <  2:1985752/93‑1 (MQ=255)
                                                       cTTAGCTGTTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAACATAACCTGct                               >  2:4260114/1‑92 (MQ=255)
                                                       cTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGcc                               >  2:2009985/1‑93 (MQ=255)
                                                        ttAGCGGGGTTTAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCa                              >  2:4226478/1‑93 (MQ=255)
                                                         tAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCAc                             <  2:3597761/93‑1 (MQ=255)
                                                          aGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCAc                             <  2:2140974/93‑1 (MQ=255)
                                                          aGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACt                            <  1:2371539/93‑1 (MQ=255)
                                                           gCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTg                           <  1:1641922/93‑1 (MQ=255)
                                                            cTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGa                          <  1:4019804/93‑1 (MQ=255)
                                                             tGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGAc                         >  1:1338690/1‑93 (MQ=255)
                                                             tGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGAc                         <  2:3695191/93‑1 (MQ=255)
                                                               gTTGAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGc                       <  2:3098107/93‑1 (MQ=255)
                                                                tttAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGcc                      >  1:2197005/1‑93 (MQ=255)
                                                                 ttAAAGGGACCGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCca                     <  2:521527/93‑1 (MQ=255)
                                                                   aaaGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCcagc                   >  2:4040299/1‑93 (MQ=255)
                                                                    aaGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCcagca                  <  1:3343225/93‑1 (MQ=255)
                                                                    aaGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCcagca                  <  1:1481011/93‑1 (MQ=255)
                                                                     aGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCcagcag                 >  2:3594952/1‑93 (MQ=255)
                                                                     aGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCcagcag                 >  1:1360728/1‑93 (MQ=255)
                                                                     aGGGAACGATCGCCCGCCCG‑‑CCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGa                <  1:3235118/93‑1 (MQ=255)
                                                                      gggAGCGACCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAAGAACCTGCCACTGACGCCAGCAGa                <  1:2403442/93‑1 (MQ=255)
                                                                      gggACCGACCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGa                <  1:2117548/93‑1 (MQ=255)
                                                                        gAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATa              >  1:3516961/1‑93 (MQ=255)
                                                                          aCGAGCGCCCGCCCCC‑‑CCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATaaa            <  2:4241632/93‑1 (MQ=255)
                                                                          aCGACCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATaaa            <  2:918789/93‑1 (MQ=255)
                                                                           cGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAg           <  1:1749710/93‑1 (MQ=255)
                                                                            gATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGt          >  1:2659282/1‑93 (MQ=255)
                                                                            gATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGt          <  2:236319/93‑1 (MQ=255)
                                                                             aTCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGtt         >  2:1675549/1‑93 (MQ=255)
                                                                             aTCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGtt         >  2:3518201/1‑93 (MQ=255)
                                                                             aTCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGtt         >  2:1948826/1‑93 (MQ=255)
                                                                             aTCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGtt         <  2:1796140/93‑1 (MQ=255)
                                                                             aTCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGtt         >  2:1643800/1‑93 (MQ=255)
                                                                              gcccccgcccg‑‑cccCCGCCGGAAAAGCGCGGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTc       <  1:3905916/90‑1 (MQ=255)
                                                                                gcccgcccgcccCCCGCAGAAAAAGCGCGGGAGTTCAGCAAGCCTGCCGCCTGGCAAACTAACCTGCCACTGACGCCAGGAGATAAAGTTTcc      <  2:42507/93‑1 (MQ=255)
                                                                                ccccgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTcc      <  2:42506/92‑1 (MQ=255)
                                                                                ccccgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTcc      <  2:2171700/92‑1 (MQ=255)
                                                                                 cccgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCg     <  1:950768/93‑1 (MQ=255)
                                                                                  ccgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCgg    >  1:3934888/1‑93 (MQ=255)
                                                                                   cgcccg‑‑cccCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGtt  >  2:3893364/1‑93 (MQ=255)
                                                                                    cccccc‑‑cccccgccGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGtt  <  1:133336/91‑1 (MQ=255)
                                                                                          |                                                                                        
ACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCG‑‑CCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTT  >  NC_000913/2037583‑2037758

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACT  >  NC_000913/2037586‑2037767
                                                                                           |                                                                                           
GGGCCTCAGTGCAACTGCCCTTTCTTTTCCTCACTCTGCGCCTGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCcgc                                                                                          <  SRR4302544.1510208/94‑4 (MQ=60)
   CATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGG                                                                                        <  SRR4302544.548997/94‑1 (MQ=60)
    ATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGA                                                                                       >  SRR4302544.653714/1‑94 (MQ=60)
      CAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAA                                                                                     >  SRR4302544.982657/1‑94 (MQ=60)
           CAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGC                                                                                >  SRR4302544.4112469/1‑94 (MQ=60)
              CTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGG                                                                             >  SRR4302544.4428903/1‑94 (MQ=60)
                    TTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCA                                                                       <  SRR4302544.34509/94‑1 (MQ=60)
                     TTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAG                                                                      <  SRR4302544.598378/94‑1 (MQ=60)
                      TCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGC                                                                     >  SRR4302544.4635522/1‑94 (MQ=60)
                         TTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAG                                                                  <  SRR4302544.3648285/94‑1 (MQ=60)
                               CACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCC                                                            <  SRR4302544.1490135/94‑1 (MQ=60)
                               CACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCC                                                            <  SRR4302544.3899965/94‑1 (MQ=60)
                                 CGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGC                                                          >  SRR4302544.4569894/1‑94 (MQ=60)
                                 CGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGC                                                          >  SRR4302544.1773869/1‑94 (MQ=60)
                                    TGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTG                                                       >  SRR4302544.4263330/1‑94 (MQ=60)
                                      CGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGC                                                     >  SRR4302544.3834830/1‑94 (MQ=60)
                                         ATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAA                                                  <  SRR4302544.1193290/94‑1 (MQ=60)
                                            CCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATA                                               >  SRR4302544.4093598/1‑94 (MQ=60)
                                             CGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAA                                              >  SRR4302544.3856700/1‑94 (MQ=60)
                                             CGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAA                                              >  SRR4302544.4343517/1‑94 (MQ=60)
                                              GATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAAC                                             >  SRR4302544.3225833/1‑94 (MQ=60)
                                                TCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCACGCCTGCCGCCTGGCAAAATAACCT                                           >  SRR4302544.556960/1‑94 (MQ=60)
                                                TCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCT                                           >  SRR4302544.1969731/1‑94 (MQ=60)
                                                TCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCCGGCAAAATAACCT                                           >  SRR4302544.570119/1‑94 (MQ=60)
                                                       AGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTG                                    <  SRR4302544.2393980/94‑1 (MQ=60)
                                                        GCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGA                                   <  SRR4302544.1658106/94‑1 (MQ=60)
                                                         CTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGAC                                  >  SRR4302544.1351910/1‑94 (MQ=60)
                                                         CTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGAC                                  <  SRR4302544.4059808/94‑1 (MQ=60)
                                                            GTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCC                               >  SRR4302544.2217721/1‑94 (MQ=60)
                                                                 AAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAG                          <  SRR4302544.3375103/94‑1 (MQ=60)
                                                                 AAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAG                          >  SRR4302544.1374238/1‑94 (MQ=60)
                                                                 AAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAG                          <  SRR4302544.1495646/94‑1 (MQ=60)
                                                                  AGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGAT                        <  SRR4302544.3266086/94‑1 (MQ=60)
                                                                   GGGAGCGACCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAAGAACCTGCCACTGACGCCAGCAGAT                        <  SRR4302544.2426253/94‑1 (MQ=60)
                                                                   GGGACCGACCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGAT                        <  SRR4302544.2137885/94‑1 (MQ=60)
                                                                    GGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATA                       >  SRR4302544.3550539/1‑94 (MQ=60)
                                                                        CGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGT                   <  SRR4302544.1766836/94‑1 (MQ=60)
                                                                        CGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGT                   >  SRR4302544.2684786/1‑94 (MQ=60)
                                                                           gccCCCGCCCGCCCCCGCCGGAAAAGCGCGGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCC               <  SRR4302544.3944724/91‑1 (MQ=60)
                                                                            cccctcCCCGCCCCCGGCGAAAAAGCGCTGCAGTTCATCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCG              <  SRR4302544.3171050/88‑1 (MQ=60)
                                                                              CCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGG             >  SRR4302544.3973965/1‑94 (MQ=60)
                                                                              CCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGG             <  SRR4302544.960116/94‑1 (MQ=60)
                                                                                CCCCCCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTA          <  SRR4302544.134574/94‑1 (MQ=60)
                                                                                 ccccgcCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTA          <  SRR4302544.3294453/88‑1 (MQ=60)
                                                                                 GCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTA          >  SRR4302544.3366034/1‑94 (MQ=60)
                                                                                    cgcCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAA       >  SRR4302544.573520/4‑94 (MQ=60)
                                                                                     gcCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAAT      >  SRR4302544.2513204/3‑94 (MQ=60)
                                                                                     CCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATA     <  SRR4302544.1580251/94‑1 (MQ=60)
                                                                                       CCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAA    >  SRR4302544.4235027/1‑94 (MQ=60)
                                                                                       CCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAA    <  SRR4302544.3229296/94‑1 (MQ=60)
                                                                                       CCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAA    >  SRR4302544.1736339/1‑94 (MQ=60)
                                                                                       CCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAA    >  SRR4302544.1426254/1‑94 (MQ=60)
                                                                                         CCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACT  <  SRR4302544.2870689/94‑1 (MQ=60)
                                                                                         CCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACT  >  SRR4302544.2341867/1‑94 (MQ=60)
                                                                                           |                                                                                           
GGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACT  >  NC_000913/2037586‑2037767

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