Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,216,746 G→C S252S (TCC→TCG yehZ ← putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,216,7460GC99.5% 742.9 / ‑6.7 217S252S (TCC→TCGyehZputative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (100/116);  minor base T (0/1);  total (100/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 225 total aligned reads displayed.

GCTGGCATTCAGTTGCTGCAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACGGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCC  >  NC_000913/2216656‑2216835
                                                                                          |                                                                                         
gCTGGCATTCAGTTGCTGCAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTTCGGATACTCCCTCAACACCGa                                                                                         >  2:3217705/1‑93 (MQ=255)
 cTGGCATTCAGTTGCTGCAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGAc                                                                                        <  1:4952947/93‑1 (MQ=255)
 cTGGCATTCAGTTGCTGCAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGAc                                                                                        <  1:1983153/93‑1 (MQ=255)
    gCATTCAGTTGCTGCAATGTTTTTGCATCGAGGCTGGCGAATACTGGCTGTATCCACTGAGCCATTTGTCGATACTCCCTTAACACCGAcacc                                                                                     >  1:5607166/1‑90 (MQ=255)
     cATTCAGTTGCTGCAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCAc                                                                                    >  2:6431043/1‑93 (MQ=255)
      aTTCAGTTGCTGCAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACg                                                                                   >  1:1797548/1‑93 (MQ=255)
       ttCAGTTGCTGCAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGc                                                                                  >  2:4350555/1‑93 (MQ=255)
           gTTGCTGCAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCCACACCGACTCACGcacca                                                                              >  1:899796/1‑93 (MQ=255)
           gTTGCTGCAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCACCACCCACTACCGcacca                                                                              >  2:2877069/1‑93 (MQ=255)
           gTTGCTGCAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGcaccc                                                                              >  1:6020602/1‑92 (MQ=255)
               ctgcAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTgg                                                                          <  1:5005370/93‑1 (MQ=255)
               ctgcAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTgg                                                                          <  1:2314547/93‑1 (MQ=255)
                tgcAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGt                                                                         >  1:3634144/1‑93 (MQ=255)
                 gctaTGTTTTTGCATCGAGGCTGGCTAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTg                                                                        <  2:1546158/90‑1 (MQ=255)
                   aaTGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCa                                                                      <  2:2466814/93‑1 (MQ=255)
                    aTGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAg                                                                     >  1:345302/1‑93 (MQ=255)
                      gTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGgc                                                                   <  1:2732095/93‑1 (MQ=255)
                       tttttGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGgcg                                                                  >  2:1982284/1‑93 (MQ=255)
                        ttttGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGgcgc                                                                 <  1:233659/93‑1 (MQ=255)
                          ttGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGt                                                               >  1:3168265/1‑93 (MQ=255)
                          ttGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGt                                                               <  2:5433984/93‑1 (MQ=255)
                            gCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAg                                                             >  1:5418408/1‑93 (MQ=255)
                            gCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAg                                                             >  1:3646225/1‑93 (MQ=255)
                             cATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGa                                                            <  1:6458364/93‑1 (MQ=255)
                             cATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGa                                                            <  1:5278749/93‑1 (MQ=255)
                              aTCGAGGCTGGCGAAGACTGGCTGTAGCCCCTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGAt                                                           >  2:1683008/1‑93 (MQ=255)
                               tCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATa                                                          <  2:623009/93‑1 (MQ=255)
                               tCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATa                                                          >  1:1834973/1‑93 (MQ=255)
                                cGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCACGCGCGTAGATAg                                                         >  1:5928471/1‑93 (MQ=255)
                                 gAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAgg                                                        <  2:2054489/93‑1 (MQ=255)
                                  aGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGGTACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGt                                                       >  1:6393311/1‑93 (MQ=255)
                                  aGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCGACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGt                                                       >  2:4132957/1‑93 (MQ=255)
                                  aGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGt                                                       <  2:520391/93‑1 (MQ=255)
                                  aGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGt                                                       <  1:4194743/93‑1 (MQ=255)
                                   ggCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGtt                                                      >  1:4501153/1‑93 (MQ=255)
                                     cTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGc                                                    <  1:5947788/93‑1 (MQ=255)
                                     cTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGGGCAGGCGCGTAGATAGGTTGc                                                    >  2:2643561/1‑93 (MQ=255)
                                       ggCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGcac                                                  >  1:3812693/1‑93 (MQ=255)
                                        gCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGcaca                                                 <  2:1172340/93‑1 (MQ=255)
                                          gAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACAcc                                               >  2:3643965/1‑93 (MQ=255)
                                           aaGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCt                                              <  2:547998/93‑1 (MQ=255)
                                           aaGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCt                                              <  2:5911946/93‑1 (MQ=255)
                                           aaGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCt                                              <  1:2683970/93‑1 (MQ=255)
                                            aGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCtt                                             >  1:680926/1‑93 (MQ=255)
                                            aGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCtt                                             <  1:1329838/93‑1 (MQ=255)
                                            aGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCtt                                             <  1:1375163/93‑1 (MQ=255)
                                              ccTGGCTGTAGCCCCTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCCCCACTGGTGCAGGCGCGTAGATAGGTTGCGCCCCTTGc                                           <  1:785085/92‑1 (MQ=255)
                                              aCTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGc                                           <  2:578432/93‑1 (MQ=255)
                                              aCTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGc                                           <  2:1539956/93‑1 (MQ=255)
                                              aCTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGc                                           <  2:5014686/93‑1 (MQ=255)
                                               cTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCg                                          <  1:3517752/93‑1 (MQ=255)
                                                 ggCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGa                                        <  2:4897673/93‑1 (MQ=255)
                                                 ggCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGa                                        >  1:348561/1‑93 (MQ=255)
                                                 ggCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGa                                        >  2:4051861/1‑93 (MQ=255)
                                                 ggCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTCCACACCTTGCGGa                                        >  2:3731896/1‑93 (MQ=255)
                                                   cTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATc                                      >  1:257608/1‑93 (MQ=255)
                                                    tGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCg                                     >  2:3844305/1‑93 (MQ=255)
                                                    tGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCg                                     >  1:2965798/1‑93 (MQ=255)
                                                    tGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTCCACACCTTGCGGATCg                                     >  1:1480026/1‑93 (MQ=255)
                                                     gTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGc                                    >  2:5257818/1‑93 (MQ=255)
                                                         ccACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTaa                                <  1:2257383/93‑1 (MQ=255)
                                                         ccACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTaa                                >  2:1259864/1‑93 (MQ=255)
                                                           aCTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAgg                              >  2:815555/1‑93 (MQ=255)
                                                            cTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGt                             >  1:2913282/1‑93 (MQ=255)
                                                            cTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGt                             <  2:6259316/93‑1 (MQ=255)
                                                            cTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGt                             <  2:447761/93‑1 (MQ=255)
                                                            cTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCAGATCGCTTAAGGt                             >  2:3665938/1‑93 (MQ=255)
                                                             tgtgCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGtt                            <  2:2383098/93‑1 (MQ=255)
                                                             tgtgCCATGTGCGGATACTCCCTCAACGCCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAggct                            <  1:5064146/93‑3 (MQ=255)
                                                              gtgCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGttt                           >  1:61362/1‑93 (MQ=255)
                                                               tgCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTg                          <  2:2954839/93‑1 (MQ=255)
                                                                gCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGc                         <  2:660558/93‑1 (MQ=255)
                                                                gCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGc                         >  2:5063530/1‑93 (MQ=255)
                                                                 ccATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCa                        <  2:1908379/93‑1 (MQ=255)
                                                                 ccATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCa                        <  1:1687161/93‑1 (MQ=255)
                                                                   aTTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGc                      >  2:5287266/1‑93 (MQ=255)
                                                                    tttGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGcc                     <  2:1635651/93‑1 (MQ=255)
                                                                    tttGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGcc                     <  2:4715538/93‑1 (MQ=255)
                                                                    tttGCGGATACTCCCCCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGcc                     >  2:5352006/1‑93 (MQ=255)
                                                                     ttGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGccc                    <  2:5965005/93‑1 (MQ=255)
                                                                      tGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGcccc                   <  2:4173708/93‑1 (MQ=255)
                                                                      tGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGcccc                   <  1:5664222/93‑1 (MQ=255)
                                                                      tGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGcccc                   <  1:1097854/93‑1 (MQ=255)
                                                                        cGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAg                 <  2:4558175/93‑1 (MQ=255)
                                                                        cGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAg                 <  1:3723518/93‑1 (MQ=255)
                                                                          gATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAgcg               >  1:959524/1‑93 (MQ=255)
                                                                            tACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGcgcc             >  2:5750389/1‑93 (MQ=255)
                                                                            tACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGcgcc             >  2:4746351/1‑93 (MQ=255)
                                                                             ccTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGcgccg            <  2:2674961/92‑1 (MQ=255)
                                                                              cTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGcgccgc           >  1:3297784/1‑93 (MQ=255)
                                                                              cTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGcgccgc           <  2:5668894/93‑1 (MQ=255)
                                                                                cccTCACCACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGTTTTGCAGCCCCAGCGCCGCTa         <  1:4091353/93‑1 (MQ=255)
                                                                                cccTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTa         >  1:2036557/1‑93 (MQ=255)
                                                                                  cTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTAcc       >  1:1567522/1‑93 (MQ=255)
                                                                                  cTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTAcc       >  1:4657209/1‑93 (MQ=255)
                                                                                   tCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCg      <  1:2090624/93‑1 (MQ=255)
                                                                                    cAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCgg     <  2:6361356/93‑1 (MQ=255)
                                                                                      acacCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGc   >  1:5084521/1‑93 (MQ=255)
                                                                                      acacCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGc   <  2:2207089/93‑1 (MQ=255)
                                                                                      acacCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGc   >  2:3286545/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCTGGCATTCAGTTGCTGCAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACGGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCC  >  NC_000913/2216656‑2216835

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