Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,744,607	G→T	6.0%	G73G (GGC→GGA)	PP_4198 ←	OmpA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,744,607	0	G	T	6.0%	328.3 / 8.5	100	G73G (GGC→GGA)	PP_4198	OmpA family protein
Reads supporting (aligned to +/- strand):  ref base G (45/49);  new base T (4/2);  total (49/51)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

CCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGA  >  NC_002947/4744462‑4744754                                                                                                                                                  |                                                                                                                                                      ccTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGcc                                                                                                                                                     >  1:54941/1‑148 (MQ=255)  cTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGa                                                                                                                                                   >  4:358523/1‑149 (MQ=255)  cTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGa                                                                                                                                                   >  5:406398/1‑149 (MQ=255)   tCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCTTGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGAt                                                                                                                                                  >  1:140242/1‑149 (MQ=255)      cGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGCACCGGTTCCGATGtg                                                                                                                                               <  5:363680/149‑2 (MQ=255)      cGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGtg                                                                                                                                               <  3:292129/149‑2 (MQ=255)         aaGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGCACCGGTTCCAATGtgcgg                                                                                                                                            <  6:161679/149‑5 (MQ=255)         aaGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACt                                                                                                                                            <  4:23742/149‑1 (MQ=255)         aaGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACt                                                                                                                                            >  7:403291/1‑149 (MQ=255)         aaGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACt                                                                                                                                            >  2:133362/1‑149 (MQ=255)         aaGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACt                                                                                                                                            <  8:113507/149‑1 (MQ=255)           ggCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTg                                                                                                                                           >  8:338867/1‑148 (MQ=255)            gCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGa                                                                                                                                         <  7:191684/149‑1 (MQ=255)              cAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGa                                                                                                                                         >  7:123472/1‑147 (MQ=255)                gCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGTTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAgg                                                                                                                                     <  3:211581/149‑1 (MQ=255)                    cAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGt                                                                                                                                                >  7:43804/1‑134 (MQ=255)                    cAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGt                                                                                                                                                <  8:43804/134‑1 (MQ=255)                      ggcgctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCg                                                                                                                               >  2:78391/1‑149 (MQ=255)                       gcgctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCgg                                                                                                                              >  7:344492/1‑149 (MQ=255)                        cgctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCCTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGgc                                                                                                                             >  8:85699/1‑149 (MQ=255)                         gctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTgg                                                                                                                                          <  5:156577/135‑1 (MQ=255)                         gctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTgg                                                                                                                                          >  6:156577/1‑135 (MQ=255)                         gctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGgcg                                                                                                                            >  7:320725/1‑149 (MQ=255)                         gctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGgcg                                                                                                                            <  8:359276/149‑1 (MQ=255)                         gctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGgcg                                                                                                                            <  3:391349/149‑1 (MQ=255)                             gcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGa                                                                                                                         <  6:8308/148‑1 (MQ=255)                             gcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGa                                                                                                                         <  1:426091/148‑1 (MQ=255)                             gcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGa                                                                                                                         <  4:99950/148‑1 (MQ=255)                             gcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGa                                                                                                                         <  5:328527/148‑1 (MQ=255)                               gctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGc                                                                                                                       >  2:308510/1‑148 (MQ=255)                                   cAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGc                                                                                                                   <  5:64649/148‑1 (MQ=255)                                    aGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGa                                                                                                                 <  8:123472/149‑1 (MQ=255)                                     ggTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGa                                                                                                                                         >  1:237541/1‑124 (MQ=255)                                     ggTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGa                                                                                                                                         <  2:237541/124‑1 (MQ=255)                                         gAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAc                                                                                                            >  6:396463/1‑149 (MQ=255)                                         gAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAc                                                                                                            <  4:279237/149‑1 (MQ=255)                                         gAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAc                                                                                                            >  8:160765/1‑149 (MQ=255)                                         gAGCTTGGCCTGGCGCGCCTGGTGCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAc                                                                                                            <  5:263580/149‑1 (MQ=255)                                           gCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAccc                                                                                                          >  4:102676/1‑149 (MQ=255)                                           gCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAccc                                                                                                          >  4:114290/1‑149 (MQ=255)                                              tGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAcccgcc                                                                                                       >  8:273515/1‑149 (MQ=255)                                               ggCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAcccgccc                                                                                                      >  4:254981/1‑149 (MQ=255)                                                gCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAcccgccc                                                                                                      <  7:160765/148‑1 (MQ=255)                                                gCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAcccgccc                                                                                                      <  3:367604/148‑1 (MQ=255)                                                 ccTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGa                                                                                                    <  8:320725/149‑1 (MQ=255)                                                   tGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTAGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCCGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCACGAAt                                                                                                  >  4:249571/1‑149 (MQ=255)                                                     gcgcTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAAt                                                                                                  <  5:145868/147‑1 (MQ=255)                                                           cTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAAtca‑cgtcac                                                                                         <  1:133362/149‑1 (MQ=255)                                                               ttCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCAGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACgcgg                                                                                     >  7:153344/1‑149 (MQ=255)                                                                   gcgcCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGgcgcg                                                                                 >  3:359107/1‑149 (MQ=255)                                                                     gcCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCAGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑Ct                                                                              <  8:153344/149‑1 (MQ=255)                                                                      ccAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTg                                                                             >  4:130071/1‑149 (MQ=255)                                                                      ccAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTg                                                                             <  8:403291/149‑1 (MQ=255)                                                                          cTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCgca                                                                         >  2:296386/1‑149 (MQ=255)                                                                          cTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCgca                                                                         <  3:358523/149‑1 (MQ=255)                                                                             ggCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCgcag                                                                        <  4:184889/147‑1 (MQ=255)                                                                              gCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCgcagca                                                                      <  3:102676/148‑1 (MQ=255)                                                                                   cTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTACGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCg                                                                 <  7:338867/148‑1 (MQ=255)                                                                                   cTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCg                                                                 <  6:326818/148‑1 (MQ=255)                                                                                   cTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCg                                                                 <  8:344492/148‑1 (MQ=255)                                                                                      gcgGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGa                                                               <  1:346370/147‑1 (MQ=255)                                                                                      gcgGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAat                                                             >  4:386705/1‑149 (MQ=255)                                                                                       cgGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAatc                                                            <  3:130071/149‑1 (MQ=255)                                                                                       cgGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAat                                                             >  5:338886/1‑148 (MQ=255)                                                                                          gCAATCTAGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAatcatc                                                         >  4:198854/1‑149 (MQ=255)                                                                                            aaTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACCGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCt                                                        <  1:78391/149‑1 (MQ=255)                                                                                            aaTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCtt                                                       >  5:191120/1‑149 (MQ=255)                                                                                            aaTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCtt                                                       >  4:417257/1‑149 (MQ=255)                                                                                             aTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCttt                                                      <  4:359107/149‑1 (MQ=255)                                                                                                  ggtGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTTGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCGCCTGCGCAGCACATCGGAATCATCTTTGa                                                    <  1:19681/145‑1 (MQ=255)                                                                                                  gCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGa                                                    <  5:315530/146‑1 (MQ=255)                                                                                                    tGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGAccccc                                               >  6:102054/1‑149 (MQ=255)                                                                                                           ctggctgATGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAgg                                        <  3:81432/149‑1 (MQ=255)                                                                                                              gctgAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCaa                                     <  3:261447/149‑1 (MQ=255)                                                                                                                tgAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACg                                   >  1:395303/1‑149 (MQ=255)                                                                                                                    gTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGaa                                          >  2:122470/1‑138 (MQ=255)                                                                                                                    gTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGaa                                          <  1:122470/138‑1 (MQ=255)                                                                                                                     tAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGc                              <  6:191120/149‑1 (MQ=255)                                                                                                                       aGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGcc                             <  5:212412/148‑1 (MQ=255)                                                                                                                               gCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAgc                    >  1:415947/1‑149 (MQ=255)                                                                                                                               cagttacaTCGGCACCGGTTCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAgc                    >  3:368697/6‑149 (MQ=255)                                                                                                                                ccGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAgcg                   >  3:64708/1‑149 (MQ=255)                                                                                                                                agttacaTCGGCACCGGTTCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAgcg                   >  2:112684/5‑149 (MQ=255)                                                                                                                                agttacaTCGGCACCGGTTCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCAACTTTGACCCCCTCGAAGGCAACGGTGGCCTCATCCAgcg                   >  7:308142/5‑149 (MQ=255)                                                                                                                                agttacaTCGGCACCGGTTCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCACAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAgcg                   >  5:326524/5‑149 (MQ=255)                                                                                                                                 cGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCAAgcgc                  >  2:274609/1‑149 (MQ=255)                                                                                                                                   cacaTCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCgg                <  2:140242/149‑1 (MQ=255)                                                                                                                                      aTCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAgcg                   <  2:41635/143‑1 (MQ=255)                                                                                                                                      aTCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAgcg                   >  1:41635/1‑143 (MQ=255)                                                                                                                                      aTCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCAc             >  6:424853/1‑149 (MQ=255)                                                                                                                                        cGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCActg           >  3:131293/1‑149 (MQ=255)                                                                                                                                         ggAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCActg           <  5:102054/148‑1 (MQ=255)                                                                                                                                          gAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgct         >  1:408640/1‑149 (MQ=255)                                                                                                                                           aaCCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgctg        <  3:417257/149‑1 (MQ=255)                                                                                                                                           aaCCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgctg        <  4:179313/149‑1 (MQ=255)                                                                                                                                             ccGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTgcg      <  3:198854/149‑1 (MQ=255)                                                                                                                                               ggTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTACTCCAGCGCGGCACTGCTGCGCt    <  7:273515/149‑1 (MQ=255)                                                                                                                                                gTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCt    >  6:445319/1‑148 (MQ=255)                                                                                                                                                gTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCt    >  7:152009/1‑148 (MQ=255)                                                                                                                                                 tGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGa  >  3:419963/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGA  >  NC_002947/4744462‑4744754

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