Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	391,590	T→G	6.3%	W622G (TGG→GGG)	soxA →	sarcosine oxidase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	391,590	0	T	G	6.3%	247.9 / 6.8	79	W622G (TGG→GGG)	soxA	sarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (34/40);  new base G (3/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

CATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGCCTCGAC  >  NC_002947/391442‑391738                                                                                                                                                        |                                                                                                                                                     cATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCt                                                                                                                                                      >  5:80924/1‑149 (MQ=255)     gccgccCGTTCGATAGGCCTCACCATCCCCGTAATGGGTTACACCATGTTTCGCCCCAAATACACGCCGGTTACTTTTGGGCTGCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTgg                                                                                                                                                    <  1:102911/149‑1 (MQ=255)             ttCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTCCACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCACGTGCCAGGCATGCAGGGCA‑GTgaagc                                                                                                                                          <  8:90182/149‑2 (MQ=255)                  tCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGGTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGgc                                                                                                                                     >  8:304342/1‑149 (MQ=255)                    ggCATCACCATCCCGGAAAGGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGGCACCTGTTCGAGCCCGTGCGCTTCACGTGCCAGGCATGCAGGGCA‑GTgaagcgcacgg                                                                                                                                    <  1:155988/148‑8 (MQ=37)                      cATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGcc                                                                                                                                  >  5:248130/1‑148 (MQ=255)                       aTCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGa                                                                                                                                <  4:169103/149‑1 (MQ=255)                          aCCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGtt                                                                                                                             >  3:129081/1‑149 (MQ=255)                          aCCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGtt                                                                                                                             >  3:124082/1‑149 (MQ=255)                            cATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACg                                                                                                                                       >  4:289859/1‑137 (MQ=255)                            cATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACg                                                                                                                                       <  3:289859/137‑1 (MQ=255)                              tCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGaa                                                                                                                         <  3:48352/149‑1 (MQ=255)                                ccGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGa                                                                                                                       <  7:267957/149‑1 (MQ=255)                                 cGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGAc                                                                                                                      >  7:169283/1‑149 (MQ=255)                                         ggTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCa                                                                                                              <  6:133833/149‑1 (MQ=255)                                           tACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGt                                                                                                            >  3:72823/1‑149 (MQ=255)                                              caccaTGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTgg                                                                                                          >  3:335352/1‑148 (MQ=255)                                              caccaTGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTgg                                                                                                          >  3:1228/1‑148 (MQ=255)                                              caccaTGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTgg                                                                                                          >  6:182768/1‑148 (MQ=255)                                               accaTGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGGTGTCACCTGTTCGAGCCCTTGCGCTTCAC‑TGCCCTGCATGCCTGTCACGTGAAGAACGGCGCCTAGTTCGAAGACGTCGGCCAGTGGaa                                                                                                        <  6:342164/149‑1 (MQ=255)                                                 caTGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAgc                                                                                                      >  7:154099/1‑149 (MQ=255)                                                     ttCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGcc                                                                                                   <  6:248130/148‑1 (MQ=255)                                                         gCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                                              >  7:285697/1‑149 (MQ=255)                                                         gCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                                              <  1:144243/149‑1 (MQ=255)                                                         gCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                                              <  4:129081/149‑1 (MQ=255)                                                         gCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTgg                                                                                              <  4:124082/149‑1 (MQ=255)                                                           cccAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTa                                                                                            <  4:335352/149‑1 (MQ=255)                                                              aaCTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACtt                                                                                         >  8:11477/1‑149 (MQ=255)                                                                     cGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGaaa                                                                                  <  1:104474/149‑1 (MQ=255)                                                                         ggTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCg                                                                              >  6:328776/1‑149 (MQ=255)                                                                          gTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCgg                                                                             <  7:11477/149‑1 (MQ=255)                                                                          gTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCgg                                                                             <  5:328776/149‑1 (MQ=255)                                                                          gTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCgg                                                                             >  5:180515/1‑149 (MQ=255)                                                                             aC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTg                                                                           <  4:1228/148‑1 (MQ=255)                                                                                    ggC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACAt                                                                    <  6:106960/149‑1 (MQ=255)                                                                                         cgGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCAt                                                                <  1:121439/149‑1 (MQ=255)                                                                                         cgGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTACCGAAAGCCGGTGAAGACATCCAt                                                                >  3:128433/1‑149 (MQ=255)                                                                                          ggTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCAtg                                                               <  6:115138/149‑1 (MQ=255)                                                                                               cGGGGCGTCACTGTTGTCACCTGTTCGATCCCGTGCGCTTCAC‑TTCCCTGCATGCCTGGCACGTTAATAACGGCGCCGATTTCGAATACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGcc                                                          <  4:128433/149‑1 (MQ=255)                                                                                                      tCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTc                                                   >  5:272927/1‑149 (MQ=255)                                                                                                      tCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTc                                                   >  3:293793/1‑149 (MQ=255)                                                                                                        aCTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTcgc                                                 <  8:169283/149‑1 (MQ=255)                                                                                                        aCTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTcgc                                                 <  4:293793/149‑1 (MQ=255)                                                                                                            tgGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGaa                                              <  8:154099/148‑1 (MQ=255)                                                                                                            tgGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAAt                                             <  6:180515/149‑1 (MQ=255)                                                                                                            tgGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAAt                                             <  6:80924/149‑1 (MQ=255)                                                                                                              gTCACCTGTTCGAGCCCGTGCGCTTCTC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATg                                            >  1:183510/1‑148 (MQ=255)                                                                                                              gTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATg                                            >  6:53254/1‑148 (MQ=255)                                                                                                              gTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATg                                            <  2:308193/148‑1 (MQ=255)                                                                                                                 aCCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAg                                        >  2:218157/1‑149 (MQ=255)                                                                                                                  ccTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAgg                                       <  6:272927/149‑1 (MQ=255)                                                                                                                    tGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGcc                                     >  3:29586/1‑149 (MQ=255)                                                                                                                        cGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAAt                                             >  5:112885/1‑137 (MQ=255)                                                                                                                        cGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAAt                                             <  6:112885/137‑1 (MQ=255)                                                                                                                         gAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTgcg                                >  7:62235/1‑149 (MQ=255)                                                                                                                           gCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCAAGAGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACCTCCATGCGGCCGTGCCTCGCGAAAGCAAGGCGGTGcgcg                              >  4:276121/1‑149 (MQ=255)                                                                                                                               gTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGATAGCCGGTGATGACATCCATGCTGACGTGACTcgcg                                                >  6:238126/1‑127 (MQ=255)                                                                                                                               gTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTcgcg                                                <  5:238126/127‑1 (MQ=255)                                                                                                                               gTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAg                          >  6:323265/1‑149 (MQ=255)                                                                                                                                cgttcTTCACGTGCCAGGCATGCAGGGCA‑GTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGc                         >  5:12879/5‑149 (MQ=38)                                                                                                                                  ttcTTCACGTGCCAGGCATGCAGGGCA‑GTGAAGAACGGCTCCGAGTTCGAAGAAGTAGGCAAGTGGAATCGCCCTTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGATGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGt                       >  2:27966/3‑149 (MQ=255)                                                                                                                                  ttcTTCACGTGCCAGGCATGCAGGGCA‑GTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGt                       >  6:252412/3‑149 (MQ=255)                                                                                                                                  cgcTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCg                        <  5:182768/148‑1 (MQ=255)                                                                                                                                   gcTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTg                      >  4:337326/1‑149 (MQ=255)                                                                                                                                   gcTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATACAAGGCCGTGCGCGACAGCGTg                      >  5:258652/1‑149 (MQ=255)                                                                                                                                     ttCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGaa                                              <  4:230158/123‑1 (MQ=255)                                                                                                                                     ttCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGaa                                              >  3:230158/1‑123 (MQ=255)                                                                                                                                     ttCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTggg                    >  4:143466/1‑149 (MQ=255)                                                                                                                                       cAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGc                   <  8:31368/148‑1 (MQ=255)                                                                                                                                       cAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATACAAGGCCGTGCGCGACAGCGTGGGc                   <  6:258652/148‑1 (MQ=255)                                                                                                                                        aC‑TGCCCTGCATGCCTGGCACGTGAAGTACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCTACAGCGTGGGCct                 <  3:276121/149‑1 (MQ=255)                                                                                                                                         c‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTc                                                   <  2:14312/114‑1 (MQ=255)                                                                                                                                         c‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTc                                                   >  1:14312/1‑114 (MQ=255)                                                                                                                                             cccTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCctgctg             <  2:183510/149‑1 (MQ=255)                                                                                                                                             cccTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCctgct              <  1:218157/148‑1 (MQ=255)                                                                                                                                                tGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGAc          >  6:50804/1‑149 (MQ=255)                                                                                                                                                  cATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGAc          <  4:29586/147‑1 (MQ=255)                                                                                                                                                  cATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGc        <  4:72823/149‑1 (MQ=255)                                                                                                                                                   aTGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGcc       <  3:13012/149‑1 (MQ=255)                                                                                                                                                     gCCTGGCACATGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGCCTc     >  4:326690/1‑149 (MQ=255)                                                                                                                                                        tGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGCCTCGAc  >  3:10844/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                     CATCGCCGCCCGTTCGATCGGCATCACCATCCCGGAAATGGGTACCACCATGTTCCGCCCCAACTACACGCCGGTTAC‑TTTCGGC‑GCGGTAGCGGGCCGTCACTGTGGTCACCTGTTCGAGCCCGTGCGCTTCAC‑TGCCCTGCATGCCTGGCACGTGAAGAACGGCGCCGAGTTCGAAGACGTCGGCCAGTGGAAGCGCCCGTGGTACTTCCCGAAAGCCGGTGAAGACATCCATGCTGCCGTGACTCGCGAATGCAAGGCCGTGCGCGACAGCGTGGGCCTGCTGGACGCCTCGAC  >  NC_002947/391442‑391738

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