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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,018,251	0	T	G	5.6%	205.8 / 2.9	72	S17R (AGC→CGC) ‡	ycbZ	putative peptidase
Reads supporting (aligned to +/- strand):  ref base T (30/38);  new base G (2/2);  total (32/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 = 4.21e-01

CTCTTCCGGGGCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAA  >  NC_000913/1018102‑1018400                                                                                                                                                      |                                                                                                                                                      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:226437/150‑1 (MQ=255)                                                                                          cACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAc                                                              <  1:31127/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                                        >  1:1132878/1‑150 (MQ=255)                                                                                                                gTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCacac                                        <  2:515272/150‑1 (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)                                                                                                                                                     cTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTGCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTa   <  1:1123828/150‑1 (MQ=255)                                                                                                                                                      tATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCATCGACAATACTCCTGGCATGGGCGTTaa  >  1:482939/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CTCTTCCGGGGCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAA  >  NC_000913/1018102‑1018400

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