Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,018,249	G→C	5.6%	S17R (AGC→AGG) ‡	ycbZ ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,018,249	0	G	C	5.6%	202.6 / 4.7	71	S17R (AGC→AGG) ‡	ycbZ	putative peptidase
Reads supporting (aligned to +/- strand):  ref base G (29/37);  new base C (2/2);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

GACTCTTCCGGGGCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTT  >  NC_000913/1018100‑1018398                                                                                                                                                      |                                                                                                                                                      gACTCTTCCGGGGCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAg                                                                                                                                                       <  1:1079973/150‑1 (MQ=255)  aCTCTTCCGGGGCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGc                                                                                                                                                      <  2:287879/150‑1 (MQ=255)   ctctTCCGGGGCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCt                                                                                                                                                     >  1:691443/1‑150 (MQ=255)      ttCCGGGGCCTTCGCCAGCATAAAAGAGGAGGATCGTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATGACTGAATAATGGATCGTTTTCGTCTATCAAATGTGGCTGAGCTAATATTTCCTGATAGCTATc                                                                                                                                                  <  2:120507/150‑1 (MQ=255)       tCCGGGGCCTTCGCCAGCATAAAAGAGGAGGATGCGCGCGTATGCAGCAACTGCTCCAGCGCAAATAGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATCTTTCCTGATAGCTATCg                                                                                                                                                 <  1:505356/150‑1 (MQ=255)            ggCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATc                                                                                                                                            <  1:1134059/150‑1 (MQ=255)               cTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCATCGCAAATTGCAGGCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAgg                                                                                                                                         <  1:258312/150‑1 (MQ=255)                 tCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGaa                                                                                                                                       >  2:427005/1‑150 (MQ=255)                   gCCAGCATAAAAGAGGCGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGCAcc                                                                                                                                     >  2:779554/1‑150 (MQ=255)                   gCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAAcc                                                                                                                                     >  2:608535/1‑150 (MQ=255)                    ccAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCTAATATTTCCTGATAGCTATCGGTATCAGGAACCa                                                                                                                                    >  2:1080285/1‑150 (MQ=255)                    ccAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCa                                                                                                                                    <  2:809913/150‑1 (MQ=255)                    ccAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaata                                                                                                                                    <  2:568276/150‑4 (MQ=255)                             aaGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTCTCGGTATCAGGAACCAGGTCACGcc                                                                                                                           <  1:1245507/150‑1 (MQ=255)                              agcggaggaGTCTCGCGTATGCAGCAACTGCTGCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCa                                                                                                                          <  2:657938/147‑1 (MQ=255)                                   aggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCaa                                                                                                                     >  2:855433/1‑150 (MQ=255)                                   aggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCaa                                                                                                                     <  1:325934/150‑1 (MQ=255)                                      aTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGtt                                                                                                                  <  1:1299041/150‑1 (MQ=255)                                         cTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCg                                                                                                               <  1:855433/150‑1 (MQ=255)                                           cgcgTATGCAGCAACTGCTCCAGCGCAAATTGCCGCCGCGGTTGAGGATCACTGAATAAAGGATCGGTTTCGTCAAGCAAATGTGGCTGCGCAAATATTTCCTGATAGCTTTCGGTATCAGGAACCCGGTCACGCCATGCAAGTTTCGTa                                                                                                             >  2:164967/1‑150 (MQ=255)                                                aTGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGt                                                                                                        <  1:239185/150‑1 (MQ=255)                                                 tGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTc                                                                                                       <  2:416535/150‑1 (MQ=255)                                                    agcaACTGCTCCAGCGCAAATTGCAGCCGCGGGTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCTAGTTTCGTAATGGTCaaa                                                                                                    <  2:949090/150‑1 (MQ=255)                                                      caACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGt                                                                                                  >  2:1194195/1‑150 (MQ=255)                                                         cTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGa                                                                                               <  1:1299627/150‑1 (MQ=255)                                                          tGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGAt                                                                                              <  1:1038247/150‑1 (MQ=255)                                                             tCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtt                                                                                           >  1:515272/1‑150 (MQ=255)                                                              ccAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGttt                                                                                          >  2:955990/1‑150 (MQ=255)                                                               catcgcTAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCATATATTTCCTTATAGCTATCGGTATCAGGATCCAGGTCACGCCATGCAAGTTTCGTAATTGTCAAAGTTGATGtttt                                                                                         <  1:96532/147‑1 (MQ=255)                                                                 gcgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttttt                                                                                       <  2:283349/150‑1 (MQ=255)                                                                 gcgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttttt                                                                                       <  2:788199/150‑1 (MQ=255)                                                                 gcgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttttt                                                                                       >  1:1051093/1‑150 (MQ=255)                                                                  cgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTa                                                                                      >  2:698454/1‑150 (MQ=255)                                                                   gcAATTTGCAGACGAGGTTGAGTATCACTGAAAAAAGGATCGTTTTCGTCAATCAAATTTGGATGAGGAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCTAAGTTGATGTTTTTTAg                                                                                     <  2:231450/150‑1 (MQ=255)                                                                   gcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                                                                     <  2:128703/150‑1 (MQ=255)                                                                       aTTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCtg                                                                                 <  1:655548/150‑1 (MQ=255)                                                                              ccGCGGATGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCaa                                                                          <  1:761903/150‑1 (MQ=255)                                                                                     tGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAcgccgcc                                                                   >  1:1155534/1‑149 (MQ=255)                                                                                     tGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGAAAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgcg                                                                   <  1:815834/150‑1 (MQ=255)                                                                                       aGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGAt                                                                 >  2:557661/1‑150 (MQ=255)                                                                                        gTATCACTGATTAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGCTGTTTTTTAGTCTGTAGTCAAAGCCGCGAtt                                                                >  2:203995/1‑150 (MQ=255)                                                                                         tATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGAGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCACCAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTCTTGTCAAAGCCGCGATta                                                               >  2:318338/1‑150 (MQ=255)                                                                                           tCACTTAATAAAGGATCTTTTTAGTCAATCAAATGTGGCTGAGCATATATTTCCTGATAGCTATCGGTATCATGAACCAGGTCACGCCATGCAAGTTTCGTAATGGACAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtata                                                             <  1:217274/150‑1 (MQ=255)                                                                                            cACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAc                                                            <  1:31127/150‑1 (MQ=255)                                                                                            cACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAc                                                            <  1:226437/150‑1 (MQ=255)                                                                                               tGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGt                                                         <  2:778576/150‑1 (MQ=255)                                                                                               tGAATAAAGGATCGTTTTCGTCAATCAAATGT‑GCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTa                                                        <  1:198055/150‑1 (MQ=255)                                                                                                gAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACAGTa                                                        >  2:1113654/1‑150 (MQ=255)                                                                                                  aTAAAGGATCGTGATCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGGTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAAc                                                      <  2:204593/150‑1 (MQ=255)                                                                                                    aaaGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCg                                                    >  2:1178849/1‑150 (MQ=255)                                                                                                     aaGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATATCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCgg                                                   <  1:874702/150‑1 (MQ=255)                                                                                                            gTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGGTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACa                                            <  2:737832/150‑1 (MQ=255)                                                                                                              tttCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGc                                          <  2:1177703/150‑1 (MQ=255)                                                                                                                  gTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCacac                                      <  2:515272/150‑1 (MQ=255)                                                                                                                  gTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCacac                                      >  1:1132878/1‑150 (MQ=255)                                                                                                                    caatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGt                                    >  1:871000/1‑150 (MQ=255)                                                                                                                      atcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAg                                  >  1:1076145/1‑150 (MQ=255)                                                                                                                              gggCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCAcc                          <  1:291768/149‑1 (MQ=255)                                                                                                                               ggCTGAGCCAATATTTCCTGATAGCTAGCGGTATCAGGAACCATGTCACGCAATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTATAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCCCTACAGCACACGTAGAAAGCACCg                         >  2:169293/1‑150 (MQ=255)                                                                                                                                 cTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTTATGGTCAAACTTGATGTTTTTTAGTCTGTTGTCAAAGCAGCGATTATACCGAAACCGGCACTACAGCACACGTAGAAAGCACCGAc                       >  2:49757/1‑150 (MQ=255)                                                                                                                                  tGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACa                      >  1:1148942/1‑150 (MQ=255)                                                                                                                                    aGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAAt                    >  2:275504/1‑150 (MQ=255)                                                                                                                                    aGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTATTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGACTGCACCGACAAt                    >  1:801481/1‑150 (MQ=255)                                                                                                                                      cAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATAc                  <  1:608535/150‑1 (MQ=255)                                                                                                                                        aaTATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTc                >  2:1007958/1‑150 (MQ=255)                                                                                                                                          tatTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCt              >  2:330727/1‑150 (MQ=255)                                                                                                                                            tttaatGATAGCTTTCGGTATCATGAACCAGTTCACGACATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTTTTGTCAAAGCCGCGTTTATACGGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTgg            <  2:110573/145‑1 (MQ=255)                                                                                                                                            gttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTgg            >  2:1265830/2‑150 (MQ=255)                                                                                                                                             ttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAACGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGc           >  2:690724/1‑150 (MQ=255)                                                                                                                                              tCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCa          <  2:520202/150‑1 (MQ=255)                                                                                                                                               ccTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCAACAATACTCCTGGCAt         >  1:1200774/1‑150 (MQ=255)                                                                                                                                                    tAGCTATCGGTAGCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATTGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCg    <  2:688252/150‑1 (MQ=255)                                                                                                                                                      gCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGCCAATACTCCTGGCATGGGCGtt  >  1:673989/1‑150 (MQ=255)                                                                                                                                                      gCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGtt  >  1:1150207/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GACTCTTCCGGGGCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTT  >  NC_000913/1018100‑1018398

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