Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,502,924	A→C	8.4%	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	8.4%	212.1 / 7.1	69	N185H (AAT→CAT)	xanP	xanthine permease
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base C (2/4);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGC  >  NC_002947/4502779‑4503057                                                                                                                                                    |                                                                                                                                      cgtcgAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAAt                                                                                                                                     <  2:329308/148‑1 (MQ=255) cgtcgAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAAt                                                                                                                                     <  7:37711/148‑1 (MQ=255) cgtcgAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGACGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATg                                                                                                                                    >  7:10991/1‑149 (MQ=255)  gtcgAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATg                                                                                                                                    <  8:53487/148‑1 (MQ=255)  gtcgAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATg                                                                                                                                    <  4:275502/148‑1 (MQ=255)   tcgAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATgg                                                                                                                                   >  1:4817/1‑148 (MQ=255)     gagaTGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATTGGCCATGGCAc                                                                                                                                <  4:303181/149‑1 (MQ=255)         tGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGGATG‑GGCGGTGG‑CTTTGGTGCCATTGGCCATGGCAccaaa                                                                                                                            <  1:12064/149‑4 (MQ=255)         tGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATTGGCCATGGCAccaaa                                                                                                                            <  3:134356/149‑4 (MQ=255)             cGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGcc                                                                                                                         >  2:284767/1‑148 (MQ=255)               tgtCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAAt                                                                                                                      <  7:309509/149‑1 (MQ=255)               tgtCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAAt                                                                                                                      <  7:283519/149‑1 (MQ=255)                   gcgcATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGc                                                                                                                   >  7:261536/1‑148 (MQ=255)                         ccTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAAcc                                                                                                            >  8:100454/1‑149 (MQ=255)                           tGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACct                                                                                                           >  2:215629/1‑148 (MQ=255)                            gCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgc                                                                                                         >  5:242841/1‑149 (MQ=255)                             cccTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgc                                                                                                         <  4:162649/148‑1 (MQ=255)                             cccTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctg                                                                                                          <  8:261536/147‑1 (MQ=255)                               cTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGa                                                                                                      >  1:226482/1‑149 (MQ=255)                                ttCGTCAAACGCCTGATCACCCCGCTGTTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CAAAGGTGCCATTGGCCATGGCAccaaa                                                                                                                            <  8:264432/126‑4 (MQ=38)                                ttCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CAAAGGTGCCATTGGCCATGGCAccaaa                                                                                                                            >  7:264432/1‑123 (MQ=38)                                   gTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTc                                                                                                  <  8:38884/149‑1 (MQ=255)                                      aaaCGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCttt                                                                                                                            >  6:17432/1‑120 (MQ=255)                                      aaaCGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCggg                                                                                               >  2:274926/1‑149 (MQ=255)                                        aCGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCg                                                                                             <  6:242841/149‑1 (MQ=255)                                         cGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGt                                                                                            <  5:63766/149‑1 (MQ=255)                                           ccTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTgg                                                                                          <  7:181666/149‑1 (MQ=255)                                            cTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTg                                                                                           <  8:133565/147‑1 (MQ=255)                                            cTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTg                                                                                           >  7:133565/1‑147 (MQ=255)                                            cTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGCCCAATGGCACCATTGCCAAGGGCGAAAACCAGCTAAAGACGGGCGTgg                                                                                          >  8:49015/1‑148 (MQ=255)                                             tGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgc                                                                                                         <  1:208709/132‑1 (MQ=255)                                             tGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgc                                                                                                         >  2:208709/1‑132 (MQ=255)                                             tGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTa                                                                                        >  7:166097/1‑149 (MQ=255)                                                 cACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTg                                                                                     <  7:313093/148‑1 (MQ=255)                                                    cccGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGa                                                                                 >  1:253251/1‑149 (MQ=255)                                                       gCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGatca                                                                              >  7:298027/1‑149 (MQ=255)                                                                tGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGGACTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGtt                                                                      >  1:271676/1‑149 (MQ=255)                                                                     tCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAAc                                                                >  7:71287/1‑149 (MQ=255)                                                                             ttGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATccc                                                        >  3:229329/1‑149 (MQ=255)                                                                               gatgatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAAcc                                                               >  6:211906/1‑140 (MQ=255)                                                                               gatgatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAAcc                                                               <  5:211906/140‑1 (MQ=255)                                                                               gatgatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCgg                                                      >  6:249707/1‑149 (MQ=255)                                                                               gatgatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCgg                                                      >  7:199667/1‑149 (MQ=255)                                                                               gatgatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCg                                                       <  2:4817/148‑1 (MQ=255)                                                                                atgatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggt                                                     >  6:4480/1‑149 (MQ=255)                                                                                  gatCGGCCTGACGCTGATCAAGGGCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggtg                                                    <  7:49015/148‑1 (MQ=255)                                                                                      ggCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCTTCGTCGTTCTCAACCGTATCCCGGTGgtgt                                                <  1:29441/148‑1 (MQ=255)                                                                                      ggCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGgtgt                                                <  5:4480/148‑1 (MQ=255)                                                                                        ccTGACGCTGATCAAGGTCGGCCTGATCAGCATGAGGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGtgtg                                               >  1:46121/1‑148 (MQ=255)                                                                                         cTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGAt                                            <  8:183279/149‑1 (MQ=255)                                                                                            aCGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATgc                                          <  8:71287/148‑1 (MQ=255)                                                                                                  aTCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAACACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTg                                   >  1:339810/1‑149 (MQ=255)                                                                                                        gTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCAt                             >  3:54585/1‑149 (MQ=255)                                                                                                         tCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATc                            <  4:54585/149‑1 (MQ=255)                                                                                                         tCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATc                            <  2:253251/149‑1 (MQ=255)                                                                                                         tCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATc                            <  1:274926/149‑1 (MQ=255)                                                                                                           ggCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCgg                                                      <  3:254323/121‑1 (MQ=255)                                                                                                           ggCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCgg                                                      >  4:254323/1‑121 (MQ=255)                                                                                                                 aTCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGc                    >  6:203189/1‑149 (MQ=255)                                                                                                                  tCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGcc                   <  8:298027/149‑1 (MQ=255)                                                                                                                     gCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTg                <  8:166097/149‑1 (MQ=255)                                                                                                                        tG‑GGCGGTGG‑CTTTTGTGCCATGGCCAATGGCACCATTTGCAATTGCGAAAACCAGCTGATGTCGGGCGGGGTACTGGGGATCATCGTCGTTATAAACCGTATCACGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATAGCCCAGGCg             >  8:256984/1‑149 (MQ=255)                                                                                                                                tGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCt      <  8:10991/149‑1 (MQ=255)                                                                                                                                  g‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTAt    <  2:226482/149‑1 (MQ=255)                                                                                                                                    cTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATg   <  2:271676/149‑1 (MQ=255)                                                                                                                                     tttGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATgc                                          <  8:147963/109‑1 (MQ=255)                                                                                                                                     tttGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATgc                                          >  7:147963/1‑109 (MQ=255)                                                                                                                                     aaaggTGCCATTGGCCATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATgc  >  4:265539/4‑149 (MQ=255)                                                                                                                                            ccATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggt                                                     <  2:319022/91‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                      CGTCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGC  >  NC_002947/4502779‑4503057

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