Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,018,254	C→A	5.5%	D16Y (GAT→TAT)	ycbZ ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,018,254	0	C	A	5.5%	200.4 / 2.9	73	D16Y (GAT→TAT)	ycbZ	putative peptidase
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base A (2/2);  total (33/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.88e-01

TTCCGGGGCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAG  >  NC_000913/1018105‑1018402                                                                                                                                                      |                                                                                                                                                     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                                          >  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)                                                                                                                                                    aTCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGTGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTaaa   <  1:1093741/150‑1 (MQ=255)                                                                                                                                                     ccGGTATCAGGAACCAGGTCACTCCATGCAAGTTTCGTAATGGTCAAAGAAGATGTTTATTAGTCGGGTGTCAAAGCCGCGATTATACCGTAACCCGCACTACAGCACTCGTAGTCAGCACCGACAATACTCCTGGCATGGGCGTTAAAg  <  2:185081/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TTCCGGGGCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAG  >  NC_000913/1018105‑1018402

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: ‑