Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,393	C→T	L389L (CTG→TTG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,393	0	C	T	98.2%	159.7 / NA	54	L389L (CTG→TTG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/36);  minor base G (1/0);  total (18/36)
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

AGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACG‑CCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTC  >  NC_000913/4619152‑4619634                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     aGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑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‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCa                                                                                                                                                                                                          <  1:411152/251‑1 (MQ=255)                                    aaaaaGGCCGAACCGGCTTCGTTGAGAGCTCCGCGAAGTAGCTGCCGACGGCTATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGAGCTGGGGAGGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCAc                                                                                                                                                                                                         <  2:133928/251‑1 (MQ=255)                                          gCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGaaa                                                                                                                                                                                                   <  1:406726/251‑1 (MQ=255)                                             caaccgacTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTT‑TCACGCGaaaaac                                                                                                                                                                                                >  1:212836/3‑248 (MQ=255)                                              aaccgacTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACg                                                                                                                                                                                               <  1:320613/250‑1 (MQ=255)                                                  gacgtcggtGAGAACTACGTGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑GCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGaaaa                                                                                                                                                                                           <  2:250084/243‑1 (MQ=255)                                                    cTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGCTACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAA‑CCACGCCAAAGGCGaaaaa                                                                                                                                                                                          >  1:649014/1‑249 (MQ=255)                                                    cTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTAT‑CCACGCGAAAGACGgaaaa                                                                                                                                                                                          >  1:573779/1‑249 (MQ=255)                                                       cGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAaacagct                                                                                                                                                                                      <  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‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTg                                                                                                                                                                                     <  1:338281/251‑1 (MQ=255)                                                             gacagacgagAAGTATCTGCCGACAGGGATGCTGAGTAAAGCAGATC‑AGTCTGATACCGAAGGGTTTGTCAGTGGAATGGATTTCCGCGGGCTGGGGAGGGCTGTGGTTGCCAGGGGCGGCGGTCG‑TCGTCAGGCAGCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCGGAGTCCGCTGGCGGTTAT‑CCACGCGAAAGCCGAAAACAGCTGGCAgg                                                                                                                                                                                <  2:157354/242‑1 (MQ=255)                                                               aCTACGCGAAGTATCTGCCGGCAGCGATGCTGACGAAAGCAG‑TCTATGCTGATTCCGAAGGTTTTGTCAGGGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGaa                                                                                                                                                                              <  2:622657/251‑1 (MQ=255)                                                                       aaGTATCTGCCGATAGCGATGCGGACGAAAGCCG‑TCTATGCAGAGACCGAAGGTTGTGTCGGTGAAATGGATCCCCTCGCGCTTGGGATGGCAGTGGTTGCAATGGGCGGCGGGCG‑TCTTCAGGCATCTGACTGCATTGATTACAGCGTCCGCTTTACTGATAGGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCTCTGGCGGTTAT‑CCACGCGAAAGACGAAAACATCTGGAAGTAAGCGGCGaa                                                                                                                                                                      <  2:57357/251‑1 (MQ=255)                                                                           aTCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑TCCCGCGAAAAACGGAAAAAGCGGGCAGGAAGCGGGGAAAGCg                                                                                                                                                                  >  1:678698/1‑251 (MQ=255)                                                                              tGCCGACAGCGATGCTGACGAAAGCAG‑TTTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGt                                                                                                                                                                <  2:485386/250‑1 (MQ=255)                                                                               tccGACAGCGATGCTGACTAAAGCAG‑GCTATGCTGATACCGAAGGTGTTGTCAGTGAAAGGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAAGGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGTCAGGAAGCGGCGAAAGCGGt                                                                                                                                                                <  2:189242/248‑1 (MQ=255)                                                                                   aCAGCGATGCTGACGAAAGCCG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGAAGCCGGTCAGCGTCCGCTGGCGGGTTT‑CCACGCGAAAGACGGAAACAGCCGGCCGGAAGCAGCGGAAGCGGAGaaga                                                                                                                                                           >  2:39956/1‑248 (MQ=255)                                                                                                 gAAAGCAG‑TCTATGCTTATACCGAAGGGTTTGTCAGTGAAATGGATCCCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTTAGCGTCCGCTGGCGTTTAT‑CCACGCGAAAGACTAAAACAGCTGTCAGGAAgcgg                                                                                                                                                                          <  2:571573/221‑1 (MQ=255)                                                                                                    aGCAG‑TCTATGCTGATCCCGAATGGTTTGTCAGTGAAATGGGTACCCGCGCGCTGTGGATGGCAGGGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGGCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCg                                                                                                                                        <  2:252251/251‑1 (MQ=255)                                                                                                                              tttGTCAGTGAAATGGATACCCGCGCGCGGGGGATGGCAGTGGTTGCAATGGGCGGCTGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                                                <  2:48413/250‑1 (MQ=255)                                                                                                                                tgtacgTGAAATGGAGACCCGCGCGCTGGGGATGGTAGTGGTTGCAATGGGCGGCTGTTGTATGT‑AGGCATCTGACACCATTGATTACAGCGGCGGCGTTACTGATTTGGCGCGTTTGGGCGACCAGTTAGACGGGCAGCGTCCGATGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑TAATTAAACTTGCCGATAAAGCACCTGAAAGCACATCaaat                                                                                                              <  2:473201/245‑3 (MQ=255)                                                                                                                                   cagggAAATGGATCCCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCTGTCT‑TCGTCAGGCATCTGAAACCATTGATTACAGCGGCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTTCCGATAAAGCACCGGAAAGCACCCCAACTGTCt                                                                                                          <  2:649014/247‑1 (MQ=255)                                                                                                                                                       cgcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                                      <  1:530687/251‑1 (MQ=255)                                                                                                                                                        gcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAAAAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAGCGG‑CCATTAAACTTGCCGATAAAGCACCGGGAAACACACCAACAGTCTATTGCCGTATCAGCGAATaa                                                                                      >  1:338733/1‑250 (MQ=255)                                                                                                                                                              gaaggCATTGGTTGCAATGGGCGGTGCTCT‑TCGTCAGGCATCAGACCCCATTGATTACAGCGTCGGCTTTACGGATATTGCGCGTTTGGGCGACCAGGTAGACGGTTAGCGTCCGCTGGCGGTTAT‑CCACTCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAATCACACCAATTGTCTATCGCCGTATCAGCGAATAACGGtata                                                                               <  2:117209/248‑1 (MQ=255)                                                                                                                                                                 ggCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActg                                                                            <  1:137083/251‑1 (MQ=255)                                                                                                                                                                     gTGGTTGCAATGGTCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCTAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCATCTGTCTATCGCCGTATCAGCGAATAACGGTATActgatct                                                                        <  2:573779/251‑1 (MQ=255)                                                                                                                                                                       ggTTACAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGTGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCTCCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatctga                                                                      <  1:70123/251‑1 (MQ=255)                                                                                                                                                                        ggtGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGGAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatctgat                                                                     <  1:33388/249‑1 (MQ=255)                                                                                                                                                                          tGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACtgatctgatc                                                                    <  1:190777/250‑1 (MQ=255)                                                                                                                                                                            cAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTAT‑CCCCGCGAAAGGCGAAAACAGCTGGCCGGAAGCGGCGGAAGCAGTGAAAGCGGGCAAT‑AAAAGTGCCGGTCAAGCACCGGGAAGCACACCCACTGTCTATCCCCGGATCAGCGAATAACGGGAAACTGATATGATCAtt                                                                 >  1:654888/1‑251 (MQ=255)                                                                                                                                                                                gggCGGCGGTTG‑TCGTCAGGCATCTGACACAATTGAGTACAGCGTCGGTATTACTGGTATGGCTCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTTATCTGATCATTTaa                                                              <  2:665587/250‑1 (MQ=255)                                                                                                                                                                                      gcggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGa                                                        <  1:101410/243‑1 (MQ=255)                                                                                                                                                                                        ggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGGATCAGCGAATAACGGGATACTGATCTGATCATTTAAATTTGa                                                        >  1:417007/7‑249 (MQ=255)                                                                                                                                                                                                    gCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCacac                                                                                                                   <  2:482602/178‑1 (MQ=255)                                                                                                                                                                                                    gCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAACGG‑CAATTAAAATTGCCGATAAAGCACCGGGAAGCacac                                                                                                                   >  1:482602/1‑178 (MQ=255)                                                                                                                                                                                                          gACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACa                                   <  2:289460/251‑1 (MQ=255)                                                                                                                                                                                                           acacCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAAc                                                                                     <  1:102069/201‑1 (MQ=255)                                                                                                                                                                                                              ccATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATg                               <  1:24847/251‑1 (MQ=255)                                                                                                                                                                                                                  tgATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCCGGGAGCGGCGAAAGCGGGGAAAGGGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGAACACCAACTGTCTATACCCGAATCAGCGAATAACGGGATACTGATCTGATCATTTTAAATTGAAACACTGAATACGGAGA‑ACATATggaa                            >  1:42788/2‑247 (MQ=255)                                                                                                                                                                                                                   gATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACg                          <  1:122481/251‑1 (MQ=255)                                                                                                                                                                                                                       aCAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAGCGG‑CAATTAAAATTGCCGATAAAGCACCGGAAAGCACACCCACTGTCTATCGCCGTATCCGCGAATAACGGGATACTGATCTGATCATTTAAATTTTAAGCACTGGGGACGGAGA‑AAATATGAAACGTGc                       >  1:487037/1‑250 (MQ=255)                                                                                                                                                                                                                          gCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCAttt                   <  1:39956/251‑1 (MQ=255)                                                                                                                                                                                                                                           atGGCGCGCTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTCATACCA‑GCGAAAGACGAAAACAGCACGGAGGCAGCGAAGAAAGCGGTGAAAACGG‑CAATTAAAAATGCCGCTAAAAGACCGGAGAGAAACCCAACTGGCTATACCCGAATCAGCGAATCAAGGGATACCGTTATGATAATTTAACTTTTGAGCACAGAGTACGGAGACACA‑ATGAAACGTGCATttattttagagcttgactc  >  1:46902/1‑237 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     AGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACG‑CCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTC  >  NC_000913/4619152‑4619634

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