Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,502,920	G→T	7.9%	M183I (ATG→ATT)	xanP →	xanthine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,502,920	0	G	T	7.9%	235.1 / 8.7	76	M183I (ATG→ATT)	xanP	xanthine permease
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base T (2/4);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-01

GCTCGTTCGTCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCG  >  NC_002947/4502772‑4503065                                                                                                                                                        |                                                                                                                                                  gCTCGTTCGTCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATg                                                                                                                                                   <  4:295527/149‑1 (MQ=255) gCTCGTTCGTCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATg                                                                                                                                                   <  2:2215/149‑1 (MQ=255)          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:230979/1‑149 (MQ=255)                       gtCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATg                                                                                                                             >  6:228968/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                                                                                           >  2:203530/1‑149 (MQ=255)                                                         aCCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGc                                                                                           >  1:159538/1‑149 (MQ=255)                                                         aCCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGc                                                                                           <  4:240080/149‑1 (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                                                    <  8:301857/149‑1 (MQ=255)                                                                                                cTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGAt                                                    <  7:7787/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                                                                  >  6:263960/1‑127 (MQ=255)                                                                                                        gATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATc                                                                  <  5:263960/127‑1 (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                                    <  6:349819/149‑1 (MQ=255)                                                                                                                tCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATc                                    <  5:58073/149‑1 (MQ=255)                                                                                                                tCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATc                                    <  3:208538/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            <  2:7958/149‑1 (MQ=255)                                                                                                                                           gCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTAt            <  7:214161/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  >  2:141950/1‑149 (MQ=255)                                                                                                                                                     cATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCg  <  3:99609/149‑1 (MQ=255)                                                                                                                                                        ggCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGc                      >  7:34898/1‑126 (MQ=255)                                                                                                                                                        ggCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGc                      <  8:34898/126‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                  GCTCGTTCGTCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCC‑TGATC‑AGCAT‑GGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCG  >  NC_002947/4502772‑4503065

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