Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,744,600	A→C	6.3%	S76A (TCC→GCC)	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,600	0	A	C	6.3%	293.5 / 9.6	95	S76A (TCC→GCC)	PP_4198	OmpA family protein
Reads supporting (aligned to +/- strand):  ref base A (42/47);  new base C (4/2);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

TTGGCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTG  >  NC_002947/4744458‑4744748                                                                                                                                               |                                                                                                                                                       ttGGCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCgg                                                                                                                                                   <  7:168976/148‑1 (MQ=255)  tGGCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGt                                                                                                                                                  <  1:3592/148‑1 (MQ=255)     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                                                                                                                                      <  8:113507/149‑1 (MQ=255)             aaGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACt                                                                                                                                      >  7:403291/1‑149 (MQ=255)             aaGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACt                                                                                                                                      <  4:23742/149‑1 (MQ=255)             aaGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACt                                                                                                                                      >  2:133362/1‑149 (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                                                                                                                                          <  8:43804/134‑1 (MQ=255)                        cAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGt                                                                                                                                          >  7:43804/1‑134 (MQ=255)                          ggcgctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCg                                                                                                                         >  2:78391/1‑149 (MQ=255)                           gcgctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCgg                                                                                                                        >  7:344492/1‑149 (MQ=255)                            cgctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCCTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGgc                                                                                                                       >  8:85699/1‑149 (MQ=255)                             gctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTgg                                                                                                                                    >  6:156577/1‑135 (MQ=255)                             gctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTgg                                                                                                                                    <  5:156577/135‑1 (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                                                                                                                   <  5:328527/148‑1 (MQ=255)                                 gcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGa                                                                                                                   <  4:99950/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                                                                                                      >  8:160765/1‑149 (MQ=255)                                             gAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAc                                                                                                      <  4:279237/149‑1 (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                                                                                                <  3:367604/148‑1 (MQ=255)                                                    gCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAcccgccc                                                                                                <  7:160765/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                                                                   <  3:358523/149‑1 (MQ=255)                                                                              cTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCgca                                                                   >  2:296386/1‑149 (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                                    <  1:122470/138‑1 (MQ=255)                                                                                                                        gTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGaa                                    >  2:122470/1‑138 (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)                                                                                                                                               |                                                                                                                                                       TTGGCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCA‑CGTCACGCGGCGCG‑CTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTG  >  NC_002947/4744458‑4744748

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