Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	306,294	C→T	S175N (AGC→AAC)	ecpD ←	fimbrial adhesin EcpD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	306,294	0	C	T	93.3%	170.3 / ‑5.7	61	S175N (AGC→AAC)	ecpD	fimbrial adhesin EcpD
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (24/32);  minor base A (3/1);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATGCAGCGCAAGCCAGTTAAAGGAGATGAAACTGGTGAATTTTCCAGCCACATATCAAACTTCCAGTTGGTATA  >  NC_000913/306048‑306535                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   aaTCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAgttgct                                                                                                                                                                                                                                                <  1:330868/251‑1 (MQ=255)  aTCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCCGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCAGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAGGttgctg                                                                                                                                                                                                                                               >  2:561991/1‑251 (MQ=255)    ccTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGCGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGCTGGTCATTGTAAAGTTGCTGcc                                                                                                                                                                                                                                             >  2:650456/1‑251 (MQ=255)     cTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGGGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCCTTGTAAAGTTGCTGCCg                                                                                                                                                                                                                                            >  2:525340/1‑251 (MQ=255)          cGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACggg                                                                                                                                                                                                                                       <  1:349591/251‑1 (MQ=255)           ggCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGc                                                                                                                                                                                                                                      <  1:72899/251‑1 (MQ=255)                  gATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATttgtt                                                                                                                                                                                                                               >  2:189105/1‑251 (MQ=255)                     ggTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCa                                                                                                                                                                                                                            <  1:525340/251‑1 (MQ=255)                            gTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTg                                                                                                                                                                                                                     <  1:509227/251‑1 (MQ=255)                                                  cGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGCCGGGCATTTGTTGCATATCCTGGTAAAACGCGTCTGacacactc                                                                                                                                                                                               >  2:149163/1‑246 (MQ=255)                                                       agagTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCat                                                                                                                                                                                                                           >  1:62889/1‑218 (MQ=255)                                                       agagGCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCat                                                                                                                                                                                                                           <  2:62889/218‑1 (MQ=255)                                                        gagTCGGTACTCCGTCTCTGTTGATAAATACTTCCGCCAGCGTGTGGGTAGTTATCAACCGTTGATTTGCTGCTGTCGTATGGGTGACGTTGCGAACATACCAGGTGCCGGAGGCCTGATCCTTACAGCGTGCTCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACGGGTAAAACGCGTCTGACAACATGCCGTg                                                                                                                                                                                          <  1:418375/250‑1 (MQ=255)                                                             ggTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATCCCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGCTGGTCATTGTAAAGTTGCTGCCGGCGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATcc                                                                                                                                                                                     >  2:74217/1‑250 (MQ=255)                                                                aCTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACtt                                                                                                                                                                                 <  1:480110/251‑1 (MQ=255)                                                                    cgtttatGTTGATAAATGTTTCCGCCATCGAGTGTGTATTTATCAACCGTAGATTTGCTGCTTGCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTGACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCTCTTcgcg                                                                                                                                                                             <  2:450196/245‑1 (MQ=255)                                                                               aTAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGAAAAACGCGTCTTACAACATGCTGTGCATCCACTTCGCGCCGACGATGGt                                                                                                                                                                  >  2:456184/1‑251 (MQ=255)                                                                                tAAATACTTACGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGCTGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTc                                                                                                                                                                 >  1:436501/1‑251 (MQ=255)                                                                                 aaaTACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTc                                                                                                                                                                 <  1:74217/250‑1 (MQ=255)                                                                                  aaTACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCAc                                                                                                                                                               <  1:206084/251‑1 (MQ=255)                                                                                    tACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCg                                                                                                                                                             >  1:41237/1‑251 (MQ=255)                                                                                     aCTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGc                                                                                                                                                                                              <  1:450190/217‑1 (MQ=255)                                                                                                cGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATaaa                                                                                                                                                 <  1:471983/251‑1 (MQ=255)                                                                                                 gAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCAttt                                                                                                                                                                                                                                  >  1:490458/1‑169 (MQ=255)                                                                                                 gAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCAttt                                                                                                                                                                                                                                  <  2:490458/169‑1 (MQ=255)                                                                                                    ggggTATTTATCAACCGTAGATTTGCTGCTGGCGTATGGGGGACGTTGCGAACATACCAGTTGCCGGAGGCCTGAGCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATcc                                                                                                                                             <  1:323964/250‑1 (MQ=255)                                                                                                      ggTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCAc                                                                                                                                           <  1:662200/251‑1 (MQ=255)                                                                                                              tCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGCTCCTTACAGCGTGCGCCGCTGCTGGCGTCATCGTTCACCGAGGTCTGGCAGGCATTGATGGTCATTGTACAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAACCGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCCCCATAAAATCCACTGGCGTCg                                                                                                                                   >  2:584382/1‑251 (MQ=255)                                                                                                              tCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACACCATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGAGTCg                                                                                                                                   >  2:216188/1‑251 (MQ=255)                                                                                                                     ccgtccctttcccgacctATGGGTGACGTTGCGAACATTCCAGTTTCCGGAGGCCTGATCCTGTCAGCTTGCGCTGCTTCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGTTAAAACGCGTCTGACAACATGCCGTGCATCCACGTCGCGCCGCCGTTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGttt                                                                                                                            <  2:423275/234‑1 (MQ=255)                                                                                                                        aTTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCg                                                                                                                         >  1:362615/1‑251 (MQ=255)                                                                                                                           tgctgcTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCa                                                                                                                                                                                    <  1:50615/189‑1 (MQ=255)                                                                                                                               gctTTCGTATGGGTTACGTTGCGAACCTTCCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGTGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCAATGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACAGGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCCCTGGCGTCGGTGGTTTGCGGCAGGAt                                                                                                                  <  2:325668/251‑1 (MQ=255)                                                                                                                                    cGTATGGGGGACGTTGCGAACATACCAGTTGCCGGCGGCATGATCCTTACAGCGTCAGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGATGCCGCCGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACACCATGCCGGGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCACCATAAAATACACTGGCTTCGGTGGTTTGCGGCAGGATAAGAc                                                                                                             >  2:567235/1‑251 (MQ=255)                                                                                                                                    cGTATGGGCGACGTTGCGAACATACCAGTCACCGGAGGCCTGCTCCTTACAGCGTGCGCCGCTGCTGGCGTCATCGATCACAGAGGTCTGGCAGGCCTTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGAAAAACGCGTCTGACAACATGCCGTGCATCCACCTCGAGCCACCGCTGGACACCGACGCGCCATAAAATCCACTGGCGCCGGTGGTTTGAGGCAggtaacaac                                                                                                             >  2:130227/1‑244 (MQ=255)                                                                                                                                       aTGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGttt                                                                                                                            >  1:423275/1‑233 (MQ=255)                                                                                                                                         gggTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACtggt                                                                                                         <  1:650456/250‑1 (MQ=255)                                                                                                                                                gTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGAATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCAATGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑Cctat                                                                                                >  2:209532/1‑251 (MQ=255)                                                                                                                                                      aaCATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACAcc                                                                                          <  1:174625/251‑1 (MQ=255)                                                                                                                                                                 tGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATATTTCACAGAGGTCTTGCAGGCATGGATGGTCATTGTAAAGTTGCTGCCGATGGGCATTTTTTGCATATACTGGTAAAACGCGTTTGACAACAGGCCGTGCAGCCACTTCGCGCCGCCGCTGGTGACCGTCGCGCCATAAAATCCACTGGCGTCTGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATGACACCCTGCGTCCCa                                                                                <  2:565461/250‑1 (MQ=255)                                                                                                                                                                     ggAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGA‑CATATTACACCCAGCGTACCAgatg                                                                            >  1:341506/1‑247 (MQ=255)                                                                                                                                                                       aGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATg                                                                         <  1:419835/251‑1 (MQ=255)                                                                                                                                                                          ccTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGAAGGATAAGACTGGGGGT‑CCTATTACACCCAGGGTACCAGTTGATGc                                                                        >  2:171851/1‑249 (MQ=255)                                                                                                                                                                                     cAGCGTGCGCCGCTGCTGGCGCCATAGTCCACAGAGGTCTGGCAGGCATTGATGTTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTTCATCCTCTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATGCAGCGCAAGCCAGt                                                           <  1:1623/251‑1 (MQ=255)                                                                                                                                                                                         gTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATGCAGCGCAAGCCAGTTaaa                                                       <  1:539420/251‑1 (MQ=255)                                                                                                                                                                                             gccgcTGCTGGCGTCATAGTTCACAGAGGGCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATGCAGCGCAAGCCAGTTAAAGgag                                                   <  1:4475/251‑1 (MQ=255)                                                                                                                                                                                              ccgcTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTGGATGCAGCGCAAGCCAGTTAAAGgaga                                                  <  1:149163/251‑1 (MQ=255)                                                                                                                                                                                                   gctgGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATGCAGCGCAAGCCAGTTAAAGGAGATGaaa                                             <  1:189105/251‑1 (MQ=255)                                                                                                                                                                                                     tgGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTc                                                                                                                                                           >  1:371102/1‑140 (MQ=255)                                                                                                                                                                                                     tgGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTc                                                                                                                                                           <  2:371102/140‑1 (MQ=255)                                                                                                                                                                                                     tgGCGTCATAGTTCACAGAGGACTGGCAGGCATTGATGGTCATTGTAAAGATGCTGACGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGCTAAGACTGGTGGG‑ACTATTAAACCCCACGGGCCAGTTGATGCAGCGCACACCAGTTAAAGGAGATGAAAcc                                           >  2:401745/1‑250 (MQ=255)                                                                                                                                                                                                            aaagttgtcagtggTCTGGCGGGCATTGATGGTTATTGTAAAGATGGAGCGGACGGGCTGTTGTTGCATATACTGGTAAAACGCGTCCAACAACATGCAGTGCATCCACTCCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGAAGGATAAGTCTGGTGGT‑CATATTACACCCAGCGTAACAGTTGATGCGGTGCAGGCGAGTTAAAGGAGATGAAACTGGTGAAt                                    <  1:532732/239‑1 (MQ=255)                                                                                                                                                                                                              aGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATGCAGCGCAAGCCAGTTAAAGGAGATGAAACTGGTGAAttt                                  <  2:100050/251‑1 (MQ=255)                                                                                                                                                                                                               gTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATGCAGCGCAAGCCAGTTAAAGGAGATGAAACTGGTGAAtttt                                 >  1:248738/1‑251 (MQ=255)                                                                                                                                                                                                                ttCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAACGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGACGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCCGGGTAAGACTGGTGGT‑CATATTACACCCCGCGTACCAGTTGATGCAGAGAAAGCCAGTTACAGGAGGTGAAACTTGTGAATTTTc                                >  2:519528/1‑251 (MQ=255)                                                                                                                                                                                                                ttCACAGAGGTAAGGCAGGCATTGCTGGTCATTGTAAAAATTCAGCCGAAGTGCATTTGGTGCATATCCTGGTAAAACGCGACTGACAACATGTCGTGCATCCACTTCGCGCCGCCGCTGGTCAACGTCGCGCCATAAAATCCACTGGCGTCGGTGGGTTGAGGCAGGATAAGACAGATGGTTCA‑ATTACACCCAGCGTACCAGCTGTTGCAGATCTAACCAGATAAACGAGATGCAAGTGGAGAAttttt                                >  2:563603/1‑250 (MQ=255)                                                                                                                                                                                                                 tCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAAGCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATGCAGCGCAAGCCAGTTAAAGGAGATGAAACTGGTGAATTTTc                                <  1:313957/250‑1 (MQ=255)                                                                                                                                                                                                                 tCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGACGACGGGCATTTGTTGCATATACTGGTAAACCGCGTCTGACACCATAACGTGCATCCACTTCGCGCCGCCGCTGGTAACCGTCGCGCCATAAAATCCACTGGCGTCCGTGGTTTGCGGCAGGCTAAGACTGGTGGT‑CATATTAAACCCAACGTACCAGTTGGTGCAGAGCAAGCCCGTTCGAGGAGAttgaccaggtgaactttac                               >  1:45675/1‑233 (MQ=255)                                                                                                                                                                                                                     agagGTCTGGCAGGCATTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATGCAGCGCAAGCCAGTTAAAGGAGATGAAACTGGTGAATTTTCCAGcc                           <  1:559464/251‑1 (MQ=255)                                                                                                                                                                                                                             ggcaggcaTTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATGCAGCGCAAGCCAGTTAAAGGAGATGAAACTGGTGAATTTTCCAGCCACatat                     <  1:519528/249‑1 (MQ=255)                                                                                                                                                                                                                               caggcaTTGATGGTCATTGTAAAGTTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATGCAGCGCAAGCCAGTTAAAGGAGATGAAACTGGTGAATTTTCCATCCACATATCaa                  >  1:665516/1‑250 (MQ=255)                                                                                                                                                                                                                                      tgttggctagtGTAAAGTTGCTGGCGTCGGGCATTTGTTGCATCGTCTGGTAAAACGATGCTTACAACATGTCTTGCTTCCACTGCGCGCTGCCTCTGGTCACCCGCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGTATAAGACTGGTGGG‑CATATTACACCCAGCGTTCCAGTTGATGCAGCGCAAGCCAGTTAAAGGAGATGAAACTGGTGAATTTTCCAGCCACATATCAAACTTCCa           <  2:658810/240‑1 (MQ=255)                                                                                                                                                                                                                                               ttGTAAAGATGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGTCGCTGGTCACCGTAGCGCCATAAAGTCCACTGACGTCGGTGGTTTGCGGCAGGTTAAGGATAGCGGT‑CACATTACAACCAACGTACCCGTTGATGCCGGGGAAGCCCATTCAATGAGATGAAAAAGGTGATATTTCCAGCCACATATACAACTTCCAGTTGGtatt  >  2:580478/1‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   AATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTGGT‑CATATTACACCCAGCGTACCAGTTGATGCAGCGCAAGCCAGTTAAAGGAGATGAAACTGGTGAATTTTCCAGCCACATATCAAACTTCCAGTTGGTATA  >  NC_000913/306048‑306535

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