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.G98.8% 514.2 / ‑0.8 162coding (41/993 nt)kdgT2‑keto‑3‑deoxy‑D‑gluconate transport system
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (90/69);  total (91/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

Only 100 of 181 total aligned reads displayed.

CTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCA  >  NC_000913/4099211‑4099388
                                                                                            |                                                                                       
cTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑g                                                                                       >  2:5254041/1‑93 (MQ=255)
cTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑g                                                                                       <  2:250183/93‑1 (MQ=255)
  gATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggg                                                                                     >  1:4820650/1‑93 (MQ=255)
      cTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggg                                                                                 <  1:1884937/93‑1 (MQ=255)
      cTCCATTGAGAAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggg                                                                                 <  1:3231751/93‑1 (MQ=255)
       tCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑gggggggg                                                                                >  2:1828789/1‑92 (MQ=255)
         cATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggat                                                                              >  1:4626926/1‑90 (MQ=255)
         cATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggat                                                                              <  1:2129964/93‑4 (MQ=255)
         cATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggat                                                                              >  1:3161850/1‑90 (MQ=255)
         cATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggag                                                                              >  2:4128476/1‑90 (MQ=255)
          aTTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggatg                                                                             >  2:933967/1‑89 (MQ=255)
              aGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGc                                                                          >  2:6537205/1‑93 (MQ=255)
               ggAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCt                                                                         >  2:3597313/1‑93 (MQ=255)
                 aaTTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCg                                                                       >  2:5037966/1‑93 (MQ=255)
                 aaTTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCg                                                                       >  2:6150063/1‑93 (MQ=255)
                     aCGCATCAGCCCTTAAAAATATGCCGACAGGAGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTccc                                                                   >  1:2855736/1‑93 (MQ=255)
                        cATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCt                                                                >  1:2266809/1‑93 (MQ=255)
                        cATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCt                                                                >  2:1123147/1‑93 (MQ=255)
                        cATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGAAGATAAAACGCTCGATTGAGAAAACCCCGGGGGGGGATGATGCTCGTCCCGCt                                                                <  1:6303926/93‑1 (MQ=255)
                          tCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTAt                                                              <  2:1814545/93‑1 (MQ=255)
                          tCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTAt                                                              <  1:6729152/93‑1 (MQ=255)
                          tCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTAt                                                              >  1:2070197/1‑93 (MQ=255)
                          tCAGCCCTAAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTAt                                                              <  1:4004418/93‑1 (MQ=255)
                           cAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTAtt                                                             >  2:4289562/1‑93 (MQ=255)
                           cAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTAtt                                                             >  2:4183579/1‑93 (MQ=255)
                            aGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTc                                                            <  2:3067205/93‑1 (MQ=255)
                             gCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTcc                                                           >  1:5750445/1‑93 (MQ=255)
                              cccTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCt                                                          <  2:841052/93‑1 (MQ=255)
                              cccTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCt                                                          >  2:2689279/1‑93 (MQ=255)
                               ccTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCtt                                                         >  1:1298126/1‑93 (MQ=255)
                               ccTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCtt                                                         >  2:2134081/1‑93 (MQ=255)
                               ccTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCtt                                                         <  2:4029566/93‑1 (MQ=255)
                                cTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTg                                                        <  1:5145517/93‑1 (MQ=255)
                                cTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTg                                                        >  2:2580582/1‑93 (MQ=255)
                                 ttAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTgg                                                       >  2:5394192/1‑93 (MQ=255)
                                   aaaaaTATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGgcg                                                     >  1:5732501/1‑93 (MQ=255)
                                   aaaaaTATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGgcg                                                     >  2:1187064/1‑93 (MQ=255)
                                     aaaTATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGgcgcc                                                   >  1:4776327/1‑92 (MQ=255)
                                      aaTATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCAc                                                  >  2:2901700/1‑93 (MQ=255)
                                      aaTATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCAGCTATTCCTTGGCGAAc                                                  >  1:5339899/1‑93 (MQ=255)
                                       atatGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACt                                                 <  2:5127127/93‑1 (MQ=255)
                                        tatGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACtg                                                >  1:3710546/1‑93 (MQ=255)
                                         atGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACtgt                                               >  2:2715248/1‑93 (MQ=255)
                                            ccGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGcc                                            >  1:6602569/1‑93 (MQ=255)
                                             cGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCca                                           >  2:6366278/1‑93 (MQ=255)
                                               aCAGGTGATGGAAATGCAGATAAACCGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCcaca                                         <  2:4912911/93‑1 (MQ=255)
                                               aCAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCcaca                                         >  2:462588/1‑93 (MQ=255)
                                                  ggTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCt                                      >  2:3657601/1‑93 (MQ=255)
                                                  ggTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCt                                      >  2:2612350/1‑93 (MQ=255)
                                                  ggTGATGGAAATGCAGATAAAACGCTAGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCt                                      <  1:4213192/93‑1 (MQ=255)
                                                   gTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCtt                                     <  1:4235359/93‑1 (MQ=255)
                                                    tGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtc                                    <  1:6159860/93‑1 (MQ=255)
                                                     gATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtct                                   <  1:4515300/93‑1 (MQ=255)
                                                     gATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtct                                   >  1:3124642/1‑93 (MQ=255)
                                                     gATGAAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtct                                   <  1:5038252/93‑1 (MQ=255)
                                                      aTGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtctc                                  >  2:5453414/1‑93 (MQ=255)
                                                      aTGAAAATTCAGAAAAACCTCTTGATTGAAAATACCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtctc                                  <  2:3693471/93‑1 (MQ=255)
                                                        ggAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGc                                >  1:3724110/1‑93 (MQ=255)
                                                          aaaTGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACTCCTTCTCGCCg                              >  2:6037240/1‑93 (MQ=255)
                                                          aaaTGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCg                              >  2:2217807/1‑93 (MQ=255)
                                                          aaaTGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCg                              <  2:1016664/93‑1 (MQ=255)
                                                            aTGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCggg                            >  2:497026/1‑93 (MQ=255)
                                                             tGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCgggg                           >  2:4564935/1‑93 (MQ=255)
                                                                aGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCg                        >  1:3505615/1‑93 (MQ=255)
                                                                  aTAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCggg                      <  1:1771941/93‑1 (MQ=255)
                                                                   aaaaaCGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCgggg                     <  1:4975390/92‑1 (MQ=255)
                                                                    aaaaCGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGa                    <  1:3956303/93‑1 (MQ=255)
                                                                    aaaaCGCTCAATTGAGAAAACCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGa                    <  1:1597276/93‑1 (MQ=255)
                                                                     aaaCGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGaa                   >  2:5866952/1‑93 (MQ=255)
                                                                      aaCGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGaaa                  >  1:5245942/1‑93 (MQ=255)
                                                                      aaCGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGaaa                  >  2:6092440/1‑93 (MQ=255)
                                                                        cGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAata                >  1:4906527/1‑93 (MQ=255)
                                                                        cGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAata                <  2:1514463/93‑1 (MQ=255)
                                                                        cGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAata                >  2:5581841/1‑93 (MQ=255)
                                                                        cGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAata                <  1:3882703/93‑1 (MQ=255)
                                                                        cGCTCGATTGAGAAAATCCCAGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAat                 <  1:6430322/93‑1 (MQ=255)
                                                                           tCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGGAATGttt             >  1:5708450/1‑93 (MQ=255)
                                                                           tCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGGAATAttt             >  1:1184119/1‑93 (MQ=255)
                                                                           tCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATAttt             <  2:6628267/93‑1 (MQ=255)
                                                                           tCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATAttt             <  2:2744896/93‑1 (MQ=255)
                                                                           tCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATAttt             <  1:5448500/93‑1 (MQ=255)
                                                                            cGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACCCCTTCTCGCCGGGGGCGGGGAAATAtttt            <  2:6722386/93‑1 (MQ=255)
                                                                             gATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGCAATATTTTg           >  1:364619/1‑93 (MQ=255)
                                                                              agtGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTgg          >  1:4449829/3‑93 (MQ=255)
                                                                              aTTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTgg          >  2:6463592/1‑93 (MQ=255)
                                                                              aTTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTgg          >  1:3833042/1‑93 (MQ=255)
                                                                               ttGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGa         <  2:5236234/93‑1 (MQ=255)
                                                                               ttGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGa         >  1:2913348/1‑93 (MQ=255)
                                                                               ttGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGa         <  2:1899270/93‑1 (MQ=255)
                                                                               ttGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGa         >  2:4888009/1‑93 (MQ=255)
                                                                                gaaaaaaaCCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGAt        <  1:4672527/89‑1 (MQ=255)
                                                                                 gagaAAATCCCGGGGGGGGCTGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTCGAta       >  2:1995809/1‑92 (MQ=255)
                                                                                 gagaAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATc       <  2:5344789/93‑1 (MQ=255)
                                                                                 gagaAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATc       <  2:4470707/93‑1 (MQ=255)
                                                                                  agaAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCa      <  2:6211255/93‑1 (MQ=255)
                                                                                  agaAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCa      >  1:830885/1‑93 (MQ=255)
                                                                                   gaAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCAt     >  2:4070917/1‑93 (MQ=255)
                                                                                     aaaTCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTc   >  2:4883022/1‑93 (MQ=255)
                                                                                     aaaTCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTc   >  2:3124647/1‑93 (MQ=255)
                                                                                      aaTCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCa  >  1:5989810/1‑93 (MQ=255)
                                                                                            |                                                                                       
CTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCA  >  NC_000913/4099211‑4099388

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: