Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,882,207:1	+A	5.1%	coding (1290/1356 nt)	uacT →	uric acid permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,882,207	1	.	A	5.1%	329.6 / 11.0	79	coding (1290/1356 nt)	uacT	uric acid permease
Reads supporting (aligned to +/- strand):  ref base . (42/33);  new base A (2/2);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

GCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTC  >  NC_002947/4882059‑4882351                                                                                                                                                       |                                                                                                                                                 gccgcAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCGCGGGGTGCGTCCGGAGTTGTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGTAc                                                                                                                                                   <  2:103968/149‑1 (MQ=255)    gcAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTa                                                                                                                                               <  2:324529/149‑1 (MQ=255)      aaCCATGTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTGCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCGACTGTCTGTGCCCTGG‑TGTTGAAC‑GTACt                                                                                                                                             <  7:210124/149‑1 (MQ=255)      aaCCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACt                                                                                                                                             >  5:72599/1‑149 (MQ=255)          aGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTc                                                                                                                                         <  5:165523/149‑1 (MQ=255)          aGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTc                                                                                                                                         <  6:72599/149‑1 (MQ=255)             tGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCaa                                                                                                                                       >  8:211350/1‑148 (MQ=255)               ttggttgCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCAAGGATGTTGAACAGTAC‑GTTCAACa                                                                                                                                     <  2:395485/149‑1 (MQ=255)                   ttgCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCtt                                                                                                                                <  6:176394/149‑1 (MQ=255)                   ttgCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCtt                                                                                                                                >  7:341294/1‑149 (MQ=255)                     gCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTgg                                                                                                                              >  1:199245/1‑149 (MQ=255)                         tGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGt                                                                                                                          <  2:314453/149‑1 (MQ=255)                            gCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGGTGAAC‑GTACTGTTCAACATCCTTGGCGGTGcc                                                                                                                       <  6:399588/149‑1 (MQ=255)                              gTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGa                                                                                                                     >  2:8033/1‑149 (MQ=255)                                cGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGAcc                                                                                                                   <  7:276612/149‑1 (MQ=255)                                    ttCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACcgcg                                                                                                                >  3:355544/1‑148 (MQ=255)                                       gggTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAAt                                                                                                            <  8:162923/149‑1 (MQ=255)                                         gTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATg                                                                                                           >  2:186623/1‑148 (MQ=255)                                         gTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGA‑TGTTGAACAGTAC‑GTTCAACacca                                                                                                                                  >  6:48382/1‑122 (MQ=255)                                         gTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGA‑TGTTGAACAGTAC‑GTTCAACacca                                                                                                                                  <  5:48382/125‑4 (MQ=255)                                          ttGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGca                                                                                                         >  2:352633/1‑149 (MQ=255)                                              tCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAac                                                                                                     <  7:211350/149‑1 (MQ=255)                                              tCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAac                                                                                                     <  2:199245/149‑1 (MQ=255)                                                      gtgCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCt                                                                                              <  1:8033/148‑1 (MQ=255)                                                               gAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCa                                                                                    >  5:180217/1‑149 (MQ=255)                                                                    cTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAAc                                                                                                                                      >  6:317284/1‑94 (MQ=255)                                                                    cTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAAc                                                                                                                                      <  5:317284/94‑1 (MQ=255)                                                                    cTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAacgac                                                                                                  <  6:286827/130‑1 (MQ=255)                                                                    cTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAacgac                                                                                                  >  5:286827/1‑130 (MQ=255)                                                                    cTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAg                                                                               >  7:44400/1‑149 (MQ=255)                                                                       tGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCAGAACGACGTCTGTCACCAGCATTGAg                                                                               <  8:44400/146‑1 (MQ=255)                                                                        gCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGcc                                                                             >  3:241416/1‑147 (MQ=255)                                                                        gCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGcc                                                                             >  7:322090/1‑147 (MQ=255)                                                                         cccAGATGCCACAGTGGATGGAGCCCTTCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGa                                                                                                                     >  7:153003/1‑106 (MQ=255)                                                                         cccAGATGCCACAGTGGATGGAGCCCTTCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGa                                                                                                                     <  8:153003/106‑1 (MQ=255)                                                                         cccAGATGCCACAGGGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAg                                                                          <  7:344800/149‑1 (MQ=255)                                                                          ccAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgg                                                                         >  2:231962/1‑149 (MQ=255)                                                                           cAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAggg                                                                        >  4:147674/1‑149 (MQ=255)                                                                             gATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGTCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGa                                                                                >  3:221354/1‑139 (MQ=255)                                                                             gATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGTCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGa                                                                                <  4:221354/139‑1 (MQ=255)                                                                                 ccACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCg                                                                   <  6:126367/148‑1 (MQ=255)                                                                                  cacaGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGt                                                                 >  8:247048/1‑149 (MQ=255)                                                                                   acaGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgggggg                                                                     >  7:293316/1‑144 (MQ=255)                                                                                   acaGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTg                                                                >  7:381772/1‑149 (MQ=255)                                                                                   acaGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTATGTCACCAGCATTGAGCCAAgggggg                                                                     <  8:293316/144‑1 (MQ=255)                                                                                       tggatggaGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATg                                                            <  3:62495/149‑1 (MQ=255)                                                                                         gatggaGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGc                                                          >  4:123213/1‑149 (MQ=255)                                                                                         gatggaGCCAATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCACCATTGAGCCAAGGGGGGCGGTGCAGGGc                                                          >  6:279146/1‑149 (MQ=255)                                                                                                cccATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCt                                                   >  4:88030/1‑149 (MQ=255)                                                                                                  cATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCg                                                 >  3:237574/1‑149 (MQ=255)                                                                                                  cATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCg                                                 >  6:11145/1‑149 (MQ=255)                                                                                                   aTCACCCCAAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACTATGTCTGTAACCAGCATTGAGCCAAGGGCGGCGGTGACTGCCTGCCCCTCGt                                                >  5:392243/1‑149 (MQ=255)                                                                                                   aTCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGt                                                <  1:352633/149‑1 (MQ=255)                                                                                                    tCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGtt                                               <  4:355544/149‑1 (MQ=255)                                                                                                       gtcactGCGGTATGGCCTGGGCCACGGTCTGTGCCCTGG‑TGTTGATC‑GTACTGTTCAACTTCCTTGGCGGTGCCGAGCGCGCAATGCACAACGTCGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCa                                             <  6:140668/146‑1 (MQ=255)                                                                                                         cacaGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGTGGCGGTGCATGGCTGCCCCTCGTTCAcc                                           >  5:386767/1‑148 (MQ=255)                                                                                                          acaGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACctc                                         >  1:361066/1‑149 (MQ=255)                                                                                                          acaGCGGTATCGCCATCGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTCTACAACATCCTTGGCGGGGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGGTCAca                                           >  5:140668/1‑146 (MQ=255)                                                                                                           caGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGc                                                                    <  2:407680/121‑1 (MQ=255)                                                                                                           caGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGc                                                                    >  1:407680/1‑121 (MQ=255)                                                                                                             gCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGTATGGCTGCCCCTCGTTCACCTCTGc                                      <  8:322090/149‑1 (MQ=255)                                                                                                              cGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCa                                     <  3:147674/149‑1 (MQ=255)                                                                                                                gTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAg                                    >  3:5713/1‑148 (MQ=255)                                                                                                                      ccATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGc                             >  3:350574/1‑149 (MQ=255)                                                                                                                      ccATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCATCCTTCACCTGTGCAGTAAccc                              >  8:252203/1‑147 (MQ=255)                                                                                                                        aTGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgggg                                                                       >  2:49551/1‑105 (MQ=255)                                                                                                                        aTGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTa                           <  1:231962/149‑1 (MQ=255)                                                                                                                        aTGGCCACGGTCTGTGCCCTGG‑TGTTCAAC‑GTACTGTTCAACATCCTTGGCGGTGCGTACCGCGAAATGCAATACGACGTCTGTCACCAGCATTGAGCCAAgggg                                                                       <  1:49551/105‑1 (MQ=255)                                                                                                                             cACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAAc                                >  8:26/1‑139 (MQ=255)                                                                                                                             cACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAAc                                <  7:26/139‑1 (MQ=255)                                                                                                                                 gTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGa                  <  8:341294/149‑1 (MQ=255)                                                                                                                                      tgCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACccgc             >  7:44498/1‑149 (MQ=255)                                                                                                                                        cccTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACccgccg           >  7:372341/1‑149 (MQ=255)                                                                                                                                        cccTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACccgccg           >  5:172755/1‑149 (MQ=255)                                                                                                                                         ccTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACccgccg           >  8:209774/1‑148 (MQ=255)                                                                                                                                           tGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGaa        >  1:19362/1‑149 (MQ=255)                                                                                                                                               tgtTGAACAGTAC‑GTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGCCCCGCCGGAATgg     >  7:345857/1‑149 (MQ=255)                                                                                                                                                   gAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTc  <  3:123213/148‑1 (MQ=255)                                                                                                                                                    aaC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCa                                     >  7:38530/1‑112 (MQ=255)                                                                                                                                                    aaC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCa                                     <  8:38530/112‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 GCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGG‑TGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTC  >  NC_002947/4882059‑4882351

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑