Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	518,281	G→A	intergenic (‑2/+59)	cnoX ← / ← ybbO	chaperedoxin/NADP(+)‑dependent aldehyde reductase YbbO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	518,281	0	G	A	97.6%	129.7 / ‑4.3	41	intergenic (‑2/+59)	cnoX/ybbO	chaperedoxin/NADP(+)‑dependent aldehyde reductase YbbO
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/19);  minor base . (1/0);  total (22/19)
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.81e-01

ATACACGGTCGGAATCGCACGCAGACCAAACTGCGCGGCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGAT  >  NC_000913/518034‑518510                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       aTACACGGTCGGAATCGCACGCAGACCAAACTGCGCGGCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑t                                                                                                                                                                                                                                    >  1:349523/1‑251 (MQ=255)   acacGGTCGGAATCGCACGCAGACCAAACTGCGCGGCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCg                                                                                                                                                                                                                                  >  2:317603/1‑251 (MQ=255)   acacGGGCGGAATCGCACGCAGACCAAACTGCGCGGCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑Tc                                                                                                                                                                                                                                   <  1:366418/250‑1 (MQ=255)           ggAATCGCACGCAGACCAAACTGCGCGGCACTCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCtgt                                                                                                                                                                                                                          >  1:167124/1‑251 (MQ=255)            gAATCGCACGCAGACCAAACTGCGCGGCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCtgtt                                                                                                                                                                                                                         <  1:2788/251‑1 (MQ=255)            gAATCGCACGCAGACCAAACTGCGCGGCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCtgtt                                                                                                                                                                                                                         >  1:327038/1‑251 (MQ=255)                      caGACCAAACTGCGCGGCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATAttt                                                                                                                                                                                                               <  1:390485/251‑1 (MQ=255)                                     ggCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATggggg                                                                                                                                                                                                <  1:5876/251‑1 (MQ=255)                                                       gTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGAAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGATATGGGGGCTTAAGCGCGCGCTTCa                                                                                                                                                                               >  2:40828/1‑250 (MQ=255)                                                                        ccAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATAtt                                                                                                                                                             <  1:317603/251‑1 (MQ=255)                                                                                       cGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACgcgc                                                                                                                                              >  2:32334/1‑251 (MQ=255)                                                                                                       gAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCt                                                                                                                              >  1:435369/1‑251 (MQ=255)                                                                                                            tttCCAGAATTGGAGTTAACTGCCAACAGTGCTGGTTACGTTCAGCCCAAAAATAGAACAGCACCGGAGTGCTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTTCGGACGTGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGTCTTAAGCGCGCTCTTCAACTCATCCCTGCAATATTTTGTCCATCACGCGCCCTGTCGGCAGGCGCTTaat                                                                                                                          <  2:384271/250‑2 (MQ=255)                                                                                                             ttCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAg                                                                                                                          <  2:215621/249‑1 (MQ=255)                                                                                                                    aTTGGGGTTACCTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGCGTTCATCGACTGTTCAAGAACCTGATGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑CCGCTCTGTGTTGTCGCTATTTTCTTTGACATGGGTGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATCTTTTGTCCATCACTCGCCCCGCCAGCAGGCGCTTAAGCTCCATTAc                                                                                                                 >  2:320948/1‑251 (MQ=255)                                                                                                                      tggtgctAACTGCAGTCATTGCTGGCTGTGTTCAGACCAAAAATATAACAGCACCGGATTGGTCATCGACTGTTCAATAGCCTGTTTCAGGTGAGATTCGTTAATGTTGACAATATTTTCTACGTACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGTCATGGGGGCTTAACCGCGCGCTTAAACTCACCCCTGCAATATTTTGTCCATCTCGCGCACAGGTATCTGTCGCGTAAGCACCATTAcc                                                                                                                <  2:98609/244‑1 (MQ=255)                                                                                                                        gggTTAAGTGCAAACAGTGCTGGCTACGTTCAGAACAAAAATAGAACAGCAGCTGAGTGGTCATCAACTGTTCAAGAGCCTGTTGCAGGTTAGATTCGGTAATGTTGACAATATTTTCTACTGACATGAAG‑TCGCTCTCTGTTGGCGAGATTTTCTTTGACATGGGGGCTTAATCTCGCGCTTCAACTCACCCCTGCGATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACTATTTCCGcc                                                                                                             <  2:196463/251‑1 (MQ=255)                                                                                                                        gggTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGcc                                                                                                             >  1:445737/1‑251 (MQ=255)                                                                                                                            tAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAgg                                                                                                         >  1:250779/1‑251 (MQ=255)                                                                                                                                            tGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCgg                                                                                         >  1:242056/1‑251 (MQ=255)                                                                                                                                                      cAGACCAAAAATAGAACAGCACCGGAGTCGTCATCGACTGTTCAAGAACCTGTTGAAGGGTAGAATCGATAATGTTGACAATATTTTCTACGGACATGACG‑TCGCTCTATGTTGTCGATATTTTCTTTGTCATGGGGGCTTAAGCGCGCGCTTAAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCAACAGCGTCACCGGATAGCGCATc                                                                               >  2:217060/1‑251 (MQ=255)                                                                                                                                                               aaTAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTtc                                                                      <  2:250779/251‑1 (MQ=255)                                                                                                                                                               aaTAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTtc                                                                      >  2:159869/1‑251 (MQ=255)                                                                                                                                                                  aGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCg                                                                   <  1:32334/251‑1 (MQ=255)                                                                                                                                                                      cAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCta                                                                >  2:351359/1‑250 (MQ=255)                                                                                                                                                                              gAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGc                                                                  >  1:253243/1‑240 (MQ=255)                                                                                                                                                                              gAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGc                                                                  <  2:253243/240‑1 (MQ=255)                                                                                                                                                                                                   gaaacTGTTGCAGGTTAGAGTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGTTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCCCCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCtccacc                                  <  1:86546/247‑1 (MQ=255)                                                                                                                                                                                                         gTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTcc                            >  2:79824/1‑251 (MQ=255)                                                                                                                                                                                                            gCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACtt                                             <  2:12505/231‑1 (MQ=255)                                                                                                                                                                                                            gCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACtt                                             >  1:12505/1‑231 (MQ=255)                                                                                                                                                                                                             cAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTc                        >  1:161978/1‑251 (MQ=255)                                                                                                                                                                                                                 ttAGATTCGTTAATGTTGGCAATATTTTCTACGGACGTG‑AGATCGCTCACTGTTGTAGCTATTTTCTTTGACATGGGGGCTTAAACGCGCGCTTCAACTCACTCCTGCAATATTTTGTCCACCACGCGCCCCGGCAGCAGGCGCTTAAGCACCCTTACCGCCCAGGTCACCAGCATCACCGGATAGCGCATCTTCGGCTTCTCGCTAATCAAAGCATGGCGTACTTTGTAAACCACCGCTTCCGGTccaa                     >  2:205715/1‑248 (MQ=255)                                                                                                                                                                                                                  tAGATTCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCaa                    <  1:250305/250‑1 (MQ=255)                                                                                                                                                                                                                      ttCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGTCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTa                >  2:99507/1‑250 (MQ=255)                                                                                                                                                                                                                       tCGTTAATGTTGACAATATTTTCTACGGCCATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCAt                                                                                                                                                     >  2:253691/1‑116 (MQ=255)                                                                                                                                                                                                                       tCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCAt                                                                                                                                                     <  1:253691/116‑1 (MQ=255)                                                                                                                                                                                                                       tCGTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTaa               >  1:241088/1‑250 (MQ=255)                                                                                                                                                                                                                        cGTTAATGTTTACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAg             <  2:262858/251‑1 (MQ=255)                                                                                                                                                                                                                         gTTAATGTTGACAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCCGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAgc            <  2:168691/251‑1 (MQ=255)                                                                                                                                                                                                                                   aCAATATTTTCTACGGACATGAAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGAt  <  2:435369/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       ATACACGGTCGGAATCGCACGCAGACCAAACTGCGCGGCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAG‑TCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGAT  >  NC_000913/518034‑518510

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