Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,240,295	G→A	intergenic (‑240/‑30)	yjbT ← / → psiE	PF17089 family protein YjbT/putative phosphate starvation‑inducible protein PsiE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,240,295	0	G	A	97.6%	129.1 / NA	42	intergenic (‑240/‑30)	yjbT/psiE	PF17089 family protein YjbT/putative phosphate starvation‑inducible protein PsiE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/21);  minor base C (1/0);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTACTTTCTCTATTTCGAATTTATCGC  >  NC_000913/4240053‑4240539                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     cATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGAGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGc                                                                                                                                                                                                                                              >  1:187061/1‑251 (MQ=255)      gTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAACTCGCCACtg                                                                                                                                                                                                                                         >  2:54870/1‑251 (MQ=255)       ttCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGCATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGGATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACtgt                                                                                                                                                                                                                                        >  2:149779/1‑251 (MQ=255)        tCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACtgtg                                                                                                                                                                                                                                       >  2:258869/1‑251 (MQ=255)        tCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACtgt                                                                                                                                                                                                                                        <  2:59870/250‑1 (MQ=255)                                     aTTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGGTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCActctct                                                                                                                                                                                                          <  2:21176/251‑1 (MQ=255)                                     aTTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGATTTGGTGCGGTTTGTCTTGCCGGACGCGCCGCCAGGACCGGCTTATCCGGCCTACGGGCCGGTATATCCGGCTTTGCAGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑CTATAGATCTCCGTCACATTTTTGCGTTATCCAGCAAGCTCGCCACTGTGAAGGAGGTACTGATATGACGTCActctct                                                                                                                                                                                                          >  1:166949/1‑251 (MQ=255)                                     aTTGCTGGAGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAAATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGATATGTCATCActttct                                                                                                                                                                                                          >  2:299533/1‑251 (MQ=255)                                        gCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGTGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAAGCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGt                                                                                                                                                                                                       <  2:5185/251‑1 (MQ=255)                                         cTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGGGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCAATCTCTCGTc                                                                                                                                                                                                      >  2:306387/1‑251 (MQ=255)                                           ggtggCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTcc                                                                                                                                                                                                     <  2:10215/250‑1 (MQ=255)                                            gtggCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCcgc                                                                                                                                                                                                   <  1:229245/251‑1 (MQ=255)                                                  gggAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGt                                                                                                                                                                                               <  1:258869/249‑1 (MQ=255)                                                                 gggTGATTTTTGAATTTCCGGTTGTGGTGCGGTGTGTCTGTCCGGATGCGCCGCCAGGTGTGGCTGATCCGGCCTACGGGTAGGTATATCAGGTTTTGGTGGAGGCGCGCTCCAAAGCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGGCACATTTGTGCGTTATTCAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCAATCACACGTCCGCTCGTGGAGCTTAGCTccacca                                                                                                                                                                              <  2:83402/251‑1 (MQ=255)                                                                  ggTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCAt                                                                                                                                                                             >  2:86518/1‑251 (MQ=255)                                                                      aTTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGATGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTg                                                                                                                                                                         >  1:318258/1‑251 (MQ=255)                                                                         tttGAATTACCGGCTTTTGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGAGGCAGCGCGCTCCAAATCTAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACACTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCTCTCTCTCGTCCGCGCGTGTAGGTTATCTCCACCATTTTGCa                                                                                                                                                                       <  1:226657/250‑1 (MQ=255)                                                                         tttGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCCAATCAAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTAGTCCGCGCGTGGAGTTTATCTCCACCATTTTGc                                                                                                                                                                        >  2:104484/1‑249 (MQ=255)                                                                            gAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTTAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGAcc                                                                                                                                                                   <  1:81357/251‑1 (MQ=255)                                                                                  ccGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGAGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGCGGAGTTTATCTCCACCATTTTGCAGACCGTACTc                                                                                                                                                             >  2:3094/1‑251 (MQ=255)                                                                                         tGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCtt                                                                                                                                                       <  1:306387/250‑1 (MQ=255)                                                                                         tGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCGTATCCGGCCTACTGGTAGGTATATCCGGCTTTGGTGGAGGCTCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTTCTATGACGTCACTCTCTCGGCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAAGCtt                                                                                                                                                       <  1:241193/250‑1 (MQ=255)                                                                                           gTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAACTCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGc                                                                                                                                                    >  1:204613/1‑251 (MQ=255)                                                                                                            ggATGCGGCGCGAGACGCGGCTAAACCGCCCTACGGGTAGGAATATCTGGCTCTGGAGGACGCGCGCTCAAAATCCAGGCTGCACAAAAAAT‑CAAAAAAACATATAGATCTCCGTCACATTTTTGCGTTATACAGACAAAACGCCACTGTAAAGGAGGTACTGCTATGAAGCCACTCTCACGACCGCGCGTGGCGTTTATCTCCACCATTTTGCAGACCGTACTCACTCTTGGCCTGCTGTGCCTCGGcat                                                                                                                                   >  2:359183/1‑249 (MQ=255)                                                                                                                    cgccAGGTGCGGCTTATCCGGCGTACGGGTAGGTATATTCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑AGATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCCCTCTCTCGTCCGCGCGTGGAGTTTATCTCCCCCATTTTGAAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTgg                                                                                                                           <  1:215094/251‑1 (MQ=255)                                                                                                                                     ccGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGTCCTGCTGTGCCTCGGCCTGATTTTGGTTGGCTTCCTCGGCaaa                                                                                                          <  2:141402/251‑1 (MQ=255)                                                                                                                                                                          cgcgcTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTTCTGCTATGACGTCACTCTCTCGTTCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAg                                                                     <  1:393833/251‑1 (MQ=255)                                                                                                                                                                                                cAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCt                                               <  1:283639/251‑1 (MQ=255)                                                                                                                                                                                                 aaaaCATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTg                                              >  1:240822/1‑251 (MQ=255)                                                                                                                                                                                                 aaaaCATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTg                                              >  1:302004/1‑251 (MQ=255)                                                                                                                                                                                                   aaCATACACAAAAA‑ATATAGATCTCCGTAACATTTTTGCGTTAGACAGAAAGCGCGCCACGGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTAAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTACTCGGCAAAGAAACGGTGCATCTGGCGGATGTGCTGGTCGCGACAGAACAATCCAGCAAATATGAGCttgt                                            <  2:320661/251‑4 (MQ=255)                                                                                                                                                                                                           caAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGCCCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGAc                                    >  1:402798/1‑251 (MQ=255)                                                                                                                                                                                                               aa‑aTATAGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACtggt                                >  2:114827/1‑251 (MQ=255)                                                                                                                                                                                                                    ataGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGtt                                                                             >  2:253707/1‑202 (MQ=255)                                                                                                                                                                                                                    ataGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGtt                                                                             <  1:253707/202‑1 (MQ=255)                                                                                                                                                                                                                    ataGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGAACGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGtt                            <  1:118737/251‑1 (MQ=255)                                                                                                                                                                                                                      aGATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCGCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGttt                           <  2:318258/250‑1 (MQ=255)                                                                                                                                                                                                                       gATCTCCGTCACATTTTTGCGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTAc                         >  2:178289/1‑251 (MQ=255)                                                                                                                                                                                                                           tcCGTCACAGTTTTGCGTTATACAGAAAGCGCGCCACGGTGAAGGATGTACTTCTATGACGTCACTCTCTCGGCTGCGCGGGGAGTGTATCTCGACCATTTTGCAGACCGGACGCAAGCTTGGCCTGCTGTGCCGCGGCCTGATTTTGGTTGGCTTCCGCGGCAAAGAAACGGTGCATGTGTCTGATGTGTTGTGCGCGCCAGAACAAACCAGCAAATATGAGCTGGtgtaatgtctgtttgtttacttt                      <  2:31213/250‑23 (MQ=255)                                                                                                                                                                                                                                          cGTTATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTACTTTCTCTATTTCGAAttt       >  2:302661/1‑250 (MQ=255)                                                                                                                                                                                                                                            ttATACAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGGGGTTTACTTTCTCTATTTCGAATTTATc    >  2:99908/1‑251 (MQ=255)                                                                                                                                                                                                                                              ataCAGAAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTACTTTCTCTAGTTCGAATTTATcgc  <  1:72950/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     CATTTGTTCCAGAGGAACATTGTCGATTTTTCGCGCATTGCTGGTGGCTGGGAATCACCTGAATGGGTGATTTTTGAATTACCGGCTTTGGTGCGGTTTGTCTTGCCGGATGCGCCGCCAGGCGCGGCTTATCCGGCCTACGGGTAGGTATATCCGGCTTTGGTGGAGGCGCGCTCCAAATCCAGGTTGAACAAAACATACACAAAAA‑ATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTACTTTCTCTATTTCGAATTTATCGC  >  NC_000913/4240053‑4240539

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