Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,520,118	C→T	18.5%	intergenic (+54/‑144)	curA → / → mcbR	curcumin/dihydrocurcumin reductase, NADPH‑dependent/colanic acid and biofilm gene transcriptional regulator, MqsR‑controlled

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,520,118	0	C	T	18.5%	145.6 / 3.4	61	intergenic (+54/‑144)	curA/mcbR	curcumin/dihydrocurcumin reductase, NADPH‑dependent/colanic acid and biofilm gene transcriptional regulator, MqsR‑controlled
Reads supporting (aligned to +/- strand):  ref base C (14/33);  new base T (10/2);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.87e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAA‑TGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAA  >  NC_000913/1519887‑1520350                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              tttgtgtcccgcgaTCCATGTGTTGAAGGGGGGTATGGGGGAG‑GGGGTGAAAGAGAATAACCTCCCCTCCCGCGAAAAAATGACTGAGGGTTAAGGGAAGGCGCCACTGACGTTTACCGCCGTGCGGAAGGGTAAAAG‑CGTCGGCAAA‑GTGGGGACCCGCGTGGCGGGTGAAGATTAATGTTTCAACGGC‑GGCTTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTtta                                                                                                                                                                                                                                        <  2:490552/226‑1 (MQ=255) ttgtGGCCACCGCATCCATGAGTTTCAGGGGGAGAGGGGGCAA‑TGGGGAAAAGAGGAAAAAACCCACTGCCGAGAAAAAATTCCTGACGGTTTAGAGAATGCGCCACAGACTTTTATCGGCCTGCTGAGGGGAAAAAA‑CTTCGGAAAA‑GTGGTGATCCGCGTGGCGGGTGATAATTAAATTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTtta                                                                                                                                                                                                                                        <  1:990973/236‑1 (MQ=255)    tGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAA‑TGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGGGATGATTAAAGTTTCAACGGC‑GGCGGAAAC‑CCCCCTAAAAATCCAGGGAAAAAAC‑TTTTTTTCAGCTttattt                                                                                                                                                                                                                                     >  1:839257/1‑239 (MQ=255)      gcccccGCCTCCATGAGTTAGAG‑GGGGGAGGGGGGAAGTGGGGAAAAGAGAATAAAACCCACTACCGCAAAGAAATTACGGGCGTTTGAGAGAACGCCCCACAGTCGTTTCGCGGCAGGATAAAGGATAAAAA‑CTTCGGCAAA‑GGGGTGACCCGCGGGGGGGGTGATGATAAAAGTTTCAACGGC‑GGCGAAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTttaatt‑aa                                                                                                                                                                                                                                  <  2:1078798/235‑1 (MQ=255)         aCCGCATCCATGAGTTTCAGAGGGAGATGGGGCAA‑TGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAAC‑CGCCCTAAAAATCCAGGCCAAAAGC‑CTATTTTCAACTtttatt‑aaatt                                                                                                                                                                                                                               >  1:936103/1‑236 (MQ=255)         aCCGCATCCATGAGTTTCAGAGGGAGATGGGGCAA‑TGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGGAAGC‑CGGCCCAAAAATCCAGGCAAAAAGC‑ATATTTTCCACTTTAATt‑aacta                                                                                                                                                                                                                               >  1:829410/1‑239 (MQ=255)         aCCGCATCCATGAGTTTCAGAGGGAGATGGGGCAA‑TGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGGGATCCGCGTGGCGGGTGATGATTAAAGATTCAACCGCGGGGGTAAG‑‑CGCCCTCAAAAAACAAGCAAAAAACATAATTTTACGCTttaata‑aaaa                                                                                                                                                                                                                                >  2:444523/1‑237 (MQ=255)         aCCGCATCCATGAGGTTCAGAGGGAGATGGGGCAA‑TGGGTGAAAGCGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGGAAAAA‑CTGCGGCAAA‑GTGGACGTCCGCGTGGCGGGGGAGGGTTTAAGGTTCAACGGCGGGGGAAAGC‑CCACA‑AAAAAACCTGGCAAAAAAC‑CTATTTTTGCCTTTAATa‑aacaa                                                                                                                                                                                                                               >  2:824213/1‑237 (MQ=255)             tatccgttgtttagggggggAGAGGGGGTAA‑GGGGTAAAGGAGGAAAAAACCAGCGACGGGAAAAAAATTGCTAGTGGTTGGGAGAGGGCCCCACAGACGTTTATCGCCCGGCGAAAGGGAAAAAA‑CTGCGGAAAA‑GTGGTGACCCGCGTGGGGGGTTATTAATAAAGTTTCACCGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACca                                                                                                                                                                                                                           <  2:936103/222‑1 (MQ=255)               tCCATGAGTTTCAGAGGGAGATGGGGCAA‑TGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGGGATCCGCGTGGCGGGGGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTTAACTTTAATT‑AACTAACccca                                                                                                                                                                                                                         >  1:418402/1‑239 (MQ=255)                        ttaagggggggtgggggAAA‑TGGGGAAAAGGGGATAAAACCCACTCCGGCAAAAAAATTACTGACGTTTTAGAGAATGCGCCACAGAGGTTTAGGGGCCTGCTGAAGGGTAAAAA‑CGTCGGAAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGtt                                                                                                                                                                                                                 <  1:444523/225‑1 (MQ=255)                          gagggggAGAGGGGGCAA‑TGGGTAAAAGAGGAAAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACGGACTTTTATCGGCCGGCGGAAGGGAAAAAC‑CTTCGGAAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTctg                                                                                                                                                                                                              <  1:83166/235‑1 (MQ=255)                               ggAGATGGGGCAA‑TGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAAC‑CGCCCTAAAAAACCAGGCAAAAAGC‑TTATTTTAAGCTTTAATT‑AACTAACCACATCAATTTTCTg‑ct‑gct                                                                                                                                                                                                       >  1:690656/1‑239 (MQ=255)                                  gaTGGGGCAA‑TGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGGTTCAACGGC‑GGCGGAAAC‑CGCCCTAAAAAACCCAGGAAAAAAC‑TTAATTTTCGCCTTTATt‑aactaacacacacaattttatgctgttaaa                                                                                                                                                                                                      >  2:33505/1‑216 (MQ=255)                                   aTGGGGAAA‑GGGGTGAAAGAGGATAAAATCCACTCCCGGAAAAAAATTACTGACGGTTTAGAGAATGCGCCACAGCCGTTTACCGGCCTGCTGAAGGGTAAAAA‑CTTCGGAAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAATTTTAAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCtaat                                                                                                                                                                                                    <  1:544597/238‑1 (MQ=255)                                        aaaa‑TGGGGAAAAGAGGAAAAAACCCAAGCCCGGAAAGAATTTACTGAGGGTTAAGAGAATGCCCCACAGACTTTTATCGGCCTGCGGAAGGTTAAAAC‑TTTCGGCAAA‑GTGGTAATCGGCGGGGCGGGTAATGATTAAAGTTTCAAGGGC‑GGTGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATaacaa                                                                                                                                                                                               <  1:33505/236‑1 (MQ=255)                                          gagggggTAAAAGGAGAAAAACCCCACTGCCGGAAAGAATTGTCTGACGGTAAAGAAGATGCCAGACAGACTTTTGCGGGCCTGCTAAAGGTAAAAAA‑CTTCGGAAAA‑GTGATGACCGGCGTGGCGGGTGATAATAAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAacaaca                                                                                                                                                                                             <  2:790260/235‑1 (MQ=255)                                                  gtaaaagGAAAAGATCCACTGCCGAGAAGAAATTCATGAGGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCGGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGGGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAg                                                                                                                                                                                     <  1:638490/234‑1 (MQ=255)                                                   aaagggaaaaaaaACCATCACCGAGAAGAAATTCCTGACGGTGTAGAGAAGGCGCCACAGACGTTTAGCGGCTTGATGAAGGGAAAAAA‑CTTGGGCAAA‑GGGGTGATCCGCGGGGCGGGGGAGGATTAAAGTTTAAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAg                                                                                                                                                                                     <  1:1043493/229‑1 (MQ=255)                                                   aaaGGGGATAAAACCCACTCCCGGAAAGAAATTACTGACGGTTTAGAGAATGCGCAACAGACTTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGAAAA‑GTGGTGACCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAg                                                                                                                                                                                     <  1:800241/238‑1 (MQ=255)                                                   aaaGAGGATAAAACCCATTCCGGGGAAAAATTGTCTGGCGTTTAAAAGAATGCCCCACAGCCGTTTTCCGGCCTGCTGAAGGTAAAAAC‑CTTCGGAAAA‑GTGGTAATCCGGGTGGCGGTTGATGATAAAATTTCCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAg                                                                                                                                                                                     <  1:927156/238‑1 (MQ=255)                                                   aaaGAGGATAAAACCCACTCCCGAGAAGAAATTACTGCGGTTTTAGAGAAGGCCCCACAGACGTTTATCGGCCTGCGGAAGGGTAAAAA‑CTTCGGAAAA‑GTGGTGATCCGCGGGGCGGGTGATGATTAAAGTTTAAAGGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAg                                                                                                                                                                                     <  1:610795/238‑1 (MQ=255)                                                        gaaTAAACTCCACTCCCGGGAAGAAATTCCTGACGGTTTAGGAAAGTCGCCAAAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTGGGCAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGtt                                                                                                                                                                                <  1:443383/236‑1 (MQ=255)                                                         gATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGGGATGGTTAAAGTTTCAACGGC‑GGCGGAAAG‑CGCCCTAAAAAACCAAGGCAAAAGC‑GTATTTTTTCCTTTAAAT‑TAATAACCAAAAAATTTTTTTT‑GT‑TCTAATAAAAAAAAAAATGAGAGGgttt                                                                                                                                                                              >  2:800241/1‑238 (MQ=255)                                                                   aatACCGCGAAGAAATTACTGACGGTGTAGAGAATGGGCAACAGGGTTTTAGCGGCCTGCTAAAGGGTAAAAA‑GTTGGGAAAA‑GGGGTGATCGGCGGGGCGGGTGATGATTAAGGTTCCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCa                                                                                                                                                                    <  1:44765/237‑1 (MQ=255)                                                                   aCTACCGCGAAAAATTTTCTGGGGGTATAAAGAATGCGACACAGACGTTTACCGGCCTGCTGAAGGGAAAAAA‑CTTCGGAAAA‑GTGGGGATCCGGGGGGCGGGTAATGATTAAAGTTTAAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCa                                                                                                                                                                    <  2:678866/239‑1 (MQ=255)                                                                    cTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACCGCCGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GAGGTGATCCGCGTGGCGGGTGATGGTTAAAGATTCCACGGC‑GGCGGAAAA‑GCCCCTAAAAATTCACGCACAAAGC‑CTATTTTCTACTTTAATT‑AACACAACACACAACTTTTCTG‑CG‑GTTAATAaaaaaaaaaaagagagggtaataataaata                                                                                                                                                                     >  2:528469/1‑215 (MQ=255)                                                                    cTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTCATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGGTTCAACGGC‑GGCGTAAAC‑CGCCCTAAAAATCCAGGCCAAAAGG‑TTATTTTTGGCGTTTATT‑AAATAACCAACTCATTGTTCTT‑CG‑GGTAATAaaaaaaaaaagaaggggtttataataaatt                                                                                                                                                                     >  2:341841/1‑215 (MQ=255)                                                                          cgAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAAC‑TTATTTTTAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTACCCTTA‑At                                                                                                                                                            >  1:678866/1‑239 (MQ=255)                                                                          aaaggaaATTAAGGAGGGTTTAGAAAATGCGCCAAAGGCGTTTACGGGCCTGGTGAAGGGTAAAAC‑CTGCGGAAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTa‑a                                                                                                                                                             <  1:497309/234‑1 (MQ=255)                                                                             gtaaaaggagcagggTTTAGAGAAGGCGCCACAGTCTTTTAGCGGCCTGCTGAAGGGTAAAAA‑CTTCGGAAAG‑GTGGTGATCGGGGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATaca                                                                                                                                                         <  1:103238/226‑1 (MQ=255)                                                                                         cGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGGGATGATTAAAGTTTCAACGGC‑GGCGGAAAC‑CGCCCTTAAAAACCCAGCAAAAAAA‑TTATTTTTCGCTTTAATTAAACTAA‑CACAACATTTTTTTGTCT‑CCT‑ATAAAAAAAAAAGTGAGGGGGTTTTTATTATTTCCttttaaaaaaaaagggctgt                                                                                                                                               >  2:497309/1‑221 (MQ=255)                                                                                            ttgagagAAGGCGCCAGAGATGTTTACCGGCTTGCTGAAGGGTAAAAA‑TTCCGGCAAA‑GTGGTGAACCGCGTGGCGGGTGATGATTAATGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAg                                                                                                                                                                                     <  1:878662/194‑1 (MQ=255)                                                                                                    aTGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCAt                                                                                                                                  >  1:935393/1‑238 (MQ=255)                                                                                                        gccgcggactttagacggccggcTGAAGGGTAAAAA‑CTTGGGCAAA‑GGGGGGATCCGGTGGGGGGGTGATAATTAAGTTTCCAAGGGC‑GGCGAAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATgcagg                                                                                                                             <  1:81744/224‑1 (MQ=255)                                                                                                        ccccctgcctttgtcgggCCTGCT‑TAGGGTAAAAACCTGCGGAAAG‑GTGTTGCCCGGGGTGGGGGTTGATGTTAAAAGTTTCAGCGGC‑GGCGAAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATgcagg                                                                                                                             <  2:935393/223‑1 (MQ=255)                                                                                                                 tttttggggctgggggAAGGG‑AAAAA‑CTGCGGGAAAGGGGGTGACCCGGGTGGGGGGTGATGATAAAATTTTAAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTg                                                                                                                     <  1:567549/223‑1 (MQ=255)                                                                                                                   tgacgggCCTGCAGAAGGGTAAAAA‑CTGGGCAAAA‑GGTGTGCTCGGCGTGGCGGGTTGAGATTAAAGTTTCAAGGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGAc                                                                                                                  <  1:627091/234‑1 (MQ=255)                                                                                                                         gCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCCTAAAAATCCAGGCAAAAAAC‑TTATTTTCAGCTTTAATT‑AACTAACCCCATCATTGTTCTG‑CT‑GCTAATAAAAACAAAGGTGAGTGGTTTATGATTATTCACCTTA‑ATACACCAGGGG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACCGaaaa                                                                                                            >  1:16942/1‑239 (MQ=255)                                                                                                                                               tCGGCAAA‑GTGGTGATCCGCGTGGCGGGGGATGATTAAAGTTTCAAACGC‑GGGGGAAGCTCCCGAAAAAAATCAAGG‑AAAAAGC‑TTATGTTCACGCTTAATT‑TACTAAACCCAAAATTTTCTTG‑CTGGCAAATACAAAAAAAGGTATGTGG‑TTATAATTATTACTTTTATAAACACCAGGCG‑‑ATTCCTTGTTCCCTACAGAGGGGGGTGTAGAAAAAAAAAACGAAGGAAATAAAAAGc                                                                                       >  2:864356/1‑239 (MQ=255)                                                                                                                                                cGGCAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATgca                                                                                                                               >  1:3721/1‑198 (MQ=255)                                                                                                                                                     aa‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTTAGCTTTAATT‑AAATAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGGG‑AAATCCTTCTGCCCTGCAGGCAGGGTTGGACAGAAAACACCCGAAGGAATTtaaactgttcctt                                                                                 >  1:166637/1‑230 (MQ=255)                                                                                                                                                           gttacccGGGTGGCGTGTAAAGATTAAAGTTGCATCGGC‑GGCGAAAGC‑CGCCATAAAAAGCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACggg                                                                             <  1:981008/233‑1 (MQ=255)                                                                                                                                                           gtgtgccggggggggggTGAAGATTAAATTTTCAAGGGC‑GCGGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACggg                                                                             <  1:824213/224‑1 (MQ=255)                                                                                                                                                                   cgTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTGAAC‑CGCCCTAAAAATCCAGGCAAAAAGC‑TTATTTTTCGCTTTAAAT‑AACTAAACACAACATTGTTCTT‑CT‑GTTAATTAAAAAAAAAGGGAGTGGGTTTTGATTTTTTCCTTTA‑ATAAACACGGTGAAATTCCTTTGCCCATCAAGGAGGGGTTGGGAAAAAAAAACCGCAAGGAATAAAAAGCCTTGGAAGGGGaaaaaaaa                                                                     >  2:1006150/1‑236 (MQ=255)                                                                                                                                                                   cgTGGCGGGGGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCGTCAAAATTCCGGCAAAAAAC‑TTATTTTTAGCCTTAATT‑AAATAACCACATCATTGTTTTG‑CT‑GCTCATTACAAAAAAGATGAGTGGGTTTTGATTATTTCCTTTT‑TAACAACCGGGG‑AATTTCTTCTTCCCGTCAGTCAGGGGTTGAAACAAAAAAACCAAAAGCATAAAAAAACCTGGAAAAGGaaaaaaaaa                                                                    >  2:508857/1‑235 (MQ=255)                                                                                                                                                                      gtctgttattgatgaAATTTTAACCGGT‑GACGCAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTGA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAAc                                                                  <  1:197620/229‑1 (MQ=255)                                                                                                                                                                        cggttaatgttTAAATTTTCAACGTC‑GAGTTAAGC‑CGCAATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACAt                                                                <  1:109270/232‑1 (MQ=255)                                                                                                                                                                                      ttttaatccggg‑ggTGTAACC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCAAATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTAA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTg                                                 <  2:16942/231‑1 (MQ=255)                                                                                                                                                                                          tCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACCCGAAGGAATTAAAATGCCTGGAAAGGGACAAATGAAACATGGCAGTTTGACCCTGCAgg                                             >  1:692275/1‑239 (MQ=255)                                                                                                                                                                                                 c‑ggcgTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCCGGTTtaga                                       >  1:562094/1‑238 (MQ=255)                                                                                                                                                                                                   ggtgTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGaa                                      <  1:1037273/235‑1 (MQ=255)                                                                                                                                                                                                        aaGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCt                                <  1:341841/239‑1 (MQ=255)                                                                                                                                                                                                               ccATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTACTT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAAAACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACAAGCAAGTTTGACTCAGAAGGATGGAAACAACCCGaaaaa                          >  2:330714/1‑237 (MQ=255)                                                                                                                                                                                                                cttAAAAATCTAGGCAAAAAGC‑TTATTTTTAGCTTTAATT‑AACTATCCACAGCATTGTTCTG‑CT‑GATAATAACAACAAAGGTGGGTGGTTTATTATTATGCACTGTA‑ATACTCCAGGTA‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACAt                         <  2:236999/237‑1 (MQ=255)                                                                                                                                                                                                                  tAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACAGTCCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACaac                        >  1:433965/1‑236 (MQ=255)                                                                                                                                                                                                                          cAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCAt                                                                                                                                  <  2:691779/124‑1 (MQ=255)                                                                                                                                                                                                                          cAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCAt                                                                                                                                  >  1:691779/1‑124 (MQ=255)                                                                                                                                                                                                                              cAAAAAGT‑TTTTTTTTAGCTTTAATG‑AACTAACCACATCATTTTTCTG‑CT‑GCTAATAAAAACAAAGGTCAGTGGTTTATGATTATTCACTGTA‑ATACAACAGGTG‑AATACCTTCTTCCAGGCAGTCAGTGTTGGACAGAGAACACGCGAAGGAATTAAAAGTCCTGGAACGGGAAAAATGAAACATGTTAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTgg           <  2:76166/239‑1 (MQ=255)                                                                                                                                                                                                                                       ttATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTTGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTaa   <  1:962440/239‑1 (MQ=255)                                                                                                                                                                                                                                        tATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTaaa  >  1:714427/1‑239 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              TTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAA‑TGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAA‑CTTCGGCAAA‑GTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGC‑GGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGC‑TTATTTTCAGCTTTAATT‑AACTAACCACATCATTGTTCTG‑CT‑GCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTA‑ATACACCAGGTG‑AATTCCTTCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAA  >  NC_000913/1519887‑1520350

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