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	95.2%	131.0 / ‑6.2	42	H89H (CAC→CAT)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/30);  minor base C (2/0);  total (12/30)
Fisher's exact test for biased strand distribution p-value = 7.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCGGAATCCA‑CCAACGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC  >  NC_000913/3815255‑3815741                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     ccTCAACTTCCTGAACCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGa                                                                                                                                                                                                                                                   <  1:423237/250‑1 (MQ=255)         tCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGGCATGATGCTcccc                                                                                                                                                                                                                                           >  1:422325/1‑249 (MQ=255)             gAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAgc                                                                                                                                                                                                                                       <  2:146882/251‑2 (MQ=255)               aTCGCCGAAATCTCCCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAgcgc                                                                                                                                                                                                                                     <  2:191190/251‑4 (MQ=255)                        aTCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTGCTGGTGCTCCggccgcgccaccggggg                                                                                                                                                                                                                            >  1:382223/1‑235 (MQ=255)                          ctcGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCAAGATGGTCCGCCAGCGccccggcggggg                                                                                                                                                                                                                          >  2:327145/1‑241 (MQ=255)                                 cgcCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAgtggtcgc                                                                                                                                                                                                                   <  1:9585/251‑4 (MQ=255)                                    ccccgTTCCGGACGATCAAGCGAAACCAATCCGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCAAATCGCCGTTCCGGCGGTGAGCTCATCATCTCCCAGCATTTCGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCCCCGTGGTCTGTTGCCTCTTGGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAgtggtcgcgac                                                                                                                                                                                                                <  2:37360/248‑7 (MQ=255)                                          tCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGa                                                                                                                                                                                                           <  2:43009/250‑1 (MQ=255)                                                        cGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAg                                                                                                                                                                                           <  1:291027/251‑1 (MQ=255)                                                                     ccAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAAcagcag                                                                                                                                                                            <  1:151919/251‑1 (MQ=255)                                                                                                            gattCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCCTGATGCTCCGCCAGCGCCCCGGCGGGAGGGGTCGCCAACGGGGTCCCCTTT‑TAAGCAGGG‑CCCAAACA‑CAAATCGGACTTAACGCCCGAATCCA‑CCAAAGCTTccccgg                                                                                                                                     >  2:97570/4‑245 (MQ=255)                                                                                                            gattCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGGCATGATGCTCCGCCAGCGCCACGGCGGGGGTGGGCGCCAACGGGATCCCTTTG‑TTAGGAGGG‑GC‑CAACAGCAGAACGCACTCCATTCCAGAAACCC‑CCAATGCTTCAGCgg                                                                                                                                     >  1:434993/4‑249 (MQ=255)                                                                                                                                    gTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAggtc                                                                                                           <  1:203251/251‑4 (MQ=255)                                                                                                                                    gTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCATGATGCTCCGCCAGCGCCACGGCGGGAGGGGTCGCAATCGGGATCCCTTTG‑TAAACAGGG‑CC‑AAACAGGCGATAGGACTCAAAGCCCGAATCCA‑CCCATGCTTCAGCGTAAAAACGGACTAACaaggccagggc                                                                                                            >  1:363690/1‑241 (MQ=255)                                                                                                                                           catcatCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGGGGTCGCAATCGGGATCCCTTTGGTAACCGGGG‑CC‑AC‑CAGCAGGTCGGACTCCATTCCAGAATCCAACCATGGCT‑CAGCGTAAAAAAGGCCCAACAC‑TGCCAGGTCGGGcccc                                                                                                    >  1:253718/1‑250 (MQ=255)                                                                                                                                           catcatCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCATGATGCTCCGCCAGCGCCACGGCGGGAGTGGGCGCAATCGGGGTCCCCTTG‑TAAGGCGGG‑CC‑CAACAGCCGATCGGACACCATGCCAAAATCCA‑CCAATGCTTCAACGGAAAAACGGCCTAACAG‑GGCCCGGGCGCGCCCg                                                                                                    >  1:101675/1‑251 (MQ=255)                                                                                                                                              catcTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGt                                                                                                  <  1:291383/250‑1 (MQ=255)                                                                                                                                               atcTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGGCATGATTCTCCGCCAGCGCCACGGCGGTTGTGGTCGCAATACGGATCCCCTTG‑TAAGCAGGGGCC‑AACCA‑CCGATCGAAATCCATGCCAGAATCCA‑CCAATGCTTCAGCGGAAAAACCGCCTAACAG‑TGCCAggtccccccccggat                                                                                                >  1:292961/1‑237 (MQ=255)                                                                                                                                                   aCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGTCCTAACAG‑TGCCAGGTCGCGCCCGGTAt                                                                                                <  1:373636/247‑1 (MQ=255)                                                                                                                                                   aCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTAtt                                                                                               <  2:434993/248‑1 (MQ=255)                                                                                                                                                   aCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTAtt                                                                                               <  1:46016/248‑1 (MQ=255)                                                                                                                                                   aCCAGCATTACGCGTCCTGGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCATAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTAt                                                                                                <  1:88163/247‑1 (MQ=255)                                                                                                                                                       gCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGc                                                                                        <  1:282872/251‑1 (MQ=255)                                                                                                                                                           tACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCg                                                                                     <  1:243618/250‑1 (MQ=255)                                                                                                                                                                                accaaccaGATTGCCGCCTTCTCCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTCCGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCtt                                                                <  2:296541/250‑1 (MQ=255)                                                                                                                                                                                accaaccaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAATCAGTCCGGCAGTTCAAGGGCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGAGCCCTTTG‑TAAGCAGGG‑CC‑AAACATCAGAGCGAACTCAATGCCAGAATCCA‑CCAATGCTCCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGTCGTTGAAGAAATAGGGGCt                                                                 <  2:432589/249‑1 (MQ=255)                                                                                                                                                                                accaaccaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCt                                                                 <  1:335850/249‑1 (MQ=255)                                                                                                                                                                                     ccaGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTgcgc                                                          <  1:432612/251‑1 (MQ=255)                                                                                                                                                                                                   gtaCCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCATGGTCATTATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGAGCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGAGCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGTCTCTTGCGCCCGGATTTCAGc                                              <  2:382223/247‑1 (MQ=255)                                                                                                                                                                                                          ggatctttgcttatgggcgtttAAAGCAGTTCGGCAGGTCAAGGTAATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCATAACAG‑TGCCAGGTCGCGACCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAAcgc                                      <  2:177027/230‑3 (MQ=255)                                                                                                                                                                                                               ttggcctcttTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGcc                                   <  1:447506/242‑1 (MQ=255)                                                                                                                                                                                                                ttGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCATGATGCTCCGCCCGCGCCACGGCGGGAGTGGGCGCCATCGGGGTCCCCTTG‑TAAGCCGGG‑CC‑AAACCACCGATCGGAATCCATGGCAGAATCCA‑CCCACGCGTCCTCGTAAAAAACGCCTAACAGCTGCC‑GGCCGCGCCCGGGATTAAACAGCCCCGCGTTGCCGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGGAAAcccacccaac                               >  2:240599/1‑242 (MQ=255)                                                                                                                                                                                                                            ggTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgc                   <  1:327145/251‑1 (MQ=255)                                                                                                                                                                                                                            ggTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgc                   <  1:367969/251‑1 (MQ=255)                                                                                                                                                                                                                               tAAAGCAGTACGGCCGGTCAATGTCAGGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgctt                 <  2:26589/250‑1 (MQ=255)                                                                                                                                                                                                                               tAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCt                                                                 <  2:424512/202‑1 (MQ=255)                                                                                                                                                                                                                               tAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgctt                 <  1:223825/250‑1 (MQ=255)                                                                                                                                                                                                                               tAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGGCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAAAAG‑TGCCAGGGCGCGCCCGGGATTAAACAGCCCGGCGGTGAAGAAATAGGGGCt                                                                 >  1:424512/1‑202 (MQ=255)                                                                                                                                                                                                                                 aaGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgcttgct              <  1:423298/251‑1 (MQ=255)                                                                                                                                                                                                                                           gCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAt    <  1:97570/251‑1 (MQ=255)                                                                                                                                                                                                                                             aggtcaaggtcaTGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCACTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTc  <  1:240599/251‑1 (MQ=255)                                                                                                                                                                                                                                              ggtcaaggtcATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGCAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTc  >  1:436868/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     GCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCGGAATCCA‑CCAACGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC  >  NC_000913/3815255‑3815741

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