Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,781,645	0	A	C	42.4%	121.3 / 12.1	67	V361G (GTT→GGT)	ypjA	adhesin‑like autotransporter
Reads supporting (aligned to +/- strand):  ref base A (31/7);  new base C (0/28);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 2.29e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATATAATTGTTGTGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACATTGTCTGCCACG  >  NC_000913/2781497‑2781794                                                                                                                                                      |                                                                                                                                                      gaataaAATTGTTGTGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAAc                                                                                                                                                      <  1:498909/148‑1 (MQ=255)   atatAATTGTTGTGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAAcccc                                                                                                                                                   >  2:255215/1‑151 (MQ=255)    tataATTGTTGTGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCg                                                                                                                                                  >  2:108354/1‑151 (MQ=255)     ataATTGTTGTGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGc                                                                                                                                                 >  2:865735/1‑151 (MQ=255)      taATTGTTGTGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGcc                                                                                                                                                >  2:740895/1‑151 (MQ=255)           gttgtGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTAt                                                                                                                                           >  2:819356/1‑151 (MQ=255)           gttgtGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTAt                                                                                                                                           >  2:296954/1‑151 (MQ=255)                     tttccccGCCTTTTTTGATTTTTTGCTTTCCGTTTTTGATTTTGGTGCGTGTGTCTATGCCATAATTATTAATATTCTGTGTCCCCCCATTAATTATGT‑TTTTTGTCGCTTTCCCTGAAACACCCATACCCCCGCCATTATCTAtt                                                                                                                                      <  1:737635/142‑1 (MQ=255)                          ccgccgCTTTTAATTTTTTGCGTCCCGCTGTTTTTTTTGGTGCCTGGGGCTATGCAATAATTATTAATATTCTGTGTGCCACCTTTAATTTTGT‑TTTTTGCCCCTTTCCCTGCAACATCCAAACCCCCGCCATTATCTATTCGGGTCGCAt                                                                                                                            <  2:392726/151‑1 (MQ=255)                          ccgcccTTTTTGATTTTTTGCGTTCCGTTTTTAATATGGGTGCTTGTGGCTTTGCAATAATTATTAATTTTCTGGGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACTTCCAAACCCCCGCCATTATCTATTCGGGTCGCAt                                                                                                                            <  2:403023/151‑1 (MQ=255)                                 ttttGATTTTTTGCTTTCCGTTGTTGATATTGGTGCTTGTGGTTATGCCATATTTATTAATTTTCTGTGTCCCCCCATTAATTATGG‑TTTTTGTCGCTTTTCCTCCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGc                                                                                                                        <  1:266094/148‑1 (MQ=255)                                 ttttAATATTTGGGGTCCCGTTGTGGATATGGGTGCCTTTGGTTATGCCATATTTTTAAATTTCTTGTGT‑CCCCCCTTAATTTTGGTTTTTTGTCGCTTTCCCTCAAACACCCATACCCCCGCCATTATCTATTCGGGTCGCATCa                                                                                                                          <  1:494973/146‑1 (MQ=255)                                     gATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGc                                                                                                                  >  1:757387/1‑150 (MQ=255)                                      atatTTTGCGTTCCGCTTTTGATTTGGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACAATTAATTATGG‑TTTTTGTCGCTTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCa                                                                                                                 <  2:288651/150‑1 (MQ=255)                                         ttttGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTTTGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TTTTTGTCGCTTTTCCTGCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCAtt                                                                                                               <  2:534743/149‑1 (MQ=255)                                           ttGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGtt                                                                                                           >  2:381635/1‑151 (MQ=255)                                              ctttCCGTTGTTGATATTGGTGCCTGGGGTTATCCCATATTTATTAATTTTCTGTGTCCCACCTTTAATTATGG‑TTTTTGTCGCTTTTCCGGCAACACCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTaa                                                                                                         <  2:181251/148‑1 (MQ=255)                                                  ccGCTTTTGATATTGTTGCTTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTATTTATGG‑TTTTTGCCGTGTTTCCTCCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTg                                                                                                    <  1:282570/151‑1 (MQ=255)                                                   cGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTg                                                                                                    >  2:92709/1‑150 (MQ=255)                                                    gttGTTGATATGGGGGCGTTTGGCTATGCAATAATTATAAATTTTTTGTGTGCCACTATAATTTATGT‑TTTTTGTCGCTTTTCCTGAAACGTCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGa                                                                                                   <  1:442032/148‑1 (MQ=255)                                                      tGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACAt                                                                                                >  1:888553/1‑151 (MQ=255)                                                        ttGATATTGGTGCTTGTGGTTATGCCATAATTATTAATATTCTGTGTCCCACCATTATTTATGT‑TTTTTGTCGCTTTTCCTGCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACAtt                                                                                               <  1:783391/150‑1 (MQ=255)                                                           atatTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTc                                                                                           >  2:595540/1‑151 (MQ=255)                                                           atatTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTc                                                                                           >  2:447884/1‑151 (MQ=255)                                                           atatTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTc                                                                                           >  2:852126/1‑151 (MQ=255)                                                                    cttttgGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TTTTTTCCGCGTTCCCTGAAACATCCAACCCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTAt                                                                                     <  1:181632/144‑1 (MQ=255)                                                                    ccTGTGGTTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCTTTTCCTCCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATc                                                                                    <  2:926155/149‑1 (MQ=255)                                                                      tgtgGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCttt                                                                                 >  2:476366/1‑150 (MQ=255)                                                                         gcttATGCAATAATTATTAATATTTTGTGTGCCCCCTTAAATTATGG‑TTTTTGTCGCTTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCt                                                                                   <  1:522033/142‑1 (MQ=255)                                                                           cTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTaataat                                                                           >  2:810695/1‑151 (MQ=255)                                                                            tATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATc                                                                          >  2:564826/1‑151 (MQ=255)                                                                            tATGCAATAATTATTAATATTCTGTGTGCCACCATTAATTATGT‑TATTTGTCGCTTTTCCTGCAACATCAATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTtaataa                                                                            <  2:239110/149‑1 (MQ=255)                                                                              tGCTATATTTTTAAATATTCTGTGTGCCACCATTATTTATGG‑TATTTGTCGTGTTTCTTACACCATCAATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTtaataa                                                                            <  2:218581/147‑1 (MQ=255)                                                                              tGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGt                                                                        >  2:672792/1‑151 (MQ=255)                                                                              tGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGt                                                                        >  2:642165/1‑151 (MQ=255)                                                                              tGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGt                                                                        >  2:783162/1‑151 (MQ=255)                                                                                ccATAATTATTAATATTTTGTGGGCCACCATTAATTATGT‑TATTTGTCGCTTTTCCTCCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCt                                                                      <  2:540240/151‑1 (MQ=255)                                                                                ccATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCt                                                                      >  2:575015/1‑151 (MQ=255)                                                                                 cataaggtgtaATTTTCTGTGTGCCACCATTAATTGTGT‑TTTTTGTCGGTTTTCCTCCAACATCCAAACCCCCGCCAGTATCTATTCGGGTCGCATCAGCTTTAGCAGGGGTTAAAACTGACATTGTTCCGGTATCTTTAAGAATCGTc                                                                       <  1:410117/140‑1 (MQ=255)                                                                                  aTAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTTGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTc                                                                       >  1:59659/1‑148 (MQ=255)                                                                                    aaTTATTAATATTCTGTGTGCCCCAATTAATTAGGG‑TATTTGTCGCTTTTCCTGCAACATCCAAACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGt                                                                        <  2:594966/145‑1 (MQ=255)                                                                                    aaTTATAAATTTTTTGTGTCCCACAATTAATTATGG‑TTTTTGTCGCTTTTCCTGAAACATCAAAACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTg                                                                    <  1:322515/149‑1 (MQ=255)                                                                                      ttattaATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTg                                                                >  2:162104/1‑151 (MQ=255)                                                                                         ttaATATTCTGTGTCCCCCCATAAATTATGG‑TTTTTGTCGCTTTCCCTCCAACACCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGc                                                               <  2:860872/149‑1 (MQ=255)                                                                                           aaTTTTTTGTGTCCCCCATTTAATTATGT‑TTTTTGGCGCTTTTCCTGAAACTCCAATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGa                                                            <  1:683548/150‑1 (MQ=255)                                                                                             tatTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAAcc                                                         >  2:319795/1‑151 (MQ=255)                                                                                              atTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAAcc                                                         >  2:262080/1‑150 (MQ=255)                                                                                              atTCTGTGTGCCACCATAAATTATGT‑TTTTTGTCGCTTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAAcc                                                         <  2:580391/150‑1 (MQ=255)                                                                                               ttCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCTTTTCCTGCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAAcc                                                         <  2:253843/149‑1 (MQ=255)                                                                                                  tgtgtgCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATg                                                    >  2:103515/1‑151 (MQ=255)                                                                                                  tctttcccacaaTTAATTAAGA‑TATTTGTCATGTTTTCTGCACAACCCTTACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCat                                                       <  2:237924/142‑1 (MQ=255)                                                                                                   gtgtgCCACCATTAATTATGG‑TTTTTGTCGCTTTTCCTGCACCATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATg                                                    <  1:518276/150‑1 (MQ=255)                                                                                                     gtgCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGt                                                 >  2:156134/1‑151 (MQ=255)                                                                                                      tgCCACCATTAATTATGT‑TTTTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCatat                                                     <  2:413861/146‑1 (MQ=255)                                                                                                         caccaTTAATTATGG‑TATTTTTCGTGTTTCCGGCAACACCAATAGCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTAt                                             <  1:486143/151‑1 (MQ=255)                                                                                                            caTTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGt                                           >  2:292117/1‑150 (MQ=255)                                                                                                                aatttTGG‑TTTTTGTCGCTTTTCCTGCAACACCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGt                                           <  1:222603/146‑1 (MQ=255)                                                                                                                       gttttttGTGGCTTTCCCTACAACACCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGac                                  <  2:202928/145‑1 (MQ=255)                                                                                                                         tttttGTCGATTTTCCTGCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGaccac                               <  2:900316/148‑1 (MQ=255)                                                                                                                          aTTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCg                             >  1:5015/1‑151 (MQ=255)                                                                                                                               tCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTc                        >  2:708477/1‑151 (MQ=255)                                                                                                                                    tttCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGc                    >  1:126496/1‑150 (MQ=255)                                                                                                                                    tttCCTGCAACATCAAAACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAg                     <  2:872175/149‑1 (MQ=255)                                                                                                                                            accacccATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACAtt            <  1:925695/145‑1 (MQ=255)                                                                                                                                            accacccATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACAt             <  1:894886/144‑1 (MQ=255)                                                                                                                                            aaCATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACATTg           >  2:623528/1‑151 (MQ=255)                                                                                                                                                 ccATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACATTGTCt        <  1:512546/149‑1 (MQ=255)                                                                                                                                                      aCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGCCATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACATTGTCTGCCACg  >  1:177171/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GGATATAATTGTTGTGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGG‑TATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACATTGTCTGCCACG  >  NC_000913/2781497‑2781794

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