Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,589	G→A	S60S (TCC→TCT)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,589	0	G	A	95.8%	149.6 / NA	48	S60S (TCC→TCT)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/33);  minor base C (2/0);  total (15/33)
Fisher's exact test for biased strand distribution p-value = 9.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCGGAATCCA‑CCAACGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAA  >  NC_000913/3815362‑3815833                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       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)                                                                                                                                                aTGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGGACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGGAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCt                                                                                >  1:230880/1‑251 (MQ=255)                                                                                                                                                       gCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTcg                                                                                                                                                                                                      >  1:237594/1‑126 (MQ=255)                                                                                                                                                       gCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTcg                                                                                                                                                                                                      <  2:237594/126‑1 (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)                                                                                                                                                                     ggtagtggtcGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCAGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGCTATGGTTTCATGCCTTCGc                                                           >  1:403232/5‑251 (MQ=255)                                                                                                                                                                      gtagtggtcGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGct                                                          <  1:322680/248‑1 (MQ=255)                                                                                                                                                                              cgCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATc                                                   <  1:240550/249‑1 (MQ=255)                                                                                                                                                                              cgCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACGCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCTTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAATCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCatcatc                                                   <  2:138404/249‑1 (MQ=255)                                                                                                                                                                                        aTCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCt                                         >  1:382114/1‑250 (MQ=255)                                                                                                                                                                                                   aaGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAAGCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCACTCCTCATCTTACTTTTCTACAGACaaaa                               <  2:292961/250‑1 (MQ=255)                                                                                                                                                                                                       aGGG‑CC‑AAACAGCAGATCGAACTCAATGCCANAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAg                          <  2:230880/251‑1 (MQ=255)                                                                                                                                                                                                       aGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTAAGCGTAAAAACGGCCTACCAG‑TGCCAGGTCGCCCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCCGCGTAAACTCGCCAAACTTTAAAACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATCTGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGCCaaaaaaa                            >  2:365291/1‑249 (MQ=255)                                                                                                                                                                                                                           gAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCt        <  2:298604/251‑1 (MQ=255)                                                                                                                                                                                                                                cAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGAATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGGTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaa    >  2:102178/1‑250 (MQ=255)                                                                                                                                                                                                                                 aaTGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaaa  <  2:363690/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       GACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCGGAATCCA‑CCAACGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAA  >  NC_000913/3815362‑3815833

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