Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,070,814	G→T	R315S (CGT→AGT)	yihS ←	sulphoquinovose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,070,814	0	G	T	96.6%	208.6 / NA	59	R315S (CGT→AGT)	yihS	sulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (33/24);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACGTACACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCAGTGGCCTGGTGTACCTCCGAACGCGCGGAAGCGATGCGCCGGGTTATCTTTGTTGTAATCCGGCAGCGGATTCCACTGG  >  NC_000913/4070536‑4071091                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                       cagACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTAATTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTc                                                                                                                                                                                                                                                                               >  1:755552/1‑288 (MQ=255)   gACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGTCCGTACTCATTGCTTCGACGATAGGCCAACTTAAACGTTcgcg                                                                                                                                                                                                                                                                            >  1:386696/1‑289 (MQ=255)          tGCAACCGCGTGATCCAGGCCGGGCGCTAACGGGATACGCTGGATCACCAGGACATGCAGCAGGTTATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGaccacc                                                                                                                                                                                                                                                                     <  2:380805/289‑1 (MQ=255)              aCCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGtt                                                                                                                                                                                                                                                                 >  2:302849/1‑289 (MQ=255)                      aGCCATGCCGGGCGCTAACGGGATACGCGGGGTCACCAGGCAATGCAGCAGGTGATAAATACCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTttccttc                                                                                                                                                                                                                                                          <  1:59436/288‑1 (MQ=255)                           tGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTc                                                                                                                                                                                                                                                    >  1:541866/1‑289 (MQ=255)                            gCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCTCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGACCCACCGGTTTTCCTTCCCAGTCa                                                                                                                                                                                                                                                   >  1:502840/1‑289 (MQ=255)                            gCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCa                                                                                                                                                                                                                                                   <  1:655083/289‑1 (MQ=255)                            gCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCa                                                                                                                                                                                                                                                   <  1:302849/289‑1 (MQ=255)                            gCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTCGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCCGTTTTCCTTCCCAGTCa                                                                                                                                                                                                                                                   >  1:430102/1‑289 (MQ=255)                              cGGGCGCGAACGGGATACGCGGGATCACCAGGAAAGGCAGCAGGTGATAAATATCCGGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCaa                                                                                                                                                                                                                                                  <  1:800295/288‑1 (MQ=255)                                     tAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGtata                                                                                                                                                                                                                                           <  1:472922/288‑1 (MQ=255)                                                      tCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGCCCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCgcgc                                                                                                                                                                                                                          >  1:81234/1‑288 (MQ=255)                                                              cAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGt                                                                                                                                                                                                                  >  1:155910/1‑288 (MQ=255)                                                                           gATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATACGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGCCAACGGTATAAACAATCCCCTCCGCGCCATCGGGTGCCCAGGCAgcgc                                                                                                                                                                                                     >  1:474125/1‑288 (MQ=255)                                                                                atatCCTGTTTGCCGGCCCAGACTTTGGTGGTGACCTTAGTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACgg                                                                                                                                                                                               <  1:557926/289‑1 (MQ=255)                                                                                   tCCTGTTTGCCGTCCAAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTATACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACACTACAGTCCGACCCATCTGGTGCCCCTGCATCGAGCACGGttg                                                                                                                                                                                            >  1:280540/1‑287 (MQ=255)                                                                                     aacagacaggTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCATGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCg                                                                                                                                                                                          <  1:108402/280‑1 (MQ=255)                                                                                          tGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTTCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTCTAACCAATCCCAACCGCCCACTAGGGTCCCCCGGCCTCCCGCCCtgtgagctcaca                                                                                                                                                                                     >  2:635741/1‑281 (MQ=255)                                                                                          tGCCGTACGAGACTTTTGTGGTGACCTTATTGGCCGCATCCACGTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTaaa                                                                                                                                                                                     <  2:81234/289‑1 (MQ=255)                                                                                             cGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAg                                                                                                                                                                                  <  1:635741/289‑1 (MQ=255)                                                                                                     tctgtggttgTGACCTGATGGTCAGCATCCAGCTCCTGCCACCAGGACCCATTTTCATAGTCCATCAGGAACTGAATGCAGTACTCCCACTATGGTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTgg                                                                                                                                                                          <  1:779638/285‑1 (MQ=255)                                                                                                         tggtggTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTCCCCATTGCTTCGACGATCGGCCACCTTCAACGTTCGCGTCCCACCGGTTTTCCTTCCACGTCCACGCCAAACACAATCCCCTCCGCCCCATCCGGTGTCCAGGCATCGCGCACCGTCCCGCTAACCAGACCtttagaa                                                                                                                                                                        >  2:141599/1‑283 (MQ=255)                                                                                                         tggtggTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCCCCATCGGGTGCCCCGGCATCGCGCCCGGTGGCCTTacaccgcc                                                                                                                                                                                >  2:108402/1‑275 (MQ=255)                                                                                                               tgACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTa                                                                                                                                                                 >  1:697495/1‑288 (MQ=255)                                                                                                                      aTTGTCCGCATCCAGCTCCTTGCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATCCCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCAt                                                                                                                                                         <  1:148628/289‑1 (MQ=255)                                                                                                                            cGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCCGCCATGCtggt                                                                                                                                                   >  1:665873/1‑289 (MQ=255)                                                                                                                            cGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGCTCGCGAACTAAGGGTTTTCCTTCCCAGTCAACGGTCTACACAATCCCCTCCGCGCCACCGGCTGCCCATGCATCGCGCACGTCCGAGTTCAGCAGACCCTTCACACCATCTAGCATCCATGCtggt                                                                                                                                                   >  2:658265/1‑289 (MQ=255)                                                                                                                                           ccaccaGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCCCCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACgg                                                                                                                                    <  1:469431/289‑1 (MQ=255)                                                                                                                                                     cctttttCTTAGTGCATCTGGTACTTAATGCAGTCCTCCCGCCCGGTTTGATACCAGGTTTCATACTGGCGGTCGCCAGTGACGGTGATGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACGTAAACGTTCGCGGACCACCGGTTTTCCGTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAggg                                                                                                                          <  2:356716/286‑1 (MQ=255)                                                                                                                                                                gTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCCTGGAtg                                                                                                               >  1:168226/1‑289 (MQ=255)                                                                                                                                                                   aatgaGGTACTTAATGCAGTACTCCCACCATGTGTTATACGAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTGCGCGGACCAACGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGTATGTg‑c                                                                                                           <  1:141599/285‑1 (MQ=255)                                                                                                                                                                                aggaggTACTCCACCAATGTTGGATCCCGGGTGACGTACTGGCGCTGGCAGGTGAGGGTGTGGAGCGCGTAGGCCGGACCCCTTGGTCCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGc                                                                                              <  2:785298/284‑1 (MQ=255)                                                                                                                                                                                      tACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCTCGTAGGTCGCAGCAACAGCATAAAAGATAAGCCAACTTAAACGTTCGCGGACCACCGGTTGTCCTTCACAGTCAAAGGTATAAACAATCCCGTCCGCGCCATCGGGGCCCCAGGCATCGCGCTCGGCTGCGTTAAACAGCCCTTGGCCATCTTCTAGCCGCCATTCTCGTGGTTGTTCCCAACGGCCTTCCCGCGCAGATTGGGTCTC‑TCGCATTAAACGGtccact                                                                                         >  2:439273/1‑285 (MQ=255)                                                                                                                                                                                      tACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCCTCGCGCACGGTGGCGTTAAACAGCCCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCAAGGGCAGCCTGGA‑GTGTCAGCATTACACGGCCCAAt                                                                                         >  1:156250/1‑288 (MQ=255)                                                                                                                                                                                      tACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCAt                                                                                         >  1:152743/1‑288 (MQ=255)                                                                                                                                                                                      tACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCAt                                                                                         >  1:32175/1‑288 (MQ=255)                                                                                                                                                                                             acccggttcgctACCCGGGTTCATGCTGGCGATCGCCAGTGACGGTGGCGAGCCCGTAGGCCGGACCCATTGCTTCGACGGTAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTGCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCa                                                                                 <  2:541866/278‑1 (MQ=255)                                                                                                                                                                                             accaTGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCaa                                                                                                                                                                                                                                                  <  2:709239/129‑1 (MQ=255)                                                                                                                                                                                             accaTGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCaa                                                                                                                                                                                                                                                  >  1:709239/1‑129 (MQ=255)                                                                                                                                                                                             accaTGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGgcag                                                                                                                       <  2:142876/252‑1 (MQ=255)                                                                                                                                                                                             accaTGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGgcag                                                                                                                       >  1:142876/1‑252 (MQ=255)                                                                                                                                                                                             accaTGTTGGATTCCAGGTTTCATACTGGGTAGCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATCGGCCAACTTCAACGTTCGCGGACCACCGGTTATCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCa                                                                                 <  1:783106/289‑1 (MQ=255)                                                                                                                                                                                                aTGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCAGTg                                                                              >  1:409838/1‑289 (MQ=255)                                                                                                                                                                                                  gTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGACCAGGCATCGCGCCCGGTGGCGTTAAACAGCCCTTTGGCATCTCCTAGCAGCCATGCGGGTGGTTGTTCGCAACGGGCTTCCAGCTCAGCATGGCTGGT‑CAATATTACACGGCCCCGTTCACTCCACTgt                                                                             >  2:768397/1‑287 (MQ=255)                                                                                                                                                                                                     tGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCACGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGATGGTGGTTGTTCGCAACGGTCTTCCAGGGCAGCATGGATGTG‑CAGCATTAACAGGCCACATTCGATCCAGTGGCCTg                                                                         >  1:239629/1‑289 (MQ=255)                                                                                                                                                                                                             ggTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCAGTGGCCTGGTGTAcc                                                                  >  1:31430/1‑288 (MQ=255)                                                                                                                                                                                                                                  acggacgtgtgtgagCGTGGAGGCCGTACCACTTGCCTCGACTATAGGCCAACTTAACCGTTCCCGGACCACCGGTTTTCCGTCCCAGTCAACGGTATAAACAAGCCCGGCCGCGCCTTCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCAGCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAATGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCAGTGGCCTGGTGTACCTCCGAGCGCGCGGAAGCGATgc                                            <  2:755552/277‑1 (MQ=255)                                                                                                                                                                                                                                   gTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCACCATTCGATCCAGTGGCCTGGTGTACCTCCGAACGCGCGGAAGCGATgc                                            >  1:608321/1‑288 (MQ=255)                                                                                                                                                                                                                                                  gcgTAGGCCGTACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCAGTGGCCTGGTGTACCTCCGAACGCGCGGAAGCGATGCGCCGGGGTCACTttgt                            >  1:360032/1‑289 (MQ=255)                                                                                                                                                                                                                                                            tACCCATTGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCATGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCCATCCAGTGGCCTGGTGTACCTCCGAACGCGCGACAGCGATGCCCCGTGATCTCtttttttcactccc                   >  2:751631/1‑280 (MQ=255)                                                                                                                                                                                                                                                                  ttGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGt                                                                                                                                                                                                                                                     >  1:46630/1‑57 (MQ=255)                                                                                                                                                                                                                                                                  ttGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGt                                                                                                                                                                                                                                                     <  2:46630/57‑1 (MQ=255)                                                                                                                                                                                                                                                                  ttGCTTCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCAGTGGCCTGGTGTACCTCCGAACGCGCGGAATCGAGTCGCCGGGTGAGCTTTGTTGTAATCCGGccgggg            >  2:778349/1‑284 (MQ=255)                                                                                                                                                                                                                                                                       tCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCAGTGGCCTGGTGTACCTCCGAACGCGCGGAAGCGATgcg                                           >  1:583970/1‑253 (MQ=255)                                                                                                                                                                                                                                                                       tCGACGATAGGCCAACTTAAACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCAGTGGCCTGGTGTACCTCCGAACGCGCGGAAGCGATgcg                                           <  2:583970/253‑1 (MQ=255)                                                                                                                                                                                                                                                                        cgacgaGAGGCCAACTTAAACGTTCGCGGACCCCCGGTTTTCCTTCCCCGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCAGTGGCCTGGTGTACCTCCGAACGCGCGGAAGCGATGCGCCGGGTTATCTTTGTTGTAATCCGGCAGCGGATTCCa      <  2:697495/289‑1 (MQ=255)                                                                                                                                                                                                                                                                        cgacgaGAGGCCAACTTAAACGTGCACGGAACACCGGTTTTCCATCCCAGTCAACGGTATAAACAATCCCGACCGCGCCATCGGGTGCCCAGTCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCGAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCAGTGGCCTGGTGTACCTCCGAACGCGCGGAAGCGATGCGCCGGGTTATCTTTGTTGTAATCCGGCAGCGGATTCCa      <  2:360032/289‑1 (MQ=255)                                                                                                                                                                                                                                                                            taTAGGCCAACGTACACGTTCGCGCACCCCCCGTTTTCCTTCCCAGTAAACAGGATAAACAATCCCTGCCGCGCCATTGGGTGCCCAGGCAGCGCGCCCGGTGGCGTTAAACAGACCTTTGGCAGCTTCTATCAGCCATGCTGGTGGTTGTGCGCAACGGGCTTCCAGGGCAGCATGGATGGG‑CAGCATTAAACGGCCCCATTCGATCCAGTGGCCTGGTGTACCTCCGAACGCGCGGAAGCGATGCGCCGGGTTATCTTTGTTGTAATCCGGCAGCGGATTCCACTgg  <  2:608321/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                       CAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTACTTAATGCAGTACTCCCACCATGTTTGATACCAGGTTTCATACTGGCGATCGCCAGTGACGGTGTAGAGCGCGTAGGCCGTACCCATTGCTTCGACGATAGGCCAACGTACACGTTCGCGGACCACCGGTTTTCCTTCCCAGTCAACGGTATAAACAATCCCGTCCGCGCCATCGGGTGCCCAGGCATCGCGCACGGTGGCGTTAAACAGACCTTTGGCATCTTCTAGCAGCCATGCTGGTGGTTGTTCGCAACGGGCTTCCAGGGCAGCATGGATGTG‑CAGCATTAAACGGCCCCATTCGATCCAGTGGCCTGGTGTACCTCCGAACGCGCGGAAGCGATGCGCCGGGTTATCTTTGTTGTAATCCGGCAGCGGATTCCACTGG  >  NC_000913/4070536‑4071091

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑