Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,065,603	T→A	29.8%	I55F (ATC→TTC)	yihO ←	putative sulphoquinovose importer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,065,603	0	T	A	29.8%	66.4 / 7.8	45	I55F (ATC→TTC)	yihO	putative sulphoquinovose importer
Reads supporting (aligned to +/- strand):  ref base T (24/7);  new base A (4/10);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 2.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GG‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACG  >  NC_000913/4065377‑4065825                                                                                                                                                                                                                                          |                                                                                                                                                                                                                               tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GG‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAGTCCGG‑TGAG‑CATATCGGAAAAC‑GCGGAAAAAAACTTCGCCACCAG‑AAAAATa                                                                                                                                                                                                                               >  2:309267/1‑228 (MQ=255) tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GG‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATa                                                                                                                                                                                                                               >  2:517299/1‑228 (MQ=255) tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GG‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑GGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATa                                                                                                                                                                                                                               >  2:517502/1‑228 (MQ=255)       aGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAACGGTGGT‑GG‑CGCTAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAACTAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCAAGCAAAAATCCGG‑TGAG‑CAGATCGGAAAAC‑GCGGTAAAAAACTTGGCAACCAG‑AAAGATAATACCg                                                                                                                                                                                                                         >  2:583109/1‑228 (MQ=255)       aGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAACGGTGGT‑GG‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCCTATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGTGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GTGGTAAAAACCTTCCACCACAG‑AAAGAAAAAccct                                                                                                                                                                                                                         >  1:140193/1‑227 (MQ=255)         tCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GG‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGcc                                                                                                                                                                                                                       >  2:13971/1‑228 (MQ=255)           cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GG‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGGAAAAAACTTCGCCACCAG‑AAAGAAAATccctgtc                                                                                                                                                                                                                      >  2:146533/1‑223 (MQ=255)           cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GG‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGAAGAAATTTGCCTTTAGGCACAATATTTTTCCGAGATACGGGGAAAAATCCGGCTGA‑‑CATATCGGAAAAC‑GCGGTAAAAAACTTCACAAACAA‑AAAGAAAATCCCGccca                                                                                                                                                                                                                     >  1:347092/1‑226 (MQ=255)                  catAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GG‑CGATAAACCGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGAAAGc                                                                                                                                                                                                              >  2:67994/1‑228 (MQ=255)                       acgcgggggggtgaaattggtttttAACCCCGAAACAAATGGGGGT‑GGTAGAAAAAACGGGGTGTTGAAATTAC‑CCCCGCCGGAACCGCTGCATAAAAAAAAAAGGGATGAAATTGGCCTTGAGGCCAAATTTTTTCCCGCGGGGCGAGAAAAAATCGGG‑TGAG‑AATACCGGTAAAC‑GCGGAAAAAACCTCCGCCCCCAG‑AAAGATAATCCCGCCATAGTAAGCAGGc                                                                                                                                                                                                          <  1:67994/207‑1 (MQ=255)                           cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GG‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCACGAGTCGAGTAAAAATCCGG‑TGGG‑CATATCGGTAAAC‑CCGGGAAAAACCTTCGCCAACAG‑AAAGAAAATCCCGCCAAAGAAAGAAGGCAAccc                                                                                                                                                                                                     >  1:362839/1‑228 (MQ=255)                           cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GG‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAAT‑AAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATACCt                                                                                                                                                                                                    >  2:342996/1‑228 (MQ=255)                                                 aCCGGCAAACAAAAGGTGGG‑GG‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTg                                                                                                                                                                               >  2:486698/1‑228 (MQ=255)                                                    cagaaaaaaaaGGGGGTGGG‑AAATAAAGGGGAGCG‑TGCAAATAACCGCCGCGGGCACCGTTCCATAAAAAAAAAAGGAAAGAATTTCCCCTTTGGCCACAATTTTTCCCCGAGGTTCGAGTAAAAACCGGG‑AGAT‑CTTACGGGAAAAC‑GCGGAAAAAAATTCCCCCACCAG‑AAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgag                                                                                                                                                                             <  1:202392/225‑1 (MQ=255)                                                                    t‑gg‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGCAATTTACCTTTAGGCCCAATATTTTTCCGCGAGTCGCGTAACAATCCGG‑TGCG‑CATATCGGCAAAG‑CCGGTAAAAACCTTCACCACCAA‑AAAGATAACCCCGCCACAGTAACCAGGCATCCCAACTTCCTCGCAATAAAATTTAGGGAGATCCA‑AGTAcacga                                                                                                                                                           >  1:379137/1‑224 (MQ=255)                                                                      gg‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTGAGGAGGTACAAGGTACCCAt                                                                                                                                                           >  1:530861/1‑228 (MQ=255)                                                                                         tGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGCATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCg                                                                                                                                          >  2:433306/1‑227 (MQ=255)                                                                                         tGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGAAGGCATCCCTAATTCCTCGGTATAAAATTTGAGGAGATACACGGTACCGATGCAGCGAATTAAATTCg                                                                                                                                          >  1:502923/1‑227 (MQ=255)                                                                                                                                aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTCAGCAGGCATCCCTAATTCATCGGTATAAAAATTGAGGAGATACACAGTACCGTTGAAGCCCATTAAATTCGAACCGACGTCTCCCAAACCATAGGGGATTTTTTccccag                                                                                                   >  1:239209/1‑225 (MQ=255)                                                                                                                                aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGACCCAACGTCGCCCATCCCATAGGCGATTTTTTCCCGCa                                                                                                   >  2:632113/1‑228 (MQ=255)                                                                                                                                aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCcaca                                                                                                   >  1:444893/1‑228 (MQ=255)                                                                                                                                aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCACCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCa                                                                                                   >  2:176460/1‑228 (MQ=255)                                                                                                                                    agaaaaTTCGCTTTGGGGCCAAATATTTCTCCGAGAGTAGAGAAAAAATCGGG‑TGAG‑CATATCGGAAAAC‑GCGGTAAAAAACTCCACCACCAA‑AAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                                               <  2:433280/225‑1 (MQ=255)                                                                                                                                     cGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTc                                                                                              >  2:463798/1‑228 (MQ=255)                                                                                                                                                  tAGGCCTAATTTTTTGCCGGGGATCGAATAAAAAGAGGG‑TGAG‑CATGGCGGAAACC‑GGGGAAAAAAACCGCGCACCCAA‑AAAAATAACTCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCa                                                                                 <  1:286555/228‑1 (MQ=255)                                                                                                                                                  tAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga                                                                                                                                                                              >  1:182510/1‑135 (MQ=255)                                                                                                                                                  tAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGAGGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGACa                                                                                 >  2:286555/1‑228 (MQ=255)                                                                                                                                                  tAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCAAAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCa                                                                                 >  2:502747/1‑228 (MQ=255)                                                                                                                                                  tAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGCATAAAATTTGAGGAGATACAGGGTACCGATGCAGAGCATTAAATTTGAACCGACGTCGCCCATACCAAAGGCGAATTTTTCCTCCAAATCCAGTTTTaacttc                                                                                  >  1:239974/1‑223 (MQ=255)                                                                                                                                                  aaggcaaaaTGTTTTTCCGCGAGCTGA‑TAAAAAACCCG‑GGAGTCATCGCGGTAACC‑GGGGTAAAAATCTTCGCCACAAG‑AAGGAAAAACCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCa                                                                                 <  2:239209/221‑1 (MQ=255)                                                                                                                                                   aGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgag                                                                                                                                                                             <  2:182510/135‑1 (MQ=255)                                                                                                                                                    gcccaaaTTTTTTCCCGCAAACCGAAAAAAACCCCGG‑GGAG‑CATACCGTAAAGC‑GCGATAAAAAACCTCCCCACGAA‑AAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGt                                                                               <  2:475893/223‑1 (MQ=255)                                                                                                                                                         aattttttccggagaaccgaaaaaaacTCCGG‑TGAG‑AATACGGGTAAAC‑GGGGAAAAAAATTTCGCCACCAA‑AAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGAtt                                                                          <  2:347092/206‑1 (MQ=255)                                                                                                                                                              ttttCCGGAGGTCGAATAAAAATCGGG‑TGAG‑ATTACCGGAAAAC‑GGGGAAAAAAACTCCGCCCCAAA‑AAAGAAAACCCCGCCATAGTTAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                                     <  2:615454/228‑1 (MQ=255)                                                                                                                                                              tttccccggattcaaatAAAAGTCCGG‑GGAC‑AATAGCGGAAAGC‑CCGGT‑AAAAACCTCCCCAACAGAAAAGATAACGCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                                     <  2:362839/212‑1 (MQ=255)                                                                                                                                                                  cggtgcggctagaaaaacTCGGG‑TGAG‑TATATCGTTAACG‑GCAGAAAAAAACGTCGCCCACAA‑AAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGa                                                                 <  1:13971/215‑1 (MQ=255)                                                                                                                                                                          tgaatAAAACTCCGG‑TG‑G‑CCTTTCGGGAAACGGAGGAAAAAAACTCCGCCAACAG‑AAAAATATGCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                                         <  1:159657/224‑1 (MQ=255)                                                                                                                                                                                                       tAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCCGATTCAGTTTTAATGTCAGTGGATTTTGGTCACACATAATCACTCCTATTAACCCCCCTCAGGGggtgt                              >  1:150335/1‑226 (MQ=255)                                                                                                                                                                                                       tAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                                     >  2:506046/1‑189 (MQ=255)                                                                                                                                                                                                         aaC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGCGGAGATACAGAGTACCGATGCAGAGCATTAACTTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGCGGATTATTGTCAGACATAATCACTCATTATAACGCCCCGCAGGGcgccgt                             >  2:384750/1‑227 (MQ=255)                                                                                                                                                                                                             gcgtaaaaaaaCTCCGCCCGCAG‑AAAGAAAACCCCGACATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat                         <  2:502923/223‑1 (MQ=255)                                                                                                                                                                                                             gcgATAAAAACTTTCCCCCCAAA‑AAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat                         <  2:556112/228‑1 (MQ=255)                                                                                                                                                                                                             gatgTAAAAAGCTCCGCCACAAA‑AAATAAAACCCCGCCATAGTAAAGAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATCCAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat                         <  1:486698/225‑1 (MQ=255)                                                                                                                                                                                                                               acaaG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGt       <  2:150335/225‑1 (MQ=255)                                                                                                                                                                                                                                     aaaGATAACCCCGCCATAGTAAGGAGGCCTCCCTAATTCATCGGTATAACATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACg  <  2:259094/228‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                               TATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GG‑CGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAG‑AAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACG  >  NC_000913/4065377‑4065825

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