Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,250,629 (C)5→6 coding (210/1071 nt) dhaK ← dihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,250,6241.C98.5% 237.4 / ‑3.8 68coding (215/1071 nt)dhaKdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (33/34);  total (34/34)
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.98e-01

CTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑CCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATG  >  NC_000913/1250493‑1250754
                                                                                                                                    |                                                                                                                                   
cTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑ccc‑cccg                                                                                                                             <  1:1384550/139‑3 (MQ=255)
               gAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAc                                                                                                               <  1:828983/139‑1 (MQ=255)
                  tatCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGa                                                                                                            <  2:1546247/139‑1 (MQ=255)
                     cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAc                                                                                                               <  1:91536/133‑1 (MQ=255)
                     cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAc                                                                                                               >  2:91536/1‑133 (MQ=255)
                     cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGatat                                                                                                         >  2:403621/1‑139 (MQ=255)
                          gtgtAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACca                                                                                                    >  2:28147/1‑139 (MQ=255)
                                                gTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGc                                                                              >  1:504603/1‑139 (MQ=255)
                                                      cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                        <  1:28147/139‑1 (MQ=255)
                                                      cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                        <  2:504603/139‑1 (MQ=255)
                                                        tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGa                                                                                                                >  2:432795/1‑97 (MQ=255)
                                                        tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGa                                                                                                                <  1:432795/97‑1 (MQ=255)
                                                         cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGTCCGGACAGGCCCC‑CCGAAAGCATccc                                                                                                                    >  1:962993/1‑92 (MQ=255)
                                                         cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGTCCGGACAGGCCCC‑CCGAAAGCATccc                                                                                                                    <  2:962993/92‑1 (MQ=255)
                                                         cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAc                                                                     >  2:985915/1‑139 (MQ=255)
                                                           ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCAt                                                                                            <  2:1112232/114‑1 (MQ=255)
                                                           ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa                                                                      <  2:738525/136‑1 (MQ=255)
                                                           ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa                                                                      >  1:738525/1‑136 (MQ=255)
                                                           ccgccACCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCAt                                                                                            >  1:1112232/1‑114 (MQ=255)
                                                               cATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAca                                                               >  1:1381829/1‑139 (MQ=255)
                                                                aTCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAcc                                                                                                              <  1:749652/91‑1 (MQ=255)
                                                                aTCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAcc                                                                                                              <  2:735803/91‑1 (MQ=255)
                                                                aTCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAcc                                                                                                              >  1:735803/1‑91 (MQ=255)
                                                                aTCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcag                                                              >  1:30380/1‑139 (MQ=255)
                                                                aTCAACTTGCATGGCGCATTCAAAGACTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAcc                                                                                                              >  2:749652/1‑91 (MQ=255)
                                                                         cATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGa                                                                                                                <  1:558874/80‑1 (MQ=255)
                                                                         cATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGa                                                                                                                >  2:558874/1‑80 (MQ=255)
                                                                           tGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAg                                                                                                                          >  2:933398/1‑68 (MQ=255)
                                                                           tGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAg                                                                                                                          <  1:933398/68‑1 (MQ=255)
                                                                           tGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑c                                                                                                                                    >  2:302906/1‑58 (MQ=255)
                                                                           tGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑c                                                                                                                                    <  1:302906/58‑1 (MQ=255)
                                                                             gcgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAg                                                                                                                          >  1:199467/1‑66 (MQ=255)
                                                                             gcgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAg                                                                                                                          <  2:199467/66‑1 (MQ=255)
                                                                              cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGa                                                                                                            <  2:811101/79‑1 (MQ=255)
                                                                              cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGa                                                                                                            >  1:811101/1‑79 (MQ=255)
                                                                              cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTc                                                >  1:1106205/1‑139 (MQ=255)
                                                                               gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGa                                                                                                                <  2:759912/74‑1 (MQ=255)
                                                                               gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGa                                                                                                                >  1:759912/1‑74 (MQ=255)
                                                                               gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc                                               <  2:803564/139‑1 (MQ=255)
                                                                                cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCt                                              <  2:1381829/139‑1 (MQ=255)
                                                                                   tCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCa                                           >  2:904030/1‑139 (MQ=255)
                                                                                   tCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCa                                           >  2:74515/1‑139 (MQ=255)
                                                                                   tCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑c                                                                                                                                    <  1:1311605/50‑1 (MQ=255)
                                                                                   tCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑c                                                                                                                                    >  2:1311605/1‑50 (MQ=255)
                                                                                      aaGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACa                                        <  1:904030/139‑1 (MQ=255)
                                                                                       aGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATccc                                                                                                                    <  1:444505/62‑1 (MQ=255)
                                                                                       aGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATccc                                                                                                                    >  2:444505/1‑62 (MQ=255)
                                                                                           tttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccg                                                                         <  1:768439/101‑1 (MQ=255)
                                                                                           tttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccg                                                                         >  2:768439/1‑101 (MQ=255)
                                                                                             tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCAt                                                                                            <  2:1007866/80‑1 (MQ=255)
                                                                                             tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCAt                                                                                            >  1:1007866/1‑80 (MQ=255)
                                                                                             tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑ccc‑cccg                                                                                                                             >  2:1590204/1‑44 (MQ=255)
                                                                                             tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑ccc‑cccg                                                                                                                             <  1:1590204/46‑3 (MQ=255)
                                                                                                   gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATGGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                           <  2:625383/139‑1 (MQ=255)
                                                                                                   gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                        <  1:534854/94‑1 (MQ=255)
                                                                                                   gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                        >  2:534854/1‑94 (MQ=255)
                                                                                                   gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                           >  1:353913/1‑139 (MQ=255)
                                                                                                   gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                           >  2:592219/1‑139 (MQ=255)
                                                                                                    cGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc                                               <  2:469583/118‑1 (MQ=255)
                                                                                                    cGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc                                               >  1:469583/1‑118 (MQ=255)
                                                                                                     gTCGGTGAGGTGAAATTTTCGCCCGGACAGG‑CCCACCGAAAGCATCCCCTGATCGa                                                                                                            >  2:208873/1‑56 (MQ=255)
                                                                                                     gTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGa                                                                                                            <  1:208873/56‑1 (MQ=255)
                                                                                                      tCGGTGAGGTGAAAATTTCGCCTGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                           >  1:988204/1‑136 (MQ=255)
                                                                                                      tCGGTGAGGTGAAAATTTCGCCTGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                           <  2:988204/136‑1 (MQ=255)
                                                                                                      tCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCggg                        <  2:421638/139‑1 (MQ=255)
                                                                                                      tCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCggg                        <  2:353913/139‑1 (MQ=255)
                                                                                                          tgaggtgaAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGAc                    <  2:1106205/139‑1 (MQ=255)
                                                                                                             ggtgaAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGc                             >  1:332568/1‑127 (MQ=255)
                                                                                                             ggtgaAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGc                             <  2:332568/127‑1 (MQ=255)
                                                                                                              gtgaAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATac                >  1:1468170/1‑139 (MQ=255)
                                                                                                                  aaaTTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg            <  2:1468170/139‑1 (MQ=255)
                                                                                                                        cGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAg                                       >  2:1158709/1‑106 (MQ=255)
                                                                                                                        cGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAg                                       <  1:1158709/106‑1 (MQ=255)
                                                                                                                         gCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTg     >  2:1119332/1‑139 (MQ=255)
                                                                                                                            cGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATg  <  1:403621/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                   
CTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑CCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATG  >  NC_000913/1250493‑1250754

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑