Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,366,233	A→C	11.6%	V181G (GTG→GGG)	dcuA ←	C4‑dicarboxylate antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,366,233	0	A	C	11.6%	175.5 / 2.6	70	V181G (GTG→GGG)	dcuA	C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base A (49/11);  new base C (2/8);  total (51/19)
Fisher's exact test for biased strand distribution p-value = 2.55e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGC  >  NC_000913/4365999‑4366461                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     gtggtgTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGGGTTT‑GGAGTTGAAA‑AGAAAACGGACCAGG‑AAGGAAAT‑ACCAGa                                                                                                                                                                                                                                 >  2:286906/1‑239 (MQ=255)      gtTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCATAAcc                                                                                                                                                                                                                              >  1:550995/1‑238 (MQ=255)       ttCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCC‑ACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAA‑AGCAACGTGACCAGG‑AAGGACATCACCAGAACCGcc                                                                                                                                                                                                                           >  1:92586/1‑239 (MQ=255)               cGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGAACATCACCAGAACCGCCAGCAggg                                                                                                                                                                                                                    >  1:378622/1‑239 (MQ=255)               cGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTCCCAGAAGCTTGCG‑ATAAATCGGATCGTCAGAGAGGTT‑GGGGTTGAAC‑AGCAAAGTGACCAGG‑AAGGACAACCCAAGAACAGccgccggag                                                                                                                                                                                                                    >  2:1091260/1‑232 (MQ=255)                    tttCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGACCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCCGAACCCCACTCAGGGTgtca                                                                                                                                                                                                               >  1:37704/1‑236 (MQ=255)                              aCCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCAAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGGGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGACCCGCCAGCAGGGCGTACGGGATGAcc                                                                                                                                                                                                      >  1:702029/1‑238 (MQ=255)                                         gtgttttgttgTTTGCATAGATAACCCCCCCAACTACGCCCGGCAGGAACAGCGCGACGGACGTTTTTGCACC‑GGATTTGATTTAAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAACCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGGGACCAGG‑AAGGACAGCACAAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGa                                                                                                                                                                                          <  2:16409/234‑1 (MQ=255)                                          gCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGAt                                                                                                                                                                                         >  1:1063360/1‑239 (MQ=255)                                              ttGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGGGAGCAGATGGAg                                                                                                                                                                                     >  1:16409/1‑239 (MQ=255)                                              ttGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGGGGACGGGATGACCACGGAAAGcctgtccct                                                                                                                                                                                     >  1:1086531/1‑231 (MQ=255)                                              ttGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGGGACCAGG‑AAGGCCAACCCCAGAACCGCCAGCAGGGGGGACGGGATGGCCACGGACAGCaaaagggg                                                                                                                                                                                     >  2:756451/1‑232 (MQ=255)                                               taatgtttgaaaaaaTACCCACGCCAACTCCGCCCAGCAGAAACAGCCAGCCGGACGTTTTTGCCCG‑GGTTTGGATTTAAATCTTCTTTTCCCCGCGCAGTTCACCCAGGCCTTCTTCCAGACGTTTGCG‑ATAAATCGGATCGTAAGAGAGTTT‑GGAGTTGAAA‑AGCATAGTGACCAGG‑AAGGACATCCCCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAgg                                                                                                                                                                                    <  1:18389/225‑1 (MQ=255)                                                      tGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGACCCGCCATCAGGGGGGCCGGGATGACCACGGAGAGCAGAAGGAGGGAGATGa                                                                                                                                                                             >  2:941700/1‑239 (MQ=255)                                                        cATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGACC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGGGGACGGGATGACCACGGAGAGCAgcagaagtaagctgcg                                                                                                                                                                            >  1:673457/1‑223 (MQ=255)                                                        cATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCCCCAGCAGGGGGGACGGGAGGACCCCGGACAGCAAACGGAGGGAGCCGag                                                                                                                                                                            >  2:976138/1‑237 (MQ=255)                                                        cATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATAACCAGAACCGCCAGCAGGGAGGCCGGGAGGACCACGAAGAGCAGATGGAGGTacatggt                                                                                                                                                                            >  1:859114/1‑232 (MQ=255)                                                                   cccataaCTACGCCCAGGAGGAAGCGCCCGGCGGA‑GTTTTTGCCCCGGGATTTGTCTTCACTTGGCTTTCCACGGCGAGTTTCAACCGGGCCTTCTCACAGACGTTTGCG‑ATAAAGCGGACCGACAGAGTGTTT‑GGAGTTGAAC‑GGCATAGGGACCAGG‑AGGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGAcc                                                                                                                                                                   <  2:700651/231‑1 (MQ=255)                                                                    cGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGGGCAGATGGAGGTAGCTGATGCCATGACCTTcc                                                                                                                                                               >  1:605099/1‑239 (MQ=255)                                                                      caaactcggCCCAGCAGAAACAGCCAGGCGGACTTTTTTCCACG‑GGTTTTGATTCCAATCTGCTTTTCACCGCGCAGTCCAACCGGGCCCTCTTCCAGCCGCTTGCG‑AAAAATCGGATCGTAAGAGAGTTT‑GGAGTTGAAC‑GGCATAGTGCCCAGG‑AAGGACACCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCAt                                                                                                                                                             <  1:976276/231‑1 (MQ=255)                                                                       cAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGCTT‑GGCGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAAGCGGCACGGGATGACCACGGAGAGCAGCTGGAGTAGGCGGAGGCCATTACCTTCCac                                                                                                                                                             >  2:66248/1‑237 (MQ=255)                                                                            aCGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCCCGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGa                                                                                                                                                       >  1:676036/1‑239 (MQ=255)                                                                              ccccctcgtaaacccacAGACGGGTTTTTTTGCCCC‑TGTTTTTATTTCAATCTACTTTCCCCGGCGCAGTTCAACGAGCCCCTCTCCCAGACGCTTGCG‑AAAAATCGGATCGTCAGA‑AGTTTGGGGGTAAAAC‑AGGATGGTGCCGAGG‑AAGGACAACCCAAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAg                                                                                                                                                     <  2:471058/223‑1 (MQ=255)                                                                              cagcaacagGACCCGCCAGCCGGTCTTTTTTGCCCC‑GTATTTAATTTAAATCTCCTTTTCCCCGCGCAGTTCAACCGGGCCCTCTCACAGCCGTTTGCA‑ATAAATCGAAGCGTAAGAGAGTTG‑GGATTGGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGaa                                                                                                                                                      <  2:599476/232‑1 (MQ=255)                                                                                ccAGCAAGAACAGCCAGAGGGTGGTTTTAGCCCC‑GGTTTTGAGTTACATCTTCTTTTCCCCGCGGAGTTCAACCAGCCCCTTTCCCAGACGCTTGGG‑AAAAGTGGGATCGTCAGAGGGTTG‑GGAG‑TGAACAAGCTAAGTGACAGGG‑AGGGACACCCACAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGAc                                                                                                                                                   <  2:317946/239‑1 (MQ=255)                                                                                   gcagTAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATg                                                                                                                                                >  1:770430/1‑239 (MQ=255)                                                                                    cagGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGAAGCTGATGCCAGGACCTTCCATCAGGGAAGACAttg                                                                                                                                               >  1:471058/1‑237 (MQ=255)                                                                                            gCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTaaa                                                                                                                                            >  1:599476/1‑234 (MQ=255)                                                                                                 acggacgTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGGAAGCAGATGGAGGGAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACGTCCGc                                                                                                                                  >  1:754405/1‑239 (MQ=255)                                                                                                 acggacgTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACGGcc                                                                                                                                    >  1:629156/1‑237 (MQ=255)                                                                                                              cACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGCAGCTAATGCCATGACCTTCCATCACGGAAGACATGGAAACAACTGCCGCTGCGATTggg                                                                                                                        >  1:57414/1‑235 (MQ=255)                                                                                                                    tggttaaaagagacgggCTTTGCCCCGCGCAGTCACACCAGGCCCTCTTCCAGACTCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGAATAGGCACCAGG‑AAGGACTCCCCAAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCgg                                                                                                                 <  1:330447/222‑1 (MQ=255)                                                                                                                    gatttgatttCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTAGAATTGGcaagggggg                                                                                                                >  2:369475/1‑231 (MQ=255)                                                                                                                    gatttgatttCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGGAGCTGATGCCAAGACCTTCCATCACGGAAAACATGAAAAACAATCCCGCCGAGAATGGCGAGGCggg                                                                                                                >  2:1019346/1‑238 (MQ=255)                                                                                                                                      tttcCAGCGTGCAGGTCAACCAGGCCTTTCTCGAGACGCGTGCC‑ATAACGCGGTGCGTAAGAGAGTTG‑GGAGTGGAAC‑TGCATAGGGACGAGG‑AAGAAAATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGa                                                                                              <  1:422102/239‑1 (MQ=255)                                                                                                                                           ccccggcAGTTCACCCAGGCCCCGTTCAAGACGTGTCGG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTGGACC‑AGAATAGTGACTAGG‑AAGGAAATCCCAAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCa                                                                                         <  1:66248/234‑1 (MQ=255)                                                                                                                                           cccGGGCGGTTCAACCGGGCCCTC‑TCCAAGCGGTTGGGAAAAAATGGGATCGTAAGAGAGTTT‑GGAGTTGAAC‑AGAATAGGGACCAGG‑AAGGACCTCACCAGACCCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCa                                                                                         <  1:976138/238‑1 (MQ=255)                                                                                                                                           cccGCGCGGTCCAACCGGGCCCTCTTCCAGACGCTGGCG‑ATAAAGCGAGTCCTCAGAGAGTTT‑GGAGTGGAAC‑AGCAAAGGACCCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCa                                                                                         <  1:620907/238‑1 (MQ=255)                                                                                                                                             cgcgCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCCCGGCAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGGGCgggaacc                                                                                          >  2:49991/1‑231 (MQ=255)                                                                                                                                               cgcAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACCTCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGAAGATGGAGGAAGCTGATGCCATGACCTTCCATCACGGAAGAACTGTAAACCACTGCCGACGAGAGTGGCGATGCGGGGATCGCCATcctgtgtcttaacacac                                                                                      >  2:330447/1‑222 (MQ=255)                                                                                                                                               cgcAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGAGGCCATGACCTTCCATTACGGAAGACATGTAAACCACTGCCGCTTAAATTGGCGATTCGGGTATCGCAATCTGCGCGGATACCACTg                                                                                      >  1:961092/1‑238 (MQ=255)                                                                                                                                                              cccTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCa                                                                                         >  1:198565/1‑220 (MQ=255)                                                                                                                                                              cccTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGCGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGCCCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATGGGCTAGGCGGGGGCCGACATCTGCACGGCAACCACACCCATAGAACGCGGACg                                                                      >  2:349091/1‑239 (MQ=255)                                                                                                                                                                       aGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCAc                                                                                        >  1:1096404/1‑212 (MQ=255)                                                                                                                                                                       aGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTACCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATCCGGGGATCGCAATCTGGGCGGAAACCAACACAGAGAAAAGCGGAGGGCAAGGttt                                                             >  2:1048426/1‑239 (MQ=255)                                                                                                                                                                         ggcccgtgaaaTAAAGCGGATCGACAGAGTGTTG‑GGAGTTGAAC‑CGCATAGTGCCCAGG‑AAGCACATCACAAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTa                                                            <  1:286906/229‑1 (MQ=255)                                                                                                                                                                           ccgttggaaaaaaGCGGATCGTCGGAGATTTT‑GGATTGAAAC‑AGCATGGTGCCCAGG‑AGGGACACCCCCAAAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTa                                                            <  1:1048426/227‑1 (MQ=255)                                                                                                                                                                             ttGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGCCGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATCCGGGGATCGCAATCTCCGCGGATACCACTGCAGTAGACAGCGTAGGGCAAGGTTTAACg                                                         >  1:519750/1‑237 (MQ=255)                                                                                                                                                                             ttGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGCCCACGGAGAGCAGATGGAGGCAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACACCTGCCCCTGAGATTGGCGCTCCGGTAATCGCAATCTGCGCGGATCCCCCTCCAGAAGACAGCGGCCGCGAAGGTTTAAGGcc                                                       >  2:676152/1‑239 (MQ=255)                                                                                                                                                                             ttGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCAATGCAGTAGACAGCGGACGGCAAGGTTTACCg                                                         >  1:405180/1‑237 (MQ=255)                                                                                                                                                                             ttGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGAAACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACg                                                         >  1:820809/1‑237 (MQ=255)                                                                                                                                                                             ttGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCCCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACg                                                         >  1:144069/1‑237 (MQ=255)                                                                                                                                                                             ttGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACGGCCGCTGAGATTGCCGATGCGGTGAGCGCAATCTGCGCGGGTACCACTGCCGTAGACAGCGGACGGCAAGGGTTAACGcc                                                       >  2:839779/1‑239 (MQ=255)                                                                                                                                                                               gCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGGAATCGCAATCTGCGCGGATACAACTGCAGTGACCAGCGGACGGAAAGGTTTAACg                                                         >  1:782051/1‑235 (MQ=255)                                                                                                                                                                                  ggaaactgtgaggagcagaaggTTG‑GGAGTGAAAC‑GGGATGGTGACCAGG‑AGGAACACCCACAGAACCTCCCGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCGTCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGt                                                   <  2:629156/218‑1 (MQ=255)                                                                                                                                                                                                tCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTc                                     >  1:1098915/1‑239 (MQ=255)                                                                                                                                                                                                          t‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATa                              >  1:802881/1‑236 (MQ=255)                                                                                                                                                                                                                  gAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTAAGCGATAACTGGCAgt                     >  1:382550/1‑238 (MQ=255)                                                                                                                                                                                                                   aaC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTATCGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAgtg                    >  1:74697/1‑238 (MQ=255)                                                                                                                                                                                                                   aaC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAAATTCACCGACaaagggaagagt                   >  2:644322/1‑229 (MQ=255)                                                                                                                                                                                                                            aggggccAGGAAAGGACACCCCCAAAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCaga            <  1:49991/234‑1 (MQ=255)                                                                                                                                                                                                                                acaagg‑aaggaCATCACAAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCa                                                                                         <  1:239108/154‑1 (MQ=255)                                                                                                                                                                                                                                aCCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCa                                                                                         >  2:239108/1‑157 (MQ=255)                                                                                                                                                                                                                                aCCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCCGGGTCGCCAGAGAGAt       >  1:845255/1‑239 (MQ=255)                                                                                                                                                                                                                                  ccagcgcgccacacACCAGACCCCCCAGGAGGGGTGGCGGGATTACCCAGGAGGGCAGATGGAGGTTGCTTATGCCATGCCCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATTCCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGAt       <  2:505481/225‑1 (MQ=255)                                                                                                                                                                                                                                   agg‑aaggaCATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGCGATTTGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGAGGGCGGCAGGGTTTAACGCCTTGTTCCTTCTCAACTCAAGCGATAACGGGAAGTGTGGCCAGGGAGATGtt    >  2:105087/1‑239 (MQ=255)                                                                                                                                                                                                                                   agg‑aaggaCATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCCAGGGTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGt     >  1:518225/1‑238 (MQ=255)                                                                                                                                                                                                                                   agg‑aaggaCATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTTAGCGATAACTGGCAGTGTCGCCAGAGAGATGt     >  1:628886/1‑238 (MQ=255)                                                                                                                                                                                                                                   agg‑aaggaCATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGt     >  1:916681/1‑238 (MQ=255)                                                                                                                                                                                                                                       aaggaCATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGc  <  1:369475/238‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     GTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACC‑GGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCG‑ATAAATCGGATCGTCAGAGAGTTT‑GGAGTTGAAC‑AGCATAGTGACCAGG‑AAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGC  >  NC_000913/4365999‑4366461

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