Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,529	T→A	T80T (ACA→ACT)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,529	0	T	A	97.6%	131.7 / NA	42	T80T (ACA→ACT)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base A (9/31);  minor base G (1/0);  total (11/31)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCGGAATCCA‑CCAACGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTT  >  NC_000913/3815287‑3815773                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          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                                                                                                                               <  1:46016/248‑1 (MQ=255)                                                                                                                   aCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTAtt                                                                                                                               <  2:434993/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)                                                                                                                                                                                                                           aTGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGGACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGGAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCt                    >  1:230880/1‑251 (MQ=255)                                                                                                                                                                                                                                  gCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTcg                                                                                                                                          <  2:237594/126‑1 (MQ=255)                                                                                                                                                                                                                                  gCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTcg                                                                                                                                          >  1:237594/1‑126 (MQ=255)                                                                                                                                                                                                                                         ccACGGCGGTAGTGGTCGCAATCGAGATCCCTTGG‑TCAGCAGGG‑CC‑AAACAGCAGATCGACCTCAATGCCAGAATCCC‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAGTCAATAAACTGGCGCTGATATGGTTTCagg      <  2:56136/251‑3 (MQ=255)                                                                                                                                                                                                                                             ggCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCtt  <  2:422325/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          CGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCGGAATCCA‑CCAACGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTT  >  NC_000913/3815287‑3815773

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