Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,389	G→T	A387A (GCG→GCT)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,389	0	G	T	97.7%	120.2 / NA	45	A387A (GCG→GCT)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (10/33);  total (10/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTG  >  NC_000913/4619150‑4619629                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       aaaGCGGCAGAAGTCTTTGGTCGTATGGTAGCTGCACAAAAAGGCCCTACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGTTTGCAATTGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCg                                                                                                                                                                                                                                              <  1:403086/250‑1 (MQ=255) aaaGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGttcgggcgg                                                                                                                                                                                                                                             >  2:422996/1‑243 (MQ=255) aaaGCGGCAGAAGTCTTTGGACGTATGGTAGCGGCACACAAAGGACCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGAGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGcgcgacggggggg                                                                                                                                                                                                                                             >  2:164434/1‑242 (MQ=255)      ggCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTAGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGACAGCAGTCTATGCTGATACCGAAGGTTTTGACAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGGGACCAgg                                                                                                                                                                                                                                        >  1:69188/1‑251 (MQ=255)                      cGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTcc                                                                                                                                                                                                                        <  2:207836/251‑1 (MQ=255)                      cGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTcc                                                                                                                                                                                                                        <  2:120397/251‑1 (MQ=255)                       gTATGGTAGCGGCACAAAAAGGCCCGACCGAATTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTAAGTGAAATGGATACCCGCGCGCTGGGGATGGAAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACAACATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTAAGCGGCCg                                                                                                                                                                                                                       >  1:32504/1‑251 (MQ=255)                           ggTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCATTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACACGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTgg                                                                                                                                                                                                                   <  1:162461/251‑1 (MQ=255)                                 ggCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTa                                                                                                                                                                                                            <  2:221040/251‑1 (MQ=255)                                    acaAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTgg                                                                                                                                                                                                                                                 <  2:96676/212‑1 (MQ=255)                                    acaAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTgg                                                                                                                                                                                                                                                 >  1:96676/1‑212 (MQ=255)                                               cgaccgacTTCGTTGAGAACTACGCGAAGTACCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCCGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑Ac                                                                                                                                                                                            <  2:385158/251‑1 (MQ=255)                                                    gacTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATTCTGATACCGAAGGTTTTGTCAGAGAAATGGCTACCCGCGCGATGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGGTTCCAGCGTCGGCTTTACTGATATGGGTCGTCTGGGAGGCCAGGTAGACGGTCCGCGTCCGCTGGAGG‑TTaaaccagacgacgcagaaa                                                                                                                                                                                          >  2:246547/1‑232 (MQ=255)                                                    gacTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑aaaa                                                                                                                                                                                      <  1:12776/251‑1 (MQ=255)                                                             gagaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCa                                                                                                                                                                             <  1:44567/251‑1 (MQ=255)                                                             gagaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGGCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGAGTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGc                                                                                                                                                                              <  1:157030/250‑1 (MQ=255)                                                               gaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAgg                                                                                                                                                                           <  1:353504/251‑1 (MQ=255)                                                                     cgcgAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAAACGGGCAGCGGCCGCTGGCGG‑CTATTCAC‑CCCAAAGGAGG‑AAaacaaagggcgggagcgg                                                                                                                                                                      >  1:238292/1‑238 (MQ=255)                                                                      gcgAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTAGGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCGATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACGGGCAGGAAgcggc                                                                                                                                                                    <  2:25694/251‑1 (MQ=255)                                                                       cgAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAgcggcg                                                                                                                                                                   <  1:30313/251‑1 (MQ=255)                                                                                 gCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTg                                                                                                                                                          <  1:239992/250‑1 (MQ=255)                                                                                    gACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGTATGGCAGTGGTTGCAATTGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACTGCGTCGGCTTTACTGATATGTCTCGGCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAATCGGTGAAAg                                                                                                                                                      <  1:215622/251‑1 (MQ=255)                                                                                              cTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTTGGCGACCAGGGAGACGGTCAGCGTCCGCTGGCGGGTTTTCCCC‑CGGAAAA‑ACG‑AAAACAAC‑GGCAGGGAGCGGCGGAAGGGGGGGAAGCgggcaata                                                                                                                                             >  2:107416/1‑244 (MQ=255)                                                                                                              aTGCTGATACCGAAGGTTTGGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATAGGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAACCGTGCCGATAAAGCa                                                                                                                           <  2:96303/251‑1 (MQ=255)                                                                                                                       ccGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGc                                                                                                                  <  1:318168/251‑1 (MQ=255)                                                                                                                           aGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGTGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAT‑ACG‑AAAACAACTGGCATGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCacac                                                                                                              <  2:286406/251‑1 (MQ=255)                                                                                                                                             ggATCCCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCTGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGt                                                                                            <  1:164434/250‑1 (MQ=255)                                                                                                                                              gATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGa                                                                                                                      <  2:48924/224‑1 (MQ=255)                                                                                                                                                     gcgcgcTGGGGAGGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCCCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTGGCCGATAAAGCTCCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGaa                                                                                   <  2:434813/251‑1 (MQ=255)                                                                                                                                                      ccgatgtgggATGGCAGGGGGTGCAATGGTCGGCGGACGCTGTCAGGCATCGTACACCATCGAGTACGGCGTCGGCCTTACGGATTTGGCTCGGCGGGGAGTCCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCTAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAAGTAAACTTGCCGAGAAAGCACCGGAAAGCACCCC‑AACCGTCTATCGCCGTATCTGCGAAt                                                                                  <  2:403254/244‑1 (MQ=255)                                                                                                                                                       gcgcTGGTGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATa                                                                                 <  2:100743/251‑1 (MQ=255)                                                                                                                                                         gcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAAc                                                                               <  1:422996/251‑1 (MQ=255)                                                                                                                                                            ggggATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACgg                                                                             <  1:315810/250‑1 (MQ=255)                                                                                                                                                             gggATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCCTCGATTACAGCGTCGGCTTTACTGAAATGGCTCGTCTGGGCGACCAGGCAGACGGTCAGCGTCCGCAGGCGG‑TTATCCCCAGGGAAAA‑AGG‑AAACCAA‑TGGGCGGCAGCGGAGAAAGCGGGGGAAGAGGACAAT‑AAAAATGCCGATAAAGCCCCGGCAAAAACACC‑AACTGTCTAtaaccgtgttagtgaaaaaagtgg                                                                           >  1:403254/1‑228 (MQ=255)                                                                                                                                                              ggATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGtat                                                                          <  1:89550/251‑1 (MQ=255)                                                                                                                                                               gATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACGGTCTATCGCCGTATCAGCGAATAACGGtata                                                                         <  2:282151/251‑1 (MQ=255)                                                                                                                                                                  ggCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATAct                                                                       <  1:2230/250‑1 (MQ=255)                                                                                                                                                                   gCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑GTATCCCC‑GCGGAAG‑ACGGAAAAAAACTG‑CAGGGAGCGGCGAAAACGGGGAAAACGGGCATTTAACCTTGCGGTAAAAACACCGGAA‑GCACACCCAACT‑TCTATCCCCGTATCaacggaaaaccgtagaacg                                                                      >  1:77246/1‑232 (MQ=255)                                                                                                                                                                        ggTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatctg                                                                 <  2:372933/250‑1 (MQ=255)                                                                                                                                                                                            gCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGc                                               <  1:311361/248‑1 (MQ=255)                                                                                                                                                                                             ccGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGGAAA‑ACG‑AAAAAAACTGGCAAGGAGCGGGGGAAAAGGGGAAAAAGG‑CAATTAAACATGCCGATAAAAAACCGGAAAAGACACC‑AAATTGCTATCGCCGCATCCACGAAAAAAGGGAAACTGATCTGATCATTTTAATTTTAAGCACTg                                           >  1:100048/1‑251 (MQ=255)                                                                                                                                                                                                 cAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGt                                        <  1:167257/250‑1 (MQ=255)                                                                                                                                                                                                 cAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGt                                        <  1:40705/250‑1 (MQ=255)                                                                                                                                                                                                             acCATCGATTACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACatat                           <  1:407683/251‑1 (MQ=255)                                                                                                                                                                                                                      tACAGCGTCGGCTTTACTGATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTg                   <  1:331217/250‑1 (MQ=255)                                                                                                                                                                                                                                        gATATGGCTCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTg  <  1:14783/249‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       AAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGG‑TTATCCAC‑GCGAAAG‑ACG‑AAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTG  >  NC_000913/4619150‑4619629

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