Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,276,422	C→T	P318S (CCG→TCG)	pdeC →	c‑di‑GMP phosphodiesterase PdeC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,276,422	0	C	T	96.7%	181.6 / NA	61	P318S (CCG→TCG)	pdeC	c‑di‑GMP phosphodiesterase PdeC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (31/27);  minor base G (0/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTT  >  NC_000913/4276200‑4276667                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      aTTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCACACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGaaa                                                                                                                                                                                                                            >  1:408974/1‑250 (MQ=255)         gTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAgggg                                                                                                                                                                                                                    >  2:447323/1‑250 (MQ=255)            gtccgccgGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAggggctg                                                                                                                                                                                                                 >  1:452096/3‑247 (MQ=255)                  cgGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAAc                                                                                                                                                                                                          <  1:118184/251‑1 (MQ=255)                     tATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCa                                                                                                                                                                                                                                <  2:339702/226‑1 (MQ=255)                     tATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCa                                                                                                                                                                                                                                >  1:339702/1‑226 (MQ=255)                         ttGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATaa                                                                                                                                                                                                    <  1:3459/250‑1 (MQ=255)                           ggggAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAATCTAGTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACGTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGGTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTg                                                                                                                                                                                                 <  1:261255/251‑1 (MQ=255)                               aGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTa                                                                                                                                                                                             <  1:234832/251‑1 (MQ=255)                                                   ctGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTcc                                                                                                                                                                         >  2:465814/1‑251 (MQ=255)                                                     ggCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTcc                                                                                                                                                                         <  2:583025/249‑1 (MQ=255)                                                         aCGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGAGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCAGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATc                                                                                                                                                                   >  1:130481/1‑251 (MQ=255)                                                                   cAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGAtt                                                                                                                                                          >  1:428964/1‑250 (MQ=255)                                                                               cctcCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTCCCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAAcaca                                                                                                                                             <  2:531821/250‑1 (MQ=255)                                                                                  cccAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAAcaca                                                                                                                                             <  1:459167/248‑1 (MQ=255)                                                                                     aaacgggcaTTGCAACGCGCAATCGAAAAAAATCAACTTTGTATTTATTCCCAGCCAATAATCGATATAAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCCAATAATGCATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGACTGATAGATCAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGc                                                                                                                                          >  2:594088/9‑248 (MQ=255)                                                                                        cGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCg                                                                                                                                    <  1:169759/251‑1 (MQ=255)                                                                                                  aaCGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCt                                                                                                                          <  1:447323/251‑1 (MQ=255)                                                                                                  aaCGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATGGATAATGTCTTCCGCTATCTCGGCGAGTACCTGGCAACACATGCAGATCGCTATGTTTCt                                                                                                                          <  2:619791/251‑1 (MQ=255)                                                                                                      cgcCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTTAGCAGGGTCAAATAATGAATGCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTAGTGATAATGTCTTCCTCGATCTTGGCGATTGCCTGGCAACACATGCAGATCGCTATGTTCCTATTa                                                                                                                      <  2:574071/251‑1 (MQ=255)                                                                                                      cgcCCTCGAAAAACATCAACTTTGTCTTTATTAACAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTAttt                                                                                                                      <  2:165695/251‑2 (MQ=255)                                                                                                         ccTCGAAAAACATCAACTTTGTCTTTCTTACCAGCCAATAATCGATATCAAAACAGACAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAAGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAAcc                                                                                                                   >  2:481202/1‑251 (MQ=255)                                                                                                                  aCATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCt                                                                                                          <  1:533734/251‑1 (MQ=255)                                                                                                                       aaCTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGAt                                                                                                     <  2:408974/251‑1 (MQ=255)                                                                                                                             gtatttATTACCAGCCAATAATCGATCTCAAACCAGAAAAATATATCGGAACAGCAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGAAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGCTGATAGAACAGATCACTGATTATGTTATTGATAATGTCTTCCGCGCTCTGGGCGATTACCTGGCAACACATGCAGCTCGCTATGTTTCTATTACCCTGTCGGACTCCGATTTTCAt                                                                                               >  2:599212/4‑251 (MQ=255)                                                                                                                                   attaCCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGCTTTTCATACGTc                                                                                          >  1:268488/1‑250 (MQ=255)                                                                                                                                       ccAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGt                                                                                     >  2:14578/1‑251 (MQ=255)                                                                                                                                        cAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGACAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGtt                                                                                    >  1:255763/1‑251 (MQ=255)                                                                                                                                          gCCAATAATCGATATCCAAACAGGAAAATGTCTCGGCGCTGAAGCGTAGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGGTCTGGGCGATTACCTGGCACCACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTg                                                                                   >  2:40562/1‑250 (MQ=255)                                                                                                                                          gCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATGACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTgt                                                                                  >  2:57279/1‑250 (MQ=255)                                                                                                                                            cAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGAt                                                                                 >  2:42323/1‑250 (MQ=255)                                                                                                                                                    gATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGaat                                                                        >  1:189265/1‑251 (MQ=255)                                                                                                                                                     atatCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGCTGATAGCGCGaatc                                                                       >  1:3472/1‑251 (MQ=255)                                                                                                                                                       atCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGcgcgtat                                                                        <  2:26345/248‑4 (MQ=255)                                                                                                                                                       atCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGaat                                                                        >  1:26345/1‑248 (MQ=255)                                                                                                                                                        tCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGaatcaat                                                                    >  1:396218/1‑251 (MQ=255)                                                                                                                                                            aaCAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGa                                                                <  1:80281/251‑1 (MQ=255)                                                                                                                                                             aCAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGa                                                                >  1:376088/1‑250 (MQ=255)                                                                                                                                                              cAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGCGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGaa                                                               >  2:29892/1‑250 (MQ=255)                                                                                                                                                               aGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGCGCAGGGGCAAATAATGAATTAGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGaac                                                              >  1:47648/1‑249 (MQ=255)                                                                                                                                                                  aaaaTGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAAc                                                            >  1:547227/1‑249 (MQ=255)                                                                                                                                                                    aaTGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACagag                                                        <  1:293587/251‑1 (MQ=255)                                                                                                                                                                      tGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGc                                                       >  2:218433/1‑250 (MQ=255)                                                                                                                                                                       gTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCa                                                      <  1:385035/250‑1 (MQ=255)                                                                                                                                                                        tATCGGAGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAAt                                                    >  1:119062/1‑251 (MQ=255)                                                                                                                                                                         aTCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGGTATGTTTCTATTAACCTGTCGGCCTCCAATGTTCATACGTCACGGTGGATAGCGCGAAGCAATCATAAAACAGAGCAAt                                                    <  2:47648/250‑1 (MQ=255)                                                                                                                                                                               gcTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCggg                                             >  1:286954/1‑250 (MQ=255)                                                                                                                                                                               gcTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGt                                             >  1:495250/1‑251 (MQ=255)                                                                                                                                                                                  gAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATGCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGc                                           <  2:200373/250‑1 (MQ=255)                                                                                                                                                                                     gCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTcc                                       <  2:495250/251‑1 (MQ=255)                                                                                                                                                                                     gCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTcc                                       <  1:201407/251‑1 (MQ=255)                                                                                                                                                                                      cGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCg                                      <  1:29892/251‑1 (MQ=255)                                                                                                                                                                                         tgttACGGTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCgca                                    <  2:549322/250‑1 (MQ=255)                                                                                                                                                                                             aCGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCaa                                >  1:9367/1‑250 (MQ=255)                                                                                                                                                                                                   gCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGGCCGCAGCaaattaaat                         <  2:235844/251‑1 (MQ=255)                                                                                                                                                                                                     cTGGTGAGCAGGGGCAAATAATGCATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAAttt                       >  1:142697/1‑251 (MQ=255)                                                                                                                                                                                                        gTGAGCAGTGGAAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGTGGATGATAGAACAGATAACTGATTAGGTTATTGATAATGTCTTCCGCGATCTGGGCGATGACCGGGCAACACATGCAGATCGCTATGTTTCTATGAACCGGTCGGCCACCGATTTGCATACGTCACGGTTGATAGCGCGAATCAAACAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAAGTAAAttttaa                    <  2:67113/251‑4 (MQ=255)                                                                                                                                                                                                          gAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGt                  <  2:265539/251‑1 (MQ=255)                                                                                                                                                                                                             cAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGAc               >  2:451524/1‑251 (MQ=255)                                                                                                                                                                                                                 ggCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAg           <  1:596002/251‑1 (MQ=255)                                                                                                                                                                                                                           gAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCAttt  <  1:530497/250‑1 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      ATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTT  >  NC_000913/4276200‑4276667

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