Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,420,539	C→T	L59L (CTG→TTG)	ydaS →	Rac prophage; toxin YdaS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,420,539	0	C	T	98.6%	245.4 / ‑4.9	73	L59L (CTG→TTG)	ydaS	Rac prophage; toxin YdaS
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (33/38);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGTTGTCGGTGGGCAATCAGCAATGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTCTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCATTTGATCATACCTGAAACATCAAGAG  >  NC_000913/1420403‑1420672                                                                                                                                         |                                                                                                                                      tgttgtCGGTGGGCAATCAGCAATGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTtgt                                                                                                                                    <  1:37915/140‑1 (MQ=255)   ttgtCGGTGGGCAATCAGCAATGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGAGAGGTTTTGTGc                                                                                                                                  >  1:419409/1‑140 (MQ=255)          tGGGCAATCAGCAATGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGTAACAAGAGGTGAGGTTTTGTGCGAAGAAc                                                                                                                           >  1:447025/1‑140 (MQ=255)            ggCAATCAGCAATGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACtt                                                                                                                         <  1:670763/140‑1 (MQ=255)             gCAATCAGCAATGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTc                                                                                                                        >  1:493380/1‑140 (MQ=255)               aaTCAGCAATGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGt                                                                                                                      >  1:534702/1‑140 (MQ=255)                aTCAGCAATGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACtt                                                                                                                         >  2:750662/1‑136 (MQ=255)                aTCAGCAATGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACtt                                                                                                                         <  1:750662/136‑1 (MQ=255)                aTCAGCAATGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTc                                                                                                                     <  2:627311/140‑1 (MQ=255)                aTCAGCAATGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTc                                                                                                                     <  1:290744/140‑1 (MQ=255)                        tGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATAtt                                                                                                             >  2:110418/1‑140 (MQ=255)                          gCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGTAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGa                                                                                                           <  2:447025/140‑1 (MQ=255)                          gCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTc                                                                                                                     <  1:227920/130‑1 (MQ=255)                          gCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTc                                                                                                                     >  2:227920/1‑130 (MQ=255)                          gCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGa                                                                                                           <  1:616552/140‑1 (MQ=255)                           cTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGAc                                                                                                          >  2:279691/1‑140 (MQ=255)                             aGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTg                                                                                                        >  2:644824/1‑140 (MQ=255)                             aGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTg                                                                                                        >  2:95786/1‑140 (MQ=255)                               gCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGt                                                                                                      <  1:470267/140‑1 (MQ=255)                               gCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGt                                                                                                      <  2:703788/140‑1 (MQ=255)                                 ttttAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCa                                                                                                    <  1:95786/140‑1 (MQ=255)                                 ttttAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCa                                                                                                    <  1:644824/140‑1 (MQ=255)                                   ttAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCata                                                                                                  >  1:455969/1‑140 (MQ=255)                                    tAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGAc                                                                                                          >  2:566582/1‑131 (MQ=255)                                    tAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGAc                                                                                                          <  1:566582/131‑1 (MQ=255)                                    tAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCatat                                                                                                 >  2:494271/1‑140 (MQ=255)                                    tAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCatat                                                                                                 >  2:589861/1‑140 (MQ=255)                                    tAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTATGTGCGAAGAACTTCGTCCTGATATTGACTGGTCatat                                                                                                 >  1:162640/1‑140 (MQ=255)                                             cACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGt                                                                                        <  2:586755/140‑1 (MQ=255)                                                  ccAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGc                                                                                   <  1:110418/140‑1 (MQ=255)                                                  ccAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGc                                                                                   <  1:494271/140‑1 (MQ=255)                                                  ccAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGc                                                                                   >  1:409938/1‑140 (MQ=255)                                                    aaGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCAt                                                                                 >  1:52539/1‑140 (MQ=255)                                                    aaGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGGCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCAt                                                                                 >  2:607339/1‑140 (MQ=255)                                                      gCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCAtgt                                                                               >  2:43022/1‑140 (MQ=255)                                                          aaaTCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATgttgtt                                                                           >  1:535913/1‑140 (MQ=255)                                                            aTCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCg                                                                         <  1:589861/140‑1 (MQ=255)                                                            aTCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCg                                                                         <  1:665668/140‑1 (MQ=255)                                                            aTCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCg                                                                         <  1:211031/140‑1 (MQ=255)                                                            aTCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCg                                                                         <  1:421222/140‑1 (MQ=255)                                                            aTCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCg                                                                         <  2:455969/140‑1 (MQ=255)                                                            aTCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCg                                                                         <  1:43022/140‑1 (MQ=255)                                                            aTCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTAAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCg                                                                         <  2:419409/140‑1 (MQ=255)                                                                 tGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCgaagaa                                                                    >  2:739794/1‑140 (MQ=255)                                                                 tGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGaa                                                                    <  1:106452/140‑1 (MQ=255)                                                                 tGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGaa                                                                    <  2:534702/140‑1 (MQ=255)                                                                 tGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTGCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGaa                                                                    >  2:701023/1‑140 (MQ=255)                                                                 tGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGGCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGaa                                                                    <  1:607339/140‑1 (MQ=255)                                                                    aTCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAatat                                                                 <  2:409938/140‑1 (MQ=255)                                                                       aaaGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTc                                                              >  2:661824/1‑140 (MQ=255)                                                                             gTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTCTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGaa                                                        >  2:232811/1‑140 (MQ=255)                                                                              tACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAg                                                       >  2:443230/1‑140 (MQ=255)                                                                               aCGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGc                                                      <  1:452379/140‑1 (MQ=255)                                                                               aCGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGc                                                      <  2:535913/140‑1 (MQ=255)                                                                                 gTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCaa                                                    >  1:292322/1‑140 (MQ=255)                                                                                     aTTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAAATTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTaa                                                >  2:550190/1‑140 (MQ=255)                                                                                      ttGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGTAGCAACTaaa                                               >  2:494426/1‑140 (MQ=255)                                                                                      ttGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTaaa                                               >  1:259798/1‑140 (MQ=255)                                                                                           tgccGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACa                                          <  1:359531/140‑1 (MQ=255)                                                                                             ccGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGt                                        >  2:743387/1‑140 (MQ=255)                                                                                             ccGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGt                                        >  1:431146/1‑140 (MQ=255)                                                                                                     tGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATc                                <  2:259798/140‑1 (MQ=255)                                                                                                       tCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCt                              >  2:41358/1‑140 (MQ=255)                                                                                                           gCAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCAtt                          <  1:743387/140‑1 (MQ=255)                                                                                                            cAATTGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCAttt                         >  2:100660/1‑140 (MQ=255)                                                                                                               ttGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCATTTGAt                      <  2:493380/140‑1 (MQ=255)                                                                                                               ttGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCATTTGAt                      >  1:146194/1‑140 (MQ=255)                                                                                                                tGAACGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCATTTGATc                     >  1:58987/1‑140 (MQ=255)                                                                                                                    cGTGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCATTTGATCATAc                 <  1:739794/140‑1 (MQ=255)                                                                                                                      tGCAACAAGAGGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCATTTGATCATACCt               <  2:431146/140‑1 (MQ=255)                                                                                                                             agagGTGAGGTTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCATTTGATCATACCTGAAACAt        <  2:328407/140‑1 (MQ=255)                                                                                                                               aggtgaggtTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCATTTGATCATACCTGAAACATCa      <  1:701023/140‑1 (MQ=255)                                                                                                                                   gaggtTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCATTTGATCATACCTGAAACATCAagag  >  1:314312/1‑140 (MQ=255)                                                                                                                                   gaggtTTTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCATTTGATCATACCTGAAACATCAagag  <  1:224478/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      TGTTGTCGGTGGGCAATCAGCAATGGCTAGGCTTTTAGGTGTATCACCTCCAAGCGTAAATCAATGGATCAAAGGGGTACGTCAATTGCCTGCCGAGAGATGTCCAGCAATTGAACGTGCAACAAGAGGTGAGGTTCTGTGCGAAGAACTTCGTCCTGATATTGACTGGTCATATTTACGACGTTCGGCATGTTGTTCGCAGAATATGTCAGTGAAGCAACTAAATGACAGTAACAAATCCTCATTTGATCATACCTGAAACATCAAGAG  >  NC_000913/1420403‑1420672

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