Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,142,163	G→A	S18L (TCG→TTG)	dcd ←	dCTP deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,142,163	0	G	A	97.0%	212.8 / ‑5.8	66	S18L (TCG→TTG)	dcd	dCTP deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/46);  minor base C (2/0);  total (20/46)
Fisher's exact test for biased strand distribution p-value = 8.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

AGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCGACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGATGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTA  >  NC_000913/2141920‑2142400                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              ggCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACg                                                                                                                                                                                                                                        <  1:420254/250‑1 (MQ=255)    ggcACCGACTCCGGCGGCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGcc                                                                                                                                                                                                                                     <  2:161048/249‑1 (MQ=255)          gACTCCAGCGTCCCCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTtcatc                                                                                                                                                                                                                               >  2:437234/1‑251 (MQ=255)                aGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCGCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCa                                                                                                                                                                                                                         <  1:264323/251‑1 (MQ=255)                  cGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAgg                                                                                                                                                                                                                       <  1:129366/251‑1 (MQ=255)                         gccagcgccagcTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGTGAACGACCTCATCGCTCATCAAGCGGTCAAGAGCGGCGCTCACTTCATCTTTGGGCCCGCTCAGATCGCTAAGCGCTGCCGTGTGACCCCGGAAGGTACGAAATTTATTGCCCAGCCGTACCTCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGGTGAACAACAAACGGCCCtcaaccaccccgggttccat                                                                                                                                                                                                                >  2:8385/1‑234 (MQ=255)                               gccagcTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑GGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCg                                                                                                                                                                                                          <  2:32220/251‑1 (MQ=255)                                ccagcTCTCCTGGGTGAAGATAGAATGCCTCGCCCTCGTCTAGAACGAGCTCAGCGCTCATCCCGCGGTCAAGGTCGGCGCTCTCTTCATCTTGGGTCCCGCTCAGAGCGATAAACGGTGCCGTGTGCGCACGGAAGGTACGAAATTTATTGCCCAGGCGTCCAGCCACCGGCGCGCCGTTAATACGCTGCACTGGCGGACG‑TTGGTTGGTCAACAAACGGCCAGCATCTAGGGAGGCTTCAATCTCTCgg                                                                                                                                                                                                         <  2:178850/251‑1 (MQ=255)                                          tGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGAc                                                                                                                                                                                               >  1:394803/1‑251 (MQ=255)                                                aaGATAGAACGCCTCGCCCTCGTCGAGAACGAGCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAAGATCTCGGTCACACATACGCATgg                                                                                                                                                                                         <  2:371280/251‑1 (MQ=255)                                                 aGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGc                                                                                                                                                                                        <  1:200762/251‑1 (MQ=255)                                                       aaCGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTtc                                                                                                                                                                                  <  1:163438/251‑1 (MQ=255)                                                        aCGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTtct                                                                                                                                                                                 <  1:437234/251‑1 (MQ=255)                                                           ccTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCt                                                                                                                                                                              >  2:328163/1‑251 (MQ=255)                                                                   ccgtcgAGAACGATCTCATCGCTAATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCa                                                                                                                                                                      <  1:69308/250‑1 (MQ=255)                                                                    cgtcgAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCgcgc                                                                                                                                                                       <  1:85493/249‑1 (MQ=255)                                                                       cgAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCAc                                                                                                                                                                  >  2:380089/1‑251 (MQ=255)                                                                       cgAGAACGAGCGCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAACCGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGGCGTTAATACTCTCCACTTGCGGACG‑TGTGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACTCATTCGCATGGCATGTCCTCCTTGCGCtcatca                                                                                                                                                                   <  2:404240/250‑1 (MQ=255)                                                                                                 cgGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAggg                                                                                                                                        >  1:243821/1‑251 (MQ=255)                                                                                                     cAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCAt                                                                                                                                     <  1:395897/250‑1 (MQ=255)                                                                                                       aGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTGAACGCCAAGTACTTCAAGGGCAt                                                                                                                                     <  1:113027/248‑1 (MQ=255)                                                                                                        gcgcGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTa                                                                                                                                 <  1:21174/251‑1 (MQ=255)                                                                                                           cggcgCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGt                                                                                                                              <  2:70531/251‑1 (MQ=255)                                                                                                           cggcgCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGGCAAGGGCATACTAGGt                                                                                                                              <  1:94547/251‑1 (MQ=255)                                                                                                                 tCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCAc                                                                                                                        >  2:430620/1‑251 (MQ=255)                                                                                                                  cACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGTTTATCAc                                                                                                                        <  1:107426/250‑1 (MQ=255)                                                                                                                       cATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAAc                                                                                                                   <  2:145581/250‑1 (MQ=255)                                                                                                                               ggACCGCTCAGATCGATAAACGCTGCCGTGTTACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCAGTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATGCTAGGTTATCACGGAACGGTACACaa                                                                                                          <  2:308549/251‑1 (MQ=255)                                                                                                                               ggACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACaa                                                                                                          <  1:69360/251‑1 (MQ=255)                                                                                                                                            cGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCtt                                                                                              >  2:363098/1‑250 (MQ=255)                                                                                                                                            cGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCtt                                                                                              <  1:147407/250‑1 (MQ=255)                                                                                                                                             gATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTAt                                                                                            >  1:305949/1‑251 (MQ=255)                                                                                                                                              aTAAAAGCTGCCGTGTGACCAAGGAAGGTACGACATTTATTGACCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCACGCTTAtt                                                                                           >  2:81249/1‑251 (MQ=255)                                                                                                                                                aaaCGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTc                                                                                          <  1:44783/250‑1 (MQ=255)                                                                                                                                                aaaCGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTc                                                                                          >  2:160447/1‑250 (MQ=255)                                                                                                                                                 aaCGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCa                                                                                         >  1:115844/1‑250 (MQ=255)                                                                                                                                                   cGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAg                                                                                      <  1:285384/251‑1 (MQ=255)                                                                                                                                                       gCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAAc                                                                                   <  1:380089/250‑1 (MQ=255)                                                                                                                                                          gtgtGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGAc                                                                               <  1:413693/251‑1 (MQ=255)                                                                                                                                                              gACCACGGAAGGTACGAAATTTATTGCCCGTGCGTACATCCACCGGCGCGCCGTTAATACGCTCCACTGGCGGACT‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTa                                                                            <  1:273108/250‑1 (MQ=255)                                                                                                                                                                  aCGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTc                                                                       <  1:105145/251‑1 (MQ=255)                                                                                                                                                                       aGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCttt                                                                  <  1:51087/251‑1 (MQ=255)                                                                                                                                                                         gtgcGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTc                                                                                          >  1:24196/4‑225 (MQ=255)                                                                                                                                                                         gtgcGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTc                                                                                          <  2:24196/222‑1 (MQ=255)                                                                                                                                                                         gTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTc                                                                                          <  2:24174/225‑1 (MQ=255)                                                                                                                                                                         gTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTc                                                                                          >  1:24174/1‑225 (MQ=255)                                                                                                                                                                           aCGAAATTTATTGACCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAat                                                              <  2:167182/251‑1 (MQ=255)                                                                                                                                                                               aaTTTATTGCCCAGGCGGACATCCACCGTCGCTCCGTTAATTCGCTCCACTGGCGGACG‑TGGGTTGATCAACAATCGGCCTTCATCAAGCCAGTCTTCAATATCTCGGTCACACAGACGCATGGCTTATTCTCCTTTCGCTCATCACTCCCTTAACGCCATTTACGTCAAGGGCATACTAGGTTATCACGGAACGGTACAGAATTTATCAAGCGTATTCAAAGAACTGACGTATGTTCGCTTTCAatatat                                                          <  2:202264/251‑1 (MQ=255)                                                                                                                                                                                 tttATTGCCAAGGCGTACATCCACCGTCGAGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCCACAAACGGCCTTCATCAAGCCCGGCTTCAATTTCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCCTACTAGGTTATCACTGAACGGAACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCCATCTatcg                                                        >  2:361211/1‑251 (MQ=255)                                                                                                                                                                                   tATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAAGATatcga                                                       <  1:312383/250‑1 (MQ=255)                                                                                                                                                                                       gCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATcgcg                                                  <  1:152982/251‑1 (MQ=255)                                                                                                                                                                                             gCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGATGCgg                                            <  2:243821/251‑1 (MQ=255)                                                                                                                                                                                                       aCCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGCTCCACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCaa                                                               >  1:336302/1‑222 (MQ=255)                                                                                                                                                                                                       aCCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCaa                                                               <  2:336302/222‑1 (MQ=255)                                                                                                                                                                                                            cgcgCCGTTAATACGCTCCACTGGCGGACG‑TGGGTGGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGATGCGGTTTTTCCCGCcgcgc                             <  1:421655/251‑1 (MQ=255)                                                                                                                                                                                                              cgcCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGATGCGGTTTTTCCCGCCGCGCgg                           <  2:337764/251‑1 (MQ=255)                                                                                                                                                                                                              cgcCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGATGCGGTTTTTCCCGCCGCGCgg                           <  1:430620/251‑1 (MQ=255)                                                                                                                                                                                                                   ttAATACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACCCAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGATGCGGTATTTCCCGCCGCGCGGCACg                       >  1:57876/1‑250 (MQ=255)                                                                                                                                                                                                                     aatccGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCAGGGCATTTTCTCCTTGCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTAGCACTGAACGGTACACAAGTGATCAAGCTTATGCAAGGAACTGACTTATGGTCGCTTTCAATATATCGATCGCGATGCGTTTGTTCGCGCCGCGCGGCTCGATa                    <  2:33869/247‑1 (MQ=255)                                                                                                                                                                                                                     aaTACGCTCCACTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGATGCGGATTTTCCCGCCGCGCGGCACGAt                     >  2:235229/1‑250 (MQ=255)                                                                                                                                                                                                                       aaCGCTCCCCTTGCGGACGGTGGGT‑GATCACCAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGATGCGTTTTTTCCCGCCGCGCGGCACGATAAt                  <  2:56569/250‑1 (MQ=255)                                                                                                                                                                                                                              cACTGGCGGCCG‑AGGGTTGATAAACAAACGGCCTTCGTCAAGCCAGGCTTCAATATCTCGGTCAAACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCAAATTACGTTAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGGTGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGc            >  2:435739/1‑250 (MQ=255)                                                                                                                                                                                                                               aCTGGCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGATGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCat          <  1:409630/251‑1 (MQ=255)                                                                                                                                                                                                                                  ggCGGACG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGATGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATAtt       <  1:127946/251‑1 (MQ=255)                                                                                                                                                                                                                                       aCG‑TGGGTTGATCAACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGATGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTa  <  1:288143/251‑1 (MQ=255)                                                                                                                                                                                                                                       aCG‑TGGGCTAATCAACAAAAGGCGTTCATCAAGCCAGGCTTAACTATTCCGGTTACACAGACGCATGGCTTTTTCTCCTTTCGCGCATCACTCCCTTAACGGCAATTACGTCAAGGGCATACTAGGTTATCACTGCACGGTCCACAATTTATCAAGCTTATTCAAATAACTGACTTATTATAGCTTTCCCTATATCGATCGCGCTGCGGTTTTTCCCGCAGCGCGGCGCGATAATGTCCGCATATGGttg   >  2:423578/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              AGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCTCAGATCGATAAACGCTGCCGTGTGACCACGGAAGGTACGAAATTTATTGCCCAGGCGTACATCCACCGTCGCGCCGTTAATACGCTCCACTGGCGGACG‑TGGGTTGATCGACAAACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAAGCTTATTCAAAGAACTGACTTATTTTCGCTTTCAATATATCGATCGCGATGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTA  >  NC_000913/2141920‑2142400

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