Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,702	T→C	S158S (AGT→AGC)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,702	0	T	C	98.0%	147.0 / NA	51	S158S (AGT→AGC)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (27/23);  minor base G (1/0);  total (28/23)
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

CCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCAC  >  NC_000913/4618465‑4618947                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       cccGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCgg                                                                                                                                                                                                                                           <  2:344859/251‑1 (MQ=255)    gATTGTTGATAAACATTCCACGGGGGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑AGGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGcgg                                                                                                                                                                                                                                        <  2:409343/251‑3 (MQ=255)                  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)                                                                                                           cgcgGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACt                                                                                                                                                                                                                                                   <  1:329865/137‑1 (MQ=255)                                                                                                           cgcgGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACt                                                                                                                                                                                                                                                   >  2:329865/1‑137 (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)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       CCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACA‑TTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCAC  >  NC_000913/4618465‑4618947

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