Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	80,529	C→G	G112G (GGG→GGC)	leuC ←	3‑isopropylmalate dehydratase subunit LeuC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	80,529	0	C	G	95.7%	135.8 / ‑1.0	47	G112G (GGG→GGC)	leuC	3‑isopropylmalate dehydratase subunit LeuC
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base G (21/24);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

CGCTACCGGTTTTACCGATAATTGCCAGCACGATATCTT‑TTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTCATGCACCAGG  >  NC_000913/80285‑80775                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       cgcTACCGGTTTTACCGATAATTGCCAGCACGATATCTT‑TTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGCCCCAtt                                                                                                                                                                                                                                                  <  1:160533/251‑1 (MQ=255)    tACCGGTTTTACCGATAATTGCCAGCACGATATCTT‑TTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACg                                                                                                                                                                                                                                               >  1:231141/1‑251 (MQ=255)         gTTTTACCGATAAGTGCCAGCACGATATCTT‑TTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGGGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGAc                                                                                                                                                                                                                                          <  1:295204/251‑1 (MQ=255)                     agtGCCAGCACGAAATCTT‑GTGCGGTAATGCCCGGTGCGGCTTTGCCCTGGACTTCAATTTGCATGGTTTTTGCGCGGCCCTGTTTCAGGGTGTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTAATCCCCGGCAAGGTGACGATCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATa                                                                                                                                                                                                                              <  2:373405/249‑1 (MQ=255)                        gCCATCACGATA‑CTTGTTGCGGTAATGTCCGGTGCGGCTTTGCCCTGGACTTTAAGTTTCATGGTTTGTGCGCGGCCCTGTTGCAGGGTTTGCGTTGTCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCATGGCGACAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGTCGATCCCCTGAtccg                                                                                                                                                                                                                            <  2:404473/250‑4 (MQ=255)                          cAGCACGATATCTT‑TTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGGCGATCCCCTGGTACgggg                                                                                                                                                                                                                         >  2:211666/1‑250 (MQ=255)                                                   gCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACg                                                                                                                                                                                                  <  1:181223/250‑1 (MQ=255)                                                     ccGGCGCGGCTTTGCCCTGGACTGCAAGTTGCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGGGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCa                                                                                                                                                                                               <  2:346343/251‑1 (MQ=255)                                                          gcgGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTc                                                                                                                                                                                          >  1:299338/1‑251 (MQ=255)                                                            ggCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCtt                                                                                                                                                                                        <  1:344484/251‑1 (MQ=255)                                                             gCTTTGCACTGGACTTCAATTTTCATGGTTTTTGCTCGGCCCTGTTGCAGGGTTTGCTTTGCCAGTACGTGTTCAACTGCGGAAGGGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTGCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCtt                                                                                                                                                                                        <  2:444251/250‑1 (MQ=255)                                                                ttGCCCTTGACTTCAATTTTCATTGTTTTTGCGCTGCCCTGTTTCAGGGTTTGCGTTGCCAGTACTTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCTGGCCCATTACGTGGACGATCCCCTGATACGGGTTATTCAGGTCATACAGTTCGACGCCAAATTCTTTg                                                                                                                                                                                      <  1:381101/249‑1 (MQ=255)                                                                 tGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGGCATACAGTTCGACGCCAAATTCTTTGc                                                                                                                                                                                     >  2:138432/1‑249 (MQ=255)                                                                                 tttCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTcc                                                                                                                                                                    <  2:30667/250‑1 (MQ=255)                                                                                    cATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACTCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCa                                                                                                                                                                <  2:185825/251‑1 (MQ=255)                                                                                     aTGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGCCGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGCTCAGTTCCTGCAt                                                                                                                                                               >  2:111173/1‑251 (MQ=255)                                                                                          ttttGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATACGAGTCGCCGCAGCCAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCTGGTCCTACAGCTCGACGCCAAATTCTTTGCAGTTTTTTATCCGTTCCTGCATCTgg                                                                                                                                                           >  2:396326/1‑250 (MQ=255)                                                                                               cgcgGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATAcgc                                                                                                                                                     <  1:138432/251‑1 (MQ=255)                                                                                               cgcgGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATAcgc                                                                                                                                                     <  1:211666/251‑1 (MQ=255)                                                                                                  ggCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGcc                                                                                                                                                  <  1:94049/251‑1 (MQ=255)                                                                                                               ggTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGCTACCAAAGGCCAGTGCGACAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGCCAATGGTCATCCCCGGCAAGGTGACGCCATGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGATGATTCAGGTCATCCAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCCTCTGGATCCGCGCCCTTTCCCCGCAgg                                                                                                                                     >  2:397130/1‑251 (MQ=255)                                                                                                                            aGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCttt                                                                                                                        >  1:345910/1‑251 (MQ=255)                                                                                                                                 gtgtTCCACTTCGGAAGTGCCGATACCAAACGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACACTGGTCATCCCCGGCACGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACACCAAATTCTTTGCGGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTggact                                                                                                                   >  2:287584/1‑248 (MQ=255)                                                                                                                                                          ccAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCCTTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGTGACGTTGTGATCCATggg                                                                                          >  2:385261/1‑250 (MQ=255)                                                                                                                                                            aaaGGCCAGTGCCCCACACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAg                                                                                        <  2:257730/251‑1 (MQ=255)                                                                                                                                                                 ccAGTGCGCCCACCGCGCCGTGGGTGGCGGTATGCGAGTCGCCTCAGACAATGGTCATCCCCGGCAAGGTGGCGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGCGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCTTTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGCGCCGTTGTGTTCCATGGTCGCGagg                                                                                   >  2:32154/1‑249 (MQ=255)                                                                                                                                                                     tgctccAAACGCGCCGTGGGGGGCGGTATGATAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCGACTTCTTTGCAGTTTTTGATCCGTTCCTGCATCTGGATACGCGCCGTTTCACCGCGGGCTTTACTGTCTTTGGTCTGGGAAGATACGTTGTGATCCATGGTAACGAAGGttt                                                                               >  2:305151/5‑251 (MQ=255)                                                                                                                                                                        gcCAAACGCGCCGTGGGTGGCGGGATGCGAGCCGCCGCAGACAATGGGCATACACGGCAAGGTGCCGCCCTGTTCCGGGCCCATTGCGTGGAAGATCCCCTGATACGGGTGATTCAGGTAAGACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCCACGGTATAGGCACCATTTCACCGACGGCATTAATTTCTTTGGGCTGGGTAGCGACGTTGTGATCAATGGTAGAGATGGtttttg                                                                            >  2:469441/1‑249 (MQ=255)                                                                                                                                                                                 gcCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGAc                                                                   >  1:334990/1‑251 (MQ=255)                                                                                                                                                                                 gcCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGAGCCATGGTAGCGAAGGTTTTGCCCGGCTGAc                                                                   <  2:336650/251‑1 (MQ=255)                                                                                                                                                                                 gcCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTACTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGAc                                                                   >  1:377398/1‑251 (MQ=255)                                                                                                                                                                                       ggtggCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCAACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCg                                                             <  2:229409/251‑1 (MQ=255)                                                                                                                                                                                          ggCGGTATGCGTGTCGCCGCAGACAATTGTCATCCCCGGCAAGGCGACGCCCTGTTCCGGGCCCATTACGAGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCTAAATTCTTTGCAGTTTTTGATCAGTTCATGCATCTGGGTAAGACTCATTTCACCGCAGGCGTTAATGTCTTTGGTCTGGGTTGAGACGCTGTGATCCATGGGAGCAAAGGTTTTGCCGGGCTGAGGTCCCGGGc                                                          >  2:371308/1‑251 (MQ=255)                                                                                                                                                                                          ggCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTGCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCggg                                                           <  1:302600/250‑1 (MQ=255)                                                                                                                                                                                             ggTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATCCGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTCGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGa                                                        >  2:383239/1‑250 (MQ=255)                                                                                                                                                                                                aTGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGt                                                    >  1:393889/1‑251 (MQ=255)                                                                                                                                                                                                 ggCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTg                                                   <  1:77082/250‑1 (MQ=255)                                                                                                                                                                                                    gAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCGGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGTCGACCGTGGgc                                                >  2:113135/1‑251 (MQ=255)                                                                                                                                                                                                      gTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGgcgc                                              <  1:383239/251‑1 (MQ=255)                                                                                                                                                                                                                     ggTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGaccg                               >  2:457747/1‑248 (MQ=255)                                                                                                                                                                                                                     ggTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACg                               <  1:314599/251‑1 (MQ=255)                                                                                                                                                                                                                               ggCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGa                     <  1:256676/251‑1 (MQ=255)                                                                                                                                                                                                                                gcaagtggACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATGTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGTCGCGCAGACCATCGAACGCCGGCGGTGAg                    <  2:444380/244‑1 (MQ=255)                                                                                                                                                                                                                                      tGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACt              >  1:286854/1‑251 (MQ=255)                                                                                                                                                                                                                                           cccTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTCATg         >  1:132813/1‑251 (MQ=255)                                                                                                                                                                                                                                                 tCCGGGCCCATTACGTGGACGATCCCCTGATTCTGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGTTAGATACGTTTTGATCCATGGTAGCGAAGTTTTTGCTCGGCTGACGTACCGGGCTACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTCATGCACCa    <  1:270379/250‑1 (MQ=255)                                                                                                                                                                                                                                                  ccGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTCATGCACCAgg  <  1:199247/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       CGCTACCGGTTTTACCGATAATTGCCAGCACGATATCTT‑TTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTCATGCACCAGG  >  NC_000913/80285‑80775

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