Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,099,309 (G)7→8 coding (48/993 nt) kdgT → 2‑keto‑3‑deoxy‑D‑gluconate transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,099,3021.G95.7% 289.7 / 5.3 94coding (41/993 nt)kdgT2‑keto‑3‑deoxy‑D‑gluconate transport system
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base G (45/45);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 128 total aligned reads displayed.

CTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCA  >  NC_000913/4099211‑4099391
                                                                                            |                                                                                          
cTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑g                                                                                          <  2:2604752/93‑1 (MQ=255)
 tGATGCTCCGTTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGGTCGATTGAGAAACTCCC‑‑gg                                                                                         >  2:1155079/1‑93 (MQ=255)
 tGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑gg                                                                                         <  1:3363527/93‑1 (MQ=255)
 tGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑gg                                                                                         <  2:3233421/93‑1 (MQ=255)
  gATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggg                                                                                        >  2:116097/1‑93 (MQ=255)
  gATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggg                                                                                        >  2:1472027/1‑93 (MQ=255)
  gATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggg                                                                                        <  1:2935277/93‑1 (MQ=255)
  gATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggg                                                                                        >  2:2535255/1‑93 (MQ=255)
   aTGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑gggg                                                                                       >  1:3548501/1‑93 (MQ=255)
   aTGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑gggg                                                                                       <  2:395639/93‑1 (MQ=255)
   aTGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑gggg                                                                                       <  2:818511/93‑1 (MQ=255)
       tCCATTGAGGTATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑gggggggg                                                                                   >  2:2798525/1‑92 (MQ=255)
       tCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGACGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑gggggggg                                                                                   >  1:388480/1‑92 (MQ=255)
       cccATTGAGGAATTACGCATCAGCCCTAAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑gggggggg                                                                                   <  2:1995320/92‑2 (MQ=255)
        ccATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggcggggga                                                                                  >  2:3027426/1‑91 (MQ=255)
         cATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCc‑‑gggcggggag                                                                                 >  1:2288967/1‑90 (MQ=255)
         cATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggat                                                                                 <  1:503382/93‑4 (MQ=255)
         cATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggat                                                                                 >  2:2959407/1‑90 (MQ=255)
         cATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggag                                                                                 >  2:3915990/1‑90 (MQ=255)
          aTTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAACCGCTCGATTGAGAAAATCCC‑‑ggggggggatg                                                                                >  1:1350333/1‑89 (MQ=255)
          aTTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggatg                                                                                >  1:1350332/1‑89 (MQ=255)
          aTTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggatg                                                                                <  1:3701431/93‑5 (MQ=255)
          aTTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggatg                                                                                >  2:289897/1‑89 (MQ=255)
          aTTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggatg                                                                                >  1:1350334/1‑89 (MQ=255)
             gAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGatgatg                                                                              >  2:2104501/1‑93 (MQ=255)
               ggAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCt                                                                            >  2:4252545/1‑93 (MQ=255)
               ggAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCt                                                                            >  1:790525/1‑93 (MQ=255)
                  aTTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGGTGATGCTCGt                                                                         >  1:3383132/1‑93 (MQ=255)
                   ttACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGAGGCTCGTc                                                                        >  1:924335/1‑93 (MQ=255)
                   ttACACATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATCCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTc                                                                        <  2:3141502/93‑1 (MQ=255)
                     aCGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTccc                                                                      <  2:3666909/93‑1 (MQ=255)
                       gCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAGTCCCGGGGGGGGATGATGCTCGTCCCGc                                                                    <  1:4006189/93‑1 (MQ=255)
                       gCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGGATGATGCTCGTcccc                                                                     >  1:1286988/1‑92 (MQ=255)
                       gCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGc                                                                    <  2:848805/93‑1 (MQ=255)
                         aTCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGGGTGATGGTCGTCCCGCt                                                                   >  1:1052925/1‑93 (MQ=255)
                             gCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTcc                                                              >  2:2073492/1‑93 (MQ=255)
                               ccTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCtt                                                            >  1:3569587/1‑93 (MQ=255)
                               ccTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCtt                                                            >  1:89959/1‑93 (MQ=255)
                                cTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTg                                                           >  1:696436/1‑93 (MQ=255)
                                 ttAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTgg                                                          >  1:1000868/1‑93 (MQ=255)
                                 ttAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTgg                                                          >  1:3002399/1‑93 (MQ=255)
                                      aaTATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCAc                                                     >  1:3130574/1‑93 (MQ=255)
                                         atGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACtgt                                                  <  1:347960/93‑1 (MQ=255)
                                          tGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACtgtg                                                 <  1:470979/93‑1 (MQ=255)
                                               aCAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCcaca                                            >  2:1441856/1‑93 (MQ=255)
                                               aCAGGGGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑GGGGGGGTCAAAGCTCGTCCCGCTGTTCCTTCGCGCACTGTGCCccac                                           >  2:4075996/1‑93 (MQ=255)
                                                cAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCacac                                           <  2:4039564/93‑1 (MQ=255)
                                                cAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCacac                                           >  1:24670/1‑93 (MQ=255)
                                                 aGGTGATGGAAATGCAGATAAAACGCTCGATTGTGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACAcc                                          <  1:4583236/93‑1 (MQ=255)
                                                 aGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACAcc                                          >  2:778002/1‑93 (MQ=255)
                                                  ggTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCCCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCt                                         <  2:508969/93‑1 (MQ=255)
                                                    tGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtc                                       <  1:2298085/93‑1 (MQ=255)
                                                    tGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtc                                       <  1:2298086/93‑1 (MQ=255)
                                                     gATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtcc                                      >  1:633786/1‑92 (MQ=255)
                                                     gATGGAAATGCAGAAAAAACGCGAGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtct                                      <  1:1275844/93‑1 (MQ=255)
                                                       tGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCg                                    >  1:4553183/1‑93 (MQ=255)
                                                       tGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCg                                    >  1:400900/1‑93 (MQ=255)
                                                        ggAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGc                                   <  2:4598555/93‑1 (MQ=255)
                                                         gAAATACAGAAAAACCGCTCGATTGAAAAAACCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGcc                                  <  1:136033/93‑1 (MQ=255)
                                                            aTGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCggg                               <  1:718076/93‑1 (MQ=255)
                                                            aTGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCggg                               <  2:4376635/93‑1 (MQ=255)
                                                            aTGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCggg                               >  1:239419/1‑93 (MQ=255)
                                                            aTGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCgggg                              >  2:808734/1‑93 (MQ=255)
                                                            aTGCAAATAAAACCCTCGATTGAAAAAACCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCggg                               <  2:3482284/93‑1 (MQ=255)
                                                             agCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCgggg                              <  1:4365634/92‑1 (MQ=255)
                                                              tcAGATAAAACGCTCTATTGAGAAAAGCCCGGGGGGGGATTATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCggggg                             <  2:1578387/92‑1 (MQ=255)
                                                              gCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCggggg                             <  1:2352793/93‑1 (MQ=255)
                                                               cAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGc                            <  2:3919654/93‑1 (MQ=255)
                                                               cAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGc                            >  2:832130/1‑93 (MQ=255)
                                                                aGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCg                           >  2:2997419/1‑93 (MQ=255)
                                                                aGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCg                           >  2:4555686/1‑93 (MQ=255)
                                                                aGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCg                           <  1:1602555/93‑1 (MQ=255)
                                                                 gATAAAACTCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCgg                          <  2:1439632/93‑1 (MQ=255)
                                                                 gATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATTCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCgg                          <  2:1529872/93‑1 (MQ=255)
                                                                    aaaaCGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCTACACCTTCTCGCCGGGGGCGgaga                       >  2:3287091/1‑93 (MQ=255)
                                                                     aaaCGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGaa                      >  1:241557/1‑93 (MQ=255)
                                                                     aaaCGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGaa                      >  2:3431298/1‑93 (MQ=255)
                                                                      aaCGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGaaa                     >  2:4413367/1‑93 (MQ=255)
                                                                       aCGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGGGGGGAAat                    >  1:3869692/1‑93 (MQ=255)
                                                                         gCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAatat                  >  2:4405128/1‑93 (MQ=255)
                                                                          cTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATAtt                 <  2:1562219/93‑1 (MQ=255)
                                                                            cGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATAtttt               >  2:86573/1‑93 (MQ=255)
                                                                            cGATTGAGAAAATCCC‑‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTg              <  1:2047301/93‑1 (MQ=255)
                                                                             gATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTg              >  2:2623371/1‑93 (MQ=255)
                                                                             gATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTg              <  2:1185783/93‑1 (MQ=255)
                                                                               ttGAGAAAATCCCGGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTgg             >  2:2040564/1‑93 (MQ=255)
                                                                               ttGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGa            >  2:1419878/1‑93 (MQ=255)
                                                                                 gagaAAGCCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATc          <  1:1775925/93‑1 (MQ=255)
                                                                                 gagaAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATc          <  2:2312885/93‑1 (MQ=255)
                                                                                 gagaAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATc          <  1:2245029/93‑1 (MQ=255)
                                                                                  agaAAGTCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCa         <  1:7413/93‑1 (MQ=255)
                                                                                  agaAAATCCCGGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATc          <  2:3861564/93‑1 (MQ=255)
                                                                                    aaaaTCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCAtt       >  2:3214670/1‑93 (MQ=255)
                                                                                    aaaaTCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCAtt       <  1:2015614/93‑1 (MQ=255)
                                                                                    aaaaTCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCAtt       <  2:772639/93‑1 (MQ=255)
                                                                                    aaaaTCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCAtt       <  2:503287/93‑1 (MQ=255)
                                                                                      aaTCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCa     <  1:1439048/93‑1 (MQ=255)
                                                                                      aaTCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCCCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAAAATTTTGGatcaatca     >  2:2381517/1‑93 (MQ=255)
                                                                                       aTCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCAc    <  2:2397788/93‑1 (MQ=255)
                                                                                       aTCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCAc    >  1:1630696/1‑93 (MQ=255)
                                                                                            |                                                                                          
CTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCA  >  NC_000913/4099211‑4099391

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