Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,018	A→G	T264A (ACG→GCG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,018	0	A	G	98.8%	313.3 / ‑3.9	83	T264A (ACG→GCG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (44/38);  total (45/38)
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.86e-01

GAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGT  >  NC_000913/4618871‑4619163                                                                                                                                                    |                                                                                                                                                   gAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCgg                                                                                                                                                 >  1:465832/1‑151 (MQ=255) gAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCgg                                                                                                                                                 >  2:719357/1‑151 (MQ=255)  aaCTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCggg                                                                                                                                                >  2:330549/1‑151 (MQ=255)  aaCTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGACGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCggg                                                                                                                                                >  4:27095/1‑151 (MQ=255)    ctctctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCggggg                                                                                                                                              >  3:220487/1‑149 (MQ=255)    ctctctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCggggg                                                                                                                                              >  2:664562/1‑149 (MQ=255)       tctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGaa                                                                                                                                            >  3:756655/1‑150 (MQ=255)          gAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAat                                                                                                                                           <  3:1184097/148‑1 (MQ=255)                   gCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCCCCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCcg                                                                                                                                <  2:1335746/150‑1 (MQ=255)                   gCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCcgc                                                                                                                               >  2:1241094/1‑151 (MQ=255)                    ccGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCcgcg                                                                                                                              <  2:280090/151‑1 (MQ=255)                     cGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGt                                                                                                                             <  2:1128483/151‑1 (MQ=255)                     cGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAATATCGTAAccc                                                                                                                                 >  3:314698/1‑147 (MQ=255)                        aGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTg                                                                                                                          <  3:1068566/151‑1 (MQ=255)                              tgttgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGa                                                                                                                     <  4:433520/150‑1 (MQ=255)                              tgttgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGa                                                                                                                     <  2:575655/150‑1 (MQ=255)                                ttgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGt                                                                                                                  >  3:1327986/1‑151 (MQ=255)                                   gCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑c                                                                                                              >  3:1109868/1‑151 (MQ=255)                                    cGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑c                                                                                                              >  2:295752/1‑150 (MQ=255)                                    cGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑Cg                                                                                                             >  2:436601/1‑151 (MQ=255)                                    cGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGGAACCCGCGTCTGTTTGATGTCA‑c                                                                                                              >  4:169754/1‑150 (MQ=255)                                          tAACGGCGCTGGCGTGCGCCCCACCGCGCTGCTCACCGACATGAATCAGGGACTGGCCCCCAGGGCAGGGAACGCGGGTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGGAAATCGTAAACCCCGGCTGTTTGTTGTCA‑CGGTGgcg                                                                                                       >  4:229855/1‑151 (MQ=255)                                           aaCGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGgcg                                                                                                       <  3:690165/150‑1 (MQ=255)                                           aaCGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGgcg                                                                                                       <  2:820001/150‑1 (MQ=255)                                           aaCGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGgc                                                                                                        <  3:181172/149‑1 (MQ=255)                                               gcgcTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTTATGTCA‑AGATGGCGCtgt                                                                                                   >  1:442389/1‑150 (MQ=255)                                                 gcTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCtgtg                                                                                                  >  1:326356/1‑149 (MQ=255)                                                      cGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGgag                                                                                           <  2:749114/151‑1 (MQ=255)                                                       gTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTgg                                                                                             >  2:443529/1‑148 (MQ=255)                                                        tGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTgg                                                                                             >  1:706/1‑147 (MQ=255)                                                        tGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTgg                                                                                             >  3:367867/1‑147 (MQ=255)                                                             accaccGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTTAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGa                                                                                    <  3:1192219/151‑1 (MQ=255)                                                             accaccGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAAACCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTg                                                                                     >  2:279002/1‑150 (MQ=255)                                                             accaccGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTg                                                                                     >  3:293561/1‑150 (MQ=255)                                                               caccGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGAtc                                                                                  <  3:1191292/151‑1 (MQ=255)                                                                   gcgcTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCg                                                                              <  4:187540/151‑1 (MQ=255)                                                                    cgcTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCgg                                                                             <  2:499031/151‑1 (MQ=255)                                                                    cgcTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCgg                                                                             <  2:1291492/151‑1 (MQ=255)                                                                       tgctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGGAACCCGCGACCGTTTGATGTTACCGAGGGC‑CTGTGGGTGGAGGTGCTGGATTCCGGcc                                                                           >  3:839786/1‑149 (MQ=255)                                                                        gctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCaaa                                                                         <  2:314668/151‑1 (MQ=255)                                                                        gctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCaa                                                                          <  2:1184418/150‑1 (MQ=255)                                                                        gctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCa                                                                           <  4:14630/149‑1 (MQ=255)                                                                         ctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCaaa                                                                         >  3:302274/1‑150 (MQ=255)                                                                           cACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTg                                                                      >  2:833102/1‑151 (MQ=255)                                                                           cACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTg                                                                      <  3:1001809/151‑1 (MQ=255)                                                                               gACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTgg                                                                     >  2:387125/1‑148 (MQ=255)                                                                                    gAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGa                                                              <  3:817415/150‑1 (MQ=255)                                                                                    gAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGAt                                                             <  2:658213/151‑1 (MQ=255)                                                                                      aTCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGAt                                                             >  2:399146/1‑149 (MQ=255)                                                                                      aTCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATg                                                            >  1:249162/1‑150 (MQ=255)                                                                                       tCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGAc                                                          >  2:672866/1‑151 (MQ=255)                                                                                        cAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGAc                                                          <  2:1280109/150‑1 (MQ=255)                                                                                         aGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGc                                                        >  2:634279/1‑151 (MQ=255)                                                                                         aGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGc                                                        >  2:61577/1‑151 (MQ=255)                                                                                             aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGGTGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGcc                                                       >  3:256116/1‑148 (MQ=255)                                                                                             aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGa                                                     <  2:981739/150‑1 (MQ=255)                                                                                               tGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgc                                                  >  1:167114/1‑151 (MQ=255)                                                                                               tGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGAGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGAACCCCGGCAAACTGGCGAAAGATGACGCCgga                                                    >  2:706627/1‑147 (MQ=255)                                                                                                ggCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcg                                                 >  3:1114041/1‑151 (MQ=255)                                                                                                ggCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgc                                                  <  4:19530/150‑1 (MQ=255)                                                                                                 gCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgc                                                <  2:708547/151‑1 (MQ=255)                                                                                                 gCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGGAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCCGCAAACTGGCGAAAGATGACGcc                                                       >  3:250802/1‑144 (MQ=255)                                                                                                         tGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGaa                                         >  2:973235/1‑150 (MQ=255)                                                                                                           cAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCCAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAtt                                      >  3:432016/1‑151 (MQ=255)                                                                                                           cAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAtt                                      <  4:46705/151‑1 (MQ=255)                                                                                                             ggTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGc                                    <  3:832463/151‑1 (MQ=255)                                                                                                                aaCGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGGCTGTTTGATGTCC‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAgg                                 >  4:464781/1‑151 (MQ=255)                                                                                                                    cgGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGt                             <  2:173910/151‑1 (MQ=255)                                                                                                                      gTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGc                           <  3:664582/151‑1 (MQ=255)                                                                                                                       ttGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGATGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCt                          >  1:293321/1‑151 (MQ=255)                                                                                                                         gAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTgg                        <  4:413555/151‑1 (MQ=255)                                                                                                                           aGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGGGCTgg                        >  3:48444/1‑149 (MQ=255)                                                                                                                              tCGTGAAGCGGTGCAGTTCCTGACGGGAGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGaaa                     >  4:279546/1‑147 (MQ=255)                                                                                                                               cGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACg                  <  3:719655/151‑1 (MQ=255)                                                                                                                                gTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTgggtaa                    >  2:657672/1‑144 (MQ=255)                                                                                                                                      cGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGc           <  3:436635/151‑1 (MQ=255)                                                                                                                                       ggTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGGGCTGGGGAACGGGaaa             >  3:953786/1‑148 (MQ=255)                                                                                                                                        gTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGGCAACGGTAAAGctg         <  3:719384/151‑3 (MQ=255)                                                                                                                                        gTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCgg         <  3:387161/151‑1 (MQ=255)                                                                                                                                            aGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGa     <  4:197246/151‑1 (MQ=255)                                                                                                                                             gTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGaa    <  4:107886/151‑1 (MQ=255)                                                                                                                                               tCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAg   <  2:679823/150‑1 (MQ=255)                                                                                                                                                 cTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCCGGACAAAGGGAAAGGGGccagagt  >  3:100397/1‑143 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCA‑CGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGT  >  NC_000913/4618871‑4619163

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