Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,856,334	C→A	9.1%	E199* (GAG→TAG) ‡	ydjH ←	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,856,334	0	C	A	9.1%	253.4 / 10.6	94	E199* (GAG→TAG) ‡	ydjH	putative kinase
Reads supporting (aligned to +/- strand):  ref base C (44/41);  new base A (3/6);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCCGGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTAC  >  NC_000913/1856192‑1856480                                                                                                                                               |                                                                                                                                                   gccgGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCaaaa                                                                                                                                             >  2:1474638/1‑150 (MQ=255)   cgGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAAtt                                                                                                                                           <  1:690970/150‑1 (MQ=255)    cgCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTAACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGACTCATTCCCGGTGAGTAATTTTGACTCGGCAAAATTa                                                                                                                                          >  1:890036/2‑150 (MQ=255)     gCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTac                                                                                                                                         <  1:643061/150‑2 (MQ=255)     gCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTac                                                                                                                                         <  1:1684451/150‑2 (MQ=255)     gCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCCTTTTCATCACCACCGTTTTTACGCCGCACTCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTac                                                                                                                                         <  1:1668194/150‑2 (MQ=255)      cACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAgg                                                                                                                                        <  2:642235/150‑1 (MQ=255)       aCCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTactc                                                                                                                                       <  2:92326/150‑4 (MQ=255)                   aTGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCg                                                                                                                           <  1:999322/150‑1 (MQ=255)                        accacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATa                                                                                                                      <  1:948599/150‑1 (MQ=255)                         ccacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTATGAAACAGATAATCGACATaa                                                                                                                     >  2:1329152/1‑150 (MQ=255)                         ccacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAAGTAACAGAGAATCGACATaa                                                                                                                     >  2:1254108/1‑150 (MQ=255)                          cacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAAc                                                                                                                    >  2:921606/1‑150 (MQ=255)                          cacGCTTGATAAAGCAGCCGACTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGACTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAAc                                                                                                                    <  1:1111306/150‑1 (MQ=255)                           acGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACt                                                                                                                   >  2:547830/1‑150 (MQ=255)                              cTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCaa                                                                                                                >  2:561535/1‑150 (MQ=255)                                tGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATg                                                                                                              <  2:1232589/150‑1 (MQ=255)                                 gATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGc                                                                                                             >  2:558366/1‑150 (MQ=255)                                  aTAAAGCAGCCGTCTTTACCCGTTTTAATCACCCCCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCt                                                                                                            >  1:1681688/1‑150 (MQ=255)                                  aTAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCt                                                                                                            >  1:1685075/1‑150 (MQ=255)                                   tAAAGCATCCGGCTTTACCCGTTTTAATCACCACCGTATTTACGCCGCACGCAAGAACGCAGTCAGCAATATCAACCAGGTTCTCTTAACAGATTAATAATTTTGCCTCGGCAAAATTATGAAAAACATAAAAGACATAAAAAAATGCtg                                                                                                           >  1:76651/1‑149 (MQ=255)                                   tAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCtt                                                                                                           >  2:60944/1‑150 (MQ=255)                                    aaaGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTc                                                                                                          <  2:1566626/150‑1 (MQ=255)                                       gcagcCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCa                                                                                                       <  2:780430/150‑1 (MQ=255)                                        cagcCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCaa                                                                                                      >  2:1608847/1‑150 (MQ=255)                                         agcCGTCTTTACCCGTTTAAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCaaa                                                                                                     <  1:431954/150‑1 (MQ=255)                                           ccGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAata                                                                                                   <  2:994088/150‑1 (MQ=255)                                           ccGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAata                                                                                                   <  1:141726/150‑1 (MQ=255)                                             gTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatc                                                                                                 >  2:24192/1‑150 (MQ=255)                                              tcattACCCGTTTTAATCACAACCGTTTTTACGCCGCACGCAAGAAAGAAGTCAGCAATTTCATCCAGTGTCTATTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGTAAACAGATAATCGACATAACTCAATGCTTCGCAAATatca                                                                                                >  1:150367/4‑150 (MQ=255)                                               cTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatcat                                                                                               <  2:363286/150‑1 (MQ=255)                                                    cccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAg                                                                                           >  2:570296/1‑149 (MQ=255)                                                    cccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAg                                                                                           <  1:570296/149‑1 (MQ=255)                                                    cccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGc                                                                                          >  1:1120525/1‑150 (MQ=255)                                                     ccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCg                                                                                         <  2:949228/150‑1 (MQ=255)                                                         tttAATCACCACCGTTTTTACGCCGCAGGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAACTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttc                                                                                     >  1:528625/1‑150 (MQ=255)                                                         tttAATCACCACAGTTTTTCCGCCGCAAGCAAGAAAGCATTCAGGAATTTCATCAAGTGTCTCTTTAGCTGTGAGTAATTTTTCGTCGTCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCTTCCAGCGTttc                                                                                     <  1:281952/150‑1 (MQ=255)                                                           tAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcat                                                                                   <  2:663268/150‑1 (MQ=255)                                                           tAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcat                                                                                   >  1:596567/1‑150 (MQ=255)                                                                accaccGTTTTTACTCCGCACGAAAGAAAGCAGTCATCAATTTCATCCAGTGTCGATATCCCGTTGAGTAAATTTGCATCGTCAAAATTAGTAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTGCATTCAAc                                                                              <  1:1221464/150‑1 (MQ=255)                                                                accaccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAAc                                                                              >  1:1599001/1‑150 (MQ=255)                                                                accaccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAAc                                                                              >  1:1303760/1‑150 (MQ=255)                                                                 ccaacGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAAcc                                                                             <  2:1286856/150‑1 (MQ=255)                                                                   accGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACcgc                                                                           <  1:1256322/150‑1 (MQ=255)                                                                   accGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACcgc                                                                           >  2:205584/1‑150 (MQ=255)                                                                         tttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTg                                                                     <  2:1425860/150‑1 (MQ=255)                                                                          ttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGa                                                                    <  1:876006/150‑1 (MQ=255)                                                                          ttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGa                                                                    >  1:60423/1‑150 (MQ=255)                                                                          ttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGa                                                                    <  1:596119/150‑1 (MQ=255)                                                                            aCGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATc                                                                  <  1:683951/150‑1 (MQ=255)                                                                             cgccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCa                                                                 <  2:479196/150‑1 (MQ=255)                                                                              gccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCat                                                                <  2:972127/150‑1 (MQ=255)                                                                              gccgcACGCAAGAAAGCAGTCAGCAATTTAATCCAGTGTCTCTATCCCGGTGAGTAATTTTGACTCGGCAATATAAGGAAACAGAAAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCaa                                                                >  1:1041761/1‑149 (MQ=255)                                                                               ccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTTCTCatg                                                               >  1:585195/1‑149 (MQ=255)                                                                               ccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCata                                                               >  1:1093900/1‑150 (MQ=255)                                                                               ccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCata                                                               <  1:227120/150‑1 (MQ=255)                                                                                   acgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCgg                                                           <  1:1486460/150‑1 (MQ=255)                                                                                    cgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGc                                                          >  2:405754/1‑150 (MQ=255)                                                                                      caAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGcac                                                        >  2:327810/1‑150 (MQ=255)                                                                                              cAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAatc                                                <  1:1474638/150‑1 (MQ=255)                                                                                                    gCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATAATATCGGCACAGATAATCATCTGc                                          >  1:36310/1‑150 (MQ=255)                                                                                                     cAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCACCCGCGGTTTGATCATATCGGCACAGATAATCATCTGcc                                         >  2:1548482/1‑150 (MQ=255)                                                                                                      aaTTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCg                                        <  2:890723/150‑1 (MQ=255)                                                                                                       aTTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATTCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCgg                                       >  1:1101247/1‑150 (MQ=255)                                                                                                       aTTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCTGTTTGATCATATCGGCACAGATAATCATCTGCCgg                                       >  1:350900/1‑150 (MQ=255)                                                                                                           cATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCtt                                   >  1:527085/1‑150 (MQ=255)                                                                                                            aTCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATAATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCttt                                  >  2:531041/1‑150 (MQ=255)                                                                                                             tCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTa                                 <  2:1093900/150‑1 (MQ=255)                                                                                                                 ttgtCTCTTTCCAGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCtt                             <  2:836306/149‑1 (MQ=255)                                                                                                                  tgtatatttCCCGGTGAGTAATTTTGCCTCTGCAAAATTAGGAATCAGATAATCGACATAACTCACTGCTTCGAAAATATCATCCAGCGTGGCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTg                            <  1:1066389/144‑1 (MQ=255)                                                                                                                  tgtCTCTTTCCCGGTTATTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTg                            <  1:1455359/150‑1 (MQ=255)                                                                                                                  tgtCTCTTTCCCGGTGAGTAATTTTGCCTCGTCAAAATTAGGAAACAGATAATCGACATAACTCATTGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATATTCATCTGCCGGGCTTTAGCTTg                            <  2:890036/150‑1 (MQ=255)                                                                                                                       ctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTaa                       >  2:1625632/1‑150 (MQ=255)                                                                                                                       ctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGGaa                       >  2:608590/1‑150 (MQ=255)                                                                                                                           cccTGTGAGTAATTTTGCCTCGGCAAAATTAGGAAAGAGATAATCGACATAACTCATTGCTTCGCAAATATCATCCATCGTTTCATTCAACCGCGGTTTGTTCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAAt                   <  2:150367/150‑1 (MQ=255)                                                                                                                           cccGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAAt                   <  1:725517/150‑1 (MQ=255)                                                                                                                             cGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAAttt                 >  1:1622283/1‑150 (MQ=255)                                                                                                                               gTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCt               <  1:205584/150‑1 (MQ=255)                                                                                                                                 taGTAATTTTGCCTCGGCAAAATTAGGAAAAAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGt             <  2:325277/149‑1 (MQ=255)                                                                                                                                 gAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGt             <  1:725559/150‑1 (MQ=255)                                                                                                                                  cctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCgg                                                           >  1:1465728/3‑103 (MQ=255)                                                                                                                                  cctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCgg                                                           <  2:1465728/101‑1 (MQ=255)                                                                                                                                  aGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGATTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCCTGCTTTAGCTTGCGTAAAAATTTCTGtt            >  2:1073203/1‑150 (MQ=255)                                                                                                                                   ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTa           >  1:1439500/2‑150 (MQ=255)                                                                                                                                   ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTa           >  1:1256721/2‑150 (MQ=255)                                                                                                                                      aTTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCAAAGATAATCATCTGCCGTCCTTTAGCTTGCGTAAAAATTACTGTTAgcg        >  1:1041424/1‑150 (MQ=255)                                                                                                                                          tGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCttt    <  1:921606/150‑1 (MQ=255)                                                                                                                                           gCCTCGGCAAAATTAGGAAACAGATAATCGATATATCTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTa   >  1:665636/1‑150 (MQ=255)                                                                                                                                           gCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTa   <  2:34662/150‑1 (MQ=255)                                                                                                                                            ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTATCTTGCGTAAAAATTTCTGTTAGCGCTTTAc  >  2:1355648/1‑150 (MQ=255)                                                                                                                                            ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTAc  <  1:742316/150‑1 (MQ=255)                                                                                                                                            ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTAc  >  2:739540/1‑150 (MQ=255)                                                                                                                                            ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTAc  >  1:1611661/1‑150 (MQ=255)                                                                                                                                            ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTAc  >  1:1618230/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                   GCCGGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTAC  >  NC_000913/1856192‑1856480

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