Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,535	G→C	A78A (GCC→GCG)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,535	0	G	C	97.8%	150.4 / ‑6.6	46	A78A (GCC→GCG)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (6/39);  total (7/39)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGG‑TGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCG  >  NC_000913/3815293‑3815775                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 cgTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCgg                                                                                                                                                                                                                                          <  2:92318/251‑1 (MQ=255)         aCGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCtt                                                                                                                                                                                                                                  <  1:57513/251‑1 (MQ=255)               aaGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAg                                                                                                                                                                                                                            <  1:346883/250‑1 (MQ=255)                   gAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAgg                                                                                                                                                                                                                        <  1:206272/251‑1 (MQ=255)                         aatacaccAGCAAGCGGCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCaaa                                                                                                                                                                                                                  <  2:86325/245‑1 (MQ=255)                                cAGCAAGCGTCCCGCCATTAGCCTGCATAATCACCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCACCCGTATTTCCTCCTTCACCGTGGTCTTTTGCCTCTTTTCGGTTAAAGCAGTGCGGCAGGTCAACGTCATGATGCTTCGCCAGCGCCACGGCGGTAGTGT‑TCGCAATCGGTATCCCGTTGTAAGCAGGTCCAAAcagcagt                                                                                                                                                                                                           <  1:426872/251‑2 (MQ=255)                                  gCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑ACGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATc                                                                                                                                                                                                         <  1:37491/251‑1 (MQ=255)                                       cGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACt                                                                                                                                                                                                    <  1:147455/251‑1 (MQ=255)                                         tCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGACTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCATTATGCTCCGCCAGCGCCACGGCGGGAGTGGGTGGAAACTGGGA‑CCCTTTGTAAGCAGGGGCCAAcaggcaatggacctcc                                                                                                                                                                                                  >  1:210915/1‑238 (MQ=255)                                           gcgcCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGGCATGATGCTCCGCCAGCGCCACGGCGGGAGTGG‑TCGCCAACGGGGTCCCCTTTTAAGGAGGGCCAAAAAGGCGATCGAACACaa                                                                                                                                                                                                 >  1:249792/1‑250 (MQ=255)                                                       tGAATAATCTCCATCGATTCGCGAATCGCCGTTCAGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGCCAAGGGCATGATGCCCCGCCAGCGCCACGGCGGGTGGGG‑GCGCCAACGGGAACCCCTTGTAAACAGGGCCAAACAGCCGATAGAAcccacaacccgaatcccc                                                                                                                                                                                    >  1:143015/1‑234 (MQ=255)                                                            aatCTCCATCGATGCGCTAATCGCCGTTCCGGCGGTGATCACATCACCTACCATCAGTACGCGTCCTTGTAACGCTCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTTTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCTGTAGTGG‑TCGCAATCGGGATCCCTGTGGAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAAGCCCCCaatg                                                                                                                                                                               <  2:174347/251‑3 (MQ=255)                                                                                    cGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACa                                                                                                                                                       <  1:234449/251‑1 (MQ=255)                                                                                     gTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAg                                                                                                                                                      <  1:16742/251‑1 (MQ=255)                                                                                            cggTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTCCCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAgg                                                                                                                                               <  2:58531/251‑2 (MQ=255)                                                                                                gATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTcg                                                                                                                                            <  1:485528/250‑1 (MQ=255)                                                                                                             aCCAGCATTACGCGTCCTTGGTACGCGCTTCCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACTGCCTAACAGTGCCAGGTCGCGCCCGGTAtt                                                                                                                                 <  2:352733/248‑1 (MQ=255)                                                                                                                aGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAg                                                                                                                           <  1:118606/251‑1 (MQ=255)                                                                                                                   aTTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTa                                                                                                                                <  2:88227/243‑1 (MQ=255)                                                                                                                   aTTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGGGG‑TCGCAATCGGGATCCCTTTTTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGCAAAAAAGGCCCAACAGTGCCAgggcccgcccgggtttt                                                                                                                                >  1:88227/1‑227 (MQ=255)                                                                                                                     tACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCg                                                                                                                       <  1:536291/250‑1 (MQ=255)                                                                                                                        gcgTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGt                                                                                                                   <  1:333036/251‑1 (MQ=255)                                                                                                                            ccTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACTGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTgaa                                                                                                               <  2:529821/251‑1 (MQ=255)                                                                                                                             cgtGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCGTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAATGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCCCCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCTGTATTAAACAGCCCGGCGTTgaag                                                                                                              <  2:105650/249‑1 (MQ=255)                                                                                                                                gTAACGCTCTACCAACCAGCATGCCGCCTTCCCCGTGGTCTTTTGCCGCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCATTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGaaa                                                                                                           <  2:73596/251‑1 (MQ=255)                                                                                                                                    cgcgCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAgg                                                                                                       <  1:235478/251‑1 (MQ=255)                                                                                                                                      cgcTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAgggg                                                                                                     <  1:454800/251‑1 (MQ=255)                                                                                                                                           ccaaccaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCtttt                                                                                                <  1:226473/251‑1 (MQ=255)                                                                                                                                                   aTTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCgg                                                                                        <  1:462967/251‑1 (MQ=255)                                                                                                                                                      gccgccTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGAtt                                                                                     <  1:26173/251‑1 (MQ=255)                                                                                                                                                      gccgccTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGAtt                                                                                     <  1:149784/251‑1 (MQ=255)                                                                                                                                                       ccgccTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGAttt                                                                                    <  1:7720/251‑1 (MQ=255)                                                                                                                                                        cgccTTCACCGTGGTCTTTTGCCTCGTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCATGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTc                                                                                   <  1:35386/251‑1 (MQ=255)                                                                                                                                                            ttCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGt                                                                                                                                                                      <  2:326408/164‑1 (MQ=255)                                                                                                                                                             tCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCg                                                                               <  1:380996/250‑1 (MQ=255)                                                                                                                                                                    ggTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCg                                                                       <  1:52213/251‑1 (MQ=255)                                                                                                                                                                                cttTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAACCCCCCAAGGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATATGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTa                                                             <  2:338287/249‑1 (MQ=255)                                                                                                                                                                                       gTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgct                                                    <  1:159385/251‑1 (MQ=255)                                                                                                                                                                                                     gCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAt                                      >  1:440436/1‑251 (MQ=255)                                                                                                                                                                                                     gCAGGTCAAGGTCATGATGCGCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGCACGAACCCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGGCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTGTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAt                                      <  2:56262/251‑1 (MQ=255)                                                                                                                                                                                                      caggtcaaggtcATGATGCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAtt                                     <  1:46344/251‑1 (MQ=255)                                                                                                                                                                                                                       gCTCCGCCAGCGCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGGAGGGCCAAACAGCAGCTCGTACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCCGGTCGCGCCCGGGATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTTCGCCCGGGTTTCAGCGTAAACTCGCCAAACTTTAACCCCTGCTTGCTAAACGCCAATTCAATAAACTTGCGCTgt                    >  1:66546/1‑250 (MQ=255)                                                                                                                                                                                                                              cagcgCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCTTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGtt             <  1:332797/247‑1 (MQ=255)                                                                                                                                                                                                                                gcgCCACGGCGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTc           <  1:328278/249‑1 (MQ=255)                                                                                                                                                                                                                                       ggCGGTAGTGG‑TCGCACGCGGGATCCTGTTGTAAGCAGGGCCAAACAGCAGATCGACCTCAATGCCAGAAGCCACCAATGCTTCAGCGTAAACACCGCCTAACCGTGCCAGGTCACGCCCGGTATTAAACAGCCCGGCGTTGAAGATATAGGGGCTTTTGCGGCCGGATTTCAGCGTAACCTCGCCACACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATCCCtt    <  2:287320/250‑1 (MQ=255)                                                                                                                                                                                                                                         cGGTAGTGG‑TCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCg  >  1:366404/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 CGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGG‑TGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCG  >  NC_000913/3815293‑3815775

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