Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,610,649	A→G	6.7%	K22E (AAG→GAG)	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,649	0	A	G	6.7%	262.4 / 6.8	89	K22E (AAG→GAG)	emrE	small multidrug resistance protein
Reads supporting (aligned to +/- strand):  ref base A (40/43);  new base G (3/3);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

CTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGA  >  NC_002947/5610504‑5610796                                                                                                                                                  |                                                                                                                                                       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                                                                            >  3:221448/1‑149 (MQ=255)                                                                       gaggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGAt                                                                            <  5:263911/149‑1 (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                                                  <  1:368041/149‑1 (MQ=255)                                                                                                 ctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTc                                                  >  3:182425/1‑149 (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                   >  1:197263/1‑149 (MQ=255)                                                                                                                                aTCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCggg                   >  7:311022/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)                                                                                                                                                  aa‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGa  <  6:98519/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                       CTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGA  >  NC_002947/5610504‑5610796

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