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,8800TC97.5% 372.2 / ‑0.3 122Q123R (CAG→CGG) wcaHGDP‑mannose mannosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base C (52/66);  total (56/66)
Fisher's exact test for biased strand distribution p-value = 9.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 130 total aligned reads displayed.

TTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAA  >  NC_000913/2123790‑2123968
                                                                                          |                                                                                        
ttCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCAtgccgc                                                                                        >  2:3017130/1‑93 (MQ=255)
ttCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCAtgccgc                                                                                        <  2:510199/93‑1 (MQ=255)
ttCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCAtgccgc                                                                                        >  2:1848943/1‑93 (MQ=255)
 tctcAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATgccgct                                                                                       >  1:3500066/1‑93 (MQ=255)
  ctcAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTc                                                                                      >  2:854874/1‑93 (MQ=255)
    cAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATCCCGCTCAt                                                                                    <  2:2514675/93‑1 (MQ=255)
     aGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATc                                                                                   <  1:1257803/93‑1 (MQ=255)
     aGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATc                                                                                   <  2:3273724/93‑1 (MQ=255)
       cGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCg                                                                                 <  1:3597905/93‑1 (MQ=255)
        gagaAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCgg                                                                                >  2:4681397/1‑93 (MQ=255)
        gagaAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCgg                                                                                >  2:2420923/1‑93 (MQ=255)
        gagaAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCgg                                                                                >  1:2213383/1‑93 (MQ=255)
          gaAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCa                                                                              >  1:3713980/1‑93 (MQ=255)
          gaAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCa                                                                              <  2:1108442/93‑1 (MQ=255)
           aaaaGAGGCGCGGCTGTTAGCATGAACATTATCGCTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAg                                                                             <  1:4251809/93‑1 (MQ=255)
               tAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAAc                                                                         >  1:2119024/1‑93 (MQ=255)
                aGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGTGCGTCCGACGTCAGCCAGCGGTCACTGTCACGCTCCTCTCCCGACCGTAACa                                                                        >  2:3991837/1‑93 (MQ=255)
                aGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCTGCAGTAGCa                                                                        <  2:4114984/93‑1 (MQ=255)
                aGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACa                                                                        >  2:2237999/1‑93 (MQ=255)
                aGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACa                                                                        >  2:1808307/1‑93 (MQ=255)
                 ggCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGCAATCGTCATGCCGCTCATCCGGCAGTAACac                                                                       >  2:3707182/1‑92 (MQ=255)
                   cgcgGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGct                                                                     >  2:1745470/1‑93 (MQ=255)
                    gcgGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGctc                                                                    <  1:3832131/93‑1 (MQ=255)
                     cgGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGctct                                                                   <  1:683810/93‑1 (MQ=255)
                      ggCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtctt                                                                  <  2:1097859/93‑1 (MQ=255)
                      ggCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtctt                                                                  >  1:946423/1‑93 (MQ=255)
                       gCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttc                                                                 >  2:2554278/1‑93 (MQ=255)
                       gCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttc                                                                 <  2:2946966/93‑1 (MQ=255)
                        cTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttct                                                                >  1:1185410/1‑93 (MQ=255)
                        cTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttct                                                                <  2:4164373/93‑1 (MQ=255)
                        cTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttct                                                                <  2:4502876/93‑1 (MQ=255)
                        cTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttct                                                                <  1:1971886/93‑1 (MQ=255)
                        cTGTTAGCATGAACATTATCACTGGCAAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCtcttct                                                                >  1:467879/1‑93 (MQ=255)
                         tGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTcttctt                                                               >  2:1320786/1‑93 (MQ=255)
                          gTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCttcttc                                                              <  1:3603537/93‑1 (MQ=255)
                          gTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCttcttc                                                              >  2:3218035/1‑93 (MQ=255)
                           ttAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTcttctttc                                                             >  2:921827/1‑91 (MQ=255)
                           ttAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTcc                                                             <  1:2716441/93‑1 (MQ=255)
                             aGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGa                                                           <  1:1730119/93‑1 (MQ=255)
                              gCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGAt                                                          <  2:703871/93‑1 (MQ=255)
                               cATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATa                                                         <  2:1589952/93‑1 (MQ=255)
                                aTGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAc                                                        <  1:741006/93‑1 (MQ=255)
                                 tGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAct                                                       >  1:2568850/1‑93 (MQ=255)
                                 tGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAct                                                       >  2:2521217/1‑93 (MQ=255)
                                   aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATActct                                                     <  1:2152136/93‑1 (MQ=255)
                                   aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATActct                                                     <  1:961230/93‑1 (MQ=255)
                                    aCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTg                                                    >  2:2016817/1‑93 (MQ=255)
                                    aCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTg                                                    <  1:4166951/93‑1 (MQ=255)
                                     cATTATCACGGGCGAGCAGCGCGTCCGACGTCATGCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGa                                                   <  2:836640/93‑1 (MQ=255)
                                      aTTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGaa                                                  <  1:1552135/93‑1 (MQ=255)
                                        tATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGc                                                <  1:3424616/93‑1 (MQ=255)
                                         aTCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCg                                               >  1:4375196/1‑93 (MQ=255)
                                             cTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGaaaa                                           >  2:3722607/1‑93 (MQ=255)
                                             cTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGaaaa                                           <  1:3928269/93‑1 (MQ=255)
                                              tGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAAc                                          <  1:4235565/93‑1 (MQ=255)
                                              tGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAAc                                          >  2:4239433/1‑93 (MQ=255)
                                              tGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAAc                                          <  1:1477578/93‑1 (MQ=255)
                                               ggCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAAcc                                         <  1:4396895/93‑1 (MQ=255)
                                                gctagcagcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCg                                        <  1:783827/90‑1 (MQ=255)
                                                gCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCg                                        <  2:3406648/93‑1 (MQ=255)
                                                gCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCg                                        <  2:4013716/93‑1 (MQ=255)
                                                   agcagcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGc                                     >  1:2788666/1‑93 (MQ=255)
                                                      agcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGcacc                                  <  1:1463997/93‑1 (MQ=255)
                                                      agcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGcacc                                  <  2:4522488/93‑1 (MQ=255)
                                                      agcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGcacc                                  <  1:2675086/93‑1 (MQ=255)
                                                          cgTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACAt                              >  1:3893705/1‑93 (MQ=255)
                                                          cgTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACAt                              <  2:138040/93‑1 (MQ=255)
                                                          cgTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACAt                              <  1:1704986/93‑1 (MQ=255)
                                                           gTCCGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATa                             <  2:1068522/93‑1 (MQ=255)
                                                             ccGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagt                           <  1:3274376/93‑1 (MQ=255)
                                                             ccGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagt                           <  2:2708499/93‑1 (MQ=255)
                                                              cGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtg                          >  2:681727/1‑93 (MQ=255)
                                                              cGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtg                          <  2:1353722/93‑1 (MQ=255)
                                                              cGACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtg                          <  1:1754431/93‑1 (MQ=255)
                                                               gatgTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtga                         >  1:3130177/4‑93 (MQ=255)
                                                               gACGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtga                         <  1:3666170/93‑1 (MQ=255)
                                                                aCGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtgag                        <  1:1305218/93‑1 (MQ=255)
                                                                 cGTCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATagtgagt                       <  1:4078029/93‑1 (MQ=255)
                                                                   tCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAgtgg                     <  1:1973833/93‑1 (MQ=255)
                                                                   tCAGCCAGCGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAgggg                     >  2:2026163/1‑93 (MQ=255)
                                                                    cagccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAgtggt                    <  1:530663/93‑1 (MQ=255)
                                                                     agccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAgtggtg                   <  2:2308047/93‑1 (MQ=255)
                                                                      gccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGa                  >  1:2015349/1‑93 (MQ=255)
                                                                       ccagcGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGaa                 >  1:3133591/1‑93 (MQ=255)
                                                                       ccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAgtggtgca                 >  1:1986777/1‑91 (MQ=255)
                                                                       ccagcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGaa                 >  1:1986776/1‑93 (MQ=255)
                                                                         agcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAAt               <  1:782862/93‑1 (MQ=255)
                                                                          gcGGTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATc              >  2:3437201/1‑93 (MQ=255)
                                                                             gTAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGt           >  1:343672/1‑93 (MQ=255)
                                                                              tAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTg          >  1:4306111/1‑93 (MQ=255)
                                                                              tAATCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTg          <  2:1857821/93‑1 (MQ=255)
                                                                               aaTCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGc         >  1:1185702/1‑93 (MQ=255)
                                                                                aTCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTTAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGcc        >  2:4510304/1‑93 (MQ=255)
                                                                                aTCGTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGcc        <  2:2884646/93‑1 (MQ=255)
                                                                                   gTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCaga     >  1:860917/1‑93 (MQ=255)
                                                                                   gTCATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCaga     >  1:2277018/1‑93 (MQ=255)
                                                                                    tCATGCGGCTCATCCGGCAGTAACAGCTCTTCTTGCGATACTCTGAAGCGCAAACCGAGCACCACCTAGTGAGTGGTGAAATCCGTGCCcgag    >  2:2191157/1‑93 (MQ=255)
                                                                                     cATGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAgaga   <  2:3973972/93‑1 (MQ=255)
                                                                                      aTGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGGCAGAGaa  >  1:2417139/1‑93 (MQ=255)
                                                                                      aTGCCGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGaa  >  1:878528/1‑93 (MQ=255)
                                                                                          |                                                                                        
TTCTCAGCGAGAAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTCCGATACTCTGAAGCGAAAACCGAGCACCACATAGTGAGTGGTGAAATCCGTGCCAGAGAA  >  NC_000913/2123790‑2123968

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: