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,123,880 T→C Q123R (CAG→CGG)  wcaH ← GDP‑mannose mannosyl hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,123,8800TC100.0% 418.7 / NA 132Q123R (CAG→CGG) wcaHGDP‑mannose mannosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (65/67);  total (65/67)

Only 100 of 141 total aligned reads displayed.

TTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAA  >  NC_000913/2123790‑2123968
                                                                                          |                                                                                        
ttCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCAtgccgc                                                                                        <  1:4255076/93‑1 (MQ=255)
ttCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCAtgccgc                                                                                        >  2:4435077/1‑93 (MQ=255)
 tctcAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATgccgct                                                                                       >  1:2933400/1‑93 (MQ=255)
 tctcAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATgccgct                                                                                       >  1:2842439/1‑93 (MQ=255)
   tcAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCa                                                                                     >  2:3012991/1‑93 (MQ=255)
    cAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCAt                                                                                    <  1:817444/93‑1 (MQ=255)
    cAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCAt                                                                                    >  2:3312757/1‑93 (MQ=255)
        gagaAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCgg                                                                                >  1:31711/1‑93 (MQ=255)
             aaTAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTa                                                                           >  2:4100602/1‑93 (MQ=255)
              aTAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTaa                                                                          >  2:880133/1‑93 (MQ=255)
               tAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGTGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAAc                                                                         >  2:4547182/1‑93 (MQ=255)
               tAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAAc                                                                         >  2:1390559/1‑93 (MQ=255)
               tAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAAc                                                                         >  2:2151242/1‑93 (MQ=255)
               tAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAAc                                                                         <  1:162119/93‑1 (MQ=255)
               tAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAAc                                                                         <  2:2797356/93‑1 (MQ=255)
                aGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACa                                                                        <  1:2726220/93‑1 (MQ=255)
                  gcgcGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGc                                                                      <  1:3544216/93‑1 (MQ=255)
                   cgcgGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGct                                                                     >  1:2381942/1‑93 (MQ=255)
                    gcgGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGctc                                                                    <  2:1165193/93‑1 (MQ=255)
                    gcgGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGctc                                                                    <  2:3012646/93‑1 (MQ=255)
                     cgGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGctct                                                                   >  1:1738212/1‑93 (MQ=255)
                       gCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttc                                                                 >  2:4572198/1‑93 (MQ=255)
                       gCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttc                                                                 >  1:4596859/1‑93 (MQ=255)
                       gCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttc                                                                 <  2:2566859/93‑1 (MQ=255)
                       gCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttc                                                                 >  1:3503014/1‑93 (MQ=255)
                       gCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttc                                                                 <  1:2798134/93‑1 (MQ=255)
                         tGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTcttctt                                                               >  1:2348558/1‑93 (MQ=255)
                         tGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTcttctt                                                               >  2:4198954/1‑93 (MQ=255)
                         tGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTcttctt                                                               >  1:4174226/1‑93 (MQ=255)
                          gTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCttcttc                                                              >  2:3847231/1‑93 (MQ=255)
                           ttAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTcc                                                             >  2:1603345/1‑93 (MQ=255)
                             aGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGa                                                           >  1:3656215/1‑93 (MQ=255)
                             aGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGa                                                           >  2:841841/1‑93 (MQ=255)
                             aGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGa                                                           >  2:1247927/1‑93 (MQ=255)
                             aGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGa                                                           <  1:2094289/93‑1 (MQ=255)
                              gCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGAt                                                          >  2:3672986/1‑93 (MQ=255)
                              gCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGCAACAGCTCTTCCTCCGAt                                                          >  2:369922/1‑93 (MQ=255)
                                aTGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAc                                                        <  2:3205869/93‑1 (MQ=255)
                                aTGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAc                                                        >  2:1048702/1‑93 (MQ=255)
                                aTGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAc                                                        <  2:140995/93‑1 (MQ=255)
                                aTGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAc                                                        <  1:553448/93‑1 (MQ=255)
                                  gAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATActc                                                      >  2:1579835/1‑93 (MQ=255)
                                     cATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGa                                                   <  2:4400498/93‑1 (MQ=255)
                                            aatGGAGAGCAGCGCGTCCGACGTCAGCCAGCGGCAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGaca                                            >  2:2959683/3‑91 (MQ=255)
                                             cTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGaaaa                                           <  2:557052/93‑1 (MQ=255)
                                             cTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGaaaa                                           >  2:4140353/1‑93 (MQ=255)
                                             cTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGaaaa                                           >  2:2754517/1‑93 (MQ=255)
                                             cTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGaaaa                                           <  2:2223320/93‑1 (MQ=255)
                                              tGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAAc                                          <  2:2862114/93‑1 (MQ=255)
                                              tGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAAc                                          <  2:2544129/93‑1 (MQ=255)
                                               ggCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAAcc                                         >  2:1426387/1‑93 (MQ=255)
                                                gCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCg                                        <  1:2747696/93‑1 (MQ=255)
                                                gCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCg                                        <  1:4524846/93‑1 (MQ=255)
                                                  gAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAg                                      <  2:4231611/93‑1 (MQ=255)
                                                  gAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAg                                      <  1:424721/93‑1 (MQ=255)
                                                  gAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAAACGAg                                      >  1:3258270/1‑93 (MQ=255)
                                                   agcagcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGc                                     >  1:1033925/1‑93 (MQ=255)
                                                   agcagcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGc                                     <  2:347463/93‑1 (MQ=255)
                                                    gcagcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGca                                    >  1:1648821/1‑93 (MQ=255)
                                                       gcgcGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGcacca                                 >  2:2238982/1‑93 (MQ=255)
                                                       gcgcGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGcacca                                 >  1:3607721/1‑93 (MQ=255)
                                                        cgcgTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGcaccac                                <  1:1494351/93‑1 (MQ=255)
                                                         gcgTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACcaca                               <  2:4181491/93‑1 (MQ=255)
                                                         gcgTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACcaca                               >  1:451873/1‑93 (MQ=255)
                                                         ccgTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACcaca                               <  2:2361541/92‑1 (MQ=255)
                                                           gTCCGACTTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATa                             <  1:200604/93‑1 (MQ=255)
                                                           gTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATa                             >  2:2253170/1‑93 (MQ=255)
                                                            tCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATag                            >  2:475834/1‑93 (MQ=255)
                                                             ccGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagt                           >  1:4037802/1‑93 (MQ=255)
                                                              cGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtg                          >  1:981566/1‑93 (MQ=255)
                                                               gACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtga                         <  2:210255/93‑1 (MQ=255)
                                                               gACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtga                         <  2:3738367/93‑1 (MQ=255)
                                                               gACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtga                         >  2:1807964/1‑93 (MQ=255)
                                                               gACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtga                         <  2:1178462/93‑1 (MQ=255)
                                                                aCGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtgag                        >  1:539604/1‑93 (MQ=255)
                                                                   tCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAgtgg                     <  2:2946971/93‑1 (MQ=255)
                                                                   tCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAgtgg                     <  2:1220808/93‑1 (MQ=255)
                                                                    cagccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAgtggt                    <  2:1091012/93‑1 (MQ=255)
                                                                    cagccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAgtggt                    <  1:2640939/93‑1 (MQ=255)
                                                                      gccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGa                  <  1:1426124/93‑1 (MQ=255)
                                                                       ccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGaa                 >  2:2277635/1‑93 (MQ=255)
                                                                        cagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGaaa                <  1:3568389/93‑1 (MQ=255)
                                                                          gcgcTAATCGTCATGCCGCTCATCCGGCAGTACCAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATc              <  2:3092503/93‑1 (MQ=255)
                                                                          gcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATc              <  1:4100602/93‑1 (MQ=255)
                                                                          gcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATc              >  2:2239594/1‑93 (MQ=255)
                                                                           cGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATcc             <  1:3848696/93‑1 (MQ=255)
                                                                           cGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATcc             <  1:3180214/93‑1 (MQ=255)
                                                                            ggTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCg            <  2:1775753/93‑1 (MQ=255)
                                                                             gTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAACCCGAGCACCACATAGTGAGTGGTGAAATCCGt           <  1:332614/93‑1 (MQ=255)
                                                                             gTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGt           >  2:683637/1‑93 (MQ=255)
                                                                             gTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGt           <  1:4435077/93‑1 (MQ=255)
                                                                                aTCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGcc        >  1:1019017/1‑93 (MQ=255)
                                                                                aTCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGcc        >  2:3831252/1‑93 (MQ=255)
                                                                                aTCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGcc        >  2:2928499/1‑93 (MQ=255)
                                                                                aTCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGcc        >  1:3827319/1‑93 (MQ=255)
                                                                                   gTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCaga     <  2:2297357/93‑1 (MQ=255)
                                                                                   gTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCaga     <  2:3580455/93‑1 (MQ=255)
                                                                                    tCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCagag    >  1:708526/1‑93 (MQ=255)
                                                                                      aTGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGaa  >  2:2589814/1‑93 (MQ=255)
                                                                                      aTGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGaa  >  1:1183145/1‑93 (MQ=255)
                                                                                          |                                                                                        
TTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAA  >  NC_000913/2123790‑2123968

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ 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
CGCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAAGTT  >  NC_000913/2123787‑2123971
                                                                                             |                                                                                           
CGCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCC                                                                                             >  SRR4302535.1910513/1‑94 (MQ=60)
CGCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCC                                                                                             <  SRR4302535.541767/94‑1 (MQ=60)
 GCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCG                                                                                            <  SRR4302535.3432728/94‑1 (MQ=60)
 GCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCG                                                                                            >  SRR4302535.1713759/1‑94 (MQ=60)
   TTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCT                                                                                          <  SRR4302535.4284227/94‑1 (MQ=60)
   TTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCT                                                                                          >  SRR4302535.2954047/1‑94 (MQ=60)
   TTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCT                                                                                          >  SRR4302535.2862596/1‑94 (MQ=60)
      TCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCAT                                                                                       <  SRR4302535.460701/94‑1 (MQ=60)
       CAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATC                                                                                      <  SRR4302535.823180/94‑1 (MQ=60)
        AGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCC                                                                                     >  SRR4302535.3055648/1‑94 (MQ=60)
          CGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGG                                                                                   >  SRR4302535.2223994/1‑94 (MQ=60)
          CGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGG                                                                                   >  SRR4302535.31954/1‑94 (MQ=60)
               AAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTA                                                                              <  SRR4302535.887344/94‑1 (MQ=60)
                  TAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACA                                                                           <  SRR4302535.163251/94‑1 (MQ=60)
                   AGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAG                                                                          >  SRR4302535.4224386/1‑94 (MQ=60)
                   AGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAG                                                                          <  SRR4302535.2745694/94‑1 (MQ=60)
                     GCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCT                                                                        >  SRR4302535.2399045/1‑94 (MQ=60)
                     GCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCT                                                                        <  SRR4302535.3568831/94‑1 (MQ=60)
                       GCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCT                                                                      >  SRR4302535.1750775/1‑94 (MQ=60)
                         GGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTC                                                                    >  SRR4302535.4628132/1‑94 (MQ=60)
                         GGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTC                                                                    >  SRR4302535.3527362/1‑94 (MQ=60)
                          GCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCT                                                                   <  SRR4302535.2818045/94‑1 (MQ=60)
                           CTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTT                                                                  >  SRR4302535.2365395/1‑94 (MQ=60)
                           CTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTT                                                                  >  SRR4302535.4202795/1‑94 (MQ=60)
                           CTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTT                                                                  <  SRR4302535.3266302/94‑1 (MQ=60)
                               TAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGA                                                              >  SRR4302535.3681524/1‑94 (MQ=60)
                                AGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGAT                                                             <  SRR4302535.2109389/94‑1 (MQ=60)
                                  CATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAC                                                           >  SRR4302535.2677946/1‑94 (MQ=60)
                                   ATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACT                                                          <  SRR4302535.557320/94‑1 (MQ=60)
                                     gaccattacccctggcggccgcccCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCT                                                        <  SRR4302535.2490597/70‑1 (MQ=60)
                                       ACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGA                                                      >  SRR4302535.4067501/1‑94 (MQ=60)
                                          TTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGC                                                   >  SRR4302535.3755540/1‑94 (MQ=60)
                                                   GCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGA                                          <  SRR4302535.2767282/94‑1 (MQ=60)
                                                   GCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGA                                          <  SRR4302535.4555845/94‑1 (MQ=60)
                                                    CGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAAACGAG                                         >  SRR4302535.3280913/1‑94 (MQ=60)
                                                     GAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGC                                        <  SRR4302535.427687/94‑1 (MQ=60)
                                                     GAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGC                                        >  SRR4302535.1041270/1‑94 (MQ=60)
                                                      AGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCA                                       >  SRR4302535.1660748/1‑94 (MQ=60)
                                                        CAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACC                                     >  SRR4302535.2773585/1‑94 (MQ=60)
                                                        CAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACC                                     >  SRR4302535.701300/1‑94 (MQ=60)
                                                         AGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCA                                    >  SRR4302535.3632741/1‑94 (MQ=60)
                                                          GCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCAC                                   <  SRR4302535.3620094/94‑1 (MQ=60)
                                                           CGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACA                                  >  SRR4302535.455038/1‑94 (MQ=60)
                                                           CGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACA                                  <  SRR4302535.1505094/94‑1 (MQ=60)
                                                              GTCCGACTTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAG                               <  SRR4302535.201976/94‑1 (MQ=60)
                                                               TCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGT                              >  SRR4302535.4065441/1‑94 (MQ=60)
                                                                CCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTG                             >  SRR4302535.988543/1‑94 (MQ=60)
                                                                  GACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAG                           >  SRR4302535.543370/1‑94 (MQ=60)
                                                                       CAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTG                      <  SRR4302535.2659858/94‑1 (MQ=60)
                                                                         GCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAA                    <  SRR4302535.1436302/94‑1 (MQ=60)
                                                                           CAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAAT                  <  SRR4302535.3593133/94‑1 (MQ=60)
                                                                            AGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATC                 >  SRR4302535.1338836/1‑94 (MQ=60)
                                                                             GCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCC                <  SRR4302535.4128679/94‑1 (MQ=60)
                                                                              CGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCG               <  SRR4302535.3875185/94‑1 (MQ=60)
                                                                              CGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCG               <  SRR4302535.3202199/94‑1 (MQ=60)
                                                                               GGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGT              <  SRR4302535.3371470/94‑1 (MQ=60)
                                                                                GTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAACCCGAGCACCACATAGTGAGTGGTGAAATCCGTG             <  SRR4302535.334914/94‑1 (MQ=60)
                                                                                GTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTG             <  SRR4302535.4465474/94‑1 (MQ=60)
                                                                                  AATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCC           >  SRR4302535.1820361/1‑94 (MQ=60)
                                                                                  AATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCC           >  SRR4302535.3853721/1‑94 (MQ=60)
                                                                                  AATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCC           >  SRR4302535.1026256/1‑94 (MQ=60)
                                                                                      GTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAG       >  SRR4302535.713499/1‑94 (MQ=60)
                                                                                        CATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAA     <  SRR4302535.4477212/94‑1 (MQ=60)
                                                                                        CATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAA     >  SRR4302535.163250/1‑94 (MQ=60)
                                                                                        CATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAA     >  SRR4302535.1191594/1‑94 (MQ=60)
                                                                                           GCCGCTCATCCGGCAGTAACAGCTTTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAAGTT  <  SRR4302535.1936507/94‑1 (MQ=60)
                                                                                           GCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAAGTT  >  SRR4302535.4603905/1‑94 (MQ=60)
                                                                                             |                                                                                           
CGCTTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAAGTT  >  NC_000913/2123787‑2123971

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