Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,182,854	G→T	6.4%	G115V (GGC→GTC)	PP_4562 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,182,854	0	G	T	6.4%	252.3 / 5.1	78	G115V (GGC→GTC)	PP_4562	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/35);  new base T (2/3);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

GCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTG  >  NC_002947/5182708‑5182999                                                                                                                                                   |                                                                                                                                                     gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGGACGGGTGTCg                                                                                                                                                    <  3:124517/149‑1 (MQ=255)    gcagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGag                                                                                                                                                  >  1:278031/1‑148 (MQ=255)     cagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCTCGCCACCCGGACGGGTGTCgcg                                                                                                                                                  <  7:354232/147‑3 (MQ=255)     cagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGagca                                                                                                                                                <  4:212002/149‑1 (MQ=255)     cagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGagc                                                                                                                                                 <  4:363836/148‑1 (MQ=255)     cagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGag                                                                                                                                                  >  2:154715/1‑147 (MQ=255)       gcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGagcagc                                                                                                                                              >  7:261789/1‑149 (MQ=255)         aCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTg                                                                                                                                            <  1:251609/149‑1 (MQ=255)           cTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGa                                                                                                                                          >  6:282094/1‑149 (MQ=255)              aaGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTgg                                                                                                                                           >  4:212564/1‑145 (MQ=255)              aaGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTgg                                                                                                                                           <  3:212564/145‑1 (MQ=255)              aaGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGAc                                                                                                                                         <  5:356826/147‑1 (MQ=255)              aaGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACtt                                                                                                                                       <  1:154715/149‑1 (MQ=255)                ggcgCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGGCCGGGTGGCGAGCAGCTGGACTTcc                                                                                                                                     <  5:7439/146‑1 (MQ=255)                ggTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTcc                                                                                                                                     <  3:27584/149‑1 (MQ=255)                   gCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGc                                                                                                                                   >  1:220251/1‑148 (MQ=255)                     ccaccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAAc                                                                                                                                >  6:77063/1‑149 (MQ=255)                       accaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACt                                                                                                                               >  4:102024/1‑148 (MQ=255)                              tCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACa                                                                                                                        <  3:301324/148‑1 (MQ=255)                                 tCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCg                                                                                                                    >  7:272428/1‑149 (MQ=255)                                 tCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCg                                                                                                                    >  4:49179/1‑149 (MQ=255)                                      gCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAgg                                                                                                               <  1:275347/149‑1 (MQ=255)                                           ccTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAgg                                                                                                                      >  7:279612/1‑137 (MQ=255)                                           ccTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAgg                                                                                                                      <  8:279612/137‑1 (MQ=255)                                           ccTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGt                                                                                                           >  5:388101/1‑148 (MQ=255)                                            cTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCa                                                                                                         <  8:261789/149‑1 (MQ=255)                                                 gCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑Ctt                                                                                                 <  5:282094/149‑1 (MQ=255)                                                   aaaTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTtc                                                                                               <  3:49179/149‑1 (MQ=255)                                                    aaaggAAGTGCGCAAGCCTGGCAAGCGAACCTCCACCATCAGCCTGGCCCTCAGCGCCGCCCGGGACTCGCTGCTGCGCGACCCACGTACGGGTGGCGAGCCGCTGGAATTCCGCAACTTCGTACAGGCGTAAGCCCGTGAGG‑‑‑CTTact                                                                                              >  3:204705/4‑149 (MQ=255)                                                         aaGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATg                                                                                         <  3:102024/149‑1 (MQ=255)                                                          aGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGc                                                                                        >  6:161198/1‑149 (MQ=255)                                                            tGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGa                                                                                      >  5:385145/1‑149 (MQ=255)                                                                  aGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCa                                                                                <  3:120809/149‑1 (MQ=255)                                                                   gCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCAt                                                                               <  7:69094/149‑1 (MQ=255)                                                                     cTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATgg                                                                             <  6:388101/149‑1 (MQ=255)                                                                     cTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATg                                                                              >  8:350261/1‑148 (MQ=255)                                                                      tGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAg                                                                                                        >  6:390696/1‑124 (MQ=255)                                                                      tGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAg                                                                                                        <  5:390696/124‑1 (MQ=255)                                                                      tGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGa                                                                            >  2:399364/1‑149 (MQ=255)                                                                      tGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCAGTCAGT‑‑‑CTTTCTCGATGCGATCGTCATGGa                                                                            >  4:48896/1‑149 (MQ=255)                                                                        gCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGAt                                                                           <  6:411009/148‑1 (MQ=255)                                                                              gccgcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATc                                                                    <  2:220251/149‑1 (MQ=255)                                                                               ccgcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCa                                                                   <  7:344121/149‑1 (MQ=255)                                                                                    ccaccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGTGATCGTCATCCATGACGATcgcatcg                                                              <  5:77063/149‑7 (MQ=255)                                                                                      aacaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCAGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTg                                                            >  7:264172/3‑149 (MQ=255)                                                                                        caTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCt                                                          >  5:244245/1‑149 (MQ=255)                                                                                        caTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCt                                                          >  1:247655/1‑149 (MQ=255)                                                                                         aTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCtt                                                         <  5:296231/149‑1 (MQ=255)                                                                                          tCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCtt                                                         >  5:819/1‑148 (MQ=255)                                                                                             gCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCa                                                                                                         >  6:201946/1‑100 (MQ=255)                                                                                             gCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCa                                                                                                         <  5:201946/100‑1 (MQ=255)                                                                                             gCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaa                                                      >  6:201982/1‑148 (MQ=255)                                                                                              ccTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaac                                                     >  4:82308/1‑148 (MQ=255)                                                                                              ccTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaac                                                     >  5:334581/1‑148 (MQ=255)                                                                                                tGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGagc                                                                                                                                                 >  8:214983/1‑57 (MQ=255)                                                                                                tGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGagc                                                                                                                                                 <  7:214983/57‑1 (MQ=255)                                                                                                   cccTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCg                                               <  2:278031/149‑1 (MQ=255)                                                                                                                  gggACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTg                                <  8:272428/149‑1 (MQ=255)                                                                                                                  gggACTCGCTGCTGCGCGACACCAGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCTATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTg                                <  7:350261/149‑1 (MQ=255)                                                                                                                         gctgctgctCGACACCCGTCCGGGTGGCGAGCAGCTTGACTTCAGCAACTTCGTACAGGCCTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCGTTCTTGAATTTGc                         <  4:204705/149‑1 (MQ=255)                                                                                                                         gctgctGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                         >  7:167368/1‑149 (MQ=255)                                                                                                                          ctgctgCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGcc                        >  3:84452/1‑149 (MQ=255)                                                                                                                             ctgcTCGCCACCCGGACGGGTGTCGCGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCtt                     >  5:416803/1‑149 (MQ=255)                                                                                                                             ctgcGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCtt                     >  5:281606/1‑149 (MQ=255)                                                                                                                              tgcGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGGGAGTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGccc                       >  1:125721/1‑149 (MQ=255)                                                                                                                              tgcGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATTACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTg                    >  7:266350/1‑149 (MQ=255)                                                                                                                                  cgACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGa                <  1:399364/149‑1 (MQ=255)                                                                                                                                    acacCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACg              <  5:201982/149‑1 (MQ=255)                                                                                                                                       cccGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGtt            >  3:25611/1‑148 (MQ=255)                                                                                                                                        ccggacGGGTGTCGCGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGa                <  6:293476/138‑1 (MQ=255)                                                                                                                                        ccggacGGGTGTCGCGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGa                >  5:293476/6‑143 (MQ=255)                                                                                                                                         cGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGaa         <  2:247655/149‑1 (MQ=255)                                                                                                                                          gTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGaa         >  3:39567/1‑148 (MQ=255)                                                                                                                                           tCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACa       <  8:167368/149‑1 (MQ=255)                                                                                                                                             cGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGc     >  4:120335/1‑149 (MQ=255)                                                                                                                                                gTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTg  >  1:298174/1‑149 (MQ=255)                                                                                                                                                 tGGCGAGCAGCTGGACTTCCACAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaac                                                     >  7:173555/1‑97 (MQ=255)                                                                                                                                                 tGGCGAGCAGCTGGACTTCCACAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaac                                                     <  8:173555/97‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGG‑‑‑CTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTG  >  NC_002947/5182708‑5182999

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