Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	874,178	Δ1 bp	coding (1200/2349 nt)	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,178	0	A	.	98.0%	228.7 / ‑2.4	51	coding (1200/2349 nt)	pdeI	putative c‑di‑GMP phosphodiesterase PdeI
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (19/31);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

CGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGTAAACGATTTTGACAGCAGCTTTGGCTACCGCAATCTGCGCATCCATAAGTCTGC  >  NC_000913/873949‑874420                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   cGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTgg                                                                                                                                                                                                                              <  1:316340/251‑1 (MQ=255)     tttGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAg                                                                                                                                                                                                                           <  2:634887/250‑1 (MQ=255)               acAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGgt                                                                                                                                                                                                                <  2:204113/250‑1 (MQ=255)                   cGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTtg                                                                                                                                                                                                            >  1:136750/1‑251 (MQ=255)                        cAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCTGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTTCTCTCACTTATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGAGTATTGAGTCTAACAATAATCAGTTTCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTCGGTGTtgatga                                                                                                                                                                                                        <  2:294922/250‑1 (MQ=255)                           aTCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAgcgc                                                                                                                                                                                                    >  1:636940/1‑251 (MQ=255)                                      gAAAATTCGAGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCAGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGATTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATATGGCTGTGGGGATCATCAA‑CTGTGTGGGCTGGAAGCCGTGGAAGCGGGGTTTGGTGTTGATGAGCGCAATAAAATCGt                                                                                                                                                                                         >  2:27961/1‑251 (MQ=255)                                          aTTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCATAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGGCAgg                                                                                                                                                                                     >  2:512405/1‑251 (MQ=255)                                             cGTGAAGGCGATGTGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGGTTTATTAAAAATATGGCGATTATGGAGTCTAACAATAATCGGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGaaa                                                                                                                                                                                  <  2:622127/251‑1 (MQ=255)                                                    gCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATc                                                                                                                                                                                           <  1:530873/235‑1 (MQ=255)                                                    gCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGGGTTTGGGGTTGATGGGCGCAATAAAATc                                                                                                                                                                                           >  2:530873/1‑235 (MQ=255)                                                      gATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTTGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAg                                                                                                                                                                         <  2:631850/251‑1 (MQ=255)                                                            aaaGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTTAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTAGTGAGTCTAACAATAATCAGTATCTGGCTGTTGGTATCATCAA‑CTGTGTGGTCTGGAAGCCGTGTAAGCGTGGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTTTGTCAGCGAAtt                                                                                                                                                                   <  2:697242/251‑1 (MQ=255)                                                                    tATCAGCATCATACATTCCATCGCGGCCATGCAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGACAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGAAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCgaga                                                                                                                                                            >  2:478616/1‑250 (MQ=255)                                                                     aTCAGCCTCAAAACTTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGCTCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGACATTGTGTCAGCGAATTGCCGAGaa                                                                                                                                                           >  2:467251/1‑250 (MQ=255)                                                                        aGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAatat                                                                                                                                                       <  1:325392/251‑1 (MQ=255)                                                                                aaaTTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTTAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAAtg                                                                                                                                               <  1:559451/251‑1 (MQ=255)                                                                                          gcgGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATc                                                                                                                                     <  1:51989/251‑1 (MQ=255)                                                                                             gAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTg                                                                                                                                  <  2:636940/251‑1 (MQ=255)                                                                                                  gAAAGATTCCGAAATTGAACGGTCAAATAAATTGTTCTCACTGATCGCTTTCGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTTGTCTGGAAGCCGTGGAAGCGGTGTTTAGTGTGGATGAGCGCAATAAAATCGGCAGGAAATTGTGTCAGCGAATGTCCGAGAAATATGCGTAATGCTGCGATACCGTGACAtt                                                                                                                             <  2:655066/251‑1 (MQ=255)                                                                                                  gAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACAtt                                                                                                                             <  1:467251/251‑1 (MQ=255)                                                                                                   aaaGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTc                                                                                                                            <  1:289438/251‑1 (MQ=255)                                                                                                     aGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCaa                                                                                                                          <  1:594927/251‑1 (MQ=255)                                                                                                      gATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATAGCGTGACATTCAAt                                                                                                                         <  1:478616/251‑1 (MQ=255)                                                                                                      gATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGAAAAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAAt                                                                                                                         >  1:478177/1‑251 (MQ=255)                                                                                                          ccGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCg                                                                                                                     <  2:612693/251‑1 (MQ=255)                                                                                                           cGAAATTGAACGGTCAAATAAATTGCTCTCCCTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATAGGCGCAATGCTGCGATATCGTGACATTCAATGCCGa                                                                                                                    <  2:155862/251‑1 (MQ=255)                                                                                                           cGAAATTGAACGGTACAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACCATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGCATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCg                                                                                                                     >  2:453650/1‑250 (MQ=255)                                                                                                                   aaCGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTAtt                                                                                                            <  2:136750/251‑1 (MQ=255)                                                                                                                    aCGGTCAAATAAATTGCTCTCACGGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAATCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGTAAATTGGTTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTAttt                                                                                                           <  1:453650/251‑1 (MQ=255)                                                                                                                        tCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACtt                                                                                                       >  2:497872/1‑251 (MQ=255)                                                                                                                               aatTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTc                                                                                                >  1:201833/1‑251 (MQ=255)                                                                                                                                        cACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGaaa                                                                                          <  1:43401/248‑1 (MQ=255)                                                                                                                                         aCTGATCTCTGACGATCAGGAAATTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTTGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGGGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGt                                                                                      <  2:406051/251‑1 (MQ=255)                                                                                                                                                   tACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACAttt                                                                            >  2:672307/1‑251 (MQ=255)                                                                                                                                                     cGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTAc                                                                          >  1:228387/1‑251 (MQ=255)                                                                                                                                                       aTCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTAAAGACATTTAccc                                                                        >  2:78785/1‑251 (MQ=255)                                                                                                                                                           ggAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTaa                                                                    <  2:275337/251‑1 (MQ=255)                                                                                                                                                                             aaagatGGCGATTATTTAGTCTAACAATCATCAGTATCCGGATGTGAGGATCAT‑AAACCGTATGGTCTGGAAGCCGCGGAAGCGGTGTTTGGTGTTGATGCGCGCCATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGTGCAATATGAGCAATGCTGCGATATCGTGACATTAAATGCCGTACTCTATTTACTTCTGTGTCGGGGAAGTGAACAGACATTTACTCGTAACATAGCGATGGTAAAca                                                   >  2:673177/5‑249 (MQ=255)                                                                                                                                                                                   ggCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGTAAACGATTTTg                                             <  1:589656/250‑1 (MQ=255)                                                                                                                                                                                      gATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGTAAACGATTTTGAcag                                         <  2:465850/251‑1 (MQ=255)                                                                                                                                                                                         tattGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGTAAACGATTTTGAcagca                                       <  1:512405/250‑1 (MQ=255)                                                                                                                                                                                           ttGAGTCTAACAATCATAAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTAGAAGCAGTGGCAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATCTGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGTAAACGATTTTGACAGCAGCt                                    >  2:328804/1‑251 (MQ=255)                                                                                                                                                                                                         aatCAGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGGGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATTTACAGACATTTACCCGTAAAATAGCGATGGTAAACGATTTTGACAGCAGCTTTGGCTACCGCAAt                      <  1:427041/251‑1 (MQ=255)                                                                                                                                                                                                             aGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTaa                                                                    >  2:618709/1‑201 (MQ=255)                                                                                                                                                                                                             aGTATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTaa                                                                    <  1:618709/201‑1 (MQ=255)                                                                                                                                                                                                             aGTATCTGGCGGGGGTGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGGGTTGATGAGCGCAATAAAATCGTCAGGAAAATGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGAGATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGTAATCGATTTTGACAGCAGCTTTGGCTACCGCAATCTgc                  <  2:466875/251‑1 (MQ=255)                                                                                                                                                                                                               tATCTGGCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGTAAACGATTTTGACAACAGCTTTGGCTACAGCAATCTgcg                 >  1:560805/1‑250 (MQ=255)                                                                                                                                                                                                                    ggCTGTGGGGATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCACTAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGTAAACGATTTTGACAGCAGCTTTGGCTACCGAAATCTGCGCATcc            >  2:698091/1‑250 (MQ=255)                                                                                                                                                                                                                          ggggATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGAGTTTGGTGTTGATGAGCGGAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGTAAACGATTTTGACAGCAGCTTTGGCTACAGCAATCTGGGCAACCATAAGt      >  2:650883/1‑250 (MQ=255)                                                                                                                                                                                                                             gATCATCAA‑CTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGTAAACGATTTTGACAGCAGCTTTGGCTACCGCAAGCTGCGCAGCCATAAGTCTGc  <  1:171394/251‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   CGTTTTTGCCTGACTCAACGGATCAAATCGAGGCTCTGAAAATTCGTGAAGGCGATGCGAAAGAGATTATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCTGTGTCGGGAAAATGTACAGACATTTACCCGTAAAATAGCGATGGTAAACGATTTTGACAGCAGCTTTGGCTACCGCAATCTGCGCATCCATAAGTCTGC  >  NC_000913/873949‑874420

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: ‑