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.7% 262.7 / ‑3.9 74coding (215/1071 nt)dhaKdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (36/37);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑CCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAG  >  NC_000913/1250490‑1250757
                                                                                                                                        |                                                                                                                                      
tCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑cccc                                                                                                                                    >  1:702898/1‑139 (MQ=255)
   cttgttCCAAGGTAAGAATCTCCCCGGTTTAAATTTTGATAATCAACAGGTAACCTTCGCCCGCCATAACCTTGCAGGGGCCATTC‑AAGAATTTAACGGGCGGCCGTGAGGTGAAAATTTCGCCCGGACAGG‑ccccccg                                                                                                                                 <  2:1594133/135‑3 (MQ=255)
   cTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑ccccccg                                                                                                                                 <  1:85242/139‑3 (MQ=255)
     gTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑ccccccgaa                                                                                                                               >  1:935/1‑135 (MQ=255)
        tCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGc                                                                                                                            >  2:1348976/1‑139 (MQ=255)
              ttAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATcccc                                                                                                                      >  2:1565759/1‑139 (MQ=255)
                 aGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGa                                                                                                                   <  2:702898/139‑1 (MQ=255)
                 aGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGa                                                                                                                   <  2:1404579/139‑1 (MQ=255)
                   aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGATTATTTCGCCCGGACAGG‑CCCCCTGATAGCGTCCCCTAAcc                                                                                                                 <  2:1564612/139‑1 (MQ=255)
                   aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGAcc                                                                                                                 <  2:349779/139‑1 (MQ=255)
                   aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGAcc                                                                                                                 <  1:1277901/139‑1 (MQ=255)
                        cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGatat                                                                                                            >  1:1064208/1‑139 (MQ=255)
                        cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGatat                                                                                                            >  2:1325430/1‑139 (MQ=255)
                        cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGatat                                                                                                            <  2:1659630/139‑1 (MQ=255)
                          ccGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATaa                                                                                                          >  2:67675/1‑139 (MQ=255)
                               gtAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACcaca                                                                                                     <  2:1459091/139‑1 (MQ=255)
                                  aTTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTg                                                                                                  >  2:588190/1‑139 (MQ=255)
                                         aTAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGc                                                                                           <  1:1325430/139‑1 (MQ=255)
                                             tCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCgt                                                                                       >  2:643310/1‑139 (MQ=255)
                                                   gTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGc                                                                                 >  1:122659/1‑139 (MQ=255)
                                                      cacCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATccc                                                                                                                       <  1:47059/98‑1 (MQ=255)
                                                      cacCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATccc                                                                                                                       >  2:47059/1‑98 (MQ=255)
                                                      cacCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGAcc                                                                                                                 >  2:692435/1‑104 (MQ=255)
                                                      cacCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGAcc                                                                                                                 <  1:692435/104‑1 (MQ=255)
                                                      cacCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgc                                                                              <  2:634896/139‑1 (MQ=255)
                                                      cacCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgc                                                                              <  2:240276/139‑1 (MQ=255)
                                                           tCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa                                                                         <  1:67675/139‑1 (MQ=255)
                                                               ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                           <  1:259547/134‑1 (MQ=255)
                                                               ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                           >  2:259547/1‑134 (MQ=255)
                                                               ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAccc                                                                      >  1:1009867/1‑139 (MQ=255)
                                                               ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAccc                                                                      >  2:73785/1‑139 (MQ=255)
                                                                   cATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAca                                                                  >  2:1048816/1‑139 (MQ=255)
                                                                   cATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAca                                                                  >  2:1141087/1‑139 (MQ=255)
                                                                         cTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAggg                                                            <  1:1048816/139‑1 (MQ=255)
                                                                               tGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGc                                                                                                                            <  2:69659/69‑1 (MQ=255)
                                                                               tGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGc                                                                                                                            >  1:69659/1‑69 (MQ=255)
                                                                                  cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCTTGCCGGCTGCCGCCACCCGCCAGCGGGGGGACTTTTc                                                   >  2:1581288/1‑139 (MQ=255)
                                                                                  cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTc                                                   >  1:16477/1‑139 (MQ=255)
                                                                                   gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc                                                  >  1:1203337/1‑139 (MQ=255)
                                                                                    cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                           <  2:1149303/113‑1 (MQ=255)
                                                                                    cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                           >  1:1149303/1‑113 (MQ=255)
                                                                                    cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACt                                                       >  2:225889/1‑133 (MQ=255)
                                                                                    cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACt                                                       <  1:225889/133‑1 (MQ=255)
                                                                                      ttCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACcac                                                                                                      >  2:216454/1‑84 (MQ=255)
                                                                                      ttCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACcac                                                                                                      <  1:216454/84‑1 (MQ=255)
                                                                                               tttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGAcc                                                                                                                 <  1:247806/64‑1 (MQ=255)
                                                                                               tttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGAcc                                                                                                                 >  2:247806/1‑64 (MQ=255)
                                                                                               tttATCGGGCGTCGGTGAGGTGAAAAATTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGc                                      <  1:1565759/139‑1 (MQ=255)
                                                                                                 tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAATGCATCCCCTGACCGa                                                                                                               <  2:1105810/64‑1 (MQ=255)
                                                                                                 tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                    >  1:1391530/1‑139 (MQ=255)
                                                                                                 tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                    >  2:49006/1‑139 (MQ=255)
                                                                                                 tATCGGGCGTCGGTGAGGTGAAAATTTCGACCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGa                                                                                                               >  1:1105810/1‑64 (MQ=255)
                                                                                                 tATCGGGCGTCGGTGAGGTGAAAATTGCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGGGCATCGGCCCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                    >  2:421432/1‑139 (MQ=255)
                                                                                                     gggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGc                                <  2:1391530/139‑1 (MQ=255)
                                                                                                     gggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑ccccccg                                                                                                                                 <  2:651588/42‑3 (MQ=255)
                                                                                                     gggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑ccccccg                                                                                                                                 >  1:651588/1‑40 (MQ=255)
                                                                                                       gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                              <  2:1339224/139‑1 (MQ=255)
                                                                                                       gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                              >  1:1339224/1‑139 (MQ=255)
                                                                                                         gTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCgg                            <  1:1348976/139‑1 (MQ=255)
                                                                                                         gTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCACTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCgg                            <  2:935/139‑1 (MQ=255)
                                                                                                           cGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGt                          >  1:1293156/1‑139 (MQ=255)
                                                                                                                aggtgaAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGAACGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGCCAt                     <  2:1064208/139‑1 (MQ=255)
                                                                                                                 ggtgaAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTCTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATa                    >  1:593779/1‑139 (MQ=255)
                                                                                                                  gtgaAAATTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATac                   <  1:798782/139‑1 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGa                                                                                                                   <  1:874126/40‑1 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGa                                                                                                                   >  2:874126/1‑40 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGAcc                                                                                                                 >  2:1046598/1‑42 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGAcc                                                                                                                 <  1:1046598/42‑1 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGa                                                                                                               >  1:13174/1‑44 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGa                                                                                                               <  2:13174/44‑1 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACcac                                                                                                      >  1:225373/1‑53 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACcac                                                                                                      <  2:225373/53‑1 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                    >  1:587640/1‑119 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                    <  2:587640/119‑1 (MQ=255)
                                                                                                                      aaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg               <  2:1203337/139‑1 (MQ=255)
                                                                                                                      aaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg               <  2:122659/139‑1 (MQ=255)
                                                                                                                      aaaTTTCGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg               <  2:1293156/139‑1 (MQ=255)
                                                                                                                            cGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg               <  2:1004871/133‑1 (MQ=255)
                                                                                                                            cGCCCGGACAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg               >  1:1004871/1‑133 (MQ=255)
                                                                                                                                   aCAGGCCCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAg  <  1:421432/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCG‑CCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑CCCCC‑GAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAG  >  NC_000913/1250490‑1250757

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