Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,351,093	G→A	G190G (GGC→GGT)	dcuS ←	sensory histidine kinase in two‑component regulatory system with DcuR, regulator of anaerobic fumarate respiration

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,351,093	0	G	A	100.0%	277.8 / NA	77	G190G (GGC→GGT)	dcuS	sensory histidine kinase in two‑component regulatory system with DcuR, regulator of anaerobic fumarate respiration
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (41/36);  total (41/36)

ACGCCTTCTTTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGG  >  NC_000913/4350959‑4351214                                                                                                                                       |                                                                                                                          acgCCTTCTTTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCaa                                                                                                                       <  2:125779/139‑1 (MQ=255)           tGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGa                                                                                                             <  1:166075/139‑1 (MQ=255)                 aCTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGAt                                                                                                       >  1:951292/1‑139 (MQ=255)                 aCTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGAt                                                                                                       >  1:332827/1‑139 (MQ=255)                 aCTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAAGACCGACCAGAt                                                                                                       >  2:857469/1‑139 (MQ=255)                  cTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATa                                                                                                      >  1:1097504/1‑139 (MQ=255)                   tGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATaa                                                                                                     >  1:1207657/1‑139 (MQ=255)                   tGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATaa                                                                                                     >  1:226970/1‑139 (MQ=255)                   tGCAACATGGCCTGGAGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATaa                                                                                                     >  1:56498/1‑139 (MQ=255)                     cAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATg                                                                                                   <  1:492825/139‑1 (MQ=255)                           ggCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATAcc                                                                                                               <  1:1008524/121‑1 (MQ=255)                           ggCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATAcc                                                                                                               >  2:1008524/1‑121 (MQ=255)                           ggCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAg                                                                                             >  1:470670/1‑139 (MQ=255)                            gCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGc                                                                                            >  2:664752/1‑139 (MQ=255)                             ccTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGGATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCg                                                                                           >  1:943204/1‑139 (MQ=255)                              cTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGa                                                                                          <  2:647732/139‑1 (MQ=255)                               tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                         <  1:353124/139‑1 (MQ=255)                                      gCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATa                                                                                                      <  1:292838/119‑1 (MQ=255)                                      gCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATa                                                                                                      >  2:292838/1‑119 (MQ=255)                                       cTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCTGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAcc                                                                                                           <  1:346530/113‑1 (MQ=255)                                       cTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAcc                                                                                                           >  2:346530/1‑113 (MQ=255)                                       cTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGa                                                                                          >  2:8575/1‑130 (MQ=255)                                       cTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGa                                                                                          <  1:8575/130‑1 (MQ=255)                                       cTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTg                                                                                 >  2:671747/1‑139 (MQ=255)                                        tCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGa                                                                                <  2:951292/139‑1 (MQ=255)                                                cGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTgg                                                                        >  1:386162/1‑139 (MQ=255)                                                cGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTgg                                                                        >  2:492971/1‑139 (MQ=255)                                                 gTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCa                                                                                    <  1:114531/126‑1 (MQ=255)                                                 gTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCa                                                                                    >  2:114531/1‑126 (MQ=255)                                                 gTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTggg                                                                       >  2:455983/1‑139 (MQ=255)                                                    gagaTTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTca                                                                    >  1:957200/1‑139 (MQ=255)                                                      gaTTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTcaca                                                                  <  2:1097504/139‑1 (MQ=255)                                                         ttCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACgg                                                               >  1:340370/1‑139 (MQ=255)                                                           cGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCt                                                             <  2:1207657/139‑1 (MQ=255)                                                              aGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGa                                                                                          <  2:625522/107‑1 (MQ=255)                                                              aGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGa                                                                                          >  1:625522/1‑107 (MQ=255)                                                                     cAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAgg                                                   >  2:1074191/1‑139 (MQ=255)                                                                      aGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCt                                                             >  1:189980/1‑128 (MQ=255)                                                                      aGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCt                                                             <  2:189980/128‑1 (MQ=255)                                                                          cGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGAt                                              >  2:250332/1‑139 (MQ=255)                                                                            aaaaGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGa                                                                                                             >  2:254227/1‑74 (MQ=255)                                                                            aaaaGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGa                                                                                                             <  1:254227/74‑1 (MQ=255)                                                                            aaaaGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCg                                            >  2:323726/1‑139 (MQ=255)                                                                             aaaGGATTTTTTTCAGTACCTTAACCAGAATGGAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGc                                           <  2:1292747/139‑1 (MQ=255)                                                                               aGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGcca                                         <  1:1120884/139‑1 (MQ=255)                                                                               aGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGcca                                         <  1:664752/139‑1 (MQ=255)                                                                                   tttttttCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCcaccac                                     <  2:1117587/139‑1 (MQ=255)                                                                                   tttttttCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCcaccac                                     >  2:66272/1‑139 (MQ=255)                                                                                                  tAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATg                      >  2:859686/1‑139 (MQ=255)                                                                                                      cAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGAttt                  <  1:323726/139‑1 (MQ=255)                                                                                                       aGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGAtttt                 >  1:1219032/1‑139 (MQ=255)                                                                                                       aGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGAtttt                 <  2:957200/139‑1 (MQ=255)                                                                                                       aGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGAtttt                 >  1:1282074/1‑139 (MQ=255)                                                                                                        gAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTc                <  2:417635/139‑1 (MQ=255)                                                                                                         aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCt                                                             <  2:1204408/93‑1 (MQ=255)                                                                                                         aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCt                                                             >  1:1204408/1‑93 (MQ=255)                                                                                                         aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCa               <  2:973926/139‑1 (MQ=255)                                                                                                         aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCa               >  1:1305392/1‑139 (MQ=255)                                                                                                          aTGCAGGTGCCAATCAGTCCGGCCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCaa                                <  1:613888/121‑1 (MQ=255)                                                                                                          aTGCAGGTGCCAATCAGTCCGGCCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCaa                                >  2:613888/1‑121 (MQ=255)                                                                                                           tGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATc             <  2:332827/139‑1 (MQ=255)                                                                                                                gTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGCCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGa        <  1:859686/139‑1 (MQ=255)                                                                                                                 tGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGAt       <  1:671747/139‑1 (MQ=255)                                                                                                                 tGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGAt       <  1:1074191/139‑1 (MQ=255)                                                                                                                 tGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGAt       <  2:226970/139‑1 (MQ=255)                                                                                                                   ccAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCg                                                                                           <  1:80185/53‑1 (MQ=255)                                                                                                                   ccAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCg                                                                                           >  2:80185/1‑53 (MQ=255)                                                                                                                   ccAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATc                                                                              >  1:673765/1‑66 (MQ=255)                                                                                                                   ccAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATc                                                                              <  2:673765/66‑1 (MQ=255)                                                                                                                   ccAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATgg     >  1:864999/1‑139 (MQ=255)                                                                                                                   ccAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATgg     >  2:841776/1‑139 (MQ=255)                                                                                                                   ccAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATgg     <  1:625456/139‑1 (MQ=255)                                                                                                                     aaTCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACg                                    >  2:574777/1‑106 (MQ=255)                                                                                                                     aaTCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACg                                    <  1:574777/106‑1 (MQ=255)                                                                                                                      aTCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATggggg  >  2:1005373/1‑139 (MQ=255)                                                                                                                          gTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGAt       >  2:1236824/1‑130 (MQ=255)                                                                                                                          gTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGAt       <  1:1236824/130‑1 (MQ=255)                                                                                                                                       |                                                                                                                          ACGCCTTCTTTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGG  >  NC_000913/4350959‑4351214

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