Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	874,129	C→T	A384V (GCG→GTG)	pdeI →	putative c‑di‑GMP phosphodiesterase PdeI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	874,129	0	C	T	97.7%	141.2 / ‑6.7	44	A384V (GCG→GTG)	pdeI	putative c‑di‑GMP phosphodiesterase PdeI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/24);  minor base A (1/0);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ATTTAATGACCACTCGTGTGTTCCGTCGGCAATTTTCTGAAATGACAGAGCTGGTTAATACGCTGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGT  >  NC_000913/873883‑874366                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              atttaATGACCACTCGTGTGTTCCGTCGGCAATTTTCTGAAATGACAGAGCTGGTTAATACGCTGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGatt                                                                                                                                                                                                                                           <  1:580881/251‑1 (MQ=255) atttaATGACCACTCGTGTGTTCCGTCGGCAATTTTCTGAAATGACAGAGCTGGTTAATACGCTGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGat                                                                                                                                                                                                                                            <  2:378828/250‑1 (MQ=255)      aTGACCACTCGTGTGTTCCGTCGGCAATTTTTTGAAATGACAGAGCTGGTTAATACGCTGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGCTGAGAAATAGATTATCAGCATCAAAAATTCGATCGCGGAATTGAAAGATGCCGAAATTGAACGGTCAAATAAATTGTTCTCACTGCTCTCTTCCGATCAGGAAAGTGGTTTTATTAGAAATATGGTGATTATTg                                                                                                                                                                                                                                       >  1:476771/1‑250 (MQ=255)        gACCACTCGTGTGTTCCGTCGGCAATTTTCTGAAATGACAGAGCTGGTTAATACGATGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGTAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAACATGGTGATTATTGAGt                                                                                                                                                                                                                                    <  2:123656/251‑1 (MQ=255)               cGTGTGTTCCGTCGGCAATTTTCTGAAATGACAGAGCTGGTTAATACGCTGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAAc                                                                                                                                                                                                                               <  1:382231/249‑1 (MQ=255)               cGTGTGTTCCGTCGGCAATTTTCTGAAATGACAGAGCTGGTTAATACGCTGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACa                                                                                                                                                                                                                              >  2:165265/1‑250 (MQ=255)                               aaTTTTCTGAAATGACAGAGCTGGTTAATACGCTGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCt                                                                                                                                                                                                             >  2:160283/1‑251 (MQ=255)                                    tCTGAAATGACAGAGCTTGTTAATACGCTGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTgg                                                                                                                                                                                                         <  2:352679/250‑1 (MQ=255)                                       gAAATGACAGAGCTGGTTAATACGCTGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACAGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTTATTATTGAGTCTAACAATAATCAGTATCTGGCTGTgggg                                                                                                                                                                                                       <  1:577969/249‑1 (MQ=255)                                                gagCTGGTTAATACGCTGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCCCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCCTCAAAAATTAGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACt                                                                                                                                                                                            >  2:98685/1‑251 (MQ=255)                                                         aaTACGCTGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCt                                                                                                                                                                                   <  2:430390/251‑1 (MQ=255)                                                                             tGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGgtgt                                                                                                                                                               >  1:621255/1‑251 (MQ=255)                                                                                     cGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGat                                                                                                                                                                                                     >  1:647204/1‑205 (MQ=255)                                                                                     cGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGat                                                                                                                                                                                                     <  2:647204/205‑1 (MQ=255)                                                                                       gATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTTGTGTtga                                                                                                                                                     <  2:7163/251‑1 (MQ=255)                                                                                            aaTCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAgc                                                                                                                                                <  2:192778/251‑1 (MQ=255)                                                                                                gAGGCTCTGAAAAGTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCaa                                                                                                                                            <  2:163405/251‑1 (MQ=255)                                                                                                  ggCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATa                                                                                                                                          <  2:321861/251‑1 (MQ=255)                                                                                                                gTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGaa                                                                                                                             >  2:259862/1‑250 (MQ=255)                                                                                                                      gCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTgtgt                                                                                                                      >  2:4451/1‑251 (MQ=255)                                                                                                                       cGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTc                                                                                                                     <  1:219927/251‑1 (MQ=255)                                                                                                                          tGCGAAAGAGATTATCACCATCAAAAATTCGATCGCGGAACTGAACGATGCCGAAATTGAACGGACAAATAACTTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGGAGCGGTGTTTGGGGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGc                                                                                                                  >  1:532449/1‑251 (MQ=255)                                                                                                                           gCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGAGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGc                                                                                                                  <  1:460702/250‑1 (MQ=255)                                                                                                                              aaaGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAAt                                                                                                              <  1:271975/251‑1 (MQ=255)                                                                                                                                         cAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGaaa                                                                                                    <  1:160283/250‑1 (MQ=255)                                                                                                                                           gCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAatat                                                                                                 <  1:489501/251‑1 (MQ=255)                                                                                                                                                  aaaTTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAAt                                                                                          >  1:49910/1‑251 (MQ=255)                                                                                                                                                             cgGAAATGAAAGATGCCGAAAGTGAACTGTCAAATAAATTGCTCTCACTGATCTCTGTCGATCAGGAAAGTGGTTTTATTAAAAAGATGGTGATTTGTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCTTCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGTTGGTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTACCGAGAAATATGCGCGCTGCTGCGTTATc                                                                               <  2:345878/251‑1 (MQ=255)                                                                                                                                                              ggAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCg                                                                              >  2:141704/1‑251 (MQ=255)                                                                                                                                                                      aaGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGc                                                                                                                  <  1:608395/207‑1 (MQ=255)                                                                                                                                                                      aaGATGCCGAAATTGAACAGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGACAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCACACTGTGTGGTCTGGACGCCGTGGACGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGc                                                                                                                  >  2:608395/1‑207 (MQ=255)                                                                                                                                                                         aTGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAAt                                                                   <  1:314116/251‑1 (MQ=255)                                                                                                                                                                             cGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCg                                                               <  1:176067/251‑1 (MQ=255)                                                                                                                                                                                    gaagggTCAAATGAATTGCTCTCACTGATCTCTTACGACCGGGAAAGTGGTTTTATTAAAAATATGGTGAGTATTGCGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAAATGTGTGGTATGGAAGCGGTGGAAGCGGTGTTTGGTGTTGATGCGAGCAATAAAATCGGCAGGAAATTGGGTCAGCGAATTGCCGAGAAGTATGCGCATTGCTGCGATATCGTGACCTTCAATGCCGTtctat                                                         >  2:234544/5‑248 (MQ=255)                                                                                                                                                                                        ggTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTAttt                                                     >  1:560249/1‑250 (MQ=255)                                                                                                                                                                                                           aatGATCTCTTCCGCTCAGGCAAGTGGTTTTATTAAAGATATGTTGATGATTCAGTCTAACAATAATCAGTATCTGACTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGCGTTTCGTGTTGATGAGCGCACTAAAATCGTCAGGAAATTCTGTCAGCGAATTTCCGAGAAATATGCGCACTGCTGCGCTATCGTGCCATTCAATGCCGATCTCTATTTACTTCTGTGTCggtcaaatc                                 >  2:398267/3‑244 (MQ=255)                                                                                                                                                                                                              gATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGAAGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAGTATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTAc                              >  1:224263/1‑251 (MQ=255)                                                                                                                                                                                                                            aGGAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGGGTCGGTAAAATGTACAGACATTTACCCGt                <  2:224263/251‑1 (MQ=255)                                                                                                                                                                                                                              gAAAGTGGTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTa               >  2:627044/1‑250 (MQ=255)                                                                                                                                                                                                                                    ggTTTTATTAAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGc        <  2:235867/251‑1 (MQ=255)                                                                                                                                                                                                                                       tttATTAAAAATATGGTGATTAGGGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTGTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAAGTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGa      <  2:532449/250‑1 (MQ=255)                                                                                                                                                                                                                                        ttattaAAAATATGGTGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATg    >  1:290710/1‑251 (MQ=255)                                                                                                                                                                                                                                          attaAAACTATGGTGATTATTGAGTCTAAGAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGCGTTTGGTGTTGATGAGCGCACTAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGCAAATCTACAGACATGTATCCGTAAAATAGCGagggg  >  2:543847/1‑247 (MQ=255)                                                                                                                                                                                                                                           ttaAAATTATGGTGATTATTGCGTCTAACAATAATCAGTCTCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATACAATCGTCAGGACATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATTTACAGACATTTACCCGTAAAATAGCGATGGt  >  2:407952/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              ATTTAATGACCACTCGTGTGTTCCGTCGGCAATTTTCTGAAATGACAGAGCTGGTTAATACGCTGGCGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGT  >  NC_000913/873883‑874366

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