Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,502,922	C→G	8.6%	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	8.6%	207.4 / 7.6	69	A184G (GCC→GGC)	xanP	xanthine permease
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base G (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 = 9.57e-01

TCGTTCGTCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGC  >  NC_002947/4502774‑4503057                                                                                                                                                       |                                                                                                                                        tcgttcgtCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGc                                                                                                                                         <  4:91070/149‑1 (MQ=255)      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)                                                                                                                                                       |                                                                                                                                        TCGTTCGTCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATG‑GGCGGTGG‑CTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGC  >  NC_002947/4502774‑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: ‑