Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,489	C→T	T87T (ACC→ACT)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,489	0	C	T	93.3%	164.4 / NA	61	T87T (ACC→ACT)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (28/28);  minor base G (1/3);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGCG  >  NC_000913/4618246‑4618735                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        aaTTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATGGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCg                                                                                                                                                                                                                                                  <  2:361341/251‑1 (MQ=255)  attattCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGt                                                                                                                                                                                                                                                 <  1:581628/251‑1 (MQ=255)     attCGTAAAAAACGGGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCg                                                                                                                                                                                                                                               <  2:164723/250‑1 (MQ=255)      ttCGTAACAAAAGTGAAGGTCATGCACTACGCGATGAAGAAATTCGTTTATTTATCAACGGCATTCGAGCCAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCATGTCGCGCTGACCATGGCGATACGAGATTCAGGAACGGTTCTAGACTTGAAAAGCATGCATGTGAATGGCACGCATGGTGATAAATATTCCACTGGTGGCGTCGGc                                                                                                                                                                                                                                             >  2:72791/1‑251 (MQ=255)      ttCGTAAAAAACGTGATGGTCATGCCCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGAGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGGATGGCCCGCTTGTTGATAAACATTCCACTGGTGGCGTCGGc                                                                                                                                                                                                                                             >  1:492518/1‑251 (MQ=255)      ttCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGc                                                                                                                                                                                                                                             >  1:149719/1‑251 (MQ=255)       tCGTAAAAAACGGGATGGGCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGc                                                                                                                                                                                                                                             <  2:667498/250‑1 (MQ=255)            aaaaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGAtgtg                                                                                                                                                                                                                                       >  1:173055/1‑251 (MQ=255)              aaaCGTTATGGTCATGCACTGAGCGATGAAGAAATTCGGTTCTTTATCAACGGCATTCGCGACAACTCTATCTCCGAAGGGCAGATTGCCGCCCTCGCGAGGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGaa                                                                                                                                                                                                                                     <  2:466643/251‑2 (MQ=255)                  gTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCAATGGTGGCGTCGGCGATGTGACTTCg                                                                                                                                                                                                                                 >  2:567911/1‑251 (MQ=255)                  gTGATGATCATGCACTGAGCGATGAAGCAATTCGTTTCTTTATCAACGGCATTCGAGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCCTTTTCTTCCACGATATGCCAATGCCTGAGCGTGTCTCGCTGCCCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGCTTGTTGATAAACATTCCACTGGGGGCGTCGTCGATGTGACTTCg                                                                                                                                                                                                                                 >  2:98009/1‑251 (MQ=255)                                  gagagATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACCCTAGCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGcc                                                                                                                                                                                                                   <  2:192080/245‑1 (MQ=255)                                         gaagaaATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCAGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGCTGTGCCTTCGCTGCTGTTGGGGCCGATGGTCGc                                                                                                                                                                                                          >  2:128305/1‑251 (MQ=255)                                           agaaATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGc                                                                                                                                                                                                          >  1:190156/1‑249 (MQ=255)                                                 tCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGCTGGTCGCAGCCTgcg                                                                                                                                                                                                  >  2:336345/1‑251 (MQ=255)                                                             cAACGGCATACGCGACAACACTATCTCCGAAGGGCCGAGTGCCGCCCTCGCGATGACCATTTTCTACCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATGCAGGAACCGTGCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGGCGGCGATGTGACTGCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATAtt                                                                                                                                                                                       <  2:158222/250‑1 (MQ=255)                                                                 ggcaTTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCg                                                                                                                                                                                    >  1:202291/4‑249 (MQ=255)                                                                         cgACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTg                                                                                                                                                                          <  1:534487/251‑1 (MQ=255)                                                                          gACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGc                                                                                                                                                                                                                                             <  1:147884/183‑1 (MQ=255)                                                                          gACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTAGACTGGAAAAGCCTGCATCTGAATGGCCCGCTTGTTGATAAACATTCCACTGGTGGCGTCGGc                                                                                                                                                                                                                                             >  2:147884/1‑183 (MQ=255)                                                                            caacaCTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTTAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTc                                                                                                                                                                       <  1:532936/251‑1 (MQ=255)                                                                            caacaCTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTc                                                                                                                                                                       <  1:546175/251‑1 (MQ=255)                                                                                 cTATCTCCGAAGGGCAGATTGCCGCCTTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGc                                                                                                                                                                  >  1:294152/1‑251 (MQ=255)                                                                                       ccGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGACAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCt                                                                                                                                                                <  1:547070/247‑1 (MQ=255)                                                                                         gAAGGGCGGAGTGCCGCACTCGCGATGACAATTTTCTTCCACGATATGCCAATGCCTGAGCGGGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATAGCCCGATTGTTGATAAACATTCCACGCGAGGCGTCGGCGCTGTGATTTCGCTGATGTTTAGGCCGCTGGTCGCAACCTGAGGCGGCTATAGTCCgatattatctgggcgaggtctagggcc                                                                                                                                                          >  2:379108/1‑227 (MQ=255)                                                                                           aGGGCATATTTCCTCCCTCGCGATGTCCTTTTTCTTCCACTATATGCCAATGCCTGAGCTTGTCTCGCTGACCATGTCGATGCGAGATTCAGGAACCGTTCTCTACTGGAAAAGCCTGCATCTGAATGGCCCTTTTGTTGATAAACATTCCACTGGTGTCGTCGGCGATGTGACTTCGCTGATGTTGGGTCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATa                                                                                                                                                        <  1:452133/251‑1 (MQ=255)                                                                                                aGATTGCCGCCCTCGCGATTACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTTCGGTGATCTCTGGTCGGGGACTCgtgcagactggg                                                                                                                                                   >  2:129507/1‑240 (MQ=255)                                                                                                aGATTGCCGCCATCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGCACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGCTCAACATTCCACTGGCGGCGTCGGCGATGTGCCTTCGATGATGTTGGGGCCGCTGGTCGAAGCCTGCGGCGGGCATATTCCGATGATCTCTGGGCGAGGCCTCGGGCATACTggg                                                                                                                                                   >  2:385968/1‑250 (MQ=255)                                                                                                                 aTGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTAAGCTCgcccaact                                                                                                                                  >  1:433308/1‑244 (MQ=255)                                                                                                                   gACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTTCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCggcacaccgg                                                                                                                                >  1:196244/1‑242 (MQ=255)                                                                                                                      cATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAat                                                                                                                             >  1:390100/1‑251 (MQ=255)                                                                                                                               ccACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑Tg                                                                                                                   <  1:98009/251‑1 (MQ=255)                                                                                                                                   gATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCggcggc                                                                                                                                                                                             >  2:442736/1‑174 (MQ=255)                                                                                                                                   gATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCtt                                                                                                               <  2:184547/251‑1 (MQ=255)                                                                                                                                      atGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCg                                                                                                             <  1:673803/250‑1 (MQ=255)                                                                                                                                        gacaaggctggAGCGTGTCTCGCTTACCATGGCGATGCGAGATTCAGGAACCGTTCTCGTCTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACTTTCCACGGGTGGCGTCGGCGATGTGACTGCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACa                                                                                                          <  2:335582/241‑1 (MQ=255)                                                                                                                                         aCAATGCCCGAGCGTGTCTCGCTGCCCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGACTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATCCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACAt                                                                                                         <  1:713/251‑1 (MQ=255)                                                                                                                                                  gAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTCTATTCCGATGTTCTCTGGGCGCGGCCTCGGTCATACTGGCGGTACGCCCGACAAACTGGAATCACTCCC‑CGGCTTCGCAATTTTTCCggg                                                                                                >  2:411470/1‑250 (MQ=255)                                                                                                                                                    gcggttctcGCTGACCATGGCGATGCGGGATTCAGGAACCGTTCTCGACTGTAAAAGCCTGCCTCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGGCGGCGAGGGGACTTCTCTGATGTGGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATg                                                                                              <  2:166121/246‑1 (MQ=255)                                                                                                                                                        gtcgcgctgcccagtgCGATGTGAGAGTAAGGAACCGTTCTCGACTGGAAAAGCCTGTATCTGAATGGCCTGATTGTTGATAATCATTCCACTTGTGGCTGCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCTGCTATATTCCGATGATCTCTGGTCGCGGCCGCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGCCaa                                                                                          <  2:463604/236‑1 (MQ=255)                                                                                                                                                                gACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGACAACCGTTTcc                                                                                  <  1:128305/251‑1 (MQ=255)                                                                                                                                                                 aCCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGACAAcc                                                                                        <  1:290838/244‑1 (MQ=255)                                                                                                                                                                    aTGGCGATGCGAGCTTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGTATGGCCCGATTGTTGCTAAACATTCCACTGGTGGCGTCGGCGAAGTGACTTCGCTGATGTTGGGGCCGCTGGTCGCAGCCTGCGGCGGCTATCTTCCGTTGCTCTCTTGTCGCGGCCGCGGTCATACTTGAGGAACGGCCGACAAATCGGAATCCGTCCCGCGGCTGCG‑CATTTTCCCGGATTACACCCGCTTCcgcgg                                                                              >  1:3172/1‑250 (MQ=255)                                                                                                                                                                    aTGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGCCAAACTGGAATCCCACCC‑TGGCTTCGACATTTTCCCGGGTGGCAACCGGTTccccgg                                                                              >  1:462691/1‑247 (MQ=255)                                                                                                                                                                             cGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCCCTACAAACTGGAATCCATCCC‑TGGCTTCGCCATTTTCCCGGATGACAACCGTTTCCCCGAAAttttta                                                                      >  1:37343/1‑250 (MQ=255)                                                                                                                                                                                  ttCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCCTCCC‑TGGCTTCGGCATTTTCCCGGGTGACAAACGTTTCCGCGCAATTATTAAAAACg                                                                >  1:434529/1‑251 (MQ=255)                                                                                                                                                                                   tCAGTAACCGGTCCCGACTGGAAAGGCCTGCATTGGAATGTCCCGATTGTTGATAAACCTTCCCCTGGGGTCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCTGCCTGCGGCGGCTATATTCCGATTATCTCTGGTCGCGGCCTCGGTCATACTGGCTGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGTCAACCGTTTCCGCGAAATTATTAAAGACg                                                                <  2:173055/250‑1 (MQ=255)                                                                                                                                                                                   tCAGGAACCGTTCCCGACTGGAAAGCCCTGCATCTGAATGGCCCGATGTTTGATAACCATTCCACGGGTGGCGGCGGCGATGTGACTTCGCGGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATGCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACg                                                                <  2:178108/250‑1 (MQ=255)                                                                                                                                                                                     aGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCg                                                             <  1:55872/251‑1 (MQ=255)                                                                                                                                                                                         aCCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGCCCTTTTCCCGGATGGCAAACGTTTCCGCGGAATTATTAAAGACGTCGGCgg                                                         >  1:255932/1‑250 (MQ=255)                                                                                                                                                                                            gTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGGTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGc                                                      <  1:518031/251‑1 (MQ=255)                                                                                                                                                                                              tctcGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATGCTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCg                                                     <  1:567911/250‑1 (MQ=255)                                                                                                                                                                                              tctcGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCCCGGCAAAATGGAATCCATCCC‑TGGCTTCGCCATTTTCCCGGATGGCAACCGTTTCCGCGAAATTTTTTAAGAAGTAGGCGTGGCg                                                     >  1:362594/1‑250 (MQ=255)                                                                                                                                                                                                 cGACTGTAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCATCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGAt                                                   <  1:345048/249‑1 (MQ=255)                                                                                                                                                                                                        aaaaGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTAAGCTCGACAAACTGGAAACCATCCC‑TGGCCTCGACATTTTCCCGGGTGTCAAACGTTTTCGCGAAATTTTTAAAATCGAAGGCCTGGCGATTATTGGTc                                           >  1:508597/1‑250 (MQ=255)                                                                                                                                                                                                         aaaGCCTGCATCTGAAAGGCCCGATTGTTGATAAACATTCCACGGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCGCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAg                                         <  2:554177/251‑1 (MQ=255)                                                                                                                                                                                                                     ggtaTGGCCCGAGTGGTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACAGTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTg                             <  2:346321/248‑1 (MQ=255)                                                                                                                                                                                                                           gCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACCTTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGGCAGACCAGCCCACTGGCTccc                       >  1:434747/1‑250 (MQ=255)                                                                                                                                                                                                                                   gttgATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACGGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGTCATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATaa               <  1:465840/251‑1 (MQ=255)                                                                                                                                                                                                                                      gATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACg            <  1:336345/251‑1 (MQ=255)                                                                                                                                                                                                                                                ccACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGACAACCGGTTCCGCGGAAATATTAAAAACGTCGGCGTGGCGGTTTTCGGTCAGAACAACTCCCTGGCTCCGGCGGATAAAAGttttttagcg  >  1:256210/1‑244 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        AATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCC‑TGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGCG  >  NC_000913/4618246‑4618735

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