Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,530,797	C→T	Q285Q (CAG→CAA)	sgcX ←	putative endoglucanase with Zn‑dependent exopeptidase domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,530,797	0	C	T	100.0%	141.1 / NA	40	Q285Q (CAG→CAA)	sgcX	putative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/19);  total (21/19)

CTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACC  >  NC_000913/4530569‑4531034                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              cTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGc                                                                                                                                                                                                                         >  2:45931/1‑251 (MQ=255)         ggAAATGTGCTGCTGACATACCGTCCAGAGCGGTCAATAGGCGGATGCAATCAGACAAACCACGCAGTCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGGTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCCGGCgaggagtc                                                                                                                                                                                                                 <  2:53995/251‑3 (MQ=255)           aaaTGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGt                                                                                                                                                                                                               <  2:8414/251‑1 (MQ=255)                           gaccggCCAGAGCGGTCAAGAGGCGGATGCAATCATGCAAATAACGCAGGCTGGCACCCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTGGGGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGTAGGCGTGATTAATCCGGCCAGc                                                                                                                                                                                               <  2:37499/246‑1 (MQ=255)                                           cAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGt                                                                                                                                                                               <  1:202601/251‑1 (MQ=255)                                                        cAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGt                                                                                                                                                                  <  1:50892/251‑1 (MQ=255)                                                            cAGTCAAATCACGCAGGCTGGCAACCTCAGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGGGATc                                                                                                                                                              >  2:73162/1‑251 (MQ=255)                                                                 aaaTCACGCAGGCTGGCAACATCCGACGGAGAATGGGTATAGCGCCAAGGAATAGAGAGAATGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAAc                                                                                                                                                         >  2:402897/1‑251 (MQ=255)                                                                         aaGGCTGTCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATGGCCGGTTTCTGGTATCACGCCGGGATCCTCTTCTCGCTGCACCGGAATGTTGGGTTTAATTGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCTTCCAGCGTTCCCCGGCCATGGTAGTTCAGGCAGGGGATGCTAACGCCCTGAt                                                                                                                                                 <  2:204078/250‑1 (MQ=255)                                                                                  aaCCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTg                                                                                                                                         >  1:14841/1‑250 (MQ=255)                                                                                    ccTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGAcc                                                                                                                                      >  1:86712/1‑251 (MQ=255)                                                                                             gagaATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATATCCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTGCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGaataa                                                                                                                             <  2:321048/251‑1 (MQ=255)                                                                                                           cGACAAGGAATCGAGAGACTGGCGCAGGGAATACCATCATGCTCAACCTGGATATATCCGGTTTCTGTTATCCCGCCGGGAGTCACTTCTCGCCGCTCCGGAAGGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGTCATCGGAATAATCGTGTAACTCAgg                                                                                                               <  2:317437/251‑1 (MQ=255)                                                                                                               aaGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGgtg                                                                                                           >  1:192481/1‑251 (MQ=255)                                                                                                                 ggAATAGAGAGACTGGAGCAGGGAATACCATCCTGCTCAACCTGGATATAACAGGTTTATGTTATAACGCCGGGAGCCACTTCTCGCTGCACAGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGCTTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTc                                                                                                         >  2:387680/1‑251 (MQ=255)                                                                                                                  gaataccgagcccagcgcAGGGAATCACAACCTGCTCAACCTCGATCTCAACGGTTTCTGTTATCACGCCGGCAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCCTTTGTTCCAACATCCGTATAAGGCGAGGCGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCATGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTc                                                                                                         >  2:379700/15‑250 (MQ=255)                                                                                                                    aTAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGTAATAATCGTGTAAATCAGGGGTGTcac                                                                                                       >  2:338790/1‑250 (MQ=255)                                                                                                                         gagaCTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATg                                                                                                 >  2:298848/1‑251 (MQ=255)                                                                                                                          agaCTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATgg                                                                                                >  1:6217/1‑251 (MQ=255)                                                                                                                              tGGCGCAGGGAATACCATCCTGCTCAAACTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCACTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTACCCGACCATGGTCGTTCAGGCAGGTGATCCCAACGCCCTGATTCATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTg                                                                                            >  1:207359/1‑251 (MQ=255)                                                                                                                                gcgcAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGat                                                                                          <  1:298848/251‑1 (MQ=255)                                                                                                                                  gcAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTTTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGatat                                                                                        <  2:14841/251‑1 (MQ=255)                                                                                                                                  gcAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCGCACGATGGAGTGatat                                                                                        <  1:301053/251‑1 (MQ=255)                                                                                                                                         aTACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATAcc                                                                                 >  1:294353/1‑251 (MQ=255)                                                                                                                                                   gCTCAACCTGGATCTAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGAGTCACACGATGGAGTGATATCAATACCAATCGCcagg                                                                       >  2:360473/1‑251 (MQ=255)                                                                                                                                                   gCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCca                                                                         <  1:45931/249‑1 (MQ=255)                                                                                                                                                       aaCCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcaggtc                                                                     >  1:22949/1‑249 (MQ=255)                                                                                                                                                           tGGATATAACCGGTTTCTGTTATCAAGCCGTGACCCACTTATCGCTGCACCGGAATGTTGTGTTCAAGCGCCGATTGTTCCAACATCCGTATCAGGCGAGGAGACGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCACTACCAATCGCCAGTTCAGGttt                                                               >  1:118582/1‑251 (MQ=255)                                                                                                                                                              atatAACCGGTTTCTGTTATCACGCCGGGCGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTac                                                            >  2:77266/1‑251 (MQ=255)                                                                                                                                                                  aaCCGGTTTCTGTCATCCCGCCGGGAGCCACTTCTCGCTGCACCGGACTGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGAAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGACTAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACgg                                                        >  2:32228/1‑251 (MQ=255)                                                                                                                                                                    ccGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGGTTTACACGGCg                                                      <  1:338790/251‑1 (MQ=255)                                                                                                                                                                      ggTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTa                                                    <  1:379304/251‑1 (MQ=255)                                                                                                                                                                           ttGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGGGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCGCGGAATAATCGTGTAAAGCAGGGGTGTCACACGATGGAGTGATAGCAATACCAATCGCCAGGTCAGGTTTTACTCGGCGTAATAc                                                <  2:22949/249‑1 (MQ=255)                                                                                                                                                                               tATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGa                                            >  2:124542/1‑250 (MQ=255)                                                                                                                                                                                        gggAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGCGGCGTGATTAATCCGGCCAGCGTTCCACGACCATGGTAGTTCAGGCTGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATATTCGTGTAACTCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGGCATACGGGAGCAATGCcgc                                  >  1:397903/1‑251 (MQ=255)                                                                                                                                                                                                  tctcGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACAGCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTaa                         <  2:294353/250‑1 (MQ=255)                                                                                                                                                                                                         gCACGGGAATGTTGTGTGCAAGAGCCGTTTGTGCCAACATCCGTATCAGGCGATGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATACTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGTTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAttatt                   <  2:312093/249‑4 (MQ=255)                                                                                                                                                                                                         gCACCGGAATGTTGTGTTCAAGAGACGTTTGTGCCAACATCCTTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCGCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGTTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTcc                 <  2:283201/251‑1 (MQ=255)                                                                                                                                                                                                                  tgttgtGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGa         >  1:362271/1‑250 (MQ=255)                                                                                                                                                                                                                        gtTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGccacc  <  1:77266/251‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              CTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACC  >  NC_000913/4530569‑4531034

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