Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,496	C→A	L91L (CTG→CTT)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,496	0	C	A	97.7%	132.0 / NA	43	L91L (CTG→CTT)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (9/33);  total (10/33)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCC‑ATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATT  >  NC_000913/3815253‑3815740                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     ccGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGcaggtcaaggtcat                                                                                                                                                                                                                                               <  1:152678/250‑3 (MQ=255)       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)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     ACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCC‑ATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATT  >  NC_000913/3815253‑3815740

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