Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,065,579	T→A	16.0%	T63S (ACC→TCC)	yihO ←	putative sulphoquinovose importer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,065,579	0	T	A	16.0%	124.9 / 4.7	68	T63S (ACC→TCC)	yihO	putative sulphoquinovose importer
Reads supporting (aligned to +/- strand):  ref base T (43/9);  new base A (2/9);  total (49/19)
Fisher's exact test for biased strand distribution p-value = 8.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTAT  >  NC_000913/4065356‑4065802                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                  tACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑GGTAAAAAACCGGTGAT‑CATATCGTT‑AAACGCC‑GGaaaaa                                                                                                                                                                                                                              >  2:624630/1‑229 (MQ=255)       cAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCg                                                                                                                                                                                                                       >  2:166368/1‑229 (MQ=255)       cAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGGGAG‑CATATCGGT‑AAACGAG‑GTAAAAAA‑ATTCg                                                                                                                                                                                                                       >  2:121218/1‑229 (MQ=255)         ggtcaacaggatgttgtaaaGTCCAAACAATATAAACGGGGGGGTTGAAATCGGCGTTTTAACCGGCAAACAAAAGGTGGTGGCTATAAACTGGAGCGGTGAAATTAACTCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑AGTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGc‑c                                                                                                                                                                                                                    <  2:583262/212‑1 (MQ=255)              acaGACTTTATGAAAGCCCAAACAACAATAACGGCGGGGTGGCAACCGGCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAg                                                                                                                                                                                                               <  2:139625/228‑1 (MQ=255)                  aCTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCAAGTCG‑AGTAAAATTCCGGGGTG‑CATCTCGGT‑AAACGCG‑GTAAAAAA‑CTTAGC‑CACCATaaa                                                                                                                                                                                                            >  1:294231/1‑228 (MQ=255)                    tGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCAGTGGG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGaa                                                                                                                                                                                                         >  2:473877/1‑228 (MQ=255)                      tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑AACCAGAAAGATaa                                                                                                                                                                                                       >  1:376814/1‑229 (MQ=255)                      tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAACATTTTTCCGCGAGTCA‑GGTAAAAATCCGGTGGG‑CATATCGGT‑AAACGCG‑GTAAAACA‑CTTCGc‑cacacaaaaaaaaa                                                                                                                                                                                                       >  2:46373/1‑218 (MQ=255)                      tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCACATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CAAATCGGT‑AAACGCG‑GTAAAAAA‑CTTCG‑‑CACCCGAAAGAAAAt                                                                                                                                                                                                      >  2:173910/1‑229 (MQ=255)                      tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTAGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑AACCGGAAAGATaa                                                                                                                                                                                                       >  1:646987/1‑229 (MQ=255)                        tGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATc                                                                                                                                                                                                     >  1:496971/1‑229 (MQ=255)                            aGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGc                                                                                                                                                                                                 >  1:411094/1‑229 (MQ=255)                                cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGCTAAACTGAGGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCAAATATTTTTCCCCAGGTCG‑CGGAAAAATCCGGAGAG‑CAAAGCGTT‑AACCGAG‑TAAAAAAA‑CTTCGCACACCA‑AAAGAAAAACCCGGCaaa                                                                                                                                                                                             >  2:520249/1‑227 (MQ=255)                                cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCCATATTTTTCCGCGAGTCA‑AGTAAAAATACGGTGAG‑CATATCGGT‑AAACCCG‑GGACAAAA‑CTTCGC‑CACCAGAACGTTAATCCCGCCATa                                                                                                                                                                                             >  1:379959/1‑229 (MQ=255)                                cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGGG‑CATATCGGT‑AAACGCG‑GTCAAAAA‑CTTCGC‑CACCAGAAAAATAATCCAGCCATa                                                                                                                                                                                             >  1:208144/1‑229 (MQ=255)                                cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGcccaa                                                                                                                                                                                             >  1:462618/1‑226 (MQ=255)                                cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATa                                                                                                                                                                                             >  1:254981/1‑229 (MQ=255)                                cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATa                                                                                                                                                                                             >  1:237646/1‑229 (MQ=255)                                cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CAGATCGGT‑AAACGCG‑GAAAAAAA‑CTTCGC‑CACAAGAAAGAAAAAccccacaaa                                                                                                                                                                                             >  2:99974/1‑223 (MQ=255)                                               gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑ATTAAAAATCCGGTGAG‑CATATCGGG‑AAACTCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGAAAGAAGGCATccc                                                                                                                                                                              >  1:228056/1‑229 (MQ=255)                                               gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑AACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATccc                                                                                                                                                                              >  1:549315/1‑229 (MQ=255)                                               gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGG‑AAACGCG‑GTAAAAAA‑ATTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATccc                                                                                                                                                                              >  1:689746/1‑229 (MQ=255)                                               gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGCTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACAGGAAAGATAATCCCGCCTCAGTAAGAAGGAAcccc                                                                                                                                                                              >  1:727718/1‑229 (MQ=255)                                                cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATACCGCCATAGAAAGCAGGCATCCCt                                                                                                                                                                             >  1:74820/1‑229 (MQ=255)                                                cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGGGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CAACAGAAAGATAATCCGGCCATAGTAAGAAGGCAcccca                                                                                                                                                                             >  2:42123/1‑228 (MQ=255)                                                 gcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GCAATAAA‑CTTCGC‑CACCCGAAAGATAATCCCGCCATAGCAAGCATGCATCCCTa                                                                                                                                                                            >  1:191109/1‑229 (MQ=255)                                                        taagcgtcttTTAAACCGGAAAAAAAAAGGTGGTGGCAGAAAATTGGAGGGTTAAAATTACCGCCTCCGGAACCGTTTAATAAAAAATAAAGGGACAAA‑TTTTGCTTTTAGTCCCAATATTTTTCCGAGAGTCG‑AGTAAAAACCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAAATAATCCCGCTATAGTAAGCAGGCATCCCTAATTCATc                                                                                                                                                                     <  2:197675/225‑1 (MQ=255)                                                                        cGGAAAAAAAAAGGGGGGGGAGAAAAAATGGGGGGGTGCAATAAACCCCGCAGAAACCGCTGAAAATAAAAGAAAGAGAAAAA‑ATTTGATTTTGGGACAAATATTTTCCGGGGAGGGG‑AGAAAAAAGCTGGGGAG‑CATGCGGGA‑AAACGCG‑GTAAAAAA‑CTCCCC‑CACAAAAAAGTAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg                                                                                                                                                     <  2:599472/229‑1 (MQ=255)                                                                           caaacaaaAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTgaggga                                                                                                                                                   >  2:516409/1‑226 (MQ=255)                                                                                        gtggCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGAGGATACAGAGAACCGa                                                                                                                                     >  1:606870/1‑229 (MQ=255)                                                                                         tggCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCATAAGTAAGCAGGCATACCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGAt                                                                                                                                    >  1:399985/1‑229 (MQ=255)                                                                                                        gCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAACGAGTACCGATGCAGAGCATTAAAtt                                                                                                                     >  1:348087/1‑229 (MQ=255)                                                                                                           ttGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGATGAGATACAGAGTACCGATGCAGAGCATTAAATTCg                                                                                                                   >  1:438582/1‑228 (MQ=255)                                                                                                            tGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCACCCCTAATTCATCGGTATAAAATTTGAGGGAATACAGAGTACCCATCCAAAGAATAAAATTCaaa                                                                                                                 >  2:299121/1‑229 (MQ=255)                                                                                                              cAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATCCATCGGTATAAAATTTTAGGAAAAACAGAGGACCGATCAAAAGCATTAAATTACAAcc                                                                                                               >  2:100180/1‑229 (MQ=255)                                                                                                                   gcacgcaggagcactggttgataaAAAAAAAAGGGACAAAGATTGGCCTTAAGGCCAAATTTTTTCCCGAGAGCCG‑AAAAAAACGCGGGGGAG‑AGTA‑CGGTGAAACAGA‑CAAAAGAA‑CCACCC‑CACCAGAAACAAAACACCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGAc                                                                                                            <  2:191109/208‑1 (MQ=255)                                                                                                                   accgccgcCGGAACCGATGTATATAAAATAAAGGGACGAA‑TTTTGCTTTTAGGTCAAATTTTTTACCGAGAGTCG‑AGTAAAAACCGGGAGAC‑CATACGGGA‑AAAGGAG‑GTAAAAAA‑CTTCGC‑CACCAAAAATAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACg                                                                                                           <  2:720422/226‑1 (MQ=255)                                                                                                                         gccgGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCGGGTGAG‑CATATCCGT‑AAACGCG‑GTAAAAAA‑CATCGC‑CACCAGAAAGATAATCCCGACATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGCGGAGAAACACTGTCCCCATTCAGAGCATTAAATTCAAACCGACGTGGCcca                                                                                                    >  2:71198/1‑229 (MQ=255)                                                                                                                          ccgGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCCTTAAATTCGAACCGACGTCGCccat                                                                                                   >  1:396688/1‑229 (MQ=255)                                                                                                                                    tGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCACAATATTTTTCCGCGAGTTG‑AGTAACAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGACTTCCCTAATTCATCGGTCTAAAATTTGAGTCGATACAGAGTACCGATCCAGATCATTAAATTCGTACCGACGCAGTCCATACCTTCGGc                                                                                          >  1:575784/1‑228 (MQ=255)                                                                                                                                             aataaaGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTcc                                                                                >  1:457959/1‑229 (MQ=255)                                                                                                                                              ataaaGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCACTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGAGTTTTTccc                                                                               >  1:676284/1‑229 (MQ=255)                                                                                                                                                aaaGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATAACAAAGGCGATTTTTTccccc                                                                             >  2:461292/1‑227 (MQ=255)                                                                                                                                                aaaGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCA‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGAAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATCCGAACCGAAGTCGCCCATACCCTAGCGGATTTTTTCACGc                                                                             >  1:612182/1‑229 (MQ=255)                                                                                                                                                  aGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgag                                                                                                                                                      >  1:18679/1‑154 (MQ=255)                                                                                                                                                  aGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAAATTACCGATGCAAAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTccccccg                                                                           >  1:462331/1‑225 (MQ=255)                                                                                                                                                  aGGGAAGAA‑TTTCGCCTTTAGGCCAAATATTTTGCCGCGAGCCG‑AGAAAAACTCCGGGGAG‑CATAGCGAT‑AAACGCG‑ATAAAAAA‑CTCCAC‑CACCAGAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgag                                                                                                                                                      <  2:18679/154‑1 (MQ=255)                                                                                                                                                     gACAAA‑ATTTGCCTTGAGCCCAAATTTTTTCCCGAGAGTCG‑AAAAAAAACCCGGGTAC‑CATTCCGGT‑AAACGCG‑ATAAAAAA‑CTCCCC‑CCCCAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                        <  1:121218/229‑1 (MQ=255)                                                                                                                                                     gACAAA‑ATTGGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAATTCCGGTGAG‑CATATCGGT‑AAACGGG‑GAAAAAAA‑CTCCGC‑CACAAGAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                        <  1:624631/229‑1 (MQ=255)                                                                                                                                                       acga‑aTTTTCCGTTAAGGCCAATTTTTTTCCGGCAGACGCAAAAAAAACCCGGTGAT‑AATAGGGGA‑AAGCGCG‑AAAAAAAC‑CTTCGC‑CACCAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTc                                                                       <  1:473877/226‑1 (MQ=255)                                                                                                                                                                gCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGAGGTCGCCCAAACCACAGGGGATTTTTTCCCCCAGATCAAGTTTTAAt                                                              >  2:514638/1‑229 (MQ=255)                                                                                                                                                                     tAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGAGGTCGCCCATACCATAGGCGATTTTTTCCCCCAGATTCACTTTTAATGTCAg                                                         >  1:10415/1‑229 (MQ=255)                                                                                                                                                                                tttctCGGCGGGCCG‑AATAAAAACCCGAGGA‑‑CACATCTGG‑AAAAGCG‑GTAAAAAACCTTCCA‑CACAAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTAGAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                              <  2:208144/229‑1 (MQ=255)                                                                                                                                                                                  tttCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCACTACCATAGGCGCTTTTTTCCCGAAGATTCAGTTTTAATTTCCATGGATTATGGTc                                             >  2:475705/1‑228 (MQ=255)                                                                                                                                                                                   ttCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCa                                            >  1:663621/1‑228 (MQ=255)                                                                                                                                                                                        cgagtcg‑agtAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGcc                                                                                                      >  1:69591/1‑165 (MQ=255)                                                                                                                                                                                        cgagtcg‑agtAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTGAGGAGATACAGAGTACCGATGCAGCGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGAGTATggcgggctct                                       >  1:742409/1‑220 (MQ=255)                                                                                                                                                                                          agtcg‑agtAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GCAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGACGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGCCATAAt                                    >  1:661028/1‑229 (MQ=255)                                                                                                                                                                                            ttgtagaaaaaaTCCGGTGAA‑CATATCGGA‑AACCTGT‑GAAAAAAA‑CTTCG‑‑CACCAAAAAGATAACCCCGGCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                  <  1:173910/222‑1 (MQ=255)                                                                                                                                                                                            aaaaaaaaaaaaGCCGGGGCG‑CATAT‑GGT‑AAACGGGAGAAAAAAA‑CTCCGA‑CACAAGAAAGAAAATGCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAAt                                    <  1:669453/221‑1 (MQ=255)                                                                                                                                                                                             tcgaaaaaaaccccGGTGAC‑CATAGCGGT‑AAACGCA‑GAAAAAAA‑CCTCCC‑CACCAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                  <  1:520249/217‑1 (MQ=255)                                                                                                                                                                                                aaaaaaaaaCCGGGGAGACATAGCGGA‑AAACGGG‑AAAAAAAA‑CTCCGC‑CACAAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACt                                <  1:100180/226‑1 (MQ=255)                                                                                                                                                                                                 aataaaaaaacgggtaggaagtCGGA‑AAACGCG‑GGTAAAAA‑CCCCGCACCACAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACt                                <  2:74820/208‑1 (MQ=255)                                                                                                                                                                                                       tcccGGTGAT‑AATCTGGAT‑AACCGGG‑GTAAAAAC‑CTGCGC‑CACAAAAAATATATCCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTtata                         <  2:228056/226‑1 (MQ=255)                                                                                                                                                                                                         ccGGGGAT‑CATATCGGA‑AAACGGT‑GAAAAAAA‑CTCCCC‑CACCAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACg                      <  1:150850/229‑1 (MQ=255)                                                                                                                                                                                                                  cATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTGAGGAGGTACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATTGCAGCGGATTATGGTCAGTCATACTCACCCCTtatattgtctcttt               >  2:669452/1‑218 (MQ=255)                                                                                                                                                                                                                               cgcg‑GAAAAAAA‑CTTCGC‑CACCAGAAAGATAATGCCGACATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat  <  1:299121/229‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                  TACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAA‑ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGT‑AAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTAT  >  NC_000913/4065356‑4065802

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