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 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,106,960 T→C D216G (GAC→GGC)  rfbX ← putative O‑antigen transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,106,9600TC98.5% 407.8 / ‑2.5 133D216G (GAC→GGC) rfbXputative O‑antigen transporter
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (80/51);  total (81/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

Only 100 of 135 total aligned reads displayed.

GACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCC  >  NC_000913/2106870‑2107048
                                                                                          |                                                                                        
  cGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGc                                                                                      <  2:332441/93‑1 (MQ=255)
   gTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCt                                                                                     >  2:4068246/1‑93 (MQ=255)
   gTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCt                                                                                     >  1:4026254/1‑93 (MQ=255)
    ttGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTa                                                                                    <  2:2249749/93‑1 (MQ=255)
    ttGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTa                                                                                    >  1:1786588/1‑93 (MQ=255)
     tGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTaa                                                                                   <  2:4387400/93‑1 (MQ=255)
      ggTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATGTAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTaaa                                                                                  <  2:4619787/93‑1 (MQ=255)
      ggTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTaaa                                                                                  <  1:3089221/93‑1 (MQ=255)
      ggTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACACGAAAACCGCCTGCTaaa                                                                                  <  1:2543315/93‑1 (MQ=255)
       gTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAg                                                                                 >  1:4072032/1‑93 (MQ=255)
        tCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGa                                                                                <  2:597136/93‑1 (MQ=255)
         ccAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGAt                                                                               <  2:3231185/93‑1 (MQ=255)
          cAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATc                                                                              <  2:1495920/93‑1 (MQ=255)
           aGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCg                                                                             >  1:2707941/1‑93 (MQ=255)
            gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                            >  1:2921758/1‑93 (MQ=255)
            gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                            >  1:3095685/1‑93 (MQ=255)
            gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGGGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                            <  1:2468241/93‑1 (MQ=255)
              tataAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACg                                                                          >  2:2970725/1‑93 (MQ=255)
               ataAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGc                                                                         >  1:962598/1‑93 (MQ=255)
                taAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcg                                                                        >  1:898730/1‑92 (MQ=255)
                taAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcg                                                                        >  2:4027523/1‑92 (MQ=255)
                 aaaTCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcgc                                                                       >  1:4000095/1‑91 (MQ=255)
                 aaaTCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcgc                                                                       <  1:434801/93‑3 (MQ=255)
                  aaTCCCAGGATAATACCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcgca                                                                      <  1:2033720/93‑1 (MQ=255)
                      ccAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTAt                                                                  >  2:1751500/1‑93 (MQ=255)
                       cAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATg                                                                 >  2:3275547/1‑93 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTa                                                               <  2:1568440/93‑1 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTa                                                               >  1:4471384/1‑93 (MQ=255)
                           ataataACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAAt                                                             >  1:1958819/1‑93 (MQ=255)
                            taataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATg                                                            >  2:130509/1‑93 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGa                                                           <  2:4348244/93‑1 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGa                                                           >  2:2001059/1‑93 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGa                                                           <  1:710490/93‑1 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGa                                                           >  1:4500021/1‑93 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGa                                                           <  1:2870861/93‑1 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGa                                                           <  2:1565436/93‑1 (MQ=255)
                              ataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGAt                                                          >  1:2635852/1‑93 (MQ=255)
                               taaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAATATCGACGCGCATTATGTAATGATa                                                         >  2:4539698/1‑93 (MQ=255)
                               taaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATa                                                         <  1:185119/93‑1 (MQ=255)
                               taaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATa                                                         >  2:3552136/1‑93 (MQ=255)
                                aaCTATTCCCGTGGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAg                                                        >  1:528614/1‑93 (MQ=255)
                                  cTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCg                                                      >  2:2680368/1‑93 (MQ=255)
                                  cTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCg                                                      >  1:1655777/1‑93 (MQ=255)
                                  cTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCg                                                      >  2:2808981/1‑93 (MQ=255)
                                    aTTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTa                                                    >  2:2624619/1‑93 (MQ=255)
                                       cccGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACt                                                 >  1:3906053/1‑93 (MQ=255)
                                       cccGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACt                                                 <  2:2317930/93‑1 (MQ=255)
                                        ccGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACtt                                                >  1:588377/1‑93 (MQ=255)
                                           tAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTAc                                             >  2:3671648/1‑93 (MQ=255)
                                              aaTATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCaa                                          >  2:2826448/1‑93 (MQ=255)
                                              aaTATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCaa                                          >  2:854349/1‑93 (MQ=255)
                                               atatAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAAt                                         <  1:3480153/93‑1 (MQ=255)
                                                tataAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATc                                        >  1:2038745/1‑93 (MQ=255)
                                                 ataAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATcc                                       <  2:2637341/93‑1 (MQ=255)
                                                 ataAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATcc                                       >  2:2320193/1‑93 (MQ=255)
                                                  taAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCa                                      >  1:1901619/1‑93 (MQ=255)
                                                  taAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCa                                      <  1:1644025/93‑1 (MQ=255)
                                                    aaCTAATAGCAGATGTGGCAATAAAAACATGAAAACCGCCTGCTCAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAAc                                    >  2:707954/1‑93 (MQ=255)
                                                      cTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCt                                  <  2:4310940/93‑1 (MQ=255)
                                                      cTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCt                                  <  1:1843078/93‑1 (MQ=255)
                                                       tAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                 >  1:4309260/1‑93 (MQ=255)
                                                       tAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                 >  2:1209759/1‑93 (MQ=255)
                                                       tAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                 >  1:2251113/1‑93 (MQ=255)
                                                       tAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                 >  2:634151/1‑93 (MQ=255)
                                                       tAATAGCAGATGTGGAAATAAAAACATGAAAACCACCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                 >  2:4042853/1‑93 (MQ=255)
                                                         aTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCa                               <  1:2988887/93‑1 (MQ=255)
                                                         aTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCa                               >  2:2478997/1‑93 (MQ=255)
                                                           aGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATg                             >  2:3599834/1‑93 (MQ=255)
                                                            gCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGa                            >  1:3444889/1‑93 (MQ=255)
                                                             cAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaa                           <  2:3212479/93‑1 (MQ=255)
                                                              aGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaac                          >  2:4591200/1‑93 (MQ=255)
                                                              aGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaac                          <  2:954730/93‑1 (MQ=255)
                                                              aGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaac                          >  2:1998150/1‑93 (MQ=255)
                                                              aGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaac                          >  2:1580096/1‑93 (MQ=255)
                                                               gATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaaca                         >  1:1315150/1‑93 (MQ=255)
                                                               gATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaaca                         >  1:2660029/1‑93 (MQ=255)
                                                               gATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaaca                         >  1:4304209/1‑93 (MQ=255)
                                                                aTGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaacaa                        <  1:4551270/93‑1 (MQ=255)
                                                                aTGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaacaa                        <  1:835393/93‑1 (MQ=255)
                                                                 tgtgGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAATCTTCATGaacaac                       >  2:3028899/1‑93 (MQ=255)
                                                                    ggAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAt                    <  2:1829571/93‑1 (MQ=255)
                                                                    ggAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAt                    <  1:2330033/93‑1 (MQ=255)
                                                                    ggAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAt                    <  1:2647455/93‑1 (MQ=255)
                                                                     gAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATa                   <  1:2477070/93‑1 (MQ=255)
                                                                       aaTAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGc                 >  2:3672294/1‑93 (MQ=255)
                                                                       aaTAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGc                 >  1:308977/1‑93 (MQ=255)
                                                                         aaaaaaCATGAAACCCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTa               <  2:289811/92‑1 (MQ=255)
                                                                          aaaaaCATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAg              <  1:3223861/93‑1 (MQ=255)
                                                                          aaaaaCATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAg              >  1:2350715/1‑93 (MQ=255)
                                                                           aaaaCATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGt             >  2:3631466/1‑93 (MQ=255)
                                                                           aaaaCATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGt             >  1:2751081/1‑93 (MQ=255)
                                                                              aCATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCa          >  1:4290687/1‑93 (MQ=255)
                                                                               cATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaa         >  1:1418110/1‑93 (MQ=255)
                                                                                aTGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaat        >  1:1268232/1‑93 (MQ=255)
                                                                                 tGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaata       >  2:2435888/1‑93 (MQ=255)
                                                                                   aaaaCCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaataat     <  2:120582/93‑1 (MQ=255)
                                                                                   aaaaCCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaataat     <  2:3378531/93‑1 (MQ=255)
                                                                                     aaCCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc   <  2:3072561/93‑1 (MQ=255)
                                                                                     aaCCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc   >  2:2416528/1‑93 (MQ=255)
                                                                                      aCCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTcc  >  2:3908749/1‑93 (MQ=255)
                                                                                          |                                                                                        
GACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCC  >  NC_000913/2106870‑2107048

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGGACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCAGCAA  >  NC_000913/2106868‑2107053
                                                                                            |                                                                                             
CGGACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCC                                                                                              >  SRR4302535.3777218/1‑94 (MQ=60)
  GACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTG                                                                                            <  SRR4302535.2793872/94‑1 (MQ=60)
    CGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCT                                                                                          >  SRR4302535.4053802/1‑94 (MQ=60)
     GTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTA                                                                                         >  SRR4302535.1799514/1‑94 (MQ=60)
        GGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAG                                                                                      >  SRR4302535.4099898/1‑94 (MQ=60)
        GGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAG                                                                                      <  SRR4302535.3110736/94‑1 (MQ=60)
        GGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACACGAAAACCGCCTGCTAAAG                                                                                      <  SRR4302535.2561628/94‑1 (MQ=60)
           CCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATC                                                                                   >  SRR4302535.3057427/1‑94 (MQ=60)
            CAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCG                                                                                  <  SRR4302535.4337579/94‑1 (MQ=60)
            CAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCG                                                                                  >  SRR4302535.2727291/1‑94 (MQ=60)
             AGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGA                                                                                 >  SRR4302535.3117244/1‑94 (MQ=60)
             AGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGA                                                                                 >  SRR4302535.2942368/1‑94 (MQ=60)
              GATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGGGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGAC                                                                                <  SRR4302535.2486036/94‑1 (MQ=60)
                TATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGC                                                                              >  SRR4302535.969436/1‑94 (MQ=60)
                 ATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCG                                                                             >  SRR4302535.905099/1‑94 (MQ=60)
                  TAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGC                                                                            >  SRR4302535.4027461/1‑94 (MQ=60)
                   AAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCA                                                                           <  SRR4302535.437831/94‑1 (MQ=60)
                    AATCCCAGGATAATACCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCAT                                                                          <  SRR4302535.2048373/94‑1 (MQ=60)
                         CAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGT                                                                     <  SRR4302535.747315/94‑1 (MQ=60)
                          AGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTA                                                                    >  SRR4302535.4502067/1‑94 (MQ=60)
                            GATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAAT                                                                  >  SRR4302535.1972919/1‑94 (MQ=60)
                              TAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGA                                                                >  SRR4302535.4530927/1‑94 (MQ=60)
                               AATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGAT                                                               <  SRR4302535.1826622/94‑1 (MQ=60)
                               AATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGAT                                                               <  SRR4302535.2891201/94‑1 (MQ=60)
                               AATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGAT                                                               >  SRR4302535.1928526/1‑94 (MQ=60)
                               AATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGAT                                                               >  SRR4302535.2654737/1‑94 (MQ=60)
                               AATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGAT                                                               <  SRR4302535.715476/94‑1 (MQ=60)
                                 TAACTATTCCCGTGGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAG                                                             >  SRR4302535.532308/1‑94 (MQ=60)
                                 TAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAG                                                             <  SRR4302535.186395/94‑1 (MQ=60)
                                  AACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGC                                                            >  SRR4302535.2242320/1‑94 (MQ=60)
                                   ACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCG                                                           >  SRR4302535.1667754/1‑94 (MQ=60)
                                    CTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGT                                                          >  SRR4302535.2398194/1‑94 (MQ=60)
                                        TCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACT                                                      >  SRR4302535.3932815/1‑94 (MQ=60)
                                         CCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTT                                                     >  SRR4302535.592503/1‑94 (MQ=60)
                                                 ATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATC                                             >  SRR4302535.2053425/1‑94 (MQ=60)
                                                 ATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATC                                             <  SRR4302535.3504330/94‑1 (MQ=60)
                                                   ATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCA                                           >  SRR4302535.2607164/1‑94 (MQ=60)
                                                   ATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCA                                           >  SRR4302535.1915307/1‑94 (MQ=60)
                                                    TAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAA                                          <  SRR4302535.1655909/94‑1 (MQ=60)
                                                        CTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTT                                      <  SRR4302535.1856368/94‑1 (MQ=60)
                                                        CTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTT                                      >  SRR4302535.2267145/1‑94 (MQ=60)
                                                        CTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTT                                      >  SRR4302535.4338741/1‑94 (MQ=60)
                                                           ATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCAT                                   <  SRR4302535.3009842/94‑1 (MQ=60)
                                                             AGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGA                                 >  SRR4302535.3468843/1‑94 (MQ=60)
                                                               CAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAAC                               <  SRR4302535.1583793/94‑1 (MQ=60)
                                                                AGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACA                              >  SRR4302535.4333657/1‑94 (MQ=60)
                                                                AGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACA                              >  SRR4302535.1324546/1‑94 (MQ=60)
                                                                AGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACA                              >  SRR4302535.2679085/1‑94 (MQ=60)
                                                                AGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGAAATGATAGCGTAACTTTACCCATCCAACCTTCATGAACA                              >  SRR4302535.2679083/1‑94 (MQ=60)
                                                                  ATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAAC                            <  SRR4302535.3177666/94‑1 (MQ=60)
                                                                  ATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAAC                            <  SRR4302535.4582375/94‑1 (MQ=60)
                                                                  ATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAAC                            <  SRR4302535.841252/94‑1 (MQ=60)
                                                                  ATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAAC                            >  SRR4302535.3580800/1‑94 (MQ=60)
                                                                     TGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAT                         >  SRR4302535.409798/1‑94 (MQ=60)
                                                                      GGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATA                        <  SRR4302535.2346699/94‑1 (MQ=60)
                                                                      GGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATA                        <  SRR4302535.2666419/94‑1 (MQ=60)
                                                                       GAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAG                       <  SRR4302535.2494935/94‑1 (MQ=60)
                                                                        AAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGC                      >  SRR4302535.311120/1‑94 (MQ=60)
                                                                           TAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAG                   >  SRR4302535.2367569/1‑94 (MQ=60)
                                                                            AAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGT                  <  SRR4302535.3246210/94‑1 (MQ=60)
                                                                            AAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGT                  >  SRR4302535.2770692/1‑94 (MQ=60)
                                                                               AACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCA               >  SRR4302535.4320034/1‑94 (MQ=60)
                                                                                ACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAA              >  SRR4302535.4443843/1‑94 (MQ=60)
                                                                                ACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAA              >  SRR4302535.1428234/1‑94 (MQ=60)
                                                                                 CATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAAT             >  SRR4302535.1277303/1‑94 (MQ=60)
                                                                                       AACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCC       >  SRR4302535.4587758/1‑94 (MQ=60)
                                                                                           GCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCAGCA   <  SRR4302535.4395903/94‑1 (MQ=60)
                                                                                            CCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCAGCAA  <  SRR4302535.4648442/94‑1 (MQ=60)
                                                                                            |                                                                                             
CGGACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCAGCAA  >  NC_000913/2106868‑2107053

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