Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	309,985	T→A	G14G (GGA→GGT)	ecpB ←	putative fimbrial chaperone EcpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	309,985	0	T	A	100.0%	215.2 / NA	60	G14G (GGA→GGT)	ecpB	putative fimbrial chaperone EcpB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (30/30);  total (30/30)

AGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTA  >  NC_000913/309738‑310228                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    agaagaAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCg                                                                                                                                                                                                                                                  >  1:532148/1‑251 (MQ=255)  gaagaaGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCCTCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGCGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCgg                                                                                                                                                                                                                                                 >  1:55999/1‑251 (MQ=255)       agCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAAc                                                                                                                                                                                                                                            <  1:566765/251‑1 (MQ=255)          ggatcacctagccggattGGGCGGGTAGCAGCACGCTGGCGGGAGTGAGAAGCAACTCATACGGCTTGTCCATTGAGGTAACCTGCCCGTCGACAAGCGGGGAAGCGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGGAATGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCCGGGCCGGAGATATACCGGAAAAcagc                                                                                                                                                                                                                                         >  2:331962/17‑251 (MQ=255)             tCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAAcagc                                                                                                                                                                                                                                         >  1:43509/1‑248 (MQ=255)                    gCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCgag                                                                                                                                                                                                                                 >  2:98843/1‑249 (MQ=255)                                gggTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCt                                                                                                                                                                                                                   <  1:69512/251‑1 (MQ=255)                                    agcagcAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAAAAGCAGAGCGAGAGGCCGAAGGAGCttttt                                                                                                                                                                                                               >  1:210681/1‑251 (MQ=255)                                      cagcaAGCTGGCGGGAGTGCGTAGCACCTCATCCAGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGAGAAGAGAGCCGTTCGAGAAGGATATTGCTAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACCGCAGAGCGAGAGGCAGAAGGGGCTTTTTcc                                                                                                                                                                                                             >  2:354636/1‑250 (MQ=255)                                       agcaAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCCTTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCCCTGTTCATAAACGATGATTTATCGCCGACATCCAGAGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCAt                                                                                                                                                                                                            >  2:89306/1‑251 (MQ=255)                                         caAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATa                                                                                                                                                                                                           >  1:223143/1‑250 (MQ=255)                                         aaaGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTtaataa                                                                                                                                                                                                          <  2:559178/250‑1 (MQ=255)                                               ggCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAAt                                                                                                                                                                                                      <  2:200944/249‑1 (MQ=255)                                                gCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTcc                                                                                                                                                                                                   <  1:444362/251‑1 (MQ=255)                                                 cGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGCAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCg                                                                                                                                                                                                  >  2:418779/1‑251 (MQ=255)                                                   ggAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATACCCTGCCCGTCGTCAAGCGGTGCAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATACACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTc                                                                                                                                                                                                >  1:30533/1‑251 (MQ=255)                                                      gtgagtAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATGCCGTCCAg                                                                                                                                                                                             <  2:30533/251‑1 (MQ=255)                                                           tAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAAt                                                                                                                                                                                         >  1:272439/1‑250 (MQ=255)                                                            aGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATaa                                                                                                                                                                                       <  1:64981/251‑1 (MQ=255)                                                                   cATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgg                                                                                                                                                                                 <  1:272990/250‑1 (MQ=255)                                                                      ccGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgggggg                                                                                                                                                                             >  2:331649/1‑251 (MQ=255)                                                                      ccGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAggggg                                                                                                                                                                              <  1:88678/250‑1 (MQ=255)                                                                         gCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCt                                                                                                                                                                           >  2:161506/1‑250 (MQ=255)                                                                            tGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTAccc                                                                                                                                                                       >  1:555945/1‑251 (MQ=255)                                                                                 aTTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCctgct                                                                                                                                                                  <  2:186830/251‑1 (MQ=255)                                                                                   tGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTgg                                                                                                                                                                >  1:315608/1‑251 (MQ=255)                                                                                      gaTAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTAc                                                                                                                                                             <  1:62619/251‑1 (MQ=255)                                                                                         aaCCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTGTGCTCAGCGTGCTGCTGTCACTGTTCAGAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGGTTTTTCATAACAATTCCGGCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATc                                                                                                                                                          <  1:520105/251‑1 (MQ=255)                                                                                            cTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGGTCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACTTCagg                                                                                                                                                       >  2:282050/1‑250 (MQ=255)                                                                                                      ccAAGCGGTGAAGAGGGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGTTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGTGCGAGAGGCAGAAGGTTCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTGCCCCTGCTGGTGCATCAGAGAGATTAACt                                                                                                                                             <  2:476158/250‑1 (MQ=255)                                                                                                        aaGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTgg                                                                                                                                           >  2:55938/1‑251 (MQ=255)                                                                                                                aaGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTc                                                                                                                                   <  1:125697/251‑1 (MQ=255)                                                                                                                   agagCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTcgcg                                                                                                                                >  1:1117/1‑251 (MQ=255)                                                                                                                     agCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTc                                                                                                                              <  2:131796/251‑1 (MQ=255)                                                                                                                          ttCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAA‑CAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTg                                                                                                                        >  1:182080/1‑251 (MQ=255)                                                                                                                            cGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTg                                                                                                                        <  1:98843/250‑1 (MQ=255)                                                                                                                            cGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGt                                                                                                                       <  2:380064/251‑1 (MQ=255)                                                                                                                                   gATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCATGTCGCGGCGAACTGTACGCTaa                                                                                                                <  1:169338/251‑1 (MQ=255)                                                                                                                                      atTGATAAGGCGACCACTGTCGGGACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACACTTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGAACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGt                                                                                                             >  2:491662/1‑251 (MQ=255)                                                                                                                                            aaGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAAAAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTcgccgc                                                                                                       >  1:367255/1‑251 (MQ=255)                                                                                                                                                           ggTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCg                                                                                        <  2:278276/251‑1 (MQ=255)                                                                                                                                                             tACTGTTTTAGCTCGCTTTGCTCAGCGTGCTGCCGTCACTGTTCATAAACGCTGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGt                                                                                      >  2:295085/1‑251 (MQ=255)                                                                                                                                                                                   cAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACt                                                                 <  1:376099/250‑1 (MQ=255)                                                                                                                                                                                      cGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCgg                                                             >  2:27670/1‑251 (MQ=255)                                                                                                                                                                                         gcttctgTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCggtgg                                                          <  2:443312/247‑1 (MQ=255)                                                                                                                                                                                         gctgctGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCggtg                                                           <  1:433075/250‑1 (MQ=255)                                                                                                                                                                                              tgTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTAcc                                                      >  1:106744/1‑250 (MQ=255)                                                                                                                                                                                                 cACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATtg                                                  <  2:435075/251‑1 (MQ=255)                                                                                                                                                                                                  aCTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATtgg                                                 <  2:223143/251‑1 (MQ=255)                                                                                                                                                                                                  aCTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACACTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATtgg                                                 >  2:137552/1‑251 (MQ=255)                                                                                                                                                                                                    tGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATtggtg                                               >  2:177373/1‑251 (MQ=255)                                                                                                                                                                                                        cATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGAGTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTGCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTTCCATTGGTGGTAc                                           <  2:45641/251‑1 (MQ=255)                                                                                                                                                                                                               gatgatATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCtgat                                    <  1:15927/251‑1 (MQ=255)                                                                                                                                                                                                                   agatatCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGatgatg                                <  2:305067/249‑1 (MQ=255)                                                                                                                                                                                                                     atatCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATgg                               <  1:447057/250‑1 (MQ=255)                                                                                                                                                                                                                             gACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTg                       >  2:72377/1‑250 (MQ=255)                                                                                                                                                                                                                             gACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTg                       >  2:150281/1‑250 (MQ=255)                                                                                                                                                                                                                                        ccTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCtgt           >  1:125957/1‑251 (MQ=255)                                                                                                                                                                                                                                           gggCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGc         <  1:161506/250‑1 (MQ=255)                                                                                                                                                                                                                                                 gagaTATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTa  <  1:150281/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    AGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTA  >  NC_000913/309738‑310228

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