Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,502,924	A→C	5.1%	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	5.1%	253.0 / 3.4	79	N185H (AAT→CAT)	xanP	xanthine permease
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base C (2/2);  total (40/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 = 9.66e-01

CGTCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGATACC  >  NC_002947/4502779‑4503071                                                                                                                                                   |                                                                                                                                                    cgtcgAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAAt                                                                                                                                                   <  2:164826/148‑1 (MQ=255)  gtcgAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATgg                                                                                                                                                 >  5:64039/1‑149 (MQ=255)     gagaTGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCAc                                                                                                                                              <  2:75584/149‑1 (MQ=255)     gagaTGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCAc                                                                                                                                              >  1:75584/1‑149 (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                                                                                                                                         <  2:31708/147‑1 (MQ=255)            gCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTg                                                                                                                                         >  1:31708/1‑147 (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                                                                                                         <  1:181885/140‑1 (MQ=255)                                                   ccccGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTg                                                                                                         >  2:181885/1‑140 (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                                                                 >  2:270153/1‑149 (MQ=255)                                                                                  gatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggtgg                                                                 >  4:391698/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                                                                       <  1:71403/131‑1 (MQ=255)                                                                                              gCTGATCAAGGTTGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATcc                                                                       >  2:71403/1‑131 (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                           <  4:335429/149‑1 (MQ=255)                                                                                                                        tGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCg                           >  7:161822/1‑149 (MQ=255)                                                                                                                        tGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCg                           >  5:126617/1‑149 (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        >  8:339257/1‑148 (MQ=255)                                                                                                                                             aTGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCg        >  3:380876/1‑148 (MQ=255)                                                                                                                                              tGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCg        <  5:236020/147‑1 (MQ=255)                                                                                                                                              tGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCg        >  6:236020/1‑147 (MQ=255)                                                                                                                                              tGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGa      <  1:375580/149‑1 (MQ=255)                                                                                                                                                  cAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGATAcc  >  2:95479/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CGTCGAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGT‑GCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGATACC  >  NC_002947/4502779‑4503071

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