Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,065,603	T→A	17.1%	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	17.1%	135.6 / 2.0	61	I55F (ATC→TTC)	yihO	putative sulphoquinovose importer
Reads supporting (aligned to +/- strand):  ref base T (33/16);  new base A (5/7);  total (38/23)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACGGGCTTCTAT  >  NC_000913/4065375‑4065834                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           tGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCa                                                                                                                                                                                                                                 >  1:847560/1‑239 (MQ=255) tGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGA‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCa                                                                                                                                                                                                                                 >  1:93325/1‑239 (MQ=255)   tATGCAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGTTAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAACTAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGGAAGC‑ATACCGGTAAACGCGGAAAAAAA‑CT‑TCACAAACAAAAAGA‑AAAA‑CCCGCCATa                                                                                                                                                                                                                               >  2:430110/1‑239 (MQ=255)   tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACACAA‑GGTGGTGGCGATAAACCGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTCAGC‑ATACCGGTAAACGCGCTAAAAAA‑CT‑TCGCCACCCAAAACG‑AAAC‑CCCGCCaat                                                                                                                                                                                                                               >  2:67968/1‑237 (MQ=255)   tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATa                                                                                                                                                                                                                               >  1:940253/1‑239 (MQ=255)   tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATACCGGTAAACGCGGTAAAAAA‑CT‑TCGCACCCAAAAAAA‑AAAT‑CCCCCCaaa                                                                                                                                                                                                                               >  2:124644/1‑237 (MQ=255)   tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGCCGCGTAAAAATC‑CGGTGAGC‑ATATCGGTAAAAGCGGTAAAAAA‑CT‑TCACCACCAGAAAGA‑TAAT‑CCCGCCATa                                                                                                                                                                                                                               >  1:294388/1‑239 (MQ=255)   tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCAAGTCGAGTAAAAATC‑CGGTAAGC‑ATACCGGTAAACTCGGTCAAAAA‑CT‑TCGCCACCAAAAAGA‑AAAT‑ACCCCCaaa                                                                                                                                                                                                                               >  2:702359/1‑237 (MQ=255)        aaGTCAAAACAAAATAAACAGCGGGGTGGCAATGGTCGTTTAAACCGGAAAACAAAA‑GGGGGGGCGGATAAACTGG‑GGGTTGGAA‑ATAACCGCC‑GCCGAAACCGCTG‑AATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCA‑AATTTTTTCCCGCGAGTCGAGAAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑CCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAg                                                                                                                                                                                                                          <  1:124644/239‑1 (MQ=255)             cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGc                                                                                                                                                                                                                     >  1:803868/1‑239 (MQ=255)             cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGAAGGc                                                                                                                                                                                                                     >  1:1018667/1‑239 (MQ=255)             cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGGGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAAAAAGA‑AAAT‑CCCGCCAAAGAAAACAgga                                                                                                                                                                                                                     >  2:678874/1‑238 (MQ=255)               aaCATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGCAAGA‑TAAT‑CCCGCCATAGTAACCAGGCa                                                                                                                                                                                                                    >  1:330673/1‑238 (MQ=255)                           gtctcggttgcactttgttttttTACCGGACAACAAAA‑GGTGGGGGGAAAAAAGGGGGGGGTGTCAATTTAC‑CCCC‑GCCGGACCCTCTG‑AATAAAAAAAAAAGGGAAGAAA‑‑TTGCTTTTGCGCCCAAATATTTTCCGGCGAGCTGAGTAAAAATC‑CGGTGATC‑TTGGGGGAAAACGCGGAAAAAAT‑CT‑AAGCCACAAAAAAAA‑AAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCAt                                                                                                                                                                                                       <  2:60538/220‑1 (MQ=255)                            gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCg                                                                                                                                                                                                     >  1:715632/1‑239 (MQ=255)                            gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGAATCCCAAA‑TTCATCg                                                                                                                                                                                                     >  1:555565/1‑239 (MQ=255)                            gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCAGCCATAGTAAGAAGGCATCACTAA‑TTCATCg                                                                                                                                                                                                     >  1:860830/1‑239 (MQ=255)                                ggTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGtat                                                                                                                                                                                                 >  1:695653/1‑239 (MQ=255)                                         ccgttgtttAACCCGAAAAAAAAAGGGTGGTGGCAATAAACGGGGGCGTTGAAA‑TTAC‑CGCC‑CCCGAACCCGCTG‑AAAAAAAAAAAAAGGGACAAAT‑TTTGCTTTTAGGCCA‑AATATTTTCCCGCGAGTCGAATAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAAAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTg                                                                                                                                                                                         <  1:678874/233‑1 (MQ=255)                                           gtttttaAACCGGAAAAAAAAA‑GGTGG‑GGGGAAAAATCGGGGAGTTTAAA‑ATAA‑ACCC‑CGCGGCACCGCGA‑GAAATAAAAAAAAGGGAAAAAAGTTGGCTTTTCGGCCA‑AATTTTTTCCCGAGGGGTGAGAAAAACGC‑GGAGGAAA‑TTATGGTAAAACGGGGAAAAAAC‑CT‑CCCCCCCCAAAAATA‑AAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTgagg                                                                                                                                                                                      <  2:994084/236‑1 (MQ=255)                                             ttttcacccGGCAAAAAAAG‑GGTGGGGGCGAAAGAGGGGAGGGGTAAAA‑TAAA‑C‑CC‑CCCTGGACCCCTGTAATAAAAAAAAAAGGGA‑AAAA‑ATTGCCTTTTGGCCCAAATTTTTTGCGGAGAGGCGAAAAAAAATC‑GGGTGGGT‑ATAGCGGTAACCG‑GGTAAAAAA‑TTGCCGCCCCAAAAAAAA‑AAAC‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTgaggag                                                                                                                                                                                    <  2:847560/232‑1 (MQ=255)                                                tacccgtccaaacaaaA‑GTTGGTGTGTAT‑AACTGGTGTTCTGCAA‑TTCA‑CCCCAGCCGGACCCTCTA‑CATATAAAATAAAGGTGCGAAT‑TTTGCTTTTACGCTG‑TATTTTTACCCGAGATTCGAAAAAAAACC‑TGGTGACA‑TTATTGGTAACCGATATAAAAAC‑CC‑TCCCCCCCAAAAAAA‑AAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATa                                                                                                                                                                                 <  2:651588/231‑1 (MQ=255)                                                  acccGGAAAAAAAAG‑GGGGGGGGCAATAACCGGGAGCGTTGCAA‑TCAA‑CGCC‑GCCGGAACGGTTG‑ATAAAAAAATAAAGGGACGAAA‑TTTGCCTTAAGGCCA‑ATTATTTTCCCGCGAGTAGAGTAAAAACC‑CGGTGAGC‑ATATCGGTAAACGCGGAAAAAAA‑CT‑TCGCCACCAAAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACa                                                                                                                                                                               <  1:363595/237‑1 (MQ=255)                                                                 ttataagagacaaaaaGTGGAGTGTTTAAA‑TTAC‑CGCC‑GCCGGACCCGCTT‑AATAAAAAAAAAAGGGAAGAAT‑TTTGCTTTTAGGCCC‑AATTTTTTCCCGCGAGGAGAGTAAAAAAC‑CGGTGAGC‑ATATCGGTAAAGGCGGTAAAAAC‑CT‑TCCCCACAAGAAATA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAATTTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCa                                                                                                                                                                 <  2:374847/226‑1 (MQ=255)                                                                      gtggCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTTTAAAATTTGAGGAGGTACAGAGT‑ACCGATGCAGAGCAt                                                                                                                                                           >  1:635370/1‑239 (MQ=255)                                                                            aTAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑ACAATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGCGC‑ATATCGGTAAACGCGGTAAAAAC‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGAAGGCATCCCTAA‑ATCTTCGGAATAAAATTGGAGGAGAACAAAAGA‑ACAGATCCAGAaaattaaatt                                                                                                                                                     >  2:435618/1‑239 (MQ=255)                                                                                         ttGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAATTTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTaaa                                                                                                                                                       >  1:374847/1‑225 (MQ=255)                                                                                          tGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAATAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGGGGAAATACA‑AGTGAACGAACAAAAGCATTAAATTCGAACCAACGTca                                                                                                                                        >  2:484115/1‑237 (MQ=255)                                                                                            gaatatAA‑CCCC‑GCCGGACCCG‑TT‑TATAAAAAAAAAAGGGAAAAAATTTTGCCTTTAGCACA‑AATTTTGTACCGGGGTGAGAATAAAAACC‑GGGAGAGC‑ATAGCGGAAAACGAGAAAAAACT‑CT‑TCCCCAACAAAAAAA‑AAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGcc                                                                                                                                      <  2:342934/234‑1 (MQ=255)                                                                                             gaatttaccccc‑gccgGACCCAGTG‑AAAAAAAAAAAAAGGGAAGATA‑TTTGGTTTTAGGCAA‑ATTTTTTTCCCGCAGGTCGATTAAAAACC‑CGGTGGGC‑ATACCGGTAAACGCGTAAAAAAC‑CT‑TCCCCACCAGAAAAA‑TAAC‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGcc                                                                                                                                      <  2:997990/229‑1 (MQ=255)                                                                                                      gcc‑gccGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGGAATACAGAGA‑ACCAAGCAAAAGACTTAAATTCGAACCGAAGTCGCcaaaaccaaa                                                                                                                             >  2:430825/1‑237 (MQ=255)                                                                                                                   gCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACag                                                                                                                                                                              >  1:90758/1‑179 (MQ=255)                                                                                                                             aaaataaaGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGAAAGCAGGCATCCCTAC‑TGCATCGGTATAAAATTTGAGGAGCTACAGAAT‑ACCGATGCAGAGCATTAAATTCGAACAGACGTTGCCCCTAACATAGACGCATTTTTCCCGCAGAt                                                                                                         >  2:1066438/1‑238 (MQ=255)                                                                                                                                    aGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGCATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCCCAGATTCGGtttt                                                                                                 >  1:770771/1‑239 (MQ=255)                                                                                                                                    aGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCAACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATACCTAC‑TTCATCGGACTAAAATTTGAGGAGATACCGAGG‑ACCGATGCAAGGCATTAAATTCGATACGACGACGCCCATACACAAGGCGAagttttacagacgattaacaatt                                                                                                 >  2:27234/1‑222 (MQ=255)                                                                                                                                         cGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATg                                                                                             >  1:946990/1‑238 (MQ=255)                                                                                                                                         cGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGt                                                                                            >  1:19395/1‑239 (MQ=255)                                                                                                                                         cGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCAAACCAAAGGCGAATTTTTCCCCCAAATTCAGTTTTAATgg                                                                                            >  2:23360/1‑238 (MQ=255)                                                                                                                                           aaa‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑tt                                                                                                                                                                                                          >  1:914365/1‑128 (MQ=255)                                                                                                                                           aaa‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAGAt                                                                                                                                                      >  1:702193/1‑179 (MQ=255)                                                                                                                                            aa‑TTTGCCTTTTGGCCCAAATTTTTTCCCGGAAGTCGAGAAAAAACC‑CGGTGAGA‑ATATCGGTAACCGCGGAAAAAAC‑CT‑TCGCCCCCAGAAAAA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCa                                                                                          <  1:23360/239‑1 (MQ=255)                                                                                                                                                   ccTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGAtt                                                                                   >  1:836812/1‑239 (MQ=255)                                                                                                                                                       tAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGTAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTTAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGGACCGCCGTCGCCCATACCATCGGCGGTTTTTTACCGCAGTTTCAATTTTAAAGTCAGTGGATTAggg                                                                               >  2:412377/1‑239 (MQ=255)                                                                                                                                                        aGGCCC‑ATGATTTTTCCGGGAGTCGAGTAAAAATCGGGGTGA‑A‑ATATCGGAAAGCGGGAAAAAAAA‑CT‑ACCCCACCAGAAAGA‑TAAG‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑TCCGATGCAGAGCATTAAATTCGAGCCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                                              <  1:556193/239‑1 (MQ=255)                                                                                                                                                                    ttttCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGAt                                                                                                                                                                    >  1:23266/1‑142 (MQ=255)                                                                                                                                                                     ttcccgcaggtcgaaaAAAACCC‑CGGTGAGA‑ATCTCGGAAAACGGGAAAAAAAA‑CC‑TCCCCCCCAAAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                                                  <  1:449775/230‑1 (MQ=255)                                                                                                                                                                     tgtccgagacgagagtAAAACTC‑CGG‑GAGCAATCTCGGAAAACGCGGTAAAAAT‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                                                  <  1:576536/227‑1 (MQ=255)                                                                                                                                                                      ttcccgcagatagaataaaaccgcGGTGAGC‑ATATCGGTAACCG‑GGTAAAAAACCT‑TCCCCACAAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                                                  <  1:292460/216‑1 (MQ=255)                                                                                                                                                                       cccGGGAGTCGAGAAAAAACC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAC‑CT‑CCGCCACCAGAAAGA‑TAAC‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACt                                                                <  1:1097917/238‑1 (MQ=255)                                                                                                                                                                         cgcgAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTAAAGGGATTATGGtaaaaaaaaaaaaaccc                                                              >  2:449775/1‑223 (MQ=255)                                                                                                                                                                           cgagtcgagtAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAAAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCAAAGGCGATTTTTTCACCCAGATTCAGTTTTAAt                                                                                              >  2:9347/1‑205 (MQ=255)                                                                                                                                                                                 gagaAAAAGCC‑CGGTGATC‑ATGTCGTAAAACGCGGAAACAAC‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACg                                                      <  1:702359/239‑1 (MQ=255)                                                                                                                                                                                             cGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCa                                                                                                                                                                 >  1:807726/1‑121 (MQ=255)                                                                                                                                                                                                   agaatatggagaaaCGGGGAAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAAGTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGcgtcg                                       <  1:505235/227‑1 (MQ=255)                                                                                                                                                                                                       tacgggaaaccgcgGAAAAAAA‑CT‑TCCCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat                                  <  1:845416/229‑1 (MQ=255)                                                                                                                                                                                                                     gTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAAt                                                                    >  1:74252/1‑191 (MQ=255)                                                                                                                                                                                                                      ataaaaa‑CT‑TCTCCCCCAAAAAAAGTAATACCCGCCATATTAAGCAGGCATCCATAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCg                     <  2:715632/237‑1 (MQ=255)                                                                                                                                                                                                                       aaaaaa‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGAGGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCAGCAGATTCAGTTTTAATGTCAATGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCCACGGGGGTCGGTTATTATTATTAccccat                    >  2:840150/1‑235 (MQ=255)                                                                                                                                                                                                                            a‑CT‑CCGCAACAAGAAAGA‑AAAC‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTa             <  1:484115/239‑1 (MQ=255)                                                                                                                                                                                                                                       ccaAAAAGA‑AAAT‑CCCGTCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTCAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACGGGCTTCt    <  1:840150/239‑1 (MQ=255)                                                                                                                                                                                                                                            aata‑TAAC‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACGGGCTTCTAt  <  1:713534/233‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           TGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAA‑GGTGGTGGCGATAAACTGGAGCGTTGCAA‑TTAA‑CGCC‑GCCGGAACCGCTG‑CATATAAAATAAAGGGACGAAA‑TTTGCCTTTAGGCCC‑AATATTTTTCCGCGAGTCGAGTAAAAATC‑CGGTGAGC‑ATATCGGTAAACGCGGTAAAAAA‑CT‑TCGCCACCAGAAAGA‑TAAT‑CCCGCCATAGTAAGCAGGCATCCCTAA‑TTCATCGGTATAAAATTTGAGGAGATACAGAGT‑ACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACGGGCTTCTAT  >  NC_000913/4065375‑4065834

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: ‑