Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	648,555	A→G	R339R (CGT→CGC)	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	648,555	0	A	G	100.0%	220.8 / NA	62	R339R (CGT→CGC)	citF	citrate lyase, citrate‑ACP transferase (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/39);  total (23/39)

GTCAGAGCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTACGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGT  >  NC_000913/648411‑648700                                                                                                                                                 |                                                                                                                                                  gTCAGAGCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGcc                                                                                                                                             >  1:794398/1‑151 (MQ=255) gTCAGAGCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGcc                                                                                                                                             <  2:283699/151‑1 (MQ=255)       gCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGAt                                                                                                                                       >  2:59819/1‑151 (MQ=255)       gCCGGTCAGCACGTTAACGTTGAACTGGGAGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGAt                                                                                                                                       >  1:84894/1‑151 (MQ=255)           gTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTgcg                                                                                                                                   >  2:521302/1‑151 (MQ=255)                cACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGTGGGCCAGCGATTGCgcagc                                                                                                                               >  2:27549/1‑150 (MQ=255)                cACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCgcagca                                                                                                                              <  2:756879/151‑1 (MQ=255)                       aCGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGt                                                                                                                       <  2:970077/151‑1 (MQ=255)                            gAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCaaa                                                                                                                   <  2:436696/150‑1 (MQ=255)                               cTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGctc                                                                                                               >  2:518774/1‑151 (MQ=255)                                tGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGctct                                                                                                              >  1:919712/1‑151 (MQ=255)                                 gggTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTg                                                                                                             <  2:653095/151‑1 (MQ=255)                                  ggTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGc                                                                                                            <  2:965961/151‑1 (MQ=255)                                     gtCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGcaca                                                                                                         <  1:326873/151‑1 (MQ=255)                                      tCAATTTCCAGCGCGCTCAGTACCACCACGTCTAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGACGTGATTGTGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGcaca                                                                                                         <  2:469205/150‑1 (MQ=255)                                         aTTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCca                                                                                                     <  1:980150/151‑1 (MQ=255)                                         aTTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCca                                                                                                     <  1:415200/151‑1 (MQ=255)                                                gcgcgcTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGttt                                                                                              >  1:1002573/1‑151 (MQ=255)                                                 cgcgcTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTg                                                                                             <  1:989382/151‑1 (MQ=255)                                                  gcgcTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGc                                                                                            >  1:846162/1‑151 (MQ=255)                                                   cgcTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCg                                                                                           <  2:530211/151‑1 (MQ=255)                                                      tCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGAt                                                                                        >  1:318989/1‑151 (MQ=255)                                                       cAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATc                                                                                       <  2:249061/151‑1 (MQ=255)                                                       cAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATc                                                                                       <  2:50404/151‑1 (MQ=255)                                                         gTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATc                                                                                       >  1:251023/1‑149 (MQ=255)                                                         gTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATc                                                                                       >  1:367836/1‑149 (MQ=255)                                                         gTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCa                                                                                      <  2:761731/150‑1 (MQ=255)                                                           accaccACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGa                                                                                    <  2:990898/150‑1 (MQ=255)                                                              accacGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCtt                                                                                <  2:1057741/151‑1 (MQ=255)                                                               ccacGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCttt                                                                               <  2:153336/151‑1 (MQ=255)                                                                cacGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCtttt                                                                              <  2:419237/151‑1 (MQ=255)                                                                cacGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCtttt                                                                              <  2:509013/151‑1 (MQ=255)                                                                cacGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCtttt                                                                              <  1:325612/151‑1 (MQ=255)                                                                     cGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTg                                                                         >  1:362843/1‑151 (MQ=255)                                                                      gagaCGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGc                                                                        >  2:759405/1‑151 (MQ=255)                                                                        gaCGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAg                                                                      <  2:615749/151‑1 (MQ=255)                                                                          cGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGt                                                                    <  1:8391/151‑1 (MQ=255)                                                                            aTCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTc                                                                   >  1:218477/1‑150 (MQ=255)                                                                                 ccGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAt                                                             <  2:856743/151‑1 (MQ=255)                                                                                   gATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtcg                                                           <  2:685706/151‑1 (MQ=255)                                                                                      gcgcCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtcgtcg                                                        <  1:996060/151‑1 (MQ=255)                                                                                            ttCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTa                                                   >  1:441944/1‑150 (MQ=255)                                                                                             tCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAAt                                                 >  1:520896/1‑151 (MQ=255)                                                                                              cGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAAt                                                 <  1:505358/150‑1 (MQ=255)                                                                                                 aCCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCg                                             >  1:318975/1‑151 (MQ=255)                                                                                                     cAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCa                                          >  2:424201/1‑150 (MQ=255)                                                                                                          aGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCg                                    <  1:1025228/151‑1 (MQ=255)                                                                                                            cGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGaa                                  <  2:643515/151‑1 (MQ=255)                                                                                                            cGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGaa                                  <  2:830523/151‑1 (MQ=255)                                                                                                            cGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGaa                                  <  1:381738/151‑1 (MQ=255)                                                                                                                  ggTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCg                              <  2:245633/149‑1 (MQ=255)                                                                                                                  ggTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGgc                            <  1:365721/151‑1 (MQ=255)                                                                                                                   gTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGgcg                           >  2:39344/1‑151 (MQ=255)                                                                                                                     tGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGgcgc                          <  1:153309/150‑1 (MQ=255)                                                                                                                      ggCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGgcgc                          <  2:575105/149‑1 (MQ=255)                                                                                                                        cgcTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAat                      >  1:638035/1‑151 (MQ=255)                                                                                                                         gcTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGTTCGAAGTCGGCGCGAata                     <  2:103886/151‑1 (MQ=255)                                                                                                                            gATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAatat                    <  1:778832/149‑1 (MQ=255)                                                                                                                            gATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATc                   <  1:722555/150‑1 (MQ=255)                                                                                                                                          tGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCAtttt    >  2:1048512/1‑151 (MQ=255)                                                                                                                                           ggggTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCAtttt    <  1:147936/150‑1 (MQ=255)                                                                                                                                            gggTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGt  >  1:495464/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GTCAGAGCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTACGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGT  >  NC_000913/648411‑648700

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