Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,502,924	A→C	7.5%	N185H (AAT→CAT)	xanP →	xanthine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,502,924	0	A	C	7.5%	243.3 / 9.8	80	N185H (AAT→CAT)	xanP	xanthine permease
Reads supporting (aligned to +/- strand):  ref base A (36/38);  new base C (2/4);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

TCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGATA  >  NC_002947/4502781‑4503069                                                                                                                                                   |                                                                                                                                                  tcgAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATgg                                                                                                                                               >  5:349819/1‑148 (MQ=255)    agaTGGGAGTGTCGCGGATCCTGCCCTTCGTCAAACGCCTGATAACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATTGGCCATGGCAcc                                                                                                                                           <  5:306532/149‑1 (MQ=255)       tGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATTGGCCATGGCAccaaa                                                                                                                                        <  5:13946/149‑4 (MQ=255)       tGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATTGGCCATGGCAccaaa                                                                                                                                        <  6:244819/149‑4 (MQ=255)       tGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCttt                                                                                                                                        <  8:87437/149‑1 (MQ=255)        gggCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTg                                                                                                                                       >  6:335991/1‑149 (MQ=255)             tgtCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAAt                                                                                                                                  <  4:126533/149‑1 (MQ=255)             tgtCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAAt                                                                                                                                  <  6:372294/149‑1 (MQ=255)             tgtCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAAt                                                                                                                                  <  6:147385/149‑1 (MQ=255)              gtCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATg                                                                                                                                 >  6:228968/1‑149 (MQ=255)              gtCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATg                                                                                                                                 >  6:230979/1‑149 (MQ=255)                   gcATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGa                                                                                                                             >  3:42015/1‑148 (MQ=255)                     aTCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGaaaa                                                                                                                          >  7:48964/1‑149 (MQ=255)                       ccTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAAcc                                                                                                                        >  6:266252/1‑149 (MQ=255)                        cTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACct                                                                                                                       <  8:48964/149‑1 (MQ=255)                           cccTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgc                                                                                                                     <  6:39765/148‑1 (MQ=255)                           cccTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGGCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctg                                                                                                                      <  5:266252/147‑1 (MQ=255)                                  tCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCg                                                                                                             >  3:240080/1‑149 (MQ=255)                                    aaaCGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCggg                                                                                                           >  6:58073/1‑149 (MQ=255)                                     aaCGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGc                                                                                                          <  4:275127/149‑1 (MQ=255)                                          cTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCAGCTGATGTCGGGCGTgg                                                                                                      >  2:261567/1‑148 (MQ=255)                                             aTCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACt                                                                                                  <  2:285568/149‑1 (MQ=255)                                                aCCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGc                                                                                               >  1:159538/1‑149 (MQ=255)                                                aCCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGc                                                                                               <  4:240080/149‑1 (MQ=255)                                                aCCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGc                                                                                               >  2:203530/1‑149 (MQ=255)                                                 ccccGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGc                                                                                               >  4:194911/1‑148 (MQ=255)                                                  cccGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGa                                                                                             >  5:247397/1‑149 (MQ=255)                                                  cccGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGa                                                                                             >  4:208538/1‑149 (MQ=255)                                                       tGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCAAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGatcatc                                                                                        >  3:255292/1‑149 (MQ=255)                                                          tGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCAtcgtc                                                                                     <  4:42015/149‑1 (MQ=255)                                                            aCTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCAtcgtcg                                                                                    <  7:246255/148‑1 (MQ=255)                                                              tGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATCAAGCAT‑GGGCGTTGGATTTGTTGCCATGGCAAATGGCACATTTGCCAATGGCGAAAACATGCTGATGTCGGGCGTGGTACTGGCGATCTTCGTCGtt                                                                                  >  1:97122/1‑149 (MQ=255)                                                                gCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTtct                                                                                <  6:345300/148‑1 (MQ=255)                                                                gCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTtc                                                                                 <  3:329063/147‑1 (MQ=255)                                                                          gtTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATcc                                                                     >  7:301857/1‑149 (MQ=255)                                                                             gatgatCGGCCTGCCGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTTTTGCCATGGCCAATGGCACCTTTGCAAATGGCGAAAACCTGCTGATGTCGGGCGGTGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCg                                                                   <  6:343120/148‑1 (MQ=255)                                                                             gatgatCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCgg                                                                  >  4:44275/1‑149 (MQ=255)                                                                              atgatCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCAGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggt                                                                 <  1:261567/149‑1 (MQ=255)                                                                               tgatCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCATGGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggt                                                                 >  1:259919/1‑149 (MQ=255)                                                                                  tCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTg                                                                                                       <  7:199928/107‑1 (MQ=255)                                                                                  tCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTg                                                                                                       >  8:199928/1‑107 (MQ=255)                                                                                    ggCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGgtgt                                                            <  3:194911/148‑1 (MQ=255)                                                                                      gcccatgCTGATC‑AGGCCGACCTTGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGa                                                         >  4:99609/7‑149 (MQ=255)                                                                                      ccTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGa                                                         <  5:335991/149‑1 (MQ=255)                                                                                       cccatgCTGATC‑AGGCCGACCTTGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGAt                                                        >  5:323303/6‑149 (MQ=255)                                                                                       cTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGAt                                                        <  7:7787/149‑1 (MQ=255)                                                                                       cTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGAt                                                        <  8:301857/149‑1 (MQ=255)                                                                                          atgCTGATC‑AGGCCGACCTTGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATgcg                                                     >  3:264813/3‑149 (MQ=255)                                                                                           cGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATgcgc                                                    <  2:97122/149‑1 (MQ=255)                                                                                            gCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGtt                                                                                  >  7:173452/1‑118 (MQ=255)                                                                                            gCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGtt                                                                                  <  8:173452/118‑1 (MQ=255)                                                                                            gCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATgcg                                                     >  1:229815/1‑147 (MQ=255)                                                                                               gATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATc                                                                      <  5:263960/127‑1 (MQ=255)                                                                                               gATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATc                                                                      >  6:263960/1‑127 (MQ=255)                                                                                                 tCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTgc                                              >  1:73229/1‑149 (MQ=255)                                                                                                   aaGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTgcg                                             <  3:44275/148‑1 (MQ=255)                                                                                                       tCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATc                                        <  3:208538/149‑1 (MQ=255)                                                                                                       tCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATc                                        <  5:58073/149‑1 (MQ=255)                                                                                                       tCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATc                                        <  6:349819/149‑1 (MQ=255)                                                                                                         gaccttGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCg                                       <  6:323303/144‑1 (MQ=255)                                                                                                            c‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATGGTCAt                                   >  8:251584/1‑149 (MQ=255)                                                                                                            c‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCAt                                   >  1:7958/1‑149 (MQ=255)                                                                                                                aTC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATc                                  <  4:255292/147‑1 (MQ=255)                                                                                                                           ggCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCgg                        <  2:73229/148‑1 (MQ=255)                                                                                                                                  gCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTAt                <  7:214161/149‑1 (MQ=255)                                                                                                                                  gCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTAt                <  2:7958/149‑1 (MQ=255)                                                                                                                                  gCTTTGGTGCCATGGCCAATGGCACCTTTCCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTAt                >  2:320642/1‑149 (MQ=255)                                                                                                                                   cTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATg               <  2:159538/149‑1 (MQ=255)                                                                                                                                    tttGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATgc              >  2:185755/1‑149 (MQ=255)                                                                                                                                    tttGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATgc              >  8:229249/1‑149 (MQ=255)                                                                                                                                    aaaggTGCCATTGGCCATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATgc              >  5:155267/4‑149 (MQ=255)                                                                                                                                       ggTGCCATTGGCCATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCt           >  6:351377/1‑149 (MQ=255)                                                                                                                                        gTGCCATTGGCCATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTg          <  6:155267/149‑1 (MQ=255)                                                                                                                                         tGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTgg         <  5:230979/149‑1 (MQ=255)                                                                                                                                            cATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCg      <  3:99609/149‑1 (MQ=255)                                                                                                                                            cATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCg      >  2:141950/1‑149 (MQ=255)                                                                                                                                               ggCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGc                          <  8:34898/126‑1 (MQ=255)                                                                                                                                               ggCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGc                          >  7:34898/1‑126 (MQ=255)                                                                                                                                                gCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGATa  >  4:315664/1‑149 (MQ=255)                                                                                                                                                gCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGATa  <  4:264813/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGATA  >  NC_002947/4502781‑4503069

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