Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,610,648	G→C	7.6%	M21I (ATG→ATC)	emrE →	small multidrug resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,610,648	0	G	C	7.6%	277.4 / 8.8	89	M21I (ATG→ATC)	emrE	small multidrug resistance protein
Reads supporting (aligned to +/- strand):  ref base G (39/43);  new base C (4/3);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGC  >  NC_002947/5610504‑5610793                                                                                                                                                 |                                                                                                                                                     cTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑gg                                                                                                                                                 >  2:19596/1‑149 (MQ=255) cTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑gg                                                                                                                                                 <  6:184476/149‑1 (MQ=255) cTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑gg                                                                                                                                                 >  7:91179/1‑149 (MQ=255) cTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑gg                                                                                                                                                 <  5:109586/149‑1 (MQ=255)  tATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCAATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGc                                                                                                                                                >  2:79167/1‑149 (MQ=255)       cccATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTTCATCGA‑GGCGGTgg                                                                                                                                           <  7:292676/149‑2 (MQ=255)       cccATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTTCATCGA‑GGCGGTgg                                                                                                                                           <  4:10182/149‑2 (MQ=255)        ccATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCgg                                                                                                                                              <  4:238674/145‑1 (MQ=255)           tatCGATACGTAATGACTTTACCCCATCAAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAgg                                                                                                                                       <  4:277942/149‑1 (MQ=255)             tCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTTCATCGA‑GGCGGTggcgatga                                                                                                                                     <  5:363693/149‑8 (MQ=255)              cGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCt                                                                                                                                    >  2:264807/1‑149 (MQ=255)               gATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCt                                                                                                                                    <  3:195008/148‑1 (MQ=255)                aTAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGa                                                                                                                                  <  5:45552/149‑1 (MQ=255)                 tAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGa                                                                                                                                  <  7:260714/148‑1 (MQ=255)                   ggTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGc                                                                                                                                <  5:165774/148‑1 (MQ=255)                    gTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCAc                                                                                                                              <  8:55028/149‑1 (MQ=255)                     tAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCt                                                                                                                                    >  2:169902/1‑142 (MQ=255)                     tAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCt                                                                                                                                    <  1:169902/142‑1 (MQ=255)                          aCTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGtt                                                                                                                        >  5:324922/1‑149 (MQ=255)                            tttACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑Gc                                                                                                                     >  6:324196/1‑149 (MQ=255)                            tttACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑Gc                                                                                                                     >  3:24082/1‑149 (MQ=255)                            tttACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATCAA‑GGCGGTGAAAGGCCTGAGCACGCCCTT‑Gc                                                                                                                     >  8:223507/1‑149 (MQ=255)                               aCCCCATCCAGGCAAGGTCTAGAATGCCCGCTGAGTAACGGAGGATGCTGCCTTGAATGCCTATACCTACCTCGGCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCa                                                                                                                   <  7:183028/148‑1 (MQ=255)                               aCCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCa                                                                                                                   <  7:188261/148‑1 (MQ=255)                               aCCCCATCCAGCCAACGTCTAGAATGCCCTCTCAGTAACGGAGGATCCTGCCATGAATGCCTCTACCTACCTCGCCAATGCCATCTGCGCCGAAGTCAGCGCCACCGCCTCTATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCa                                                                                                                   <  6:315797/148‑1 (MQ=255)                                       cAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCActgctgct                                                                                                           <  1:340457/148‑1 (MQ=255)                                        aGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGa                                                                                                         >  1:243955/1‑149 (MQ=255)                                          gCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGAtg                                                                                                       <  4:24082/149‑1 (MQ=255)                                            aaCGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGAtggt                                                                                                     >  1:203505/1‑149 (MQ=255)                                            aaCGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGAtggt                                                                                                     >  8:178525/1‑149 (MQ=255)                                             aCGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAACGGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGAtggt                                                                                                     <  1:264807/149‑1 (MQ=255)                                              cGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGAtggtgg                                                                                                   >  8:143793/1‑149 (MQ=255)                                                tCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGt‑t                                                                                                <  7:178525/149‑1 (MQ=255)                                                  tAGAAAGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGAT‑CCCACTGCTGCTGATGGTGGT‑Tgg                                                                                              >  8:92455/1‑149 (MQ=255)                                                   aGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑Tgg                                                                                              <  7:223507/148‑1 (MQ=255)                                                    gAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑Tg                                                                                               <  7:143793/146‑1 (MQ=255)                                                    gAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGc                                                                                             >  2:368041/1‑148 (MQ=255)                                                         cccGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATgg                                                                                        <  1:19596/148‑1 (MQ=255)                                                          ccGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCa                                                                                      >  2:127504/1‑149 (MQ=255)                                                                aGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCg                                                                                 <  8:91179/148‑1 (MQ=255)                                                                  tAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGNTNNTGGT‑TGGCTATGGCATTGCGTTc                                                                              >  8:279971/1‑149 (MQ=255)                                                                  tAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTc                                                                              >  8:275030/1‑149 (MQ=255)                                                                    aCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCt                                                                             <  2:243955/148‑1 (MQ=255)                                                                       gaggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGAt                                                                         <  5:263911/149‑1 (MQ=255)                                                                       gaggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGAt                                                                         >  3:221448/1‑149 (MQ=255)                                                                            tCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCAGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGa                                                                    <  4:282008/149‑1 (MQ=255)                                                                             ccTGCCATGAATGCCTATACCTACCTCGCCATTGTCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGAc                                                                   >  5:82214/1‑149 (MQ=255)                                                                             ccTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGAc                                                                   >  2:76707/1‑149 (MQ=255)                                                                               tGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGc                                                                 >  6:286002/1‑149 (MQ=255)                                                                               tGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGc                                                                 >  7:252023/1‑149 (MQ=255)                                                                                gCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCt                                                                >  1:301143/1‑149 (MQ=255)                                                                                 ccATGAATGCCTATACCTACCTCGCCATTGTCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGc                                                                 <  6:82214/147‑1 (MQ=255)                                                                                  cATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCtgg                                                              <  8:252023/149‑1 (MQ=255)                                                                                   aTGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCtggt                                                             <  6:324922/149‑1 (MQ=255)                                                                                    tGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtg                                                            <  5:343783/149‑1 (MQ=255)                                                                                    tGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtg                                                            <  1:79167/149‑1 (MQ=255)                                                                                     gAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtgg                                                           >  3:112055/1‑149 (MQ=255)                                                                                      caTGCCTATCCCTACCTCGCCATTCCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTggtggt                                                          <  4:45060/148‑1 (MQ=255)                                                                                           cTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgca                                                     <  2:131246/149‑1 (MQ=255)                                                                                           cTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCg                                                       >  8:81439/1‑147 (MQ=255)                                                                                              tacctaccTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGTCTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                                                    <  1:76707/148‑1 (MQ=255)                                                                                                 ctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTc                                               >  3:182425/1‑149 (MQ=255)                                                                                                 ctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTc                                               <  1:368041/149‑1 (MQ=255)                                                                                                    cctCGCCATTGCCCTCTGCGCCTAAGTGCTCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCg                                             <  6:18374/148‑1 (MQ=255)                                                                                                       cGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgg                                         <  7:298742/149‑1 (MQ=255)                                                                                                        gCCATTGCCATCTTCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgg                                         <  5:324196/148‑1 (MQ=255)                                                                                                        gCCATTGCCATCTGCGCCGAAGTCCTCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTggg                                        <  7:279971/149‑1 (MQ=255)                                                                                                        gCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTggg                                        <  7:275030/149‑1 (MQ=255)                                                                                                                       gcCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTTAGCCAATGCTGCTGAGGGTGGT‑TGGCTATGGAATTGCGTGCTGGATGATGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCa                          >  1:289618/1‑149 (MQ=255)                                                                                                                        ccGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATc                        >  1:246179/1‑149 (MQ=255)                                                                                                                             gTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTc                   >  3:151244/1‑149 (MQ=255)                                                                                                                                aTCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCggg                >  7:311022/1‑149 (MQ=255)                                                                                                                                aTCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCggg                >  1:197263/1‑149 (MQ=255)                                                                                                                                aTCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCggg                <  4:112055/149‑1 (MQ=255)                                                                                                                                 tCGCCACCGGCTCGATGAA‑GGCGGTGAAAGGCCTGCGCACGCCGTT‑GCCACTGCTGGTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGGATTCCGGTGGGGATTGCCTACGCCATCTGGTCgggg               <  7:92455/149‑1 (MQ=255)                                                                                                                                 tCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTggg                                        >  6:50357/1‑124 (MQ=255)                                                                                                                                 ctttcACCGCCTTCATCGA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCgggg               >  2:6062/5‑149 (MQ=255)                                                                                                                                  tttcACCGCCTTCATCGA‑GGCGGTGGAAGGCATGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGc              >  1:190789/4‑149 (MQ=255)                                                                                                                                  ttccACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGc              >  7:225840/3‑149 (MQ=255)                                                                                                                                  cGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCt                                >  4:186610/1‑131 (MQ=255)                                                                                                                                  cGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCt                                <  3:186610/131‑1 (MQ=255)                                                                                                                                    tcACCGCCTTCATCGA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTg            >  7:14390/2‑149 (MQ=255)                                                                                                                                    ccACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTg            >  6:264271/1‑149 (MQ=255)                                                                                                                                    ccACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTg            >  2:69735/1‑149 (MQ=255)                                                                                                                                       ccGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTggg          >  5:253992/1‑148 (MQ=255)                                                                                                                                          ccTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGAtt      <  7:84716/149‑1 (MQ=255)                                                                                                                                          ccTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGAtt      <  1:127504/149‑1 (MQ=255)                                                                                                                                              gATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGc  >  7:57328/1‑149 (MQ=255)                                                                                                                                              gATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGc  >  4:83561/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGC  >  NC_002947/5610504‑5610793

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