Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,181,292	T→A	8.1%	L66H (CTC→CAC)	yegV →	putative sugar kinase YegV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,181,292	0	T	A	8.1%	206.5 / 8.8	74	L66H (CTC→CAC)	yegV	putative sugar kinase YegV
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base A (2/4);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCC  >  NC_000913/2181144‑2181436                                                                                                                                                     |                                                                                                                                                 cAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTc                                                                                                                                                 >  1:945489/1‑150 (MQ=255)    tCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGAGGCGCTTTAACGCCACGGc                                                                                                                                              <  1:393400/150‑1 (MQ=255)    tCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGAGGCGCTTTAACGCCACGGc                                                                                                                                              <  2:1652281/150‑1 (MQ=255)    tCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGAGGCGCTTTAACGCCACGGc                                                                                                                                              <  2:1415749/150‑1 (MQ=255)    tCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGAGGCGCTTTAACGCCACGGc                                                                                                                                              <  2:1163487/150‑1 (MQ=255)      ggTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCAt                                                                                                                                            >  1:1949371/1‑150 (MQ=255)      ggTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCAt                                                                                                                                            >  2:367428/1‑150 (MQ=255)       gTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATc                                                                                                                                           <  1:395633/150‑1 (MQ=255)           tGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAg                                                                                                                                       >  1:773591/1‑150 (MQ=255)             gtcgtcGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCg                                                                                                                                     >  2:756042/1‑150 (MQ=255)             gtcgtcGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCg                                                                                                                                     >  1:1937253/1‑150 (MQ=255)             gtcgtcGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCg                                                                                                                                     >  2:1266074/1‑150 (MQ=255)               cgtcgGCGCGGCGGTCATTGTCGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCggg                                                                                                                                   <  2:1350011/150‑1 (MQ=255)               cgtcgGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCggg                                                                                                                                   <  1:711352/150‑1 (MQ=255)               cgtcgGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCggg                                                                                                                                   <  1:1354024/150‑1 (MQ=255)                   ggcgcggcgGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCg                                                                                                                                                <  1:1990792/133‑1 (MQ=255)                   ggcgcggcgGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCg                                                                                                                                                >  2:1990792/1‑133 (MQ=255)                           ggTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGcc                                                                                                                       <  2:374275/150‑1 (MQ=255)                           ggTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGcc                                                                                                                       >  2:672076/1‑150 (MQ=255)                            gTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCg                                                                                                                      <  2:1163466/150‑1 (MQ=255)                                ttGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGTATGCTCACTTAATATTGCCGTGGCGTTAAATCGTATAGGCATCGAAGCGGGTAATGTCTTGCCggttg                                                                                                                  <  2:991835/150‑5 (MQ=255)                                       gATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAggg                                                                                                           >  1:592824/1‑150 (MQ=255)                                       gATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAggg                                                                                                           >  2:931751/1‑150 (MQ=255)                                         tCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCg                                                                                                         >  1:1357617/1‑150 (MQ=255)                                           gCCGACACTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCgtg                                                                                                       >  1:1060885/1‑150 (MQ=255)                                            ccGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCgtgt                                                                                                      >  1:134827/1‑150 (MQ=255)                                               aCGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTggg                                                                                                   <  1:316022/150‑1 (MQ=255)                                                cGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGc                                                                                                  >  2:1075201/1‑150 (MQ=255)                                                   ttATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGga                                                                                               <  1:1658959/150‑1 (MQ=255)                                                       gCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAgatg                                                                                           >  2:1507189/1‑150 (MQ=255)                                                       gCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAgatg                                                                                           >  2:1924570/1‑150 (MQ=255)                                                        cccTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAgatga                                                                                          >  2:1621580/1‑150 (MQ=255)                                                            ccccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAAAGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCg                                                                                      >  2:1203474/1‑150 (MQ=255)                                                            actgacgCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCg                                                                                      <  1:1307651/144‑1 (MQ=255)                                                                  gCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCg                                                                                <  1:32366/150‑1 (MQ=255)                                                                   cGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCgg                                                                               <  2:1889650/150‑1 (MQ=255)                                                                   cGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCgg                                                                               <  1:918970/150‑1 (MQ=255)                                                                   cGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCgg                                                                               <  1:1507189/150‑1 (MQ=255)                                                                        gTGCGATTTCGTACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCTGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGc                                                                          <  1:1798069/150‑1 (MQ=255)                                                                               atCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAagag                                                                   >  1:1708723/1‑150 (MQ=255)                                                                                   aaCTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCGCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCt                                                               <  2:396948/150‑1 (MQ=255)                                                                                     cTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTa                                                             <  2:1714145/150‑1 (MQ=255)                                                                                       gAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAAt                                                           >  1:1812777/1‑150 (MQ=255)                                                                                       gAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAAt                                                           >  1:463387/1‑150 (MQ=255)                                                                                       gAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAAt                                                           >  2:830487/1‑150 (MQ=255)                                                                                         aaCAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCa                                                         >  1:1045585/1‑150 (MQ=255)                                                                                            agcagAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGcc                                                      <  2:870706/150‑1 (MQ=255)                                                                                            agcagAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGcc                                                      <  2:1019484/150‑1 (MQ=255)                                                                                             gcagAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCt                                                     >  2:18501/1‑150 (MQ=255)                                                                                              cagAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTg                                                    <  2:4534/150‑1 (MQ=255)                                                                                                 agCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATc                                                 <  1:1176791/150‑1 (MQ=255)                                                                                                  gCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCg                                                <  1:1266074/150‑1 (MQ=255)                                                                                                   cGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGa                                               >  1:1947412/1‑150 (MQ=255)                                                                                                     ttAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATa                                             <  2:945489/150‑1 (MQ=255)                                                                                                          gTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGcc                                        <  1:830487/150‑1 (MQ=255)                                                                                                          gTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGcc                                        >  1:1837532/1‑150 (MQ=255)                                                                                                            tGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGa                                      <  1:1727996/150‑1 (MQ=255)                                                                                                            tGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGa                                      <  2:1937253/150‑1 (MQ=255)                                                                                                            tGGCGGATGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGa                                      <  2:786036/150‑1 (MQ=255)                                                                                                                      gcACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAAc                            <  2:1949371/150‑1 (MQ=255)                                                                                                                       cACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACg                           >  1:1218362/1‑150 (MQ=255)                                                                                                                        aCTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACgg                          <  1:756042/150‑1 (MQ=255)                                                                                                                           gAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTg                       >  1:1886311/1‑150 (MQ=255)                                                                                                                               gatGCCGAGGCGCTTTAACGCCACGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGgtgt                   >  2:1460062/3‑150 (MQ=255)                                                                                                                                atGCCGAGGCGCTTTAACGCCACGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTc                  >  1:1711963/2‑150 (MQ=255)                                                                                                                                  gCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTg                >  1:783311/1‑150 (MQ=255)                                                                                                                                   ccGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTgg               <  1:931751/150‑1 (MQ=255)                                                                                                                                   ccGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTgg               >  1:134849/1‑150 (MQ=255)                                                                                                                                   ccGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTgg               >  1:1150436/1‑150 (MQ=255)                                                                                                                                     gTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGgcg             <  2:1708723/150‑1 (MQ=255)                                                                                                                                         cGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTgg         >  2:1288135/1‑150 (MQ=255)                                                                                                                                            tAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTg      <  2:463387/150‑1 (MQ=255)                                                                                                                                               aGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTAGCGCTGGTTGAGc   <  2:1357617/150‑1 (MQ=255)                                                                                                                                                gcgcCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGcc  <  1:1924570/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCC  >  NC_000913/2181144‑2181436

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