Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,106,773 A→G V120V (GTA→GTG sbp → periplasmic sulfate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,106,7730AG99.8% 1558.3 / ‑5.4 479V120V (GTA→GTGsbpperiplasmic sulfate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (220/258);  total (220/259)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

Only 100 of 489 total aligned reads displayed.

AATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGG  >  NC_000913/4106683‑4106862
                                                                                          |                                                                                         
 gttGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGt                                                                                        >  1:2943465/2‑93 (MQ=255)
  ttGCGGAGCGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTa                                                                                       <  2:8569296/93‑1 (MQ=255)
  ttGCGAAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTa                                                                                       <  1:7524325/93‑1 (MQ=255)
   tGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTaa                                                                                      <  1:1653485/93‑1 (MQ=255)
   tGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTaa                                                                                      <  2:1926986/93‑1 (MQ=255)
     cGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAgg                                                                                    <  1:5135801/93‑1 (MQ=255)
      ggAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAggg                                                                                   <  2:7716210/93‑1 (MQ=255)
          cgcgGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAAt                                                                               <  1:9438227/93‑1 (MQ=255)
           gcgggcggATTGATAAAGAGTGGATCAAGCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATc                                                                              <  1:6963960/93‑1 (MQ=255)
           gcgggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATc                                                                              <  2:5880542/93‑1 (MQ=255)
              ggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGa                                                                           >  2:4545420/1‑93 (MQ=255)
              ggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGa                                                                           >  2:2131212/1‑93 (MQ=255)
               gcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCgta                                                                          >  1:4445843/1‑91 (MQ=255)
               gcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGaa                                                                          <  2:3813114/93‑1 (MQ=255)
               gcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGaa                                                                          <  1:5958197/93‑1 (MQ=255)
               gcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGaa                                                                          >  1:2191122/1‑93 (MQ=255)
                cggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAg                                                                         <  1:2023168/93‑1 (MQ=255)
                   aTTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAg                                                                      >  2:567162/1‑93 (MQ=255)
                    ttGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGa                                                                     >  1:1684890/1‑93 (MQ=255)
                    ttGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGa                                                                     <  2:838348/93‑1 (MQ=255)
                    ttGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGa                                                                     >  1:9612904/1‑93 (MQ=255)
                    ttGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGa                                                                     >  2:7915943/1‑93 (MQ=255)
                     tGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGAt                                                                    >  1:4641849/1‑93 (MQ=255)
                        tAAAGAGTGGATCAAACGTCTGCCGGATACCTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCa                                                                 <  2:2562784/93‑1 (MQ=255)
                         aaaGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCAt                                                                <  1:6454437/93‑1 (MQ=255)
                           agagTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGa                                                              >  1:7171420/1‑93 (MQ=255)
                           agagTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGa                                                              >  2:6306640/1‑93 (MQ=255)
                           agagTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGa                                                              <  2:9442830/93‑1 (MQ=255)
                            gagTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCTACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGAc                                                             >  2:4455095/1‑93 (MQ=255)
                            gagTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGAc                                                             <  1:4166286/93‑1 (MQ=255)
                             agTGGATCAAACGGCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACt                                                            <  1:2431666/93‑1 (MQ=255)
                              gTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTg                                                           >  2:2880478/1‑93 (MQ=255)
                                ggATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCGATGACTGGa                                                         >  1:3227372/1‑93 (MQ=255)
                                ggATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGa                                                         <  2:1698941/93‑1 (MQ=255)
                                  aTCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAAc                                                       <  1:14268/93‑1 (MQ=255)
                                    cAAACGTCTGCCGTATAACTCCGCACCGTACACTTCCACAATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGa                                                     <  1:10003331/93‑1 (MQ=255)
                                    cAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGa                                                     >  2:7133092/1‑93 (MQ=255)
                                    cAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGa                                                     >  1:3995160/1‑93 (MQ=255)
                                     aaaCGTCTGCGGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGTTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGAt                                                    <  1:5893905/93‑1 (MQ=255)
                                     aaaCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGAt                                                    <  1:8688097/93‑1 (MQ=255)
                                      aaCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATc                                                   <  1:4535402/93‑1 (MQ=255)
                                      aaCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATc                                                   >  1:1289882/1‑93 (MQ=255)
                                        cGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTg                                                 >  1:6663271/1‑93 (MQ=255)
                                         gTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCTGATCCATGACTGGAACGATCTGa                                                >  2:7229492/1‑93 (MQ=255)
                                         gTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGa                                                >  2:101171/1‑93 (MQ=255)
                                         gTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGa                                                >  1:2622686/1‑93 (MQ=255)
                                         gTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGa                                                <  2:5084424/93‑1 (MQ=255)
                                          tCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGAt                                               >  1:6344164/1‑93 (MQ=255)
                                           cTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGAtt                                              <  1:2806253/93‑1 (MQ=255)
                                            tGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTa                                             >  1:3782487/1‑93 (MQ=255)
                                             gCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTaa                                            <  2:2332432/93‑1 (MQ=255)
                                               cgaaTACCTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAAc                                          <  1:7871584/90‑1 (MQ=255)
                                               cGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAAc                                          <  2:7474025/93‑1 (MQ=255)
                                                ggATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAAcc                                         >  1:7987155/1‑93 (MQ=255)
                                                 gATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATTACTGGAACGATCTGATTAAACCg                                        >  2:5584824/1‑93 (MQ=255)
                                                 gATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCg                                        >  1:5724047/1‑93 (MQ=255)
                                                  aTAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCgg                                       >  2:4933048/1‑93 (MQ=255)
                                                  aTAACTCAGCACCGTGCCCTTCCACCATTGTTTTCCTGGTGAGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCgg                                       <  1:8295100/93‑1 (MQ=255)
                                                   tAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCggg                                      <  1:8526969/93‑1 (MQ=255)
                                                   tAACTCCGCACCGTACACTTCAACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCggg                                      <  2:7049266/93‑1 (MQ=255)
                                                       tCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGtt                                  <  1:9719218/93‑1 (MQ=255)
                                                        ccGCACCGTCCACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGttt                                 <  2:8174469/93‑1 (MQ=255)
                                                        ccGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGttt                                 >  2:406717/1‑93 (MQ=255)
                                                        ccGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGttt                                 <  2:984874/93‑1 (MQ=255)
                                                         cGCACCGTACCCTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAACCCGGGTGTTTc                                <  1:1751366/93‑1 (MQ=255)
                                                         cGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGTAACGATCTGATTAAACCGGGTGTTTc                                <  2:3194883/93‑1 (MQ=255)
                                                         cGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTc                                >  1:7325663/1‑93 (MQ=255)
                                                          gCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCg                               >  2:6812958/1‑93 (MQ=255)
                                                          gCACCGTACACTTCCACCATTGTTTTCCTGGGGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCg                               >  2:7975781/1‑93 (MQ=255)
                                                            aCCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGt                             >  1:6223172/1‑93 (MQ=255)
                                                               gTACACTTCCCCCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGAt                          <  2:4793228/93‑1 (MQ=255)
                                                                 acacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCa                        <  1:7092423/93‑1 (MQ=255)
                                                                 acacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCa                        <  2:332445/93‑1 (MQ=255)
                                                                 acacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCa                        <  2:8801186/93‑1 (MQ=255)
                                                                 acacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCa                        <  2:1311238/93‑1 (MQ=255)
                                                                 acacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCa                        <  1:9208484/93‑1 (MQ=255)
                                                                    cTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGc                     <  1:3634979/93‑1 (MQ=255)
                                                                     ttCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGcc                    >  1:9073057/1‑93 (MQ=255)
                                                                      tCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCt                   >  1:8897328/1‑93 (MQ=255)
                                                                        caccaTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTaa                 >  1:3845541/1‑93 (MQ=255)
                                                                          ccaTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATc               >  1:6526326/1‑93 (MQ=255)
                                                                          ccaTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATc               >  2:6566308/1‑93 (MQ=255)
                                                                           caTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATcc              <  2:9001696/93‑1 (MQ=255)
                                                                           caTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATcc              >  1:3536543/1‑93 (MQ=255)
                                                                           caTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATcc              <  2:2222502/93‑1 (MQ=255)
                                                                           caTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGATTCGGTGATCCCGCCTAATcc              >  2:3558161/1‑93 (MQ=255)
                                                                               gTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaaa          <  1:5641803/93‑1 (MQ=255)
                                                                               gTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaaa          <  1:5279933/93‑1 (MQ=255)
                                                                                ttttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGCTCACGCCTAATTCGcaag         >  1:9312312/1‑92 (MQ=38)
                                                                                ttttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaaaa         <  2:1156970/93‑1 (MQ=255)
                                                                                ttttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaaaa         <  1:4889725/93‑1 (MQ=255)
                                                                                ttttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaaaa         <  2:1413997/93‑1 (MQ=255)
                                                                                  ttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGc       >  2:2014260/1‑93 (MQ=255)
                                                                                  ttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGc       >  2:6113501/1‑93 (MQ=255)
                                                                                    ccTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGctc     >  1:3300224/1‑93 (MQ=255)
                                                                                    ccTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGctc     <  2:9678404/93‑1 (MQ=255)
                                                                                     cTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGctct    <  1:8120295/93‑1 (MQ=255)
                                                                                      tGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtg   <  1:7742064/93‑1 (MQ=255)
                                                                                       ggTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtgg  >  2:6440118/1‑93 (MQ=255)
                                                                                       ggTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtgg  <  1:3159924/93‑1 (MQ=255)
                                                                                          |                                                                                         
AATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGG  >  NC_000913/4106683‑4106862

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: