Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,856,513	A→T	6.1%	E147V (GAA→GTA) ‡	PP_5133 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,856,513	0	A	T	6.1%	197.8 / 2.8	65	E147V (GAA→GTA) ‡	PP_5133	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base T (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

GGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCATGCTA  >  NC_002947/5856371‑5856658                                                                                                                                               |                                                                                                                                                  gggCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCtt                                                                                                                                             <  3:259013/149‑1 (MQ=255)    cAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCAAAGCGTTCGCTTAGc                                                                                                                                          <  2:254686/149‑1 (MQ=255)         cccGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCACTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCcgc                                                                                                                                      >  2:238607/1‑148 (MQ=255)          ccGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGa                                                                                                                                    >  4:350359/1‑149 (MQ=255)              tGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGcc                                                                                                                                <  1:382285/149‑1 (MQ=255)              tGAAGGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGACTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGcc                                                                                                                                <  4:325853/149‑1 (MQ=255)               gAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGcc                                                                                                                                >  7:52064/1‑148 (MQ=255)                 aTGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAAc                                                                                                                                                 >  5:216539/1‑129 (MQ=255)                 aTGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAAc                                                                                                                                                 <  6:216539/129‑1 (MQ=255)                  tgtgTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCAGGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAgg                                                                                                                            <  2:255522/149‑1 (MQ=255)                     gtTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGcc                                                                                                                                         <  4:70355/133‑1 (MQ=255)                     gtTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGcc                                                                                                                                         >  3:70355/1‑133 (MQ=255)                        cTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCcgcg                                                                                                                       >  7:53199/1‑148 (MQ=255)                         tGGGCGCTGGGCAGCCTGCTCGGGGTGAACCTGCTGATGCTTCTTGGCTTGGGCTTGGGTTTTGCCCTGAGCGGCGAGCATGGCGGAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCa                                                                                                                     <  4:69987/149‑1 (MQ=255)                           ggCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCAcc                                                                                                                   >  6:342797/1‑149 (MQ=255)                            gcgcTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCt                                                                                                                  >  1:326500/1‑149 (MQ=255)                                    cAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAg                                                                                                          <  4:5293/149‑1 (MQ=255)                                          gCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCgctgc                                                                                                    <  8:52064/149‑1 (MQ=255)                                              gggCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGt                                                                                                >  7:95588/1‑149 (MQ=255)                                                       ctgctgATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGtgctgctgc                                                                                        >  4:211697/1‑148 (MQ=255)                                                         gctgATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAAc                                                                                     >  3:381002/1‑149 (MQ=255)                                                          ctgATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGt                                                                                                >  2:349795/1‑137 (MQ=255)                                                            gATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAAGGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGgcgc                                                                                                             >  2:93860/1‑122 (MQ=255)                                                            gATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAAGGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGgcgc                                                                                                             <  1:93860/122‑1 (MQ=255)                                                               gCTACATGGCTGGGGCATGGGTATTGCCCTGACCGGCGAGCATGGCGCAGCGCTGGGCCTCTAGTGGCTGTGGCTAAGCGAACGCTATGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTCCTGGTGCTGCTGCAACGTCAAc                                                                               >  6:114209/1‑149 (MQ=255)                                                                  aCTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGc                                                                             >  8:115769/1‑148 (MQ=255)                                                                   cTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGcc                                                                            <  8:95588/148‑1 (MQ=255)                                                                       gCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAAc                                                                       <  3:350359/149‑1 (MQ=255)                                                                               aTGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGgct                                                               >  4:42387/1‑149 (MQ=255)                                                                                  ggTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCCACCGGTGgctgct                                                            >  5:322269/1‑149 (MQ=255)                                                                                       tGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAATGTCAACGCCTCAACCGGTGGCTGCTCGGc                                                        <  1:237043/148‑1 (MQ=255)                                                                                          ccTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCa                                                                                                           >  1:333750/1‑94 (MQ=255)                                                                                          ccTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCa                                                                                                           <  2:333750/94‑1 (MQ=255)                                                                                          ccTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGt                                                    >  1:16004/1‑149 (MQ=255)                                                                                           cTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGtt                                                   >  5:33027/1‑149 (MQ=255)                                                                                               gcggcgAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTa                                               <  6:33027/149‑1 (MQ=255)                                                                                                cggggAGCATGGGGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCTAGGCGGCGCACCTGGCGCCAGCGCTGGTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTac                                              <  5:114209/149‑1 (MQ=255)                                                                                                cggcgAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTac                                              <  7:115769/149‑1 (MQ=255)                                                                                                 ggcgAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTaca                                             >  6:200483/1‑149 (MQ=255)                                                                                                   cgAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACg                                           <  4:381002/149‑1 (MQ=255)                                                                                                         tGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCtgt                                     <  8:262271/149‑1 (MQ=255)                                                                                                           gcgcAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCtgtg                                    <  2:130723/148‑1 (MQ=255)                                                                                                            cgcAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACgg                                          <  7:286037/141‑1 (MQ=255)                                                                                                            cgcAGGGCTGGGCCGCTTGTGGCAGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACgg                                          >  8:286037/1‑141 (MQ=255)                                                                                                                 ggCTGGGCCGCTTGTAGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGctgctg                             >  7:241391/1‑149 (MQ=255)                                                                                                                  gCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTgg                            <  2:326500/149‑1 (MQ=255)                                                                                                                    tGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGctgctg                             <  5:311478/146‑1 (MQ=255)                                                                                                                    tGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGctgctg                             >  6:311478/1‑146 (MQ=255)                                                                                                                    tGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGc                           >  5:23004/1‑148 (MQ=255)                                                                                                                       gCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATg                       <  3:34384/149‑1 (MQ=255)                                                                                                                       gCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATg                       <  3:42387/149‑1 (MQ=255)                                                                                                                       gCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGATGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGATGGTACACGGCCTGTGGCTGCTGGCGAt                        >  3:396895/1‑148 (MQ=255)                                                                                                                        ccGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGc                      >  2:277795/1‑149 (MQ=255)                                                                                                                            ttGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCAc                  >  5:113851/1‑149 (MQ=255)                                                                                                                                 tcgcgggcaAAGCGTTCGCTTAGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCgcgc             >  6:13293/6‑149 (MQ=255)                                                                                                                                 gctgtggctAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCgcgc             >  3:102156/1‑149 (MQ=255)                                                                                                                                  ctgtggctAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAGCGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCt            >  3:43190/1‑149 (MQ=255)                                                                                                                                    cgggcaAAGCGTTCGCTTAGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACACCTCAACCGGTGGCTGCTCGGCCTGTTGGTa                                               >  6:276773/3‑112 (MQ=255)                                                                                                                                    cgggcaAAGCGTTCGCTTAGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACACCTCAACCGGTGGCTGCTCGGCCTGTTGGTa                                               <  5:276773/110‑1 (MQ=255)                                                                                                                                       gctAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGc        <  6:113851/148‑1 (MQ=255)                                                                                                                                         tAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCAt      >  5:299273/1‑148 (MQ=255)                                                                                                                                           aGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCAt      <  6:299273/146‑1 (MQ=255)                                                                                                                                           aGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCATGc    >  6:233722/1‑148 (MQ=255)                                                                                                                                            gCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAGCGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCATGCTa  <  4:43190/149‑1 (MQ=255)                                                                                                                                            gCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACgg                                          >  8:324388/1‑109 (MQ=255)                                                                                                                                               |                                                                                                                                                  GGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCATGCTA  >  NC_002947/5856371‑5856658

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