Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,338	C→T	R370R (CGC→CGT)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,338	0	C	T	98.3%	176.0 / ‑4.9	61	R370R (CGC→CGT)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/40);  total (20/41)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACG‑CCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAA  >  NC_000913/4619095‑4619575                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               ggCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCTCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAgg                                                                                                                                                                                                                                         <  1:521550/251‑1 (MQ=255)    gAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCAt                                                                                                                                                                                                                                      <  1:300743/251‑1 (MQ=255)          ggACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGaca                                                                                                                                                                                                                                <  2:176690/250‑1 (MQ=255)               ccGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCAtt                                                                                                                                                                                                                           <  2:171774/251‑2 (MQ=255)               ccGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGACC‑ATGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGCGAACTACACGGCGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTAAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGCCG‑TCGTCAGGCATCTGACACCAtt                                                                                                                                                                                                                           >  1:313234/1‑250 (MQ=255)                         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                                                                                                                                                                       <  1:275533/205‑1 (MQ=255)                                                                                                                  gTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACg                                                                                                                                                                       >  2:275533/1‑205 (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)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               GGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACG‑CCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAA  >  NC_000913/4619095‑4619575

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