Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,065,604	T→A	27.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,604	0	T	A	27.8%	65.5 / 4.9	41	I55F (ATC→TTC)	yihO	putative sulphoquinovose importer
Reads supporting (aligned to +/- strand):  ref base T (24/5);  new base A (4/8);  total (28/13)
Fisher's exact test for biased strand distribution p-value = 3.65e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GGCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACG  >  NC_000913/4065378‑4065826                                                                                                                                                                                                                                        |                                                                                                                                                                                                                               tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GGCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAGTCCGG‑TGAG‑CATATCGGAAAAC‑GCGGAAAAAAACTTCGCCACCAGAAAAATa                                                                                                                                                                                                                               >  1:309267/1‑228 (MQ=255) tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GGCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATa                                                                                                                                                                                                                               >  1:517299/1‑228 (MQ=255) tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GGCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑GGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATa                                                                                                                                                                                                                               >  1:517502/1‑228 (MQ=255)       aGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAACGGTGGT‑GGCGCTAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAACTAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCAAGCAAAAATCCGG‑TGAG‑CAGATCGGAAAAC‑GCGGTAAAAAACTTGGCAACCAGAAAGATAATACCg                                                                                                                                                                                                                         >  1:583109/1‑228 (MQ=255)       aGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAACGGTGGT‑GGCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCCTATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGTGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GTGGTAAAAACCTTCCACCACAGAAAGAAAAAccct                                                                                                                                                                                                                         >  2:140193/1‑227 (MQ=255)         tCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GGCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGcc                                                                                                                                                                                                                       >  1:13971/1‑228 (MQ=255)           cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GGCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGGAAAAAACTTCGCCACCAGAAAGAAAATccctgtc                                                                                                                                                                                                                      >  1:146533/1‑223 (MQ=255)           cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GGCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGAAGAAATTTGCCTTTAGGCACAATATTTTTCCGAGATACGGGGAAAAATCCGGCTGA‑‑CATATCGGAAAAC‑GCGGTAAAAAACTTCACAAACAAAAAGAAAATCCCGccca                                                                                                                                                                                                                     >  2:347092/1‑226 (MQ=255)                  catAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GGCGATAAACCGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGAAAGc                                                                                                                                                                                                              >  1:67994/1‑228 (MQ=255)                           cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GGCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCACGAGTCGAGTAAAAATCCGG‑TGGG‑CATATCGGTAAAC‑CCGGGAAAAACCTTCGCCAACAGAAAGAAAATCCCGCCAAAGAAAGAAGGCAAccc                                                                                                                                                                                                     >  2:362839/1‑228 (MQ=255)                           cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GGCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAAT‑AAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATACCt                                                                                                                                                                                                    >  1:342996/1‑228 (MQ=255)                                                 aCCGGCAAACAAAAGGTGGG‑GGCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTg                                                                                                                                                                               >  1:486698/1‑228 (MQ=255)                                                    cagaaaaaaaaGGGGGTGGGAAATAAAGGGGAGCG‑TGCAAATAACCGCCGCGGGCACCGTTCCATAAAAAAAAAAGGAAAGAATTTCCCCTTTGGCCACAATTTTTCCCCGAGGTTCGAGTAAAAACCGGG‑AGAT‑CTTACGGGAAAAC‑GCGGAAAAAAATTCCCCCACCAGAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgag                                                                                                                                                                             <  2:202392/225‑1 (MQ=255)                                                                    t‑ggCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGCAATTTACCTTTAGGCCCAATATTTTTCCGCGAGTCGCGTAACAATCCGG‑TGCG‑CATATCGGCAAAG‑CCGGTAAAAACCTTCACCACCAAAAAGATAACCCCGCCACAGTAACCAGGCATCCCAACTTCCTCGCAATAAAATTTAGGGAGATCCA‑AGTAcacga                                                                                                                                                           >  2:379137/1‑224 (MQ=255)                                                                      ggCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTGAGGAGGTACAAGGTACCCAt                                                                                                                                                           >  2:530861/1‑228 (MQ=255)                                                                                        tGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGCATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCg                                                                                                                                          >  1:433306/1‑227 (MQ=255)                                                                                        tGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGAAGGCATCCCTAATTCCTCGGTATAAAATTTGAGGAGATACACGGTACCGATGCAGCGAATTAAATTCg                                                                                                                                          >  2:502923/1‑227 (MQ=255)                                                                                                                               aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTCAGCAGGCATCCCTAATTCATCGGTATAAAAATTGAGGAGATACACAGTACCGTTGAAGCCCATTAAATTCGAACCGACGTCTCCCAAACCATAGGGGATTTTTTccccag                                                                                                   >  2:239209/1‑225 (MQ=255)                                                                                                                               aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGACCCAACGTCGCCCATCCCATAGGCGATTTTTTCCCGCa                                                                                                   >  1:632113/1‑228 (MQ=255)                                                                                                                               aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCcaca                                                                                                   >  2:444893/1‑228 (MQ=255)                                                                                                                               aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCACCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCa                                                                                                   >  1:176460/1‑228 (MQ=255)                                                                                                                                   agaaaaTTCGCTTTGGGGCCAAATATTTCTCCGAGAGTAGAGAAAAAATCGGG‑TGAG‑CATATCGGAAAAC‑GCGGTAAAAAACTCCACCACCAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                                               <  1:433280/225‑1 (MQ=255)                                                                                                                                    cGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTc                                                                                              >  1:463798/1‑228 (MQ=255)                                                                                                                                                 tAGGCCTAATTTTTTGCCGGGGATCGAATAAAAAGAGGG‑TGAG‑CATGGCGGAAACC‑GGGGAAAAAAACCGCGCACCCAAAAAAATAACTCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCa                                                                                 <  2:286555/228‑1 (MQ=255)                                                                                                                                                 tAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga                                                                                                                                                                              >  2:182510/1‑135 (MQ=255)                                                                                                                                                 tAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGAGGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGACa                                                                                 >  1:286555/1‑228 (MQ=255)                                                                                                                                                 tAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCAAAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCa                                                                                 >  1:502747/1‑228 (MQ=255)                                                                                                                                                 tAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGCATAAAATTTGAGGAGATACAGGGTACCGATGCAGAGCATTAAATTTGAACCGACGTCGCCCATACCAAAGGCGAATTTTTCCTCCAAATCCAGTTTTaacttc                                                                                  >  2:239974/1‑223 (MQ=255)                                                                                                                                                 aaggcaaaaTGTTTTTCCGCGAGCTGA‑TAAAAAACCCG‑GGAGTCATCGCGGTAACC‑GGGGTAAAAATCTTCGCCACAAGAAGGAAAAACCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCa                                                                                 <  1:239209/221‑1 (MQ=255)                                                                                                                                                  aGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgag                                                                                                                                                                             <  1:182510/135‑1 (MQ=255)                                                                                                                                                   gcccaaaTTTTTTCCCGCAAACCGAAAAAAACCCCGG‑GGAG‑CATACCGTAAAGC‑GCGATAAAAAACCTCCCCACGAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGt                                                                               <  1:475893/223‑1 (MQ=255)                                                                                                                                                             ttttCCGGAGGTCGAATAAAAATCGGG‑TGAG‑ATTACCGGAAAAC‑GGGGAAAAAAACTCCGCCCCAAAAAAGAAAACCCCGCCATAGTTAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                                     <  1:615454/228‑1 (MQ=255)                                                                                                                                                                         tgaatAAAACTCCGG‑TG‑G‑CCTTTCGGGAAACGGAGGAAAAAAACTCCGCCAACAGAAAAATATGCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                                         <  2:159657/224‑1 (MQ=255)                                                                                                                                                                                                      tAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCCGATTCAGTTTTAATGTCAGTGGATTTTGGTCACACATAATCACTCCTATTAACCCCCCTCAGGGggtgt                              >  2:150335/1‑226 (MQ=255)                                                                                                                                                                                                      tAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                                     >  1:506046/1‑189 (MQ=255)                                                                                                                                                                                                        aaC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGCGGAGATACAGAGTACCGATGCAGAGCATTAACTTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGCGGATTATTGTCAGACATAATCACTCATTATAACGCCCCGCAGGGcgccgt                             >  1:384750/1‑227 (MQ=255)                                                                                                                                                                                                            gcgtaaaaaaaCTCCGCCCGCAGAAAGAAAACCCCGACATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat                         <  1:502923/223‑1 (MQ=255)                                                                                                                                                                                                            gcgATAAAAACTTTCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat                         <  1:556112/228‑1 (MQ=255)                                                                                                                                                                                                            gatgTAAAAAGCTCCGCCACAAAAAATAAAACCCCGCCATAGTAAAGAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATCCAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat                         <  2:486698/225‑1 (MQ=255)                                                                                                                                                                                                                              acaaGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGt       <  1:150335/225‑1 (MQ=255)                                                                                                                                                                                                                                   aaaGATAACCCCGCCATAGTAAGGAGGCCTCCCTAATTCATCGGTATAACATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACg  <  1:259094/228‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                               TATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGT‑GGCGATAAACTGGAGCGTTGCAATTAA‑CGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGG‑TGAG‑CATATCGGTAAAC‑GCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACG  >  NC_000913/4065378‑4065826

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