Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,717	T→G	A163A (GCT→GCG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,717	0	T	G	98.0%	136.7 / NA	51	A163A (GCT→GCG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (29/20);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATC  >  NC_000913/4618482‑4618958                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   atccgccggggtgggtgggCGATGTGACTGCGCTGATTTTGGGGGCGGTGGGCGCAGCCTGTGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTAATACTGGCGGTATGCTTGACAAACGTG‑AATCCAGCCCTGGCTTCGACA‑TTGTCCCGG‑ATGACAACCTTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTAGCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTAc                                                                                                                                                                                                                                    <  2:211589/234‑1 (MQ=255)    ccACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGGTCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAAATGG‑AATCCATCCCTCGCCTCG‑CACTTTTTTCGG‑GTGGCAACCCTTTTCACGAA‑ATTTTTTAAGACGTCGGCGGGGGGGTTTTCGGTTAGACCAACTCACCGGGTCCGGCCGGGAAACGTTT‑TTAcgcc                                                                                                                                                                                                                                 >  2:162307/1‑250 (MQ=255)        tggTGGCGTCGTCGATGTGACTTCGCTGATGTTGTGGCCGATGGTCGCAGCCTGCGGCGGCGATATTCCGATGATCTCTGGGCTCGGCCTCGGGCATACTGTCGGTACGCTCGATTAATTGG‑AATCCGTCCCTGGCTTCGACA‑TTTTTCCGG‑ATGACAACCGTTTCCGCTAA‑ATGATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGAcc                                                                                                                                                                                                                              <  1:291433/249‑1 (MQ=255)         ggTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGACTCGGTCATACTGGCGGTACGCTCCACAAAAAGG‑ACTACATCCATGGCATCGACA‑TTTTCCCGG‑GTGAAAACCGTTTCCGCGAA‑ATTATTAAAGACGCCGGCGGGGCGATTATCGGTCAAACCAGGTCACTGGCCCCGGCGGATAAACGTTT‑TTAcacgcaccg                                                                                                                                                                                                                            >  1:300131/1‑243 (MQ=255)                  ggCGATGTGACTTCGCTGATGTTGGGGCCGATGGTAGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTAAGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑AATATTAAAGACGTCGGCGTGGGGATTTTCGGTTCGACCCACTCACTGGCTCCGGCGGGAAAACGTTT‑CTTCGCGACCCGTGATCTTAc                                                                                                                                                                                                                   >  1:424693/1‑251 (MQ=255)                   gCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTAcc                                                                                                                                                                                                                  <  1:51560/251‑1 (MQ=255)                                                     cGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGcc                                                                                                                                                                                 <  1:331936/250‑1 (MQ=255)                                                     cGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGcc                                                                                                                                                                                 <  1:368634/250‑1 (MQ=255)                                                          ccTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCCTGACAAACTGG‑GATTCATCCCTGGCTTCGTCC‑TTTTCCCGG‑ATGCCAACCGTTTCCCCGAA‑ATTATTAAAAACGGCGGCGTGGGGATTTTCGGGCCGACCAGCCCACCGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGAACTTACCGCCACCGTGGAATCCCTCCCGCTCATCACCGccccaatt                                                                                                                                                                           >  1:305665/1‑245 (MQ=255)                                                              cggcggCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACC‑TTTTCCCGG‑GTGACAACCGTTTCCGCGGA‑ATTATTAAAGACGTCCGCGTGGCGATTATCGGTCAAAACAGCTCACAGGCTCCGGCGGATAACAGGTT‑CTAAGCGGCCCGTGAAATTACCGCACACGTGGACCACATCCCGCTACTCCCCCCCTCTATTCTgg                                                                                                                                                                       >  1:406221/1‑251 (MQ=255)                                                                            cgattagctctGGCCGCGCCCTCGGTCATACTGGCGGTCAGCTCGAAAACCTGG‑AAACCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACACCCTTTCCCGAGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGAAAAACGTTT‑CTCCGCGACCCGTGATATTCCCGCAACCGGGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCg                                                                                                                                                         <  1:94667/244‑1 (MQ=255)                                                                                   ctctGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGGTTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGGCCCGTGATTTTACCGCAACCGTGGACCCCATCCCGCTGGTCCCCGCCCCTCTTCTTGGGCAGAAACCGGCGGAAgggc                                                                                                                                                  >  2:393980/1‑249 (MQ=255)                                                                                     ctGGTCGCGGCCTCGGTCATACTGGCGGTAAGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGGTT‑CTACGCGACCCGTGATATTACCGCAACCGAGGACTCCATCCCGCTGATCACCGCCTCTCTTCTGGCGAAGaaaatggcgggagggctg                                                                                                                                                >  1:111905/1‑236 (MQ=255)                                                                                       ggTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑GTGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGGAGAAACTGGCGGAAGGTCTgg                                                                                                                                               >  1:196691/1‑250 (MQ=255)                                                                                           ccgccccggtccatccgggCGGTACGCTGGAAAAACGGA‑AACCCACCCTTGGCTTCGAC‑‑TTTTCCCGGGATGACACCCTTTTCCGAGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGGGAATAAACGTTT‑CTACGCGCCCCGTGATATTCCCGCACCCGTGGACTCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGAcgcg                                                                                                                                          <  1:424468/234‑1 (MQ=255)                                                                                              gcctcggcaatacgggcgttccccTCGACAA‑CTGGGAACCCATCCCTGGCTCCGTGT‑TTTTCCCGT‑ATGAAACCCGTTTCCGCGAA‑ATTATTAAAGCCGGCGGTGTGGCGATTATGGGCCAGCCCAGCTCATTGGCGCGGGCGGATAACCGTTT‑CTACGGGACCCGTGTTATTACCGCACCCGGGGACTCCACCCGGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGGCTGGACGCGCTg                                                                                                                                       <  1:279764/228‑1 (MQ=255)                                                                                                    gTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGGGATATTACCGCAACCGTGGACTCCATCCCGCTGATCCCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAAGACTGGACGCGCTGGTGATg                                                                                                                                 >  1:476031/1‑251 (MQ=255)                                                                                                        aacgggCGGTACGCTCGAAAATCTTG‑AACCCCTCCCTGTCTTCAACT‑TTTCCCCGG‑TAGACACCCGTTCCCACGAA‑ATTATAAAAGACGTCGGCGGGGCTATTATCGGTCAGCCCAGCTCACTGGCTCGGGCGAATAAACGTTT‑CTACGCGCCCCGGGATATTCCCGAACCCGTGGACTCCATCCCGCTGATCTCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACg                                                                                                                             <  2:305665/247‑1 (MQ=255)                                                                                                            ggCGGTCCGCTCGACAAACTGG‑AATCCATCCCTGGCTTCAAC‑‑TTTTCCCGGGATGACAACCGTTTCCGCAAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTCCGCGCCCCGTGATATTCCCGCAACCGGGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGa                                                                                                                          <  1:365809/250‑1 (MQ=255)                                                                                                               ggTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGACACTGGCGGAAGGTCTGGCCGCGCTGGTgcgggcaggggaagt                                                                                                                      >  2:176533/1‑237 (MQ=255)                                                                                                                  aCGCGCGACAACCTGG‑AACCCCCCCCTGGCTTCGACA‑TTTCCCCGG‑AAGACAACGGTTTCCGCAAA‑ATCATAAAAGACGGCGGCGTGGCGATTATCGGTAAGACCAGCTCACGGGCTCCGGCGGAAAATCGTTT‑CTACGCGCCCTGTGATATTCCCGCAACCGTGGCCCCCACCCCGCTAATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTggg                                                                                                                   <  2:207518/251‑1 (MQ=255)                                                                                                                        gACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTTGACCCGCTGGGGATGGACGTGAAAGTGGGTAgggg                                                                                                             >  2:185994/1‑248 (MQ=255)                                                                                                                        gACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATAACAACCGTTTCCGCGAA‑ATTATTAAAGACGTAGGCGTGGCGATTATCGGTCAGACAAGCTCAATGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTCCCGCAACCGTGGACTCCATCCCGCTGGTCACCGCCTCTATTCTGGCGAAGCAACTGGCGGAGGGCCTGGACGCGCTGGGGATGGACGTggacgggggaagcgg                                                                                                             >  1:71509/1‑237 (MQ=255)                                                                                                                                    atccatccCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGc                                                                                                  >  1:59938/1‑251 (MQ=255)                                                                                                                                      ccatccCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCg                                                                                                >  1:459529/1‑251 (MQ=255)                                                                                                                                                 ttCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAAct                                                                                     >  1:469309/1‑251 (MQ=255)                                                                                                                                                   cGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTg                                                                                                                                                                        >  1:91486/1‑166 (MQ=255)                                                                                                                                                    gACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGGCCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGACGGTCTGGACGCGCTGGTGATGGACGTGAAAGGGGGTAGCGGCGCGTTTTTTCCGACCTACGAActctc                                                                                  >  2:369973/1‑251 (MQ=255)                                                                                                                                                     aCA‑TTTCCCCGG‑AGGCCAACC‑TTTCCCGGAATATTATTAAAGTCGGCGCCGTGGCATTTATCGGTCAGACCAGCTCACTGGCTCCGGCGAAAAAACTTTT‑CTACGCGCCCCGGGATATTCCCGAAACCGGGGCCTCCACCCCGCGGATCGCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAActct                                                                                   <  2:476031/249‑1 (MQ=255)                                                                                                                                                          tttCCCGG‑ATGACACCCGTTCCCGCAAA‑ATTATAAAAGACGTGGGCGTGGCGATTATCGGTCAGACCAGCTCACGGGCTCCGGCGGATAAACTTTT‑CTACGCGCCCCGTGATATTACCGCACCCGTGGCCTCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGaa                                                                              <  1:169502/250‑1 (MQ=255)                                                                                                                                                             cccGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGGCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGCCCTACgggatcgggagagcgg                                                                          >  2:62713/1‑236 (MQ=255)                                                                                                                                                             cccGA‑ATGACAACCGTTTCCGAGAA‑ATTATTAAAGCCGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTGCGCGCCCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACg                                                                                         <  1:62713/236‑1 (MQ=255)                                                                                                                                                                 ggaagacacccttttccggaaa‑aTTATTAAGGACGCGGGCGGGGCGATTATCGGTCAGACCAGCTCACTGGCCCCGGCGGAAAAACGTTT‑CTACGCGCCCCGGGATATTCCCGCACCCGGGGACTCCACCCCGCTGAGCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTg                                                                       <  1:472320/231‑1 (MQ=255)                                                                                                                                                                           cGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGAt                                                             >  1:84738/1‑251 (MQ=255)                                                                                                                                                                                  gcaaa‑aTTATTAAAGCCGTCGGCGTGGCGTTTATCGGCCAGCCCAGCTCACTGGCTCCGGCGGATAAACTTTT‑CTACGCGCCCCGTGATATTCCCGCACCCGTGGACCCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGc                                                      <  1:369973/248‑1 (MQ=255)                                                                                                                                                                                  gaaaa‑aTTATAAAAGACGTCGGCGTGGCGATTACCGGTCAGACCAGCTCACTGGCTCGGGGGGATAAACGTTT‑CTCCGCGCCCCGTGATATTCCCGCACCCGGGGACTCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGc                                                      <  1:403568/248‑1 (MQ=255)                                                                                                                                                                                   cgAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGAGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGAATCCATCCCGCTGATCACCGCCTCTCTTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGGCCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGGTGGCg                                                     >  1:69636/1‑251 (MQ=255)                                                                                                                                                                                     aa‑aTTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTg                                                   >  1:319490/1‑251 (MQ=255)                                                                                                                                                                                      aaattattaagtccggcGGCGTGGCGATTAGCGGTCAGCCCAGCCCACTGCCCCCGGCGGATAA‑CGTTTTCTCCCCGCCCCGTTTTATTCCCGCACCCGTGGACTCCTCCCCGTTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTg                                                   <  1:43123/235‑1 (MQ=255)                                                                                                                                                                                         ttattaAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAACGCTTGTTGGCGTGGCt                                                >  1:359745/1‑251 (MQ=255)                                                                                                                                                                                             taAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTTACg                                            >  1:399139/1‑251 (MQ=255)                                                                                                                                                                                                 gACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGCCCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCCCTATTCTGGCGAAGGAACTGGCGGAAGGTCTGGACGCGCTGGTGCTGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTAAGAACTCTCTGAAGCCCTTGCCGAAGCGGTTGTTGGCGTGGCCAACGgccc                                        >  1:254799/1‑249 (MQ=255)                                                                                                                                                                                                      ggggggggagATTGTCGGTCAGCCCAGCTCACTGGCTCGGGGGGAAAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGGGGACCCCATCCCGCTGATCTCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCg                                   <  1:438208/242‑1 (MQ=255)                                                                                                                                                                                                                       gTCAGACCAGCTCACTGGCTACGGCGGATAACCGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGAGGGTAGCGGCGCGTTTATGGCGAGCTACGAACTCTCTGAGGCCCTTGCCGAAGCGGTTGTAGGCGAGGCTAACGGGGCTGGCGTGGGCACCACcgcgttg                  >  2:93898/1‑248 (MQ=255)                                                                                                                                                                                                                         cAGAACAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCAACGCgctgct                >  1:155417/1‑251 (MQ=255)                                                                                                                                                                                                                          aGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGAGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGGCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGGTTGTTGGCGTGGCTAACGGCGCGGGCGGGCGCACCCCCGCgctgccc               >  2:320578/1‑249 (MQ=255)                                                                                                                                                                                                                            cccAGCTCACGGGCTCCGGCTGATAAACGTTC‑CTCCGCGCCCCTTGATAGTCCCGCAACCGTGTACCCCACCCCGCTGACCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCAc             <  1:35148/250‑1 (MQ=255)                                                                                                                                                                                                                              cagctCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTccccc           >  1:396823/5‑247 (MQ=255)                                                                                                                                                                                                                                  tCACTGGCTCCGGCGGATAACCGTTT‑CTACGCGCCCCGTGATATTCCCGCACCCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTggg                                                                                                                   <  1:334178/143‑1 (MQ=255)                                                                                                                                                                                                                                  tCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTggg                                                                                                                   >  2:334178/1‑143 (MQ=255)                                                                                                                                                                                                                                         cTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGCAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATc  >  1:70522/1‑249 (MQ=255)                                                                                                                                                                                                                                            cggcggATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGCAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTAc                                                                                          >  2:381819/1‑158 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   ACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATC  >  NC_000913/4618482‑4618958

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