Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,502,922	C→G	5.1%	A184G (GCC→GGC)	xanP →	xanthine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,502,922	0	C	G	5.1%	252.1 / 3.3	78	A184G (GCC→GGC)	xanP	xanthine permease
Reads supporting (aligned to +/- strand):  ref base C (37/37);  new base G (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

GTTCGTCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGA  >  NC_002947/4502776‑4503067                                                                                                                                                    |                                                                                                                                                  gttcgtCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGcc                                                                                                                                                  >  4:395603/1‑148 (MQ=255) gttcgtCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGcc                                                                                                                                                  >  1:199186/1‑148 (MQ=255) gttcgtCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACACCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGcc                                                                                                                                                  >  7:352169/1‑148 (MQ=255)    cgtcgAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAAt                                                                                                                                               <  2:164826/148‑1 (MQ=255)     gtcgAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATgg                                                                                                                                             >  5:64039/1‑149 (MQ=255)        gagaTGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCAc                                                                                                                                          >  1:75584/1‑149 (MQ=255)        gagaTGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCAc                                                                                                                                          <  2:75584/149‑1 (MQ=255)          gaTGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCt                                                                                                                                        >  3:335429/1‑149 (MQ=255)            tGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATTGGCCATGGCAccaaa                                                                                                                                      <  1:134252/149‑4 (MQ=255)             gggCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCttt                                                                                                                                      <  8:239273/148‑1 (MQ=255)              ggCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTg                                                                                                                                     <  8:349999/148‑1 (MQ=255)              ggCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGc                                                                                                                                    <  2:199186/149‑1 (MQ=255)               gCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTg                                                                                                                                     >  1:31708/1‑147 (MQ=255)               gCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTg                                                                                                                                     <  2:31708/147‑1 (MQ=255)                  tgtCGCGCATCCTGCCCTTCGTCAAACGCCTGATAACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAAt                                                                                                                                >  5:215678/1‑149 (MQ=255)                   gtCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATg                                                                                                                               >  4:334822/1‑149 (MQ=255)                    tCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATgg                                                                                                                              <  5:33924/149‑1 (MQ=255)                    tCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGG‑GCCATGGCCAATGGCACCTTTGCCAATgg                                                                                                                              <  2:315434/149‑1 (MQ=255)                      gcgcATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGc                                                                                                                             >  6:81816/1‑148 (MQ=255)                           tCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAAc                                                                                                                       <  3:334822/149‑1 (MQ=255)                             cTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACct                                                                                                                     <  3:395603/149‑1 (MQ=255)                               gCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgc                                                                                                                   >  6:372553/1‑149 (MQ=255)                                cccTTCGTCACACGCCTGATCACCCCGCTGGGGACTGGCATCGTGGTGGTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgct                                                                                                                  <  8:259839/149‑1 (MQ=255)                                 ccTTCGGCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTTGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgctg                                                                                                                 <  8:220919/149‑1 (MQ=255)                                       tCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCg                                                                                                           >  1:375120/1‑149 (MQ=255)                                          aaCGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGc                                                                                                        >  2:375580/1‑149 (MQ=255)                                             gCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTg                                                                                                     >  4:383088/1‑149 (MQ=255)                                               cTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTgg                                                                                                    >  6:69304/1‑148 (MQ=255)                                                    cACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTgg                                                                                              >  7:349612/1‑149 (MQ=255)                                                     aCCCGGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTgg                                                                                              <  7:142983/148‑1 (MQ=255)                                                      ccccGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTg                                                                                                     >  2:181885/1‑140 (MQ=255)                                                      ccccGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTg                                                                                                     <  1:181885/140‑1 (MQ=255)                                                                 aCTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATTGGCCATGGCAccaaa                                                                                                                                      <  8:42379/96‑4 (MQ=255)                                                                 aCTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATTGGCCATGGCAccaaa                                                                                                                                      >  7:42379/1‑93 (MQ=255)                                                                   tGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTtc                                                                               >  8:62166/1‑149 (MQ=255)                                                                   tGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTtc                                                                               >  6:230004/1‑149 (MQ=255)                                                                     gCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTtc                                                                               <  5:372553/147‑1 (MQ=255)                                                                       aTCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCaa                                                                           >  5:406376/1‑149 (MQ=255)                                                                        tCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAAc                                                                          <  6:215678/149‑1 (MQ=255)                                                                         cGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAAcc                                                                         >  3:376879/1‑149 (MQ=255)                                                                               gtTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCTTCGTTATCAACCGTATca                                                                   >  3:272788/1‑148 (MQ=255)                                                                                ttGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATcc                                                                   <  4:319959/148‑1 (MQ=255)                                                                                  gatgatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCg                                                                 <  5:69304/148‑1 (MQ=255)                                                                                     gatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTg                                                                                                     <  7:107947/109‑1 (MQ=255)                                                                                     gatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTg                                                                                                     >  8:107947/1‑109 (MQ=255)                                                                                     gatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggtgg                                                             >  4:391698/1‑149 (MQ=255)                                                                                     gatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggtgg                                                             >  2:270153/1‑149 (MQ=255)                                                                                      atCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggtggt                                                            >  4:106356/1‑149 (MQ=255)                                                                                      atCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggtgg                                                             <  3:383088/148‑1 (MQ=255)                                                                                        cgtccTGTCGCTGATGTTGGTCGTCCTGATGAGCATGGGCGGTGTCTTTGGT‑GCCTTGGCCTATGGCACCTTTGCCAATGGCTAAAACCTGCTGATGTCGGGCGTGGTTCTTGCGATCATCGTCGTTCTCAACCGTATCCCGgtggtg                                                           <  6:64039/145‑1 (MQ=255)                                                                                              gACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTAGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATg                                                     >  1:293094/1‑149 (MQ=255)                                                                                               aCGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATgcg                                                   <  5:230004/149‑1 (MQ=255)                                                                                               aCGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATgcg                                                   <  2:275553/149‑1 (MQ=255)                                                                                                cGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATgcgc                                                  <  3:221293/149‑1 (MQ=255)                                                                                                cGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATgcgc                                                  <  2:375120/149‑1 (MQ=255)                                                                                                 gCTGATCAAGGTTGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATcc                                                                   >  2:71403/1‑131 (MQ=255)                                                                                                 gCTGATCAAGGTTGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATcc                                                                   <  1:71403/131‑1 (MQ=255)                                                                                                  cTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTCGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTCATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAg                                                >  5:86466/1‑149 (MQ=255)                                                                                                      tCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgctg                                                                                                                 >  3:182250/1‑80 (MQ=255)                                                                                                      tCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgctg                                                                                                                 <  4:182250/80‑1 (MQ=255)                                                                                                            tCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATc                                      <  1:405793/149‑1 (MQ=255)                                                                                                               gCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGt                                    >  8:371635/1‑148 (MQ=255)                                                                                                                  tGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATc                                >  5:346167/1‑149 (MQ=255)                                                                                                                    aTCAGCATGGGCTGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCTTCATc                                <  8:349612/147‑1 (MQ=255)                                                                                                                       aGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCt                           <  1:270153/149‑1 (MQ=255)                                                                                                                        gCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTg                          >  6:61180/1‑149 (MQ=255)                                                                                                                           tGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCg                       >  7:161822/1‑149 (MQ=255)                                                                                                                           tGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCg                       >  5:126617/1‑149 (MQ=255)                                                                                                                           tGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCg                       <  4:335429/149‑1 (MQ=255)                                                                                                                                ggTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCgg                  <  5:81816/149‑1 (MQ=255)                                                                                                                                 gTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGc                 >  7:55371/1‑149 (MQ=255)                                                                                                                                    gcaaaggT‑GCCATTGGCCATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTAt              >  6:387642/6‑149 (MQ=255)                                                                                                                                    gCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTAt              <  3:106356/149‑1 (MQ=255)                                                                                                                                      tttGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCCGTCGGCTATgc            >  3:46872/1‑149 (MQ=255)                                                                                                                                          gT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCt         <  6:406376/148‑1 (MQ=255)                                                                                                                                              ccATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGcc     <  3:125998/149‑1 (MQ=255)                                                                                                                                                aTGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCg    >  3:380876/1‑148 (MQ=255)                                                                                                                                                aTGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCg    >  8:339257/1‑148 (MQ=255)                                                                                                                                                 tGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCg    >  6:236020/1‑147 (MQ=255)                                                                                                                                                 tGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCg    <  5:236020/147‑1 (MQ=255)                                                                                                                                                 tGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGa  <  1:375580/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GTTCGTCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGA  >  NC_002947/4502776‑4503067

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