Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	185,367	A→C	6.9%	R82R (CGT→CGG)	PP_0161 ←	putative transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	185,367	0	A	C	6.9%	211.2 / 3.3	69	R82R (CGT→CGG)	PP_0161	putative transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base C (3/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGCCATCGGGCAGGGTGAAGACGCGTCGCTCGCCCGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCGTGGTGGCTGCTGCGCCATTGCGCGCAC  >  NC_002947/185220‑185514                                                                                                                                                     |                                                                                                                                                     cGAGCCATCGGGCAGGGTGAAGACGCGTCGCTCGCCCGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGAc                                                                                                                                                     >  4:35096/1‑149 (MQ=255)            gCAGGGTGAAGACGCGTCGCTCGCCCGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGgcgc                                                                                                                                         <  8:98782/149‑1 (MQ=255)                ggTGAAGACGCGTCGCTCGCCCGTGCCGGTGGCATTAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGCCGGCG‑CAGGCgctgc                                                                                                                                      >  1:242684/1‑149 (MQ=255)                ggTGAAGACGCGTCGCTCGCCCGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTG‑CGCCGTCGACAGGCgctgaa                                                                                                                                     <  7:188325/149‑3 (MQ=255)                  tGAAGACGCGTCGCTCGCCCGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTgctgg                                                                                                                                   >  5:83075/1‑149 (MQ=255)                     aGACGCGTCGCTCGCCCGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTgctggctg                                                                                                                                <  4:359357/149‑1 (MQ=255)                     aGACGCGTCGCTCGCCCGTGCCGGTGGCA‑TAGGCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTgctggctg                                                                                                                                <  3:23197/149‑1 (MQ=255)                      gACGCGTCGCTCGCCCGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGctggctgg                                                                                                                               >  3:353395/1‑149 (MQ=255)                            tcgctcgcCCGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCa                                                                                                                         <  4:353395/149‑1 (MQ=255)                            tcgctcgcCCGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCa                                                                                                                         <  5:82750/149‑1 (MQ=255)                            tcgctcacCCGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCa                                                                                                                         <  3:294058/149‑1 (MQ=255)                             cgctcgcCCGTGCCGGTTTCA‑TAGTCTGACGTCCAGGCAGTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCTTCCCTTCCATGACAACGCCACCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAg                                                                                                                        <  3:268315/149‑1 (MQ=255)                              gctcgcCCGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAg                                                                                                                        <  2:242684/148‑1 (MQ=255)                                    ccGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTg                                                                                                                  >  3:332972/1‑148 (MQ=255)                                         ccGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGc                                                                                                            <  4:109466/149‑1 (MQ=255)                                           ggTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCa                                                                                                          <  6:41130/149‑1 (MQ=255)                                            gTGGCA‑TAGTCTGACGTCCAGGCACTCAGAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGc                                                                                                            >  1:214040/1‑146 (MQ=255)                                            gTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGc                                                                                                            <  2:214040/146‑1 (MQ=255)                                                        tGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCgctgc                                                                                                                                      >  2:114519/1‑109 (MQ=255)                                                           cGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGa                                                                                           <  7:249832/149‑1 (MQ=255)                                                                    aCTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCa                                                                                  <  3:35096/149‑1 (MQ=255)                                                                    aCTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCa                                                                                  <  4:332972/149‑1 (MQ=255)                                                                    aCTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCa                                                                                  <  4:98155/149‑1 (MQ=255)                                                                      tCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAgg                                                                                >  7:170046/1‑149 (MQ=255)                                                                       cacaGTATCGAGGTCTTTGGCCCGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGt                                                                               >  2:242602/1‑149 (MQ=255)                                                                                gAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCt                                                                      >  6:187072/1‑149 (MQ=255)                                                                                  ggTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCa                                                                                  >  2:204413/1‑135 (MQ=255)                                                                                  ggTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCa                                                                                  <  1:204413/135‑1 (MQ=255)                                                                                   gTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTAt                                                                   >  6:159986/1‑149 (MQ=255)                                                                                      tttGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATg                                                                  >  4:70637/1‑147 (MQ=255)                                                                                      tttGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCA‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATg                                                                  >  8:234093/1‑147 (MQ=255)                                                                                          gccagccaGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGt                                                             >  1:187352/1‑148 (MQ=255)                                                                                          gccagccaGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGt                                                             <  4:129927/148‑1 (MQ=255)                                                                                            cagccagcTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGaa                                                          >  7:203646/1‑149 (MQ=255)                                                                                            cagccagcTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGGGAAGGTCCTTATGAAGGTgca                                                          >  8:311370/1‑147 (MQ=255)                                                                                                cagcTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATa                                                                                         >  4:242901/1‑114 (MQ=255)                                                                                                cagcTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATa                                                                                         <  3:242901/114‑1 (MQ=255)                                                                                                cagcTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGaa                                                                >  1:161683/1‑139 (MQ=255)                                                                                                cagcTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGaa                                                                <  2:161683/139‑1 (MQ=255)                                                                                                  gcTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAAc                                                         <  4:252347/144‑1 (MQ=255)                                                                                                  gcTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAAc                                                         >  3:252347/1‑144 (MQ=255)                                                                                                    tCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTa                                                                    <  1:80938/131‑1 (MQ=255)                                                                                                    tCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTa                                                                    >  2:80938/1‑131 (MQ=255)                                                                                                    tCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGc                                                  >  3:360304/1‑149 (MQ=255)                                                                                                         gACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGt                                                             <  6:272066/133‑1 (MQ=255)                                                                                                         gACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGt                                                             >  5:272066/1‑133 (MQ=255)                                                                                                            ccTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAgg                                                              >  2:35394/1‑129 (MQ=255)                                                                                                            ccTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGc                                          <  2:369691/149‑1 (MQ=255)                                                                                                            ccTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGCATAGCGCGCAGGTCAAGGTCCTTATGAAgg                                                              <  1:35394/129‑1 (MQ=255)                                                                                                             cTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGCACCATACGCTGCCAGGcc                                         >  4:41694/1‑149 (MQ=255)                                                                                                                       aaCGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGcc                                              >  6:212179/1‑134 (MQ=255)                                                                                                                       aaCGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGcc                                              <  5:212179/134‑1 (MQ=255)                                                                                                                       aaCGCCGCCAAGCAGCTTAAGCGCCTTTCGACGGCG‑CAGGCGCTGCTTGCTCGTATCAAGTGACTGCAATGCCAGCCCCGCACCCGGGATAGCGCGAAGGTGAAGGTCCTTATGAATGTGAACCATACGCTGCCAGGCCTGTTCGttt                                >  5:259313/1‑148 (MQ=255)                                                                                                                        aCGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTc                              <  3:70637/149‑1 (MQ=255)                                                                                                                          gccgccAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCgt                            <  6:83075/149‑1 (MQ=255)                                                                                                                             gccAAGCAGCTTCAGCGCCTGTCGACGGCA‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCgtggt                         <  7:234093/149‑1 (MQ=255)                                                                                                                              ccAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCgtggtg                        >  7:279912/1‑149 (MQ=255)                                                                                                                               cAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAg                                            <  4:173378/128‑1 (MQ=255)                                                                                                                                     ccagcAGCGCCTG‑CGCCGTCGACAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCGTGGTGGctgctg                 >  6:20027/5‑149 (MQ=255)                                                                                                                                      cTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCGTGGTGGCtgctgc                >  1:365367/1‑149 (MQ=255)                                                                                                                                       agcAGCGCCTG‑CGCCGTCGACAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCgt                            <  1:241788/134‑1 (MQ=255)                                                                                                                                       agcAGCGCCTG‑CGCCGTCGACAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGCTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCgt                            >  2:241788/3‑136 (MQ=255)                                                                                                                                        tCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCGTGGTGGCTGCTgcgc              >  3:163349/1‑149 (MQ=255)                                                                                                                                            cgcCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCTGGCCTGTTCGTTTTCGTGGTGGCTGCTGCGCCAtt          <  4:163349/149‑1 (MQ=255)                                                                                                                                             gcCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCGTGGTGGCTGCTGCGCCATTg         >  6:129098/1‑149 (MQ=255)                                                                                                                                                  tCGACGTCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCGTGGTGGCTGCTGCGCCATTgcgcgc    <  8:170046/149‑1 (MQ=255)                                                                                                                                                  tCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCGTGGTGGCTGCTGCGCCATTgcgcg     >  4:313555/1‑148 (MQ=255)                                                                                                                                                  tCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCGTGGTGGCTGCTGCGCCATTgcgc      <  7:311370/147‑1 (MQ=255)                                                                                                                                                    gACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGAGCGCAGGTCAATGTCCTTATTAATGTGAAACATACGCTGCCAGGCCTGTTCGTTTTCGTGGTGGCTGCTGCGCCATTGCGCGCAc  <  6:259313/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CGAGCCATCGGGCAGGGTGAAGACGCGTCGCTCGCCCGTGCCGGTGGCA‑TAGTCTGACGTCCAGGCACTCACAGTATCGAGGTCTTTGGCCAGCCAGCTCGCCGACCCTGCCATGACAACGCCGCCAAGCAGCTTCAGCGCCTGTCGACGGCG‑CAGGCGCTGCTGGCTGGTTTCCAGTGTCTGCAATGCCAGCCCGGCACCCGGGATAGCGCGCAGGTCAAGGTCCTTATGAAGGTGAACCATACGCTGCCAGGCCTGTTCGTTTTCGTGGTGGCTGCTGCGCCATTGCGCGCAC  >  NC_002947/185220‑185514

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