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

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

ATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCGGAATCCA‑CCAACGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC  >  NC_000913/3815250‑3815741                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           atcccGCTTAACTTCCTTAATCGCCGAAATCTCGCCGCGCCCCCGGTCCTGACTATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGTCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCTGTTAAAGCAGTACGGcaggtcaaggt                                                                                                                                                                                                                                                        <  1:35227/247‑1 (MQ=255)      ccTCAACTTCCTGAACCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGa                                                                                                                                                                                                                                                   <  1:423237/250‑1 (MQ=255)              tCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGGCATGATGCTcccc                                                                                                                                                                                                                                           >  1:422325/1‑249 (MQ=255)                  gAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAgc                                                                                                                                                                                                                                       <  2:146882/251‑2 (MQ=255)                   aaTCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGAGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGcagggccagggtatggtggcccgcccccg                                                                                                                                                                                                                                      >  2:375105/1‑232 (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:367969/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCAGAATCCA‑CCAATGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCtgc                   <  1:327145/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)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           ATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTG‑TAAGCAGGG‑CC‑AAACAGCAGATCGAACTCAATGCCGGAATCCA‑CCAACGCTTCAGCGTAAAAACGGCCTAACAG‑TGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC  >  NC_000913/3815250‑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: ‑