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,8800TC99.5% 670.3 / ‑4.8 214Q123R (CAG→CGG) wcaHGDP‑mannose mannosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (79/134);  total (79/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

Only 100 of 223 total aligned reads displayed.

TTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAAG  >  NC_000913/2123790‑2123969
                                                                                          |                                                                                         
ttCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCAtgccgc                                                                                         <  2:3152588/93‑1 (MQ=255)
 tctcAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATAGTCATgccgct                                                                                        >  2:2451227/1‑93 (MQ=255)
   tcAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCATCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCa                                                                                      >  1:2530412/1‑93 (MQ=255)
   tcAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCa                                                                                      >  1:7700850/1‑93 (MQ=255)
   tcAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCa                                                                                      <  2:2337789/93‑1 (MQ=255)
     aGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATc                                                                                    <  1:4002383/93‑1 (MQ=255)
      gCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATcc                                                                                   <  2:5152644/93‑1 (MQ=255)
       cGAGAAAATAGGCGCGGTTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCg                                                                                  <  1:2085286/93‑1 (MQ=255)
       cGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCg                                                                                  <  1:4261329/93‑1 (MQ=255)
         agaAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGc                                                                                >  1:803288/1‑93 (MQ=255)
             aaTAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTa                                                                            >  1:2499955/1‑93 (MQ=255)
                 ggCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAg                                                                        >  1:5077287/1‑93 (MQ=255)
                  gcgcGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGc                                                                       <  2:1893184/93‑1 (MQ=255)
                  gcgcGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGc                                                                       >  1:2347838/1‑93 (MQ=255)
                     cgGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGctct                                                                    >  2:4887552/1‑93 (MQ=255)
                      ggCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtctt                                                                   <  2:500904/93‑1 (MQ=255)
                      ggCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtctt                                                                   <  1:1331502/93‑1 (MQ=255)
                       gCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttc                                                                  <  2:1340858/93‑1 (MQ=255)
                       gCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttc                                                                  <  1:4842524/93‑1 (MQ=255)
                       gCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttc                                                                  <  2:4175389/93‑1 (MQ=255)
                        cTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttct                                                                 >  1:6827981/1‑93 (MQ=255)
                         tGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTcttctt                                                                <  1:3217097/93‑1 (MQ=255)
                         tGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTcttctt                                                                >  2:253326/1‑93 (MQ=255)
                         tGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTcttctt                                                                >  2:7163222/1‑93 (MQ=255)
                          gTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCttcttc                                                               <  2:7578736/93‑1 (MQ=255)
                           ttAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTcc                                                              <  2:7219132/93‑1 (MQ=255)
                           ttAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTcc                                                              <  1:7402537/93‑1 (MQ=255)
                           ttAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTcc                                                              >  2:5366903/1‑93 (MQ=255)
                           ttAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTcc                                                              <  1:2966437/93‑1 (MQ=255)
                             aGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGa                                                            <  1:423040/93‑1 (MQ=255)
                             aGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGa                                                            <  2:2223236/93‑1 (MQ=255)
                             aGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGa                                                            <  2:1872138/93‑1 (MQ=255)
                             aGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGa                                                            <  2:7324887/93‑1 (MQ=255)
                               cATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATa                                                          <  2:4818538/93‑1 (MQ=255)
                                aTGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAc                                                         <  1:101558/93‑1 (MQ=255)
                                aTGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAc                                                         >  1:5952928/1‑93 (MQ=255)
                                 tGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAct                                                        >  1:1878426/1‑93 (MQ=255)
                                 tGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAct                                                        >  1:3139742/1‑93 (MQ=255)
                                 tGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAct                                                        <  1:5005163/93‑1 (MQ=255)
                                  gAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATActc                                                       >  1:6494860/1‑93 (MQ=255)
                                   aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATActct                                                      <  1:5121877/93‑1 (MQ=255)
                                   aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATActct                                                      >  1:1549574/1‑93 (MQ=255)
                                     cATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGa                                                    <  1:7802913/93‑1 (MQ=255)
                                     cATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGa                                                    <  1:2870712/93‑1 (MQ=255)
                                      aTTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGaa                                                   <  2:1414229/93‑1 (MQ=255)
                                      aTTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGaa                                                   >  1:3076877/1‑93 (MQ=255)
                                        tATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGc                                                 <  1:175258/93‑1 (MQ=255)
                                         ttCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCg                                                <  2:2233490/92‑1 (MQ=255)
                                            ccTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGaaa                                             <  1:3902776/92‑1 (MQ=255)
                                             cgggCGAGCACCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGaaaa                                            <  2:1566543/91‑1 (MQ=255)
                                             cTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGaaaa                                            <  1:3437943/93‑1 (MQ=255)
                                                gCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCg                                         >  1:3869442/1‑93 (MQ=255)
                                                 cGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGa                                        <  2:7666003/93‑1 (MQ=255)
                                                   agcagcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGc                                      >  1:6191051/1‑93 (MQ=255)
                                                    gcagcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGca                                     <  1:3358413/93‑1 (MQ=255)
                                                     cagcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGcac                                    <  1:6022621/93‑1 (MQ=255)
                                                      agcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGcacc                                   <  2:3389124/93‑1 (MQ=255)
                                                           gTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATa                              <  1:5482245/93‑1 (MQ=255)
                                                           gTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATa                              >  1:1849943/1‑93 (MQ=255)
                                                            tCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACCGCTCTTTTTCCGATACTCTGAACCAACAACAGAGCCCAACATag                             >  1:3727201/1‑93 (MQ=255)
                                                            tCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATag                             <  2:7006130/93‑1 (MQ=255)
                                                            tCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATag                             >  2:3959362/1‑93 (MQ=255)
                                                              cGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtg                           <  2:601328/93‑1 (MQ=255)
                                                              cGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtg                           <  2:7725099/93‑1 (MQ=255)
                                                                 cGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtgagt                        <  2:711557/93‑1 (MQ=255)
                                                                 cGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtgagg                        >  1:4476095/1‑92 (MQ=255)
                                                                  gTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtgagtg                       <  1:6605256/93‑1 (MQ=255)
                                                                  gTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtgagtg                       >  1:6421169/1‑93 (MQ=255)
                                                                  ctCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtgagtg                       <  1:2589517/92‑1 (MQ=255)
                                                                   tCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAgtgg                      >  2:2967188/1‑93 (MQ=255)
                                                                     agccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAgtggtg                    <  2:822094/93‑1 (MQ=255)
                                                                     agccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAgtggtg                    <  2:1701048/93‑1 (MQ=255)
                                                                      gccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGa                   >  1:1462223/1‑93 (MQ=255)
                                                                      gccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGa                   >  2:1217691/1‑93 (MQ=255)
                                                                      gccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGa                   >  2:3874359/1‑93 (MQ=255)
                                                                       ccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGaa                  >  2:4223461/1‑93 (MQ=255)
                                                                        cagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGaaa                 <  2:4571131/93‑1 (MQ=255)
                                                                         agcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATGCTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAAt                <  1:4528416/93‑1 (MQ=255)
                                                                         agcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAAt                >  1:7061999/1‑93 (MQ=255)
                                                                          gcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATc               <  1:1829376/93‑1 (MQ=255)
                                                                          gcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATc               <  2:6169839/93‑1 (MQ=255)
                                                                          gcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATc               <  2:5529222/93‑1 (MQ=255)
                                                                           cgttAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATcc              <  2:636956/90‑1 (MQ=255)
                                                                            ggTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCg             >  1:949615/1‑93 (MQ=255)
                                                                             gTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGt            >  1:6285007/1‑93 (MQ=255)
                                                                             gTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGt            <  1:5546807/93‑1 (MQ=255)
                                                                             gTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGt            <  2:2998851/93‑1 (MQ=255)
                                                                               aaTCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGc          <  1:5638370/93‑1 (MQ=255)
                                                                                aTCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGcc         <  1:3567948/93‑1 (MQ=255)
                                                                                aTCGTCATGCCGATCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTAcc         <  2:7023641/93‑1 (MQ=255)
                                                                                 tCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCa        <  1:1100356/93‑1 (MQ=255)
                                                                                  cGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCag       <  1:3804162/93‑1 (MQ=255)
                                                                                   gTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCaga      <  2:3628969/93‑1 (MQ=255)
                                                                                   gTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCaga      >  2:3706314/1‑93 (MQ=255)
                                                                                    tCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCagag     <  2:3795477/93‑1 (MQ=255)
                                                                                     cATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAgaga    >  2:6871440/1‑93 (MQ=255)
                                                                                     cATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAgaga    <  2:4093946/93‑1 (MQ=255)
                                                                                      aTGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGaa   <  1:7618315/93‑1 (MQ=255)
                                                                                      aTGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGaa   <  2:1951455/93‑1 (MQ=255)
                                                                                       tgccgcTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAAg  >  2:6668195/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAAG  >  NC_000913/2123790‑2123969

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