Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,341,475	C→T	A899T (GCG→ACG)	yfaL ←	putative autotransporter adhesin YfaL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,341,475	0	C	T	97.8%	159.9 / NA	46	A899T (GCG→ACG)	yfaL	putative autotransporter adhesin YfaL
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/17);  total (29/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCGTGATCGCGTCGCTCCATCATAAACGCCTGATTTGCCGCCCGCAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGCTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATA  >  NC_000913/2341231‑2341715                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 gCGTGCTCGCGTCGATACATCATACACGACTGATTTGCCGCTCGCAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTACGAAGGCGTAGGTTCCGGGTCGGGCGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGCTCAGGTGGCGATGGCGGCGGTACGTCTTGCGATCGCAGATAACCGTCGTTGCTATCTTCCACCAGCGTGTAGTCATACGTTCCCat                                                                                                                                                                                                                                            >  2:169290/1‑251 (MQ=255)     gATCGCGTCGCTCCATCATAAACGCCTGATTTGCCGCCCGCAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTg                                                                                                                                                                                                                                        <  2:71358/251‑1 (MQ=255)         gcgTCGCTCCATCATAAACGCCTGATTTGCCGCCCGCAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGt                                                                                                                                                                                                                                    >  1:141360/1‑251 (MQ=255)           gTCGCTCCATCATAAACGCCTGATTTGCCGCCCGCAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAg                                                                                                                                                                                                                                  <  1:154964/251‑1 (MQ=255)                    tcatAAACGCCTGATTTGCCGCCCGCAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGctgccc                                                                                                                                                                                                                         >  2:46525/1‑250 (MQ=255)                              cTGATTTGCCGCCCGCAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTg                                                                                                                                                                                                                >  1:347497/1‑250 (MQ=255)                                gATTTGCCGCCCGCAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCAGATTGACGTAGCCGCTGCCGGCCAGACTGa                                                                                                                                                                                                               <  1:304715/249‑1 (MQ=255)                                       cgcccgcAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCAt                                                                                                                                                                                                      <  1:52286/251‑1 (MQ=255)                                            gcAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGATGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGttt                                                                                                                                                                                                 <  2:263479/251‑1 (MQ=255)                                                            aCCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGAt                                                                                                                                                                                 >  1:216565/1‑251 (MQ=255)                                                            aCCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATAAGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGAt                                                                                                                                                                                 >  1:216588/1‑251 (MQ=255)                                                             ccGCAAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATc                                                                                                                                                                                <  2:139416/251‑1 (MQ=255)                                                                 cAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCt                                                                                                                                                                            >  1:206644/1‑251 (MQ=255)                                                                  aaCTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTg                                                                                                                                                                           >  1:358019/1‑251 (MQ=255)                                                                      ttGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCAGACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGa                                                                                                                                                                        <  1:248572/250‑1 (MQ=255)                                                                        ggCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACAAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATAGGTTCCCATATTGAAGTAGCCGCTGCCTGCCAGCCTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGaaa                                                                                                                                                                      >  1:206860/1‑250 (MQ=255)                                                                              cAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAAcc                                                                                                                                                                <  1:33225/250‑1 (MQ=255)                                                                                  acCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACttt                                                                                                                                                           <  1:46525/251‑1 (MQ=255)                                                                                               cgggCGTAGGGTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCTCTTTAATGCCTGTCGAt                                                                                                                                              <  2:370441/249‑1 (MQ=255)                                                                                                gtgcttAGGTTCCGGGGCGGGTGGGGGATCAGGATCCGGCGTGGGATCAGGGTCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACGTCTTGCGGTCGCAGATACCAGTCGTTGTTATCTTCCCCCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGCGATCAGCTGCGAATTCACCCACTTTAATCCCTTTCGATg                                                                                                                                             <  2:73299/246‑1 (MQ=255)                                                                                                                          gATCAGGATCCGGCGTGGGATCAGGATCGGGCGTAGGGTCTGGATCAGGTGGCGCCAGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATgg                                                                                                                     >  1:32326/1‑249 (MQ=255)                                                                                                                           aTCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTTCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGt                                                                                                                  <  1:27118/251‑1 (MQ=255)                                                                                                                                 aTCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCAGATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCTCTTTAATTCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAacca                                                                                                            <  1:218025/251‑1 (MQ=255)                                                                                                                                    cGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAaccacc                                                                                                          >  1:43901/1‑250 (MQ=255)                                                                                                                                              tCAGGATCGGGAGTCGGGTCTGGATCCGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCCCAGTTGTGtt                                                                                               >  2:167267/1‑251 (MQ=255)                                                                                                                                                         aGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTAt                                                                                    <  2:97899/251‑1 (MQ=255)                                                                                                                                                              ggTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGGCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCg                                                                               <  2:25688/251‑1 (MQ=255)                                                                                                                                                                          gTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAAt                                                                   >  1:363233/1‑251 (MQ=255)                                                                                                                                                                               gATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATc                                                              >  2:119604/1‑251 (MQ=255)                                                                                                                                                                                       cgTTACTTCTTGCGATCGCAGATCCCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGt                                                                                                                                                                                                   >  1:321498/1‑110 (MQ=255)                                                                                                                                                                                       cgTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGt                                                                                                                                                                                                   <  2:321498/110‑1 (MQ=255)                                                                                                                                                                                       cgTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCg                                                      >  1:171123/1‑251 (MQ=255)                                                                                                                                                                                                  gCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCg                                           <  1:155552/251‑1 (MQ=255)                                                                                                                                                                                                          aGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTCGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTc                                   >  1:160288/1‑251 (MQ=255)                                                                                                                                                                                                          aGATACCAGTCGTTGTTATCTTCCACCAGCGGGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATAGCTTACCGAGTCATCGCCGTTTAATTc                                   >  2:210428/1‑251 (MQ=255)                                                                                                                                                                                                                     gttgttATCCTCCCCCAGCGAGTACACATACGCTCCCCTATTGACGTGGTCGCTGACCGCCTGACTGAACTACGCATTGTTTTACAACCGCCTGGGACCAGCTACGAACTCAACCACTTTACTGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGCGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGCCATCGCCGTTTAATTCGCTATCGAGc                         >  2:221483/1‑250 (MQ=255)                                                                                                                                                                                                                       tgttATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGt                      >  1:118616/1‑251 (MQ=255)                                                                                                                                                                                                                          tATCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGc                   >  2:281150/1‑251 (MQ=255)                                                                                                                                                                                                                           aTCTTCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTACTTCGCTATCGAGCAGTAGCg                  >  2:261268/1‑251 (MQ=255)                                                                                                                                                                                                                               tCCACCAGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTAcc              >  2:91063/1‑251 (MQ=255)                                                                                                                                                                                                                                   ccaGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTAccgccg          >  1:125468/1‑251 (MQ=255)                                                                                                                                                                                                                                    caGCGTGTAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTAccgccgc         >  2:109707/1‑251 (MQ=255)                                                                                                                                                                                                                                        gtgtAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGt     >  1:10193/1‑251 (MQ=255)                                                                                                                                                                                                                                         tgtAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGt     >  1:348352/1‑250 (MQ=255)                                                                                                                                                                                                                                         tgtAGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGt     >  2:360838/1‑250 (MQ=255)                                                                                                                                                                                                                                            aGTCATACGTTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCCGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGtata  >  2:190336/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 GCGTGATCGCGTCGCTCCATCATAAACGCCTGATTTGCCGCCCGCAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGATGGCGGCGTTACTTCTTGCGATCGCAGATACCAGTCGTTGTTATCTTCCACCAGCGTGTAGTCATACGCTCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATA  >  NC_000913/2341231‑2341715

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