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

TCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAG  >  NC_000913/1250490‑1250757
                                                                                                                                        |                                                                                                                                      
tCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGcccc                                                                                                                                    >  1:702898/1‑139 (MQ=255)
   cttgttCCAAGGTAAGAATCTCCCCGGTTTAAATTTTGATAATCAACAGGTAACCTTCGCCCGCCATAACCTTGCAGGGGCCATTCAAGAATTTAACGGGCGGCCGTGAGGTGAAAATTTCGCCCGGACAGGccccccg                                                                                                                                 <  2:1594133/135‑3 (MQ=255)
   cTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGccccccg                                                                                                                                 <  1:85242/139‑3 (MQ=255)
     gTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGccccccgaa                                                                                                                               >  1:935/1‑135 (MQ=255)
        tCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGc                                                                                                                            >  2:1348976/1‑139 (MQ=255)
              ttAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATcccc                                                                                                                      >  2:1565759/1‑139 (MQ=255)
                 aGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGa                                                                                                                   <  2:702898/139‑1 (MQ=255)
                 aGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGa                                                                                                                   <  2:1404579/139‑1 (MQ=255)
                   aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGATTATTTCGCCCGGACAGGCCCCCTGATAGCGTCCCCTAAcc                                                                                                                 <  2:1564612/139‑1 (MQ=255)
                   aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                                 <  2:349779/139‑1 (MQ=255)
                   aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                                 <  1:1277901/139‑1 (MQ=255)
                        cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGatat                                                                                                            >  1:1064208/1‑139 (MQ=255)
                        cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGatat                                                                                                            >  2:1325430/1‑139 (MQ=255)
                        cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGatat                                                                                                            <  2:1659630/139‑1 (MQ=255)
                          ccGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATaa                                                                                                          >  2:67675/1‑139 (MQ=255)
                               gtAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACcaca                                                                                                     <  2:1459091/139‑1 (MQ=255)
                                  aTTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTg                                                                                                  >  2:588190/1‑139 (MQ=255)
                                         aTAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGc                                                                                           <  1:1325430/139‑1 (MQ=255)
                                             tCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCgt                                                                                       >  2:643310/1‑139 (MQ=255)
                                                   gTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGc                                                                                 >  1:122659/1‑139 (MQ=255)
                                                      cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATccc                                                                                                                       <  1:47059/98‑1 (MQ=255)
                                                      cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATccc                                                                                                                       >  2:47059/1‑98 (MQ=255)
                                                      cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                                 >  2:692435/1‑104 (MQ=255)
                                                      cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                                 <  1:692435/104‑1 (MQ=255)
                                                      cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgc                                                                              <  2:634896/139‑1 (MQ=255)
                                                      cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgc                                                                              <  2:240276/139‑1 (MQ=255)
                                                           tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa                                                                         <  1:67675/139‑1 (MQ=255)
                                                               ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                           <  1:259547/134‑1 (MQ=255)
                                                               ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                           >  2:259547/1‑134 (MQ=255)
                                                               ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAccc                                                                      >  1:1009867/1‑139 (MQ=255)
                                                               ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAccc                                                                      >  2:73785/1‑139 (MQ=255)
                                                                   cATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAca                                                                  >  2:1048816/1‑139 (MQ=255)
                                                                   cATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAca                                                                  >  2:1141087/1‑139 (MQ=255)
                                                                         cTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAggg                                                            <  1:1048816/139‑1 (MQ=255)
                                                                               tGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGc                                                                                                                            <  2:69659/69‑1 (MQ=255)
                                                                               tGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGc                                                                                                                            >  1:69659/1‑69 (MQ=255)
                                                                                  cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCTTGCCGGCTGCCGCCACCCGCCAGCGGGGGGACTTTTc                                                   >  2:1581288/1‑139 (MQ=255)
                                                                                  cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTc                                                   >  1:16477/1‑139 (MQ=255)
                                                                                   gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc                                                  >  1:1203337/1‑139 (MQ=255)
                                                                                    cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                           <  2:1149303/113‑1 (MQ=255)
                                                                                    cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                           >  1:1149303/1‑113 (MQ=255)
                                                                                    cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACt                                                       >  2:225889/1‑133 (MQ=255)
                                                                                    cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACt                                                       <  1:225889/133‑1 (MQ=255)
                                                                                      ttCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACcac                                                                                                      >  2:216454/1‑84 (MQ=255)
                                                                                      ttCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACcac                                                                                                      <  1:216454/84‑1 (MQ=255)
                                                                                               tttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                                 <  1:247806/64‑1 (MQ=255)
                                                                                               tttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                                 >  2:247806/1‑64 (MQ=255)
                                                                                               tttATCGGGCGTCGGTGAGGTGAAAAATTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGc                                      <  1:1565759/139‑1 (MQ=255)
                                                                                                 tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAATGCATCCCCTGACCGa                                                                                                               <  2:1105810/64‑1 (MQ=255)
                                                                                                 tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                    >  1:1391530/1‑139 (MQ=255)
                                                                                                 tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                    >  2:49006/1‑139 (MQ=255)
                                                                                                 tATCGGGCGTCGGTGAGGTGAAAATTTCGACCGGACAGGCCCCCCGAAAGCATCCCCTGACCGa                                                                                                               >  1:1105810/1‑64 (MQ=255)
                                                                                                 tATCGGGCGTCGGTGAGGTGAAAATTGCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGGGCATCGGCCCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                    >  2:421432/1‑139 (MQ=255)
                                                                                                     gggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGc                                <  2:1391530/139‑1 (MQ=255)
                                                                                                     gggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGccccccg                                                                                                                                 <  2:651588/42‑3 (MQ=255)
                                                                                                     gggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGccccccg                                                                                                                                 >  1:651588/1‑40 (MQ=255)
                                                                                                       gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                              <  2:1339224/139‑1 (MQ=255)
                                                                                                       gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                              >  1:1339224/1‑139 (MQ=255)
                                                                                                         gTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCgg                            <  1:1348976/139‑1 (MQ=255)
                                                                                                         gTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCACTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCgg                            <  2:935/139‑1 (MQ=255)
                                                                                                           cGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGt                          >  1:1293156/1‑139 (MQ=255)
                                                                                                                aggtgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAACGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGCCAt                     <  2:1064208/139‑1 (MQ=255)
                                                                                                                 ggtgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTCTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATa                    >  1:593779/1‑139 (MQ=255)
                                                                                                                  gtgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATac                   <  1:798782/139‑1 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGa                                                                                                                   <  1:874126/40‑1 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGa                                                                                                                   >  2:874126/1‑40 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                                 >  2:1046598/1‑42 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                                 <  1:1046598/42‑1 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGa                                                                                                               >  1:13174/1‑44 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGa                                                                                                               <  2:13174/44‑1 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACcac                                                                                                      >  1:225373/1‑53 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACcac                                                                                                      <  2:225373/53‑1 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                    >  1:587640/1‑119 (MQ=255)
                                                                                                                     aaaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                    <  2:587640/119‑1 (MQ=255)
                                                                                                                      aaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg               <  2:1203337/139‑1 (MQ=255)
                                                                                                                      aaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg               <  2:122659/139‑1 (MQ=255)
                                                                                                                      aaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg               <  2:1293156/139‑1 (MQ=255)
                                                                                                                            cGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg               <  2:1004871/133‑1 (MQ=255)
                                                                                                                            cGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg               >  1:1004871/1‑133 (MQ=255)
                                                                                                                                   aCAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAg  <  1:421432/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAG  >  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: