Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,359	C→T	I377I (ATC→ATT)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,359	0	C	T	98.4%	186.6 / NA	63	I377I (ATC→ATT)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (20/42);  total (21/42)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACG‑CCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATG  >  NC_000913/4619119‑4619605                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        cgcgAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCg                                                                                                                                                                                                                                               <  1:6286/251‑1 (MQ=255)    gAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCg                                                                                                                                                                                                                                            >  1:84922/1‑251 (MQ=255)     aaaTTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCgg                                                                                                                                                                                                                                           <  1:74842/251‑1 (MQ=255)      aaTTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGc                                                                                                                                                                                                                                          <  1:96923/251‑1 (MQ=255)            aagcagTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCtttcct                                                                                                                                                                                                                                    >  2:190777/6‑248 (MQ=255)                tgTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAACAAGACCCGACCGACTTCGTTGAGCACTACGCGAAGGATCTGCAGACAGCGATGCTGACGAAAGAAG‑TCTATGCTAATACCGAAGGGTTTGTCAGTGAAATGGATACCCGCGAGCTGGGGCAGGCAGTGGTTGCAATGGACGTCGGGCG‑CAGTCAGCCATCTCACACCAGCGGTGACATCGTCGGCTTTACTg                                                                                                                                                                                                                                   >  2:14032/2‑248 (MQ=255)                agTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGata                                                                                                                                                                                                                                <  1:415265/250‑1 (MQ=255)                     tGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGgcg                                                                                                                                                                                                                           >  1:484570/1‑251 (MQ=255)                      ggACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGgcg                                                                                                                                                                                                                           <  2:135213/250‑1 (MQ=255)                      ggACAACGGGAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGgc                                                                                                                                                                                                                            <  1:404594/249‑1 (MQ=255)                         aaaCGGTAAAGCGGCAGAAGTCT‑TTGGTCGTTTGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGtt                                                                                                                                                                                                                       <  2:592256/250‑2 (MQ=255)                           aCGGTAAGGCGGCAGAAGTCTGTTGG‑CGTATGGTAGCGGCACAAAAAGGCCCAGCCGACTTCGTTGAGAACTACGCGAAGTATCTGTCGACAGCGATGCTGACGAAAGCAG‑GCGATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGGCACCATTGATTACAGCGTCGGTTTTACTGATATGGGGCGtttg                                                                                                                                                                                                                     <  2:84922/251‑4 (MQ=255)                           aCGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGtttg                                                                                                                                                                                                                     >  1:485386/1‑248 (MQ=255)                                 aaGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGGGGTTGCAATGGGCGGCGGACG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACGGATATGGCGCGttttggagg                                                                                                                                                                                                                >  2:137083/1‑242 (MQ=255)                                  aGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGAcc                                                                                                                                                                                                              >  1:675663/1‑251 (MQ=255)                                                    gTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCgg                                                                                                                                                                                             >  1:651206/1‑250 (MQ=255)                                                      cGTATGGTAGCGGCACAAAAAGGCCCAACCGTGTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCAGATACCGAAGGTTTTGGCAGTGAAATGGATACCCGCGCGCTGGGGATGTCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTc                                                                                                                                                                                            <  2:416079/250‑1 (MQ=255)                                                         aTGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAAACCCGCGCGCTGGGGATGGCAGGGGTTGCAATGGGCGGCGGTCG‑TCGTCCGGCATCTGACACCATTGTTTAACGAGTCGGCTTTACTGTTATGGCGCGTGTGGGCGTCCAGGAAGAAGGTCAGCgccagct                                                                                                                                                                                        >  2:487495/1‑245 (MQ=255)                                                           ggTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGAAACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTgg                                                                                                                                                                                      >  1:250084/1‑251 (MQ=255)                                                            gTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGGTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGc                                                                                                                                                                                     <  2:480116/251‑1 (MQ=255)                                                                 ggCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCgggtt                                                                                                                                                                                >  1:289460/1‑248 (MQ=255)                                                                  gCACAAAAAGGCCCAACCGACTTCGATTAGAACTACGCGAAGTATCTGCCGACAGCTATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGAGATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAt                                                                                                                                                                               <  2:330476/251‑1 (MQ=255)                                                                     caAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCa                                                                                                                                                                            <  1:411152/251‑1 (MQ=255)                                                                      aaaaaGGCCGAACCGGCTTCGTTGAGAGCTCCGCGAAGTAGCTGCCGACGGCTATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGAGCTGGGGAGGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCAc                                                                                                                                                                           <  2:133928/251‑1 (MQ=255)                                                                            gCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGaaa                                                                                                                                                                     <  1:406726/251‑1 (MQ=255)                                                                               caaccgacTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTTCACGCGaaaaac                                                                                                                                                                  >  1:212836/3‑248 (MQ=255)                                                                                aaccgacTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACg                                                                                                                                                                 <  1:320613/250‑1 (MQ=255)                                                                                    gacgtcggtGAGAACTACGTGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑GCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGaaaa                                                                                                                                                             <  2:250084/243‑1 (MQ=255)                                                                                      cTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGCTACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAACCACGCCAAAGGCGaaaaa                                                                                                                                                            >  1:649014/1‑249 (MQ=255)                                                                                      cTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTATCCACGCGAAAGACGgaaaa                                                                                                                                                            >  1:573779/1‑249 (MQ=255)                                                                                         cGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAaacagct                                                                                                                                                        <  1:596344/251‑4 (MQ=255)                                                                                          gTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACg                                                                                                                                                                                                     >  2:275533/1‑205 (MQ=255)                                                                                          gTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACg                                                                                                                                                                                                     <  1:275533/205‑1 (MQ=255)                                                                                          gTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTg                                                                                                                                                       <  1:338281/251‑1 (MQ=255)                                                                                               gacagacgagAAGTATCTGCCGACAGGGATGCTGAGTAAAGCAGATC‑AGTCTGATACCGAAGGGTTTGTCAGTGGAATGGATTTCCGCGGGCTGGGGAGGGCTGTGGTTGCCAGGGGCGGCGGTCG‑TCGTCAGGCAGCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCGGAGTCCGCTGGCGGTTATCCACGCGAAAGCCGAAAACAGCTGGCAgg                                                                                                                                                  <  2:157354/242‑1 (MQ=255)                                                                                                 aCTACGCGAAGTATCTGCCGGCAGCGATGCTGACGAAAGCAG‑TCTATGCTGATTCCGAAGGTTTTGTCAGGGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGaa                                                                                                                                                <  2:622657/251‑1 (MQ=255)                                                                                                         aaGTATCTGCCGATAGCGATGCGGACGAAAGCCG‑TCTATGCAGAGACCGAAGGTTGTGTCGGTGAAATGGATCCCCTCGCGCTTGGGATGGCAGTGGTTGCAATGGGCGGCGGGCG‑TCTTCAGGCATCTGACTGCATTGATTACAGCGTCCGCTTTACTGATAGGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCTCTGGCGGTTATCCACGCGAAAGACGAAAACATCTGGAAGTAAGCGGCGaa                                                                                                                                        <  2:57357/251‑1 (MQ=255)                                                                                                             aTCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATTCCCGCGAAAAACGGAAAAAGCGGGCAGGAAGCGGGGAAAGCg                                                                                                                                    >  1:678698/1‑251 (MQ=255)                                                                                                                tGCCGACAGCGATGCTGACGAAAGCAG‑TTTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGt                                                                                                                                  <  2:485386/250‑1 (MQ=255)                                                                                                                 tccGACAGCGATGCTGACTAAAGCAG‑GCTATGCTGATACCGAAGGTGTTGTCAGTGAAAGGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAAGGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGTCAGGAAGCGGCGAAAGCGGt                                                                                                                                  <  2:189242/248‑1 (MQ=255)                                                                                                                     aCAGCGATGCTGACGAAAGCCG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGAAGCCGGTCAGCGTCCGCTGGCGGGTTTCCACGCGAAAGACGGAAACAGCCGGCCGGAAGCAGCGGAAGCGGAGaaga                                                                                                                             >  2:39956/1‑248 (MQ=255)                                                                                                                                   gAAAGCAG‑TCTATGCTTATACCGAAGGGTTTGTCAGTGAAATGGATCCCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTTAGCGTCCGCTGGCGTTTATCCACGCGAAAGACTAAAACAGCTGTCAGGAAgcgg                                                                                                                                            <  2:571573/221‑1 (MQ=255)                                                                                                                                      aGCAG‑TCTATGCTGATCCCGAATGGTTTGTCAGTGAAATGGGTACCCGCGCGCTGTGGATGGCAGGGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGGCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCg                                                                                                          <  2:252251/251‑1 (MQ=255)                                                                                                                                                                tttGTCAGTGAAATGGATACCCGCGCGCGGGGGATGGCAGTGGTTGCAATGGGCGGCTGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                  <  2:48413/250‑1 (MQ=255)                                                                                                                                                                  tgtacgTGAAATGGAGACCCGCGCGCTGGGGATGGTAGTGGTTGCAATGGGCGGCTGTTGTATGT‑AGGCATCTGACACCATTGATTACAGCGGCGGCGTTACTGATTTGGCGCGTTTGGGCGACCAGTTAGACGGGCAGCGTCCGATGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑TAATTAAACTTGCCGATAAAGCACCTGAAAGCACATCaaat                                                                                <  2:473201/245‑3 (MQ=255)                                                                                                                                                                     cagggAAATGGATCCCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCTGTCT‑TCGTCAGGCATCTGAAACCATTGATTACAGCGGCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTTCCGATAAAGCACCGGAAAGCACCCCAACTGTCt                                                                            <  2:649014/247‑1 (MQ=255)                                                                                                                                                                                         cgcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                        <  1:530687/251‑1 (MQ=255)                                                                                                                                                                                          gcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAAAAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAGCGG‑CCATTAAACTTGCCGATAAAGCACCGGGAAACACACCAACAGTCTATTGCCGTATCAGCGAATaa                                                        >  1:338733/1‑250 (MQ=255)                                                                                                                                                                                                gaaggCATTGGTTGCAATGGGCGGTGCTCT‑TCGTCAGGCATCAGACCCCATTGATTACAGCGTCGGCTTTACGGATATTGCGCGTTTGGGCGACCAGGTAGACGGTTAGCGTCCGCTGGCGGTTATCCACTCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAATCACACCAATTGTCTATCGCCGTATCAGCGAATAACGGtata                                                 <  2:117209/248‑1 (MQ=255)                                                                                                                                                                                                   ggCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActg                                              <  1:137083/251‑1 (MQ=255)                                                                                                                                                                                                       gTGGTTGCAATGGTCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCTAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCATCTGTCTATCGCCGTATCAGCGAATAACGGTATActgatct                                          <  2:573779/251‑1 (MQ=255)                                                                                                                                                                                                         ggTTACAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGTGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCTCCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatctga                                        <  1:70123/251‑1 (MQ=255)                                                                                                                                                                                                          ggtGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGGAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatctgat                                       <  1:33388/249‑1 (MQ=255)                                                                                                                                                                                                            tGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACtgatctgatc                                      <  1:190777/250‑1 (MQ=255)                                                                                                                                                                                                              cAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTATCCCCGCGAAAGGCGAAAACAGCTGGCCGGAAGCGGCGGAAGCAGTGAAAGCGGGCAAT‑AAAAGTGCCGGTCAAGCACCGGGAAGCACACCCACTGTCTATCCCCGGATCAGCGAATAACGGGAAACTGATATGATCAtt                                   >  1:654888/1‑251 (MQ=255)                                                                                                                                                                                                                  gggCGGCGGTTG‑TCGTCAGGCATCTGACACAATTGAGTACAGCGTCGGTATTACTGGTATGGCTCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTTATCTGATCATTTaa                                <  2:665587/250‑1 (MQ=255)                                                                                                                                                                                                                        gcggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGa                          <  1:101410/243‑1 (MQ=255)                                                                                                                                                                                                                          ggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGGATCAGCGAATAACGGGATACTGATCTGATCATTTAAATTTGa                          >  1:417007/7‑249 (MQ=255)                                                                                                                                                                                                                                      gCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCacac                                                                                     <  2:482602/178‑1 (MQ=255)                                                                                                                                                                                                                                      gCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAACGG‑CAATTAAAATTGCCGATAAAGCACCGGGAAGCacac                                                                                     >  1:482602/1‑178 (MQ=255)                                                                                                                                                                                                                                            gACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACa      <  2:289460/251‑1 (MQ=255)                                                                                                                                                                                                                                             acacCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAAc                                                       <  1:102069/201‑1 (MQ=255)                                                                                                                                                                                                                                                ccATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATg  <  1:24847/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        CGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACG‑CCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATG  >  NC_000913/4619119‑4619605

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