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

GGATATAATTGTTGTGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACATTGTCTGCCACG  >  NC_000913/2781497‑2781794                                                                                                                                                     |                                                                                                                                                      ggaTATAATTGTTGTGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAAcc                                                                                                                                                     <  1:377535/151‑1 (MQ=255)           gttgtGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTAt                                                                                                                                           >  1:87252/1‑151 (MQ=255)                 tCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTc                                                                                                                                     >  1:776025/1‑151 (MQ=255)                     tttccccGCCGTTTTTGATATTTGGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTCCCACCATTAATTATGGTATTTGTCGCGTTTCCGGCAACACCCATACCCCCGCCAt                                                                                                                                              <  1:314562/134‑1 (MQ=255)                       ttCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCg                                                                                                                               >  1:1021433/1‑151 (MQ=255)                        ccccGCCCTTTTTGATATTTTGCTTCCCGTTTTTAATATTGGTGCTGGTGGTTATGCAATATTTATAAATATTCTTTGTCCCACCATTAATTTTGTTTTTTGTCGCTTTTCCTGAAACACCAATACCCCCGCCATTATCTATTCGGGTCGc                                                                                                                              <  1:957293/150‑1 (MQ=255)                             cggtttttGATATTTTGCGTCCCGTTGTTGATATTGGTCCCTGTGGCAATGCAATAATTATTAATATTCTGTGTGCCACAATTAATTTTGGTTTTTGTCGATTTTCTGGAAACATCTAGACCCCCGCCATTATCTATTCGGGTCGCATCAg                                                                                                                         <  2:761724/147‑1 (MQ=255)                                   ttaattttttGCGTCCCGTTGTTGTTTTTGGTGCCTGTGGTTATGCCATATTTATTAATATTCTGTGTGCCCCCTTTAATTATGGTTTTTGTCGGTTTTCCGGCAACATCCAAACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTa                                                                                                                    <  2:459917/144‑1 (MQ=255)                                       tatTTTGCTTTCCCCTGTTGATTTTGGTGCCTGTGGCTATGCCATATTTTTTAATATTCTGTGTCCCCCCATAAATTATGTTTTTTGTCGCTTTTCCTCCACCACCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCAt                                                                                                                <  2:469180/150‑1 (MQ=255)                                         ttttGCGTCCCGTTTTTGATATTGGTCCTTGTGGTTATGCAATAATTATTAATATTCTGTGTGCCACCATTATTTATGTTTTTTGTCGCGTTCCCTCCAACATCCAAACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTg                                                                                                              <  1:895635/150‑1 (MQ=255)                                            tGCTTTCCGTTTTTGATTTTGTTGCTTTTGGCTTTGCCTTAATTATTAATTTTCTGTGTCCCCCATTTATTTATGTTTTTTGTCGCTTTTCCTGCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGtt                                                                                                           <  1:524060/150‑1 (MQ=255)                                             gCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTa                                                                                                          >  2:228996/1‑150 (MQ=255)                                                  ccgctttctatatGGGCGCGTGTGGCTGTGCTATAATTCTTAATATTCTGGCGCCCCCCATGAATTATGCAATTTGTCGAGTTCCCTCAAACCTCCTTAGCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACt                                                                                                     <  1:302328/141‑1 (MQ=255)                                                  ccGCTGTTGATTTGGGTCCCTGGGGCTTTCCAATAATTATTAATATTCTGTGTGCCCCATTTAATTATGGTTTTTTTCGCTTTTCCTCCAACACCAATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTg                                                                                                    <  1:638127/151‑1 (MQ=255)                                                    gCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGAc                                                                                                  >  2:1080204/1‑151 (MQ=255)                                                        ttGATATTGTTGCCTGTGGCTATGCAATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTTTTTGTCGCTTTTCCTCAACCATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACAtt                                                                                               <  1:803196/150‑1 (MQ=255)                                                        ttGATATTGGTGCCTGTGTCTATGCCATAATTATTAATTTTCTGTGTCCCACCATTAATTATGTTATTTGTCGCTTTTCCTGCAACACCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACAtt                                                                                               <  2:538907/150‑1 (MQ=255)                                                         tGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGt                                                                                             >  2:627750/1‑151 (MQ=255)                                                              ttGGTGCCTGTGGTTCTGCCATAATTATTAATTTTCTGTGTCCCCCTATAAATTTTGTTTTTTGTCGCTTTTCCTGAAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGtt                                                                                            <  1:592962/147‑1 (MQ=255)                                                              tgggTGCTTGTGGCTTTCCAATATTTATAAATTTTCTGTGTGCCACCATTAATTATGGTTTTTGTCGCTTTCCCTGCACCATCCAAACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCt                                                                                         <  1:63101/148‑1 (MQ=255)                                                                ggTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTAGGGTTTTTGTCGCTTTTCCTGCAACATCCAAAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCtt                                                                                        <  2:618892/149‑1 (MQ=255)                                                                    cttgtgGTTATGCCATATTTATAAATATTCTGTGTCCCACTTTTATTTATGTTTTTTTTCGCTTTCCCTGCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTAt                                                                                     <  1:833057/146‑1 (MQ=255)                                                                    ccTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCtt                                                                                  >  2:1011317/1‑151 (MQ=255)                                                                    ccTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCtt                                                                                  >  2:777718/1‑151 (MQ=255)                                                                    ccTGTGGCTATGCCATAATTATTAATATTCTGTGTCCCACCATTAATTATGGTTTTTGTCGCTTTTCCTGCACCATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCtt                                                                                  <  2:990864/151‑1 (MQ=255)                                                                    ccTGTGGCTATCCCATAATTATTAATATTCTGTGTGCCACAATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCAAAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATc                                                                                    <  2:417795/149‑1 (MQ=255)                                                                              tGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGt                                                                        >  2:738953/1‑151 (MQ=255)                                                                                      ttattaATTTTCTGTGTCCCCCCATTATTTATGTTTTTTGTCGGGTTTCCTCCAAAATCCAAACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTg                                                                    <  1:224830/147‑1 (MQ=255)                                                                                      ttattaATATTTTGTGTCCCCCCTTTAATTATGTTTTTTGTCGCTTTCCCTGCAACTTCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGt                                                                        <  1:618769/143‑1 (MQ=255)                                                                                      ttattaATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTg                                                                >  1:508874/1‑151 (MQ=255)                                                                                         ttaATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCg                                                             >  1:91501/1‑151 (MQ=255)                                                                                            atatTCTGTGTGCCACCATTAATTATGTTATTTGTCGGTTTTCCTGCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAAc                                                          <  2:80862/151‑1 (MQ=255)                                                                                                tCTGTGTGCCACCATTAATTATGTTATTTGTCGCTTTTCCTGGAACATCAAAAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCat                                                       <  2:808446/150‑1 (MQ=255)                                                                                                tCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCat                                                       >  2:569035/1‑150 (MQ=255)                                                                                                tCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCat                                                       >  1:179718/1‑150 (MQ=255)                                                                                                     gtgCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGt                                                 >  2:1057576/1‑151 (MQ=255)                                                                                                      tgCCACCATTATTTATGTTTTTTTTCGCGTTTCCTGCAACACCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCata                                                      <  1:913594/145‑1 (MQ=255)                                                                                                      tgCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGtt                                                >  2:1011727/1‑151 (MQ=255)                                                                                                        ccaccaTTAATTATGTTATTTGTCGCTTTTCCTGCAACTTCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGttt                                               <  2:194463/150‑1 (MQ=255)                                                                                                        ccacaaTTATTTATGTTTTTTGTCGCTTTTCCCGCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGtt                                                <  1:141946/149‑1 (MQ=255)                                                                                                         caccaTTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTAt                                             >  1:841586/1‑151 (MQ=255)                                                                                                          accaTTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATg                                            >  1:1017617/1‑151 (MQ=255)                                                                                                           ccaTTAATTATGTTATTTGTCGCTTTTCCTGCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTAt                                             <  2:136060/149‑1 (MQ=255)                                                                                                           ccaTCAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGt                                           >  1:308243/1‑151 (MQ=255)                                                                                                           catttaattaTGGTTTTTGTCGCTTTTCCTCCAACACCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATg                                            <  2:921939/147‑1 (MQ=255)                                                                                                           caattatttaTGTTTTTTTTCGCTTTTCCTAAAACATCTAAAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTAt                                             <  1:33515/147‑1 (MQ=255)                                                                                                             attaattaTGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCt                                         >  2:693653/1‑151 (MQ=255)                                                                                                                 attaTGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAAt                                     >  2:931536/1‑151 (MQ=255)                                                                                                                  ttaTGTTTTTTGTCGCTTTTCCTCCAACTCCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTaaa                                      <  2:570101/149‑1 (MQ=255)                                                                                                                     tGGTATTTGTCGTTTTCCCTGCAACACCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGac                                  <  2:101019/150‑1 (MQ=255)                                                                                                                      ggTTTTTGTCGCTTTTCCTGCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGacc                                 <  2:361039/150‑1 (MQ=255)                                                                                                                        tttttGTCGCTTTTCCGCCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGac                                  <  1:561690/145‑1 (MQ=255)                                                                                                                        tATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGACCCATATGCGTTTATGTCTAAATGaccacc                              >  2:373577/1‑151 (MQ=255)                                                                                                                         ttttttgggttttCCGCCATCATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGa                                   <  1:686595/135‑1 (MQ=255)                                                                                                                         aTTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCg                             >  1:942312/1‑151 (MQ=255)                                                                                                                          tttGTCGCTTTTCCTGCAACATCAAAACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCg                             <  2:106788/150‑1 (MQ=255)                                                                                                                           ttGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGtt                           >  2:85606/1‑151 (MQ=255)                                                                                                                               cgcgTTTCCTGCAACACCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGttt                          <  2:88014/148‑1 (MQ=255)                                                                                                                                gcgTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTcc                       >  1:400743/1‑150 (MQ=255)                                                                                                                                  gTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCa                      >  1:1096844/1‑149 (MQ=255)                                                                                                                                   tttCCTGCAACATCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAg                     <  2:695003/149‑1 (MQ=255)                                                                                                                                   tttCCTGAAACTCCCATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCa                      <  1:780286/148‑1 (MQ=255)                                                                                                                                       cTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAacac               >  2:138248/1‑151 (MQ=255)                                                                                                                                        tGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAacac               >  2:917482/1‑150 (MQ=255)                                                                                                                                           aaCATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACATTg           >  2:672792/1‑151 (MQ=255)                                                                                                                                                ccATACCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACATTGTCt        <  1:762206/149‑1 (MQ=255)                                                                                                                                                ccATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACATTGTCTGc      >  2:390100/1‑151 (MQ=255)                                                                                                                                                 cATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACATTGTCTGcc     >  2:1029667/1‑151 (MQ=255)                                                                                                                                                    aaCCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACATTGTCTGCCACg  >  1:167591/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GGATATAATTGTTGTGTCAGCTTTCCCGCCGCTTTTGATATTTTGCGTTCCGCTGTTGATATTGGTGCCTGTGGCTATGCCATAATTATTAATATTCTGTGTGCCACCATTAATTATGGTATTTGTCGCGTTTCCTGCAACATCCATAACCCCGCCATTATCTATTCGGGTCGCATCAGCTTTAGCATTGGTTAAAACTGACATTGTTCCTTTATCTTTAATAATCGTCTTGTTTGCCGAACCATATGCGTTTATGTCTAAATGACCACCGTTTTCCAGCAACACATTGTCTGCCACG  >  NC_000913/2781497‑2781794

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