Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,502	G→A	H89H (CAC→CAT)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,502	0	G	A	97.7%	135.8 / NA	44	H89H (CAC→CAT)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/34);  minor base C (1/0);  total (10/34)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

AACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCC‑ATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAAC  >  NC_000913/3815259‑3815748                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       aaCTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCt                                                                                                                                                                                                                                                 <  1:17763/251‑1 (MQ=255)     tCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCACGCGCCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGCTCACATCATCTACCCGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTTAAGCAGTACGGCAGGTCAAGGTCATGATCCTCCGc                                                                                                                                                                                                                                             >  1:439511/1‑251 (MQ=255)          aaTCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCATGATGCTCCGCCAgcg                                                                                                                                                                                                                                        >  1:349876/1‑249 (MQ=255)                  aaaTCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGGCAGGTCATGATGCTCCGCCCGCGCcacgggg                                                                                                                                                                                                                                >  1:228067/1‑247 (MQ=255)                  aaaTCTCGCCGCGCCCCCGGTCCTGACGAGCAAGGGAAATCAATACGCCAGCAAGGGTCGCGCGATGGGGCTGAGTAGTCGCC‑AGCGATGCGCGAATCGCCGTTCGGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAGCG‑CGCTACCAGCCAGATGGCCGCCTTCACCGGGGTCTTGTGCCTCGTTGCGGTTAAAGCAGTACGGCAGGTCAAGgtggtggtggtcggggggggcggcgggg                                                                                                                                                                                                                                <  2:158919/251‑23 (MQ=255)                                            cGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCttt                                                                                                                                                                                                      <  1:100923/251‑1 (MQ=255)                                              atcaaGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGt                                                                                                                                                                                                    <  1:36864/247‑1 (MQ=255)                                                    cGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAg                                                                                                                                                                                              <  1:61057/251‑1 (MQ=255)                                                    cGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTCTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCGCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAg                                                                                                                                                                                              <  1:397270/251‑1 (MQ=255)                                                       aatcaatccGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCGCGTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGc                                                                                                                                                                                           <  2:317140/243‑1 (MQ=255)                                                          caatacGCCAGCAAGCGTCGCGCCATTAGCCTGAATAACCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCaa                                                                                                                                                                                        <  2:177370/247‑1 (MQ=255)                                                              ccGCCAGCAAGCGTCGCGCCATGAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGTTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTTAGCAGGGCCAAAca                                                                                                                                                                                     <  2:429597/249‑1 (MQ=255)                                                                       aGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTGCCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTGTGTAAGCAGGGCCAAACAGCAGATCGaa                                                                                                                                                                           <  1:396111/251‑1 (MQ=255)                                                                            cgcgCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCaa                                                                                                                                                                      <  1:149029/251‑1 (MQ=255)                                                                             ccgccatgagCCTGAATAA‑CTCCCATCGATTCTCGAATCTCCGTGCCGGCGGTGGTCACACCATCTACCAGCATTCCGCGTCCTT‑TTACGTCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCGTTGCGGTTAAAGCAGTACGTCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATGGGGATCCCTTTGTAAGCATGGCCAAACTGCAGAGCGAACTCAAt                                                                                                                                                                     <  2:475419/242‑1 (MQ=255)                                                                                        cTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGGCATGATGCCCCGCCAGCGCCACGGCGGGAGGGGGCGCCAACGGGGTCCCTTTGTAAGCAGGGGCAAAACGCAGATAGAACTccattgcagaatccc                                                                                                                                                          >  1:171302/1‑237 (MQ=255)                                                                                         tGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGaa                                                                                                                                                                           <  2:356827/233‑1 (MQ=255)                                                                                          gAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCGCTTTGCGGTTAAAGCAGTACGTCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCTGTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATTCCAGAATcctcc                                                                                                                                                        <  2:257230/251‑4 (MQ=255)                                                                                                   cc‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCATGATGCTCCGCCAGCGCCACGGCGGGAGTGGGCGCCATCGGGATCCCCTTTTAAACAGGGCCAAACAGCAGAACGAACTCAACTCCAGAATTCACCAAAGCTTc                                                                                                                                                >  1:377798/1‑250 (MQ=255)                                                                                                            ccgcgAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCGTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGAGGCGCCGCGAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTGTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAATGCTTCAGCGTaaaa                                                                                                                                       <  2:439511/250‑1 (MQ=255)                                                                                                                gAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCaa                                                                                                                                                                                        <  2:75943/198‑1 (MQ=255)                                                                                                                    cGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCt                                                                                                                                <  1:258214/250‑1 (MQ=255)                                                                                                                          tCCGTCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGt                                                                                                                          <  1:74729/250‑1 (MQ=255)                                                                                                                              gcgttgttcatctcagctcccAGCATGACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCTTGGTCTTTTGCCTCGTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGAGGCTCCGCCAGCGCCACGGCTGTAGTGTTCGCAATCGGGATCCCTTTATAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCGTCAGCGTAAATACTGCCTAACAGGGCCAg                                                                                                                     <  2:402945/232‑1 (MQ=255)                                                                                                                                 gTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAggtc                                                                                                                  <  1:60538/251‑4 (MQ=255)                                                                                                                                       ccatcatCTACCAGCATTACGAGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTTAGGATGCTCCGCCAGCGCCACGGCAGTAGTGGGCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGccc                                                                                                            <  2:356393/250‑1 (MQ=255)                                                                                                                                              cTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGGTCCCTTTGTAAGCAGGGCCAAACAGCAGATCCAACTCAATGCCAGAATCCACCAATGCTTCAACGTAAAAAAGGCCTAACAAGGCCAGGGCGCGCCCGGGAtt                                                                                                      >  1:262865/1‑250 (MQ=255)                                                                                                                                                tcccgCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGACTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTa                                                                                                     <  2:283666/245‑1 (MQ=255)                                                                                                                                                aCCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTAtt                                                                                                      <  1:456851/248‑1 (MQ=255)                                                                                                                                                aCCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTaa                                                                                                    <  1:274004/250‑1 (MQ=255)                                                                                                                                                                  tttagcg‑cgCTCCCAACCAGAGTTCCGCCTTAAACGTGGTCTGTTGCCTATTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGTTGCGCCGGCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCGTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGCAATCCACAAATGCTACAGCGTAAAAACGGCCTAACAGTCCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTgaag                                                                                   <  2:358261/244‑1 (MQ=255)                                                                                                                                                                  tGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGGAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCGTCAGCTTAAAAACGGCCTAACAGTGCCAGGTCGCGCCGGGTATTAAACAGCCCGGCGTTgaag                                                                                   <  2:435376/249‑1 (MQ=255)                                                                                                                                                                       cg‑cgCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGGAGGGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCCGAAACCACCCATGCTTCCGCGTAAAAACGGCCTAACAGTGCCAGGGCGCGCCCGGGATTAAACAGCCCGGCGTTGGAGAAATAgg                                                                            >  1:183821/1‑251 (MQ=255)                                                                                                                                                                            cTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTACAGCAGTACGGCAGGTAACGGTCAGGATGCTCCGCCAGCGCCCCGGCGGTAGTTGGCGCCATCCGGATCCCTTTTTAAGCAGGGGCAAAAAACAAATAGAAATCAAATCCAAAAACCACCAAAGCTTCAACGTCAAAAAGACCCAAAAATGCCAGGTCGCGCCCCGGATTAAACAACCCCGCCGTGAAGAAAtcagggct                                                                        >  1:248260/1‑244 (MQ=255)                                                                                                                                                                              accaaccaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCaatg                                                                                                                                                    <  1:472985/173‑3 (MQ=255)                                                                                                                                                                               ccaaccaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCCGCGCCACGGCGGGAGTGGTCGCAATTGGGATCCCCTTGTAAGCAGGGCCAAACAGCCAAACGGAATCAAAGCCCGAATCCACCAATGCTTTAGCGGAAAAAAGGCCCAACAAGGCCAGGTCGCCCCCCGTATTAAACCACCCCGCGGTGAAGAAAAAGGGGGtttt                                                                     >  2:407299/1‑251 (MQ=255)                                                                                                                                                                                  accaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGTCCTGATGCTCCGCCAGCGCCACGGCGGGAGGGGGCGCAATCCGGATCCCTTTGTAAGCAGGGGCAAAAAGCCGAACGGACTCAAAGCCAGAATCCCCCAATGGTTCAACGTAAAAACGGCCTAACAATGCCCGGTCGCCCCCCGGGTTAAACAGCCCCGCGtttgagaaaaaaggggtttttcc                                                                  >  1:74925/1‑230 (MQ=255)                                                                                                                                                                                    caGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGcc                                                                <  1:192732/251‑1 (MQ=255)                                                                                                                                                                                          gccgccgcctccGTGGTCTTTTCCCTCTTTGCGGTTAAAGCAGTACTGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGTTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACACAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTgaa                                                                                    <  2:197342/215‑1 (MQ=255)                                                                                                                                                                                                                 cctcttTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATTATGCACCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACGCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGGCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACttt                                   <  2:278356/249‑1 (MQ=255)                                                                                                                                                                                                                         cGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtg                           <  1:89547/251‑1 (MQ=255)                                                                                                                                                                                                                         cGGTTAAAGCAGTACGGCAGGTCAAGGTCAGGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtg                           <  2:345057/251‑1 (MQ=255)                                                                                                                                                                                                                             tAAAGGAGTACGTCAGGTCAAGGTGATGATGCTCCGCCAGCGCCACGGCGGTATTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCATATCGAACTCAATGCCAGAATCCACCAAGGCTTCAGCGTAAAAACGGCCTAACAGTTCCAGGTCGCGCCCGGTATTAAACAGCCCGGCTTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgctg                        <  1:276761/250‑2 (MQ=255)                                                                                                                                                                                                                                          caggtcaaggtcATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAtt          <  2:183821/251‑1 (MQ=255)                                                                                                                                                                                                                                                  ggtcATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACGTTAACACCTGCTTGCTAAGCGCAAATTCAATAAAc  <  1:355585/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       AACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCC‑ATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAAC  >  NC_000913/3815259‑3815748

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