Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	649,372	A→G	F67S (TTC→TCC)	citF ←	citrate lyase, citrate‑ACP transferase (alpha) subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	649,372	0	A	G	100.0%	273.8 / NA	76	F67S (TTC→TCC)	citF	citrate lyase, citrate‑ACP transferase (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/39);  total (37/39)

TGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGA  >  NC_000913/649241‑649500                                                                                                                                     |                                                                                                                                 tgtgTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAAc                                                                                                                           >  2:1112546/1‑139 (MQ=255)      tCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCa                                                                                                                      >  2:740452/1‑139 (MQ=255)      tCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCa                                                                                                                      >  1:447266/1‑139 (MQ=255)         aCCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGc                                                                                                                   >  2:342994/1‑139 (MQ=255)           cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                 <  1:740452/139‑1 (MQ=255)           cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                 <  1:651812/139‑1 (MQ=255)           cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                 <  2:198496/139‑1 (MQ=255)                gcgcATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTg                                                                                                            >  2:472342/1‑139 (MQ=255)                   cATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTaa                                                                                                         >  1:799028/1‑139 (MQ=255)                                  gggAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAAc                                                                                                       >  2:14836/1‑126 (MQ=255)                                  gggAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAAc                                                                                                       <  1:14836/126‑1 (MQ=255)                                     aGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCACCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                 <  2:382649/113‑1 (MQ=255)                                     aGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCACCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                 >  1:382649/1‑113 (MQ=255)                                       cTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGctt                                                                                     >  1:1624579/1‑139 (MQ=255)                                         ggACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACGGCCAC‑GGAAAGCGTGATgggagga                                                                                                                              >  1:959067/1‑96 (MQ=255)                                         ggACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATgggagga                                                                                                                              <  2:959067/96‑1 (MQ=255)                                                     cAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGgcg                                                                       <  2:447266/139‑1 (MQ=255)                                                       ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     <  1:884526/139‑1 (MQ=255)                                                       ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     <  2:799028/139‑1 (MQ=255)                                                       ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     <  1:342994/139‑1 (MQ=255)                                                       ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     <  1:345479/139‑1 (MQ=255)                                                                  cccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCg                                                                                        >  1:105030/1‑109 (MQ=255)                                                                  cccATCTTCGCGATGACGTCCATCACCATAATGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCg                                                                                        <  2:105030/109‑1 (MQ=255)                                                                  atggacgtCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGAt                                                          >  1:1651505/8‑139 (MQ=255)                                                                        ttCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcgcg                                                    <  1:587168/139‑1 (MQ=255)                                                                          cgcgATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATgggagga                                                                                                                              >  1:948214/1‑63 (MQ=255)                                                                          cgcgATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATgggagga                                                                                                                              <  2:948214/63‑1 (MQ=255)                                                                           gcgATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     <  1:261642/119‑1 (MQ=255)                                                                           gcgATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     >  2:261642/1‑119 (MQ=255)                                                                               tGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                       >  2:794476/1‑97 (MQ=255)                                                                               tGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                       <  1:794476/97‑1 (MQ=255)                                                                                 aCGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                 >  1:284768/1‑69 (MQ=255)                                                                                 aCGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                 <  2:284768/69‑1 (MQ=255)                                                                                 aCGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTa                                                                                                          <  1:414286/76‑1 (MQ=255)                                                                                 aCGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTa                                                                                                          >  2:414286/1‑76 (MQ=255)                                                                                 aCGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTa                                                                                                          <  1:262722/76‑1 (MQ=255)                                                                                 aCGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTa                                                                                                          >  2:262722/1‑76 (MQ=255)                                                                                    tCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATgggagga                                                                                                                              <  2:68690/53‑1 (MQ=255)                                                                                    tCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATgggagga                                                                                                                              >  1:68690/1‑53 (MQ=255)                                                                                     ccATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGcc                                         <  1:616106/137‑1 (MQ=255)                                                                                     ccATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGcc                                         >  2:616106/1‑137 (MQ=255)                                                                                              tatTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCttcttc                                                                                 >  1:569489/1‑88 (MQ=255)                                                                                              tatTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCttcttc                                                                                 <  2:569489/88‑1 (MQ=255)                                                                                              tatTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTg                              >  1:642627/1‑139 (MQ=255)                                                                                                ttGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCt                            <  2:1651505/139‑1 (MQ=255)                                                                                                   aCGGTCAGGTCACCGCCACGGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     <  1:1546870/96‑1 (MQ=255)                                                                                                   aCGGTCAGGTCACCGCCACGGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     >  2:1546870/1‑96 (MQ=255)                                                                                                   aCGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     <  1:1224704/95‑1 (MQ=255)                                                                                                   aCGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     >  2:1224704/1‑95 (MQ=255)                                                                                                     ggtcaggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     >  2:240634/1‑93 (MQ=255)                                                                                                     ggtcaggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     <  1:240634/93‑1 (MQ=255)                                                                                                     ggtcaggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                       <  1:472342/139‑1 (MQ=255)                                                                                                     ggtcaggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                       <  2:1624579/139‑1 (MQ=255)                                                                                                      gtcaggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                       <  1:1107765/74‑1 (MQ=255)                                                                                                      gtcaggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                       >  2:1107765/1‑74 (MQ=255)                                                                                                      gtcaggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGgcg                                                                       >  1:204739/1‑90 (MQ=255)                                                                                                      gtcaggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGgcg                                                                       <  2:204739/90‑1 (MQ=255)                                                                                                       tcaggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     >  2:1135431/1‑91 (MQ=255)                                                                                                       tcaggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     <  1:1135431/91‑1 (MQ=255)                                                                                                        caggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGct                                                                                      >  2:1208912/1‑73 (MQ=255)                                                                                                        caggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGct                                                                                      <  1:1208912/73‑1 (MQ=255)                                                                                                          ggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGcc                                                                                                                  >  2:232509/1‑43 (MQ=255)                                                                                                          ggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGcc                                                                                                                  <  1:232509/43‑1 (MQ=255)                                                                                                          ggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                       >  1:578531/1‑134 (MQ=255)                                                                                                          ggtcaCCGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                       <  2:578531/134‑1 (MQ=255)                                                                                                                cGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGct                                                                                      <  1:531087/65‑1 (MQ=255)                                                                                                                cGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGct                                                                                      >  2:531087/1‑65 (MQ=255)                                                                                                                cGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAgg                                                                             >  2:180386/1‑74 (MQ=255)                                                                                                                cGCCAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAgg                                                                             <  1:180386/74‑1 (MQ=255)                                                                                                                   cAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAgg                                                                             >  1:1551762/1‑71 (MQ=255)                                                                                                                   cAC‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAgg                                                                             <  2:1551762/71‑1 (MQ=255)                                                                                                                     c‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGct                                                                                      <  1:1352308/60‑1 (MQ=255)                                                                                                                     c‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGct                                                                                      >  2:1352308/1‑60 (MQ=255)                                                                                                                     c‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTCCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                       <  2:991623/123‑1 (MQ=255)                                                                                                                     c‑GGAAAGCGTGATGGGAGGAAACCGTCATGCCATCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                       >  1:991623/1‑123 (MQ=255)                                                                                                                           aGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                       <  2:73454/54‑1 (MQ=255)                                                                                                                           aGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                       >  1:73454/1‑54 (MQ=255)                                                                                                                           aGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGa  >  2:599335/1‑139 (MQ=255)                                                                                                                              gTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     <  2:457485/69‑1 (MQ=255)                                                                                                                              gTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                     >  1:457485/1‑69 (MQ=255)                                                                                                                              gTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa     <  2:774231/133‑1 (MQ=255)                                                                                                                              gTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa     >  1:774231/1‑133 (MQ=255)                                                                                                                                     |                                                                                                                                 TGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCAC‑GGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGA  >  NC_000913/649241‑649500

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