Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,804	T→G	L192L (CTT→CTG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,804	0	T	G	98.6%	221.9 / ‑3.8	69	L192L (CTT→CTG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (38/30);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

CGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACATTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTT  >  NC_000913/4618557‑4619046                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      cgattagctctGGCCGCGCCCTCGGTCATACTGGCGGTCAGCTCGAAAACCTGG‑AAACCATCCCTGGCTTCGACATTTTCCCGG‑ATGACACCCTTTCCCGAGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGAAAAACGTTT‑CTCCGCGACCCGTGATATTCCCGCAACCGGGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCg                                                                                                                                                                                                                                                    <  1:94667/244‑1 (MQ=255)        ctctGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACATTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGGTTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGGCCCGTGATTTTACCGCAACCGTGGACCCCATCCCGCTGGTCCCCGCCCCTCTTCTTGGGCAGAAACCGGCGGAAgggc                                                                                                                                                                                                                                             >  2:393980/1‑249 (MQ=255)            ggTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACATTTTCCCGG‑GTGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGGAGAAACTGGCGGAAGGTCTgg                                                                                                                                                                                                                                          >  1:196691/1‑250 (MQ=255)                ccgccccggtccatccgggCGGTACGCTGGAAAAACGGA‑AACCCACCCTTGGCTTCGAC‑TTTTCCCGGGATGACACCCTTTTCCGAGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGGGAATAAACGTTT‑CTACGCGCCCCGTGATATTCCCGCACCCGTGGACTCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGAcgcg                                                                                                                                                                                                                                     <  1:424468/234‑1 (MQ=255)                   gcctcggcaatacgggcgttccccTCGACAA‑CTGGGAACCCATCCCTGGCTCCGTGTTTTTCCCGT‑ATGAAACCCGTTTCCGCGAA‑ATTATTAAAGCCGGCGGTGTGGCGATTATGGGCCAGCCCAGCTCATTGGCGCGGGCGGATAACCGTTT‑CTACGGGACCCGTGTTATTACCGCACCCGGGGACTCCACCCGGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGGCTGGACGCGCTg                                                                                                                                                                                                                                  <  1:279764/228‑1 (MQ=255)                         gTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACATTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGGGATATTACCGCAACCGTGGACTCCATCCCGCTGATCCCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAAGACTGGACGCGCTGGTGATg                                                                                                                                                                                                                            >  1:476031/1‑251 (MQ=255)                             aacgggCGGTACGCTCGAAAATCTTG‑AACCCCTCCCTGTCTTCAACTTTTCCCCGG‑TAGACACCCGTTCCCACGAA‑ATTATAAAAGACGTCGGCGGGGCTATTATCGGTCAGCCCAGCTCACTGGCTCGGGCGAATAAACGTTT‑CTACGCGCCCCGGGATATTCCCGAACCCGTGGACTCCATCCCGCTGATCTCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACg                                                                                                                                                                                                                        <  2:305665/247‑1 (MQ=255)                                 ggCGGTCCGCTCGACAAACTGG‑AATCCATCCCTGGCTTCAAC‑TTTTCCCGGGATGACAACCGTTTCCGCAAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTCCGCGCCCCGTGATATTCCCGCAACCGGGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGa                                                                                                                                                                                                                     <  1:365809/250‑1 (MQ=255)                                    ggTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACATTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGACACTGGCGGAAGGTCTGGCCGCGCTGGTgcgggcaggggaagt                                                                                                                                                                                                                 >  2:176533/1‑237 (MQ=255)                                       aCGCGCGACAACCTGG‑AACCCCCCCCTGGCTTCGACATTTCCCCGG‑AAGACAACGGTTTCCGCAAA‑ATCATAAAAGACGGCGGCGTGGCGATTATCGGTAAGACCAGCTCACGGGCTCCGGCGGAAAATCGTTT‑CTACGCGCCCTGTGATATTCCCGCAACCGTGGCCCCCACCCCGCTAATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTggg                                                                                                                                                                                                              <  2:207518/251‑1 (MQ=255)                                             gACAAACTGG‑AATCCATCCCTGGCTTCGACATTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTTGACCCGCTGGGGATGGACGTGAAAGTGGGTAgggg                                                                                                                                                                                                        >  2:185994/1‑248 (MQ=255)                                             gACAAACTGG‑AATCCATCCCTGGCTTCGACATTTTCCCGG‑ATAACAACCGTTTCCGCGAA‑ATTATTAAAGACGTAGGCGTGGCGATTATCGGTCAGACAAGCTCAATGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTCCCGCAACCGTGGACTCCATCCCGCTGGTCACCGCCTCTATTCTGGCGAAGCAACTGGCGGAGGGCCTGGACGCGCTGGGGATGGACGTggacgggggaagcgg                                                                                                                                                                                                        >  1:71509/1‑237 (MQ=255)                                                         atccatccCTGGCTTCGACATTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGc                                                                                                                                                                                             >  1:59938/1‑251 (MQ=255)                                                           ccatccCTGGCTTCGACATTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCg                                                                                                                                                                                           >  1:459529/1‑251 (MQ=255)                                                                      ttCGACATTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAAct                                                                                                                                                                                >  1:469309/1‑251 (MQ=255)                                                                         gACATTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTT‑CTACGCGGCCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGACGGTCTGGACGCGCTGGTGATGGACGTGAAAGGGGGTAGCGGCGCGTTTTTTCCGACCTACGAActctc                                                                                                                                                                             >  2:369973/1‑251 (MQ=255)                                                                          aCATTTCCCCGG‑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‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGCAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑c                                                                                            >  1:70522/1‑249 (MQ=255)                                                                                                                                                                cggcggATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGCAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTAc                                                                                                                                                                                     >  2:381819/1‑158 (MQ=255)                                                                                                                                                                        aaaCGTTT‑CTACGCGCCCCGAGATCTTACCGCAACAGGGGACTCCCTCCCGCTGATCACCGCCTCTATTCTGGCGAAGAACCTGGCGGACGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGAAGCGGCGGGTTTATGCCGCCCACCGCACTCTCTGCAACCCTGGCCGAAGCCGTTGTTTGCGTGGCTAACAGCATTGGCGGTCGGCCCACCGCGCTGGTCACCCACTCGACT‑C‑AGGTACTGGACTc                                                                              >  2:269028/1‑251 (MQ=255)                                                                                                                                                                           cGTTT‑CTACGCGCCCTGTGATTTTCCCGCACCCGGGGACTCCCTCCCGCTGATCCCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCa                                                                            <  1:406816/250‑1 (MQ=255)                                                                                                                                                                               gtccccgccacctgtataatacccGCACCCGTGGACCCCACCCCGCTGACCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTg                                                                         <  2:84738/228‑1 (MQ=255)                                                                                                                                                                                      cgACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTa                                                                                                                                          >  2:314844/1‑180 (MQ=255)                                                                                                                                                                                      ccccccGTAATATTACCCCACCCGGGGCCTCCACCCCTCTGATCACCGCCTCTATACTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAAc                                                                <  1:119662/248‑1 (MQ=255)                                                                                                                                                                                          ccGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGCCCTGAAT‑C‑AGGTACTGGCCTCCCGTGCAGG‑AAACGCgg                                                            >  2:355513/1‑251 (MQ=255)                                                                                                                                                                                              gATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTAc                                                                                                                                                                                     >  1:150917/1‑129 (MQ=255)                                                                                                                                                                                                      cccacccgtgatcccctccccGCTGACCTCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCTCTGGCGTGCGCGCCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTg                                                <  1:284778/233‑1 (MQ=255)                                                                                                                                                                                                        cAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGAGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGACCTCTCTGAAGCCCTTACCGAAACGGTTGTTGGCGTGGCTCACGGCGCTGGCGTGCGCACCCCCGCGCAGCTCACCCACCTGGATAC‑AG‑TACAGGCCCCCCGTGCAGG‑TAACGCGGTTGagcgtcggta                                               >  2:105605/1‑241 (MQ=255)                                                                                                                                                                                                         atcgggtgACTCCACCCCACTGACCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAg                                             <  1:76775/244‑1 (MQ=255)                                                                                                                                                                                                         aaCCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAg                                                                                                                                                                                                           >  1:366266/1‑96 (MQ=255)                                                                                                                                                                                                             gTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTTAAGCGGt                                         >  1:388590/1‑251 (MQ=255)                                                                                                                                                                                                                    ccATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGGGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCCCACCGACATGAAT‑CCAGG‑ACTGGCACCCCGGGCAGG‑GACCGCGGGGGAAGTTCGGGAAGCGGTGGAGTTc                                  >  2:446711/1‑251 (MQ=255)                                                                                                                                                                                                                     ccccccGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTcc                                 <  1:454073/248‑1 (MQ=255)                                                                                                                                                                                                                         ccGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGAc                             >  1:122097/1‑251 (MQ=255)                                                                                                                                                                                                                          cGCTGATCCCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACg                            <  1:320578/251‑1 (MQ=255)                                                                                                                                                                                                                            cTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGCAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCACTGCAGGGTAA‑GCGGTTGAAGGTCGGGAAGGGGTGCAGTTCCTGGCggg                          >  2:365441/1‑251 (MQ=255)                                                                                                                                                                                                                              gattgccgccgCTATTCGGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTg                        <  2:274578/246‑1 (MQ=255)                                                                                                                                                                                                                              agccgccgactctATTCTGGCGAAGAAACTGGCGGAATGTCTGGACGCGCTGGTTATGGACGTGAAAGTGGTGAGCGTTGCTTTTATGCCGACCTACGATCTCTCTGAAGCCCTTGCCGAATCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCCCCGCGCTGCTCTCCGACATGAAT‑C‑AGGTACTGTCCTCTATTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTTCAGTTTCTGACGGGt                         <  2:468739/241‑1 (MQ=255)                                                                                                                                                                                                                                   ccGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGCAT‑C‑AGGTACTGGCCTCCCGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAgtat                   >  1:363682/1‑251 (MQ=255)                                                                                                                                                                                                                                   ccGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAgtat                   <  1:262023/251‑1 (MQ=255)                                                                                                                                                                                                                                    cgccgcTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATc                  <  2:389216/251‑1 (MQ=255)                                                                                                                                                                                                                                          tATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCCGG‑TAACGCGGTTGAAGTTCGTGAAGCGGGGCAGTTCCTGACGGGTGGGTATCGTAAcc            >  1:112292/1‑251 (MQ=255)                                                                                                                                                                                                                                           atgcTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTaa              <  1:176533/245‑1 (MQ=255)                                                                                                                                                                                                                                           aTTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCTCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAG‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAAc             <  2:396823/249‑1 (MQ=255)                                                                                                                                                                                                                                              cTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCCCCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCCGG‑GAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGTATCGTAACCcgcg        >  2:45486/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCt     <  1:393980/251‑1 (MQ=255)                                                                                                                                                                                                                                                 gCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCACTTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCt     <  1:284241/251‑1 (MQ=255)                                                                                                                                                                                                                                                    aagaaACTGGCGGAAGGTCTGGACGCGCTGGTGTTGGACGTGAAAGTGGGTAGCCGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGTAT‑C‑AGGGACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGtt  >  1:398955/1‑251 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      CGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGG‑AATCCATCCCTGGCTTCGACATTTTCCCGG‑ATGACAACCGTTTCCGCGAA‑ATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTT‑CTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAT‑C‑AGGTACTGGCCTCCAGTGCAGG‑TAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTT  >  NC_000913/4618557‑4619046

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