Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,274,920	T→C	5.3%	R379R (CGT→CGC)	ychO →	putative invasin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,274,920	0	T	C	5.3%	226.8 / 3.2	75	R379R (CGT→CGC)	ychO	putative invasin
Reads supporting (aligned to +/- strand):  ref base T (34/37);  new base C (2/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

GTAGTCGCTATGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTC  >  NC_000913/1274775‑1275066                                                                                                                                                  |                                                                                                                                                   gTAGTCGCTATGACAATCCGCAGCGAAATTCGCTAGCGACTCTTGAGTACCGAGAGCGTAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAAACTGGCGAAACAGTGCCGCTTAAATTACAAATCCGCATTCGTTa                                                                                                                                                  <  6:216689/148‑1 (MQ=255)   aGTCGCTATGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACgg                                                                                                                                               >  2:197584/1‑149 (MQ=255)    gTCGCTATGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAATAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAACCCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACgg                                                                                                                                               <  3:101146/148‑1 (MQ=255)    gTCGCTATGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACgg                                                                                                                                               <  5:362785/148‑1 (MQ=255)         tATGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGGTAACGGTGTTTCTGGCGCCACCGCCGTGGGATCTAAAAACTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCg                                                                                                                                         >  8:392555/1‑149 (MQ=255)          aTGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCgg                                                                                                                                        >  7:258618/1‑149 (MQ=255)          aTGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCgg                                                                                                                                        >  4:88553/1‑149 (MQ=255)          aTGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCCCAGTCGTTACGGTATTCgg                                                                                                                                        >  4:83632/1‑149 (MQ=255)           tGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGc                                                                                                                                       <  4:86/149‑1 (MQ=255)           tGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAAATACAAATCCGCAGTCGTTACGGTATTCgg                                                                                                                                        >  2:301921/1‑148 (MQ=255)               aaTCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACt                                                                                                                                   >  8:264927/1‑149 (MQ=255)                 tCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGa                                                                                                                                 <  3:82071/149‑1 (MQ=255)                    gCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGAttt                                                                                                                              <  7:175014/149‑1 (MQ=255)                    gCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGAttt                                                                                                                              >  4:345981/1‑149 (MQ=255)                      aGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTTGTATCTAAAACCTGGCGAAACAGTGCCGCTGATATTACAAATAAGCAGTCGTTACTGTATTCGGCAACTGATTTgg                                                                                                                            <  2:151167/149‑1 (MQ=255)                      aGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTgg                                                                                                                            <  8:258618/149‑1 (MQ=255)                      aGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTgg                                                                                                                            <  1:197584/149‑1 (MQ=255)                        cGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGGGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCa                                                                                                                          <  4:247705/149‑1 (MQ=255)                          aaaTAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAgg                                                                                                                        >  5:321321/1‑149 (MQ=255)                               atCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGa                                                                                                                    <  7:335494/148‑1 (MQ=255)                               atCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGa                                                                                                                    >  3:31976/1‑148 (MQ=255)                               atCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGAt                                                                                                                   <  6:5651/149‑1 (MQ=255)                                     agACTATTGAGTACCGACAGCGAAAAACGTTAACGGTTTTTCTGGCGAAACCGACGTGTGATATAAAACCTGGCGAAACAGTGCCTCTGAAATTACAAATACGCAGTCGTTACGGTCTTCGGCAACTGATTTGGCAGGGCGATACGCAg                                                                                                             <  6:273464/148‑1 (MQ=255)                                       aCTCTTGAGTACCGCCAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGGGGGATCTAAAACCTGGCGAAACCGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGa                                                                                                            >  6:122479/1‑148 (MQ=255)                                             gAGTACCGACAGCGAAAAACGTTAACGGTGTATCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAg                                                                                                     >  6:33858/1‑149 (MQ=255)                                                  ccGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGa                                                                                                <  8:23760/149‑1 (MQ=255)                                                   cGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGCTACGGTATTCgg                                                                                                                                        >  3:291680/1‑108 (MQ=255)                                                   cGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGCTACGGTATTCgg                                                                                                                                        >  3:296685/1‑108 (MQ=255)                                                   cGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGCTACGGTATTCgg                                                                                                                                        <  4:296685/108‑1 (MQ=255)                                                   cGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGCTACGGTATTCgg                                                                                                                                        <  4:291680/108‑1 (MQ=255)                                                        gCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCgg                                                                                          <  8:342885/149‑1 (MQ=255)                                                         cGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATTCGCAGATATTAAGTTTGACGCCggg                                                                                         >  7:394312/1‑149 (MQ=255)                                                         cGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCggg                                                                                         >  8:293678/1‑149 (MQ=255)                                                          gAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGgc                                                                                        <  8:244552/149‑1 (MQ=255)                                                            aaaaCGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATCCGGCAACTGAGTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGgcgc                                                                                      >  6:72358/1‑149 (MQ=255)                                                             aaaCGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGt                                                                                                    >  1:175576/1‑134 (MQ=255)                                                             aaaCGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGt                                                                                                    <  2:175576/134‑1 (MQ=255)                                                               aCGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGcac                                                                                    >  6:187518/1‑148 (MQ=255)                                                                cGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACaa                                                                                  <  4:31976/149‑1 (MQ=255)                                                                  ttAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGAGATTAAGTTTGACGCCGGGCGCACAAGc                                                                                >  3:5929/1‑149 (MQ=255)                                                                     aCGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCa                                                                              <  7:89107/148‑1 (MQ=255)                                                                       ggTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAGACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCaa                                                                             >  7:143006/1‑147 (MQ=255)                                                                       ggTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTATGGCATTc                                                                                                                                          <  7:102163/86‑1 (MQ=255)                                                                       ggTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTAGGGAATTc                                                                                                                                          >  8:102163/1‑86 (MQ=255)                                                                        gtgtTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAg                                                                          >  1:9843/1‑149 (MQ=255)                                                                        gtgtTTCTGGCGACACCGCCGGGGGATCTAAAACCTGGCGAAACAGGGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAg                                                                          <  5:309806/149‑1 (MQ=255)                                                                               tGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAg                                                                   >  8:108907/1‑149 (MQ=255)                                                                                ggCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAgg                                                                  >  4:251387/1‑149 (MQ=255)                                                                                 gCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGa                                                                    <  4:115343/146‑1 (MQ=255)                                                                                 gCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGa                                                                    <  4:5929/146‑1 (MQ=255)                                                                                 gCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAggg                                                                 >  3:326877/1‑149 (MQ=255)                                                                                   gACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCt                                                               <  7:293678/149‑1 (MQ=255)                                                                                    acacCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTg                                                              >  8:84044/1‑149 (MQ=255)                                                                                            gTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGa                                                      <  3:88553/149‑1 (MQ=255)                                                                                            gTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGa                                                      <  3:83632/149‑1 (MQ=255)                                                                                               ggATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATc                                                    <  1:301921/148‑1 (MQ=255)                                                                                                gATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCAt                                                  >  4:356674/1‑149 (MQ=255)                                                                                                 aTCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCATCTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATg                                                 >  3:316733/1‑149 (MQ=255)                                                                                                       aaCCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGAc                                           <  7:264927/149‑1 (MQ=255)                                                                                                         ccTGGCGAGACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTg                                         <  8:143006/149‑1 (MQ=255)                                                                                                          cTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTgg                                        <  3:345981/149‑1 (MQ=255)                                                                                                          cTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTgg                                        >  7:369583/1‑149 (MQ=255)                                                                                                          cTGGCGAAACAGTGCCGCGGAAATTACAAATCCGCAGTCGTTACGGTATTCCGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTgg                                        <  6:283170/149‑1 (MQ=255)                                                                                                            ggCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTgg                                        >  8:58210/1‑147 (MQ=255)                                                                                                               gAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGaa                                   <  7:58210/149‑1 (MQ=255)                                                                                                                aaaCAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGAGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGAGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAAc                                  >  1:84457/1‑149 (MQ=255)                                                                                                                aaaCAGTGCCGCTGAAATTACAAATCCGCAGGCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAAc                                  <  7:108907/149‑1 (MQ=255)                                                                                                                      tGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGaa                            <  5:187518/149‑1 (MQ=255)                                                                                                                        ccGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAg                           <  2:9843/148‑1 (MQ=255)                                                                                                                         cGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAggg                         <  5:33858/149‑1 (MQ=255)                                                                                                                                 ttACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGgcaagcaa                 >  3:127832/1‑149 (MQ=255)                                                                                                                                  tACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAAt                >  5:25794/1‑149 (MQ=255)                                                                                                                                           gCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGCACGGGGAAGGGGCAAGCAATCACTGGCGa       >  6:269851/1‑149 (MQ=255)                                                                                                                                           gCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGa       >  2:136958/1‑149 (MQ=255)                                                                                                                                                 gTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTc  <  6:321321/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GTAGTCGCTATGACAATCCGCAGCGAAATAATCTACCGACTCTTGAGTACCGACAGCGAAAAACGTTAACGGTGTTTCTGGCGACACCGCCGTGGGATCTAAAACCTGGCGAAACAGTGCCGCTGAAATTACAAATCCGCAGTCGTTACGGTATTCGGCAACTGATTTGGCAGGGCGATACGCAGATATTAAGTTTGACGCCGGGCGCACAAGCCAACAGCGCGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTC  >  NC_000913/1274775‑1275066

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