Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,386,775	G→A	5.7%	S251L (TCG→TTG)	PP_3863 ←	tail protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,386,775	0	G	A	5.7%	213.7 / 2.0	70	S251L (TCG→TTG)	PP_3863	tail protein
Reads supporting (aligned to +/- strand):  ref base G (33/33);  new base A (2/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

CCGCACGATCATGTTGTGCCGCCATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATC  >  NC_002947/4386630‑4386920                                                                                                                                                  |                                                                                                                                                  ccGCACGATCATGTTGTGCCGCCATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACt                                                                                                                                                >  7:104622/1‑149 (MQ=255) ccGCACGATCATGTTGTGCCGCCATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACt                                                                                                                                                <  7:49488/149‑1 (MQ=255)  cGCACGATCATGTTGTGCCGCCATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCCCCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTg                                                                                                                                               >  8:42654/1‑149 (MQ=255)   gCACGATCATGTTGTGCCGCCATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTgg                                                                                                                                              <  8:30620/149‑1 (MQ=255)       gATCATGTTGTGCCGCCATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGa                                                                                                                                          >  6:240295/1‑149 (MQ=255)          cATGTTGTGCCGCCATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATc                                                                                                                                        <  4:210425/148‑1 (MQ=255)           aTGTTGTGCCGCCATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCAc                                                                                                                                      <  4:112715/149‑1 (MQ=255)              ttgtGCCGCCATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACC‑‑CGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGa                                                                                                                                 <  6:401503/149‑1 (MQ=255)              ttgtGCCGCCATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGGCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCt                                                                                                                                    <  5:260249/148‑1 (MQ=255)                gtgCCGCCATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGCGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTg                                                                                                                                  <  2:283932/148‑1 (MQ=255)                    cgccATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACccg                                                                                                                             >  6:191193/1‑149 (MQ=255)                                ccATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCAc                                                                                                                 <  7:151771/149‑1 (MQ=255)                                       tctGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCAcc                                                                                                                                     <  1:356352/122‑1 (MQ=255)                                       tctGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCAcc                                                                                                                                     >  2:356352/1‑122 (MQ=255)                                       tctGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTc                                                                                                           >  6:5746/1‑148 (MQ=255)                                            cgccACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGc                                                                                                     <  6:315030/149‑1 (MQ=255)                                                    cTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTgg                                                                                             >  8:380613/1‑149 (MQ=255)                                                    cTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCAACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTgg                                                                                             >  5:30877/1‑149 (MQ=255)                                                    cTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCAACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTgg                                                                                             >  3:47829/1‑149 (MQ=255)                                                       cACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGcc                                                                                           >  1:352548/1‑148 (MQ=255)                                                       cACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGcc                                                                                           <  2:352548/148‑1 (MQ=255)                                                            tgtAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGt                                                                                                                         >  1:51410/1‑113 (MQ=255)                                                            tgtAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGt                                                                                                                         <  2:51410/113‑1 (MQ=255)                                                            tgtAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGAcc                                                                                     >  6:323101/1‑149 (MQ=255)                                                               aGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCa                                                                                  >  5:84066/1‑149 (MQ=255)                                                                     gggTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGgcc                                                                            <  5:323101/149‑1 (MQ=255)                                                                        tCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGgccgcc                                                                         <  6:84066/149‑1 (MQ=255)                                                                            ggCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGc                                                                     >  2:115087/1‑149 (MQ=255)                                                                                 gACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTc                                                                                         <  8:210966/124‑1 (MQ=255)                                                                                 gACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGaa                                                                 >  1:399412/1‑148 (MQ=255)                                                                                       tCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCgtg                                                          <  7:380613/149‑1 (MQ=255)                                                                                       tCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCgtg                                                          >  6:261343/1‑149 (MQ=255)                                                                                          cgTTGGCCCGCTCCTACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGt                                                                       <  3:185547/133‑1 (MQ=255)                                                                                          cgTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCg                                                       <  5:261343/149‑1 (MQ=255)                                                                                          cgTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCAACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCg                                                       <  4:47829/149‑1 (MQ=255)                                                                                          cgTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCAACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCg                                                       <  6:30877/149‑1 (MQ=255)                                                                                          cgTTGGCCCGCTCCAACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGt                                                                       >  4:185547/1‑133 (MQ=255)                                                                                                cccGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCcg                                                 >  6:390908/1‑149 (MQ=255)                                                                                                 ccGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCTGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACAt                                                              <  1:8902/135‑1 (MQ=255)                                                                                                 ccGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCg                                                                                        >  6:301097/1‑109 (MQ=255)                                                                                                 ccGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCAGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACAt                                                              >  2:8902/1‑135 (MQ=255)                                                                                                 ccGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGAGCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCg                                                                                        <  5:301097/109‑1 (MQ=255)                                                                                                       ccACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCcgctgcgct                                          >  6:49081/1‑149 (MQ=255)                                                                                                       ccACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCcgctgcgct                                          >  4:75494/1‑149 (MQ=255)                                                                                                       ccACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCcgctgcgc                                           <  5:390908/148‑1 (MQ=255)                                                                                                             cACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTc                                                                                         <  6:156906/96‑1 (MQ=255)                                                                                                             cACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTc                                                                                         >  5:156906/1‑96 (MQ=255)                                                                                                               cccGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCgg                                  <  1:115087/149‑1 (MQ=255)                                                                                                                ccGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGc                                 >  6:248824/1‑149 (MQ=255)                                                                                                                  gcgcgcgcgcCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTg                               <  2:399412/149‑1 (MQ=255)                                                                                                                    gcgcgcgcCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGAt                             <  5:49081/149‑1 (MQ=255)                                                                                                                    gcgcgcgcCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGAt                             <  5:240295/149‑1 (MQ=255)                                                                                                                      gcgcgcCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATg                            >  5:282238/1‑148 (MQ=255)                                                                                                                           ccATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCg                      >  7:50082/1‑149 (MQ=255)                                                                                                                            cATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCAGGCCGTCATCCGCCAGTTTTGCCTCGGCCTCCCTGGCCGCCCTGCCGAACCTAGTGTACGTGCCGCTGCGCTGGGCCCGGCGGATGTACtag                      >  1:318771/1‑146 (MQ=255)                                                                                                                            cATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCg                      >  2:314768/1‑148 (MQ=255)                                                                                                                             aTGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCgg                     >  3:97431/1‑148 (MQ=255)                                                                                                                                gCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGaaa                  <  5:248824/148‑1 (MQ=255)                                                                                                                                    gggTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGAc              <  5:5746/148‑1 (MQ=255)                                                                                                                                    gggTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACg             >  8:246996/1‑149 (MQ=255)                                                                                                                                      gTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGtt           <  5:191193/149‑1 (MQ=255)                                                                                                                                        cATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTg          >  5:220868/1‑148 (MQ=255)                                                                                                                                         aTCTGACCCGACTGGTTGATCACCTAGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGga        <  4:61949/149‑1 (MQ=255)                                                                                                                                          tCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTg          >  1:12842/1‑146 (MQ=255)                                                                                                                                           cTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGgag       >  3:363381/1‑148 (MQ=255)                                                                                                                                              aCCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTg          <  8:188990/142‑1 (MQ=255)                                                                                                                                              aCCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTg          >  7:188990/1‑142 (MQ=255)                                                                                                                                              aCCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGaaa    <  4:144501/148‑1 (MQ=255)                                                                                                                                                ccGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATc  >  5:209857/1‑148 (MQ=255)                                                                                                                                                ccGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATc  <  6:209857/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CCGCACGATCATGTTGTGCCGCCATAGCGCGCCATTGCTCTGCCGCCACCCCTGCACCGTGTAGTTCAGGGTCATGGCCCGACTGATCGCGTTGGCCCGCTCCCACTCCACCCGCGCGCGCGCCATGGCTGGGGTCATCTGACCCGACTGGTTGATCACCTTGACCCGCCGTCGAGCCACGCGGTCATCCGCCAGTTTGGCCTCGACCTCACTGGCCGCCGTGCCGAACATCGTGTCGGTGCCGCTGCGCTGGCCCCGGCTGATGTACTCGGAAAAGACGTTGGAGAAATC  >  NC_002947/4386630‑4386920

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