Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,906,609	G→C	8.3%	intergenic (+26/+34)	ygcG → / ← eno	TPM domain protein, putative phosphatase/enolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,906,609	0	G	C	8.3%	292.1 / 15.8	96	intergenic (+26/+34)	ygcG/eno	TPM domain protein, putative phosphatase/enolase
Reads supporting (aligned to +/- strand):  ref base G (38/49);  new base C (5/3);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.07e-01

AGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGA  >  NC_000913/2906466‑2906756                                                                                                                                                |                                                                                                                                                     agCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTg                                                                                                                                                <  1:1498236/150‑1 (MQ=255)      cTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑tt                                                                                                                                          <  1:630536/150‑1 (MQ=255)        ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCA‑tttt                                                                                                                                        <  2:1256882/150‑1 (MQ=255)        ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑tttt                                                                                                                                        <  1:785749/150‑1 (MQ=255)        ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑tttt                                                                                                                                        >  2:472353/1‑150 (MQ=255)         gTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑ttttt                                                                                                                                       <  2:1585960/150‑1 (MQ=255)         gTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGCTGGCGGCGCATCCGGCCGCTGTTAACTCAACCATAAAAAATGACAGCCCGGAGGCGGGCA‑ttttt                                                                                                                                       >  2:1100594/1‑150 (MQ=255)          tAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCA‑tttttt                                                                                                                                      <  2:1204387/150‑2 (MQ=255)          tAGAAGCTATCATTCAGACAACAGCGGTTCCGCTGGAGGCGCAGATCGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCA‑tttttt                                                                                                                                      <  1:94548/150‑2 (MQ=255)               gTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCa                                                                                                                                 >  1:638099/1‑150 (MQ=255)                   tCATTCAGACAACAGCGGTTCCGGTTGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCTGCCCGGAGGCTGGCA‑TTTTTAAATCAGATa                                                                                                                             <  2:341106/150‑1 (MQ=255)                    cATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATaa                                                                                                                            >  1:1461535/1‑150 (MQ=255)                       tCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAagt                                                                                                                         <  1:122743/150‑1 (MQ=255)                        cAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAagtc                                                                                                                        >  1:17257/1‑150 (MQ=255)                         aGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAagtca                                                                                                                       >  2:1103618/1‑150 (MQ=255)                            caacaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAagtcagtc                                                                                                                    >  1:1469172/1‑150 (MQ=255)                             aacaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCt                                                                                                                   >  1:482442/1‑150 (MQ=255)                             aacaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCt                                                                                                                   >  1:1176337/1‑150 (MQ=255)                              acaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCAGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCACTCtt                                                                                                                  >  1:1448848/1‑150 (MQ=255)                                 gCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATg                                                                                                               <  2:503460/150‑1 (MQ=255)                                     gtCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTg                                                                                                           <  2:701627/149‑1 (MQ=255)                                      tCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTgg                                                                                                          <  1:1570927/150‑1 (MQ=255)                                        cGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGcc                                                                                                        >  1:225567/1‑150 (MQ=255)                                         gCTGGAGGCTGAGATAGCTGCGGATTTAGCGTCGGGGGCGGTGATTCTGGAGGTGGCGGCGCATCCGGCGGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAATTCAGATAAAGTCAGTCTTATGCCTGGCCt                                                                                                       <  1:387153/150‑1 (MQ=255)                                          cTGGTGGCTCAGATAGCGGCGGTTTTAGCGGCGGGGGCTGTTCTTCTTGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATATAGTCAGTCTTATGCCTGGCCtt                                                                                                      <  1:1090520/150‑1 (MQ=255)                                          cTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATTCCTGGCCtt                                                                                                      >  1:1456449/1‑150 (MQ=255)                                            ggAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTg                                                                                                    <  2:534652/150‑1 (MQ=255)                                            ggAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTg                                                                                                    <  2:609415/150‑1 (MQ=255)                                            ggAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTg                                                                                                    <  2:541362/150‑1 (MQ=255)                                             gAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGa                                                                                                   <  1:898565/150‑1 (MQ=255)                                               ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGTGGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAt                                                                                                  <  2:1529684/150‑1 (MQ=255)                                               ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtc                                                                                                 <  1:1365007/150‑1 (MQ=255)                                               ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtc                                                                                                 <  2:17257/150‑1 (MQ=255)                                               ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtc                                                                                                 <  1:1709955/150‑1 (MQ=255)                                               ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCAT‑AAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtct                                                                                                >  2:1204339/1‑150 (MQ=255)                                                gCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCACTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtct                                                                                                <  2:139404/150‑1 (MQ=255)                                                 cTCAGATAGCGGCGGATTTAGGGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTAAGTCTTATGCCTGGCCTTTGAtctc                                                                                               <  2:182426/150‑1 (MQ=255)                                                 cTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtctc                                                                                               <  2:236035/150‑1 (MQ=255)                                                  tCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATctct                                                                                              <  2:301904/150‑1 (MQ=255)                                                  tCAGAACGCGGCGGAGTACGCGGCGGGGGCGGCTCTGCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATctct                                                                                              <  2:832669/150‑1 (MQ=255)                                                        aGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGa                                                                                        <  1:476063/150‑1 (MQ=255)                                                         gcggcgGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGAc                                                                                       <  1:1123357/150‑1 (MQ=255)                                                           ggcggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGGAACTCAACCATAAAAAATGCCAGACCGGAGCCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCGGGCCTTTGATCTCTTTACGACCg                                                                                     >  2:1103261/1‑150 (MQ=255)                                                               gATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGt                                                                                 >  1:1160034/1‑150 (MQ=255)                                                               gATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGt                                                                                 >  2:1281613/1‑150 (MQ=255)                                                                        gcggGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCtt                                                                        >  1:541582/1‑150 (MQ=255)                                                                           ggggCGGTTCTTCTTGAGGTGGAGGCTCATCCGGACGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCttttt                                                                     <  2:169882/150‑1 (MQ=255)                                                                            gggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTc                                                                    >  1:1398285/1‑150 (MQ=255)                                                                             ggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTACCTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCGCTTTACGACCGTTGTACGGTGCTTTTTCg                                                                   >  1:1136279/1‑150 (MQ=255)                                                                             ggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGGCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCg                                                                   <  1:203010/150‑1 (MQ=255)                                                                              gCGGTTCTTCTGGAGGTGGCGGCACATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGc                                                                  <  1:96487/150‑1 (MQ=255)                                                                               cGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGACA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGcc                                                                 >  2:1051140/1‑150 (MQ=255)                                                                                ggTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGccc                                                                <  2:1054534/150‑1 (MQ=255)                                                                                 gTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCa                                                               >  1:319264/1‑150 (MQ=255)                                                                                     ttctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGc                                                           >  2:986537/1‑150 (MQ=255)                                                                                         ggAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttct                                                       >  1:1284776/1‑150 (MQ=255)                                                                                            ggTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCg                                                    <  1:1025620/150‑1 (MQ=255)                                                                                              tGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGAt                                                  <  1:1062257/150‑1 (MQ=255)                                                                                                gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATAc                                                <  2:245577/150‑1 (MQ=255)                                                                                                gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATAc                                                <  2:1461535/150‑1 (MQ=255)                                                                                                 gggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACg                                               <  2:429277/149‑1 (MQ=255)                                                                                                 cggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACg                                               >  1:1002520/1‑150 (MQ=255)                                                                                                   gcgcATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGaa                                             <  1:1620722/150‑1 (MQ=255)                                                                                                      cATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGATTTTTCGCCCAGAGCTTCTTCGAAACGAATCa                                          >  2:1043140/1‑150 (MQ=255)                                                                                                       aTCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAg                                         <  1:929373/150‑1 (MQ=255)                                                                                                            gCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGt                                    >  1:1183723/1‑150 (MQ=255)                                                                                                            gCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGt                                    <  2:611000/150‑1 (MQ=255)                                                                                                            gCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGt                                    <  1:188324/150‑1 (MQ=255)                                                                                                              cGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTg                                  <  2:482442/150‑1 (MQ=255)                                                                                                              cGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTg                                  <  2:961340/150‑1 (MQ=255)                                                                                                              cGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTg                                  <  2:1176337/150‑1 (MQ=255)                                                                                                               gCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGt                                 >  2:686953/1‑150 (MQ=255)                                                                                                                 tGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTAt                               >  1:436307/1‑150 (MQ=255)                                                                                                                 tGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTAt                               <  1:244083/150‑1 (MQ=255)                                                                                                                  ggTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTAtt                              <  1:898555/150‑1 (MQ=255)                                                                                                                   gTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTAttt                             <  2:319264/150‑1 (MQ=255)                                                                                                                     aaCTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAg                           <  2:156491/150‑1 (MQ=255)                                                                                                                       cTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGGCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCa                         <  1:1103618/150‑1 (MQ=255)                                                                                                                        tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCaa                        >  1:1173512/1‑150 (MQ=255)                                                                                                                           aCCATAAAAAATGCCAGCCCGGAGGCTGTCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAacac                     <  1:1423109/150‑1 (MQ=255)                                                                                                                            gatttaaaaaTGCCAGCCTCCGGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTTGTTGTATTTAGCAACACg                    >  2:272826/6‑150 (MQ=255)                                                                                                                              aTAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGt                  >  2:396221/1‑150 (MQ=255)                                                                                                                               ttaaaaaTGCCAGCCTCCGGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTc                 >  2:108347/3‑150 (MQ=255)                                                                                                                                aaaaaaTGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCa                >  1:462774/1‑150 (MQ=255)                                                                                                                                 aaaaaTGCCAGCCTCCGGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTAAGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGCTCAg               >  2:848749/1‑150 (MQ=255)                                                                                                                                  aaaaTGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAAGCAGCTGGTTGTTTTTAGCAACACGGTCAGa              >  1:742577/1‑150 (MQ=255)                                                                                                                                    aaTGCCAGCCTCCGGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCACCACGGTCAGAAc            >  2:540971/1‑150 (MQ=255)                                                                                                                                     aTGCCAGCCTCCGGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACg           >  1:1667541/1‑150 (MQ=255)                                                                                                                                          aGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCGTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTGTTTCGCCCAGAGCTTAGTCGATACGAAGCAGGTGGTGGTTTTTAGCAACCGGGTCAGAACGGCTCa      >  2:285065/1‑150 (MQ=255)                                                                                                                                          aGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGc         >  1:1560777/1‑147 (MQ=255)                                                                                                                                          aGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGc         <  2:1560777/147‑1 (MQ=255)                                                                                                                                          aGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCa      >  1:1490932/1‑150 (MQ=255)                                                                                                                                            cccGGAGGCTGGCA‑TTTTTATATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATa    >  1:431753/1‑150 (MQ=255)                                                                                                                                             ccGTAGGCTGGAA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAg   <  2:1448848/150‑1 (MQ=255)                                                                                                                                             ccGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAg   >  1:198067/1‑150 (MQ=255)                                                                                                                                              cGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGAGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACTGTCAGAACGGCTCATAGa  >  2:1293035/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                     AGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGA  >  NC_000913/2906466‑2906756

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