Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,104,208	A→G	17.9%	G86G (GGA→GGG)	csgB →	curlin nucleator protein, minor subunit in curli complex

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,104,208	0	A	G	17.9%	166.6 / 8.2	69	G86G (GGA→GGG)	csgB	curlin nucleator protein, minor subunit in curli complex
Reads supporting (aligned to +/- strand):  ref base A (17/38);  new base G (11/2);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 9.58e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCGACTTTTAAATCAATCCGATG  >  NC_000913/1103970‑1104436                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          ttatGATTTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGga                                                                                                                                                                                                                                           <  1:718908/239‑1 (MQ=255)        ttttaaaaTTCTGGTTTCGCCTGGGATTGAAGCCGCAGGAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATtat                                                                                                                                                                                                                                     <  2:1006801/232‑1 (MQ=255)             cAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCtt                                                                                                                                                                                                                              <  1:331421/239‑1 (MQ=255)               attcTGGTTGCGCCTGGGATTGCAGCCGCAGCAGGTTTTGATTTTGCTAATTCAGAATATAATCTCGCGGTAAATGAATTGAGTAAGTTTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGc                                                                                                                                                                                                                            <  2:534910/236‑1 (MQ=255)               attcTGGGTGCTCCTGGGATTGCCGGCGCAGAAGGTTTTGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTCTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGc                                                                                                                                                                                                                            <  2:435583/236‑1 (MQ=255)               attatgggttgcgcTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGc                                                                                                                                                                                                                            <  2:248422/226‑1 (MQ=255)               aTACTGGGGGCGCCTGGGATTGCAACCGCAGCAAGTTATGATTTAGCTAAATCATAATTGAACTTCGAGGTAAATGAATTGAGTAAGGCTTCATTTAATCAGGAAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGATCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGc                                                                                                                                                                                                                            <  1:502070/239‑1 (MQ=255)                       tGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑GTTTTTCGCAAGAAGGGAGTAGCAA‑CCAGGCAAAAATTGACCAGGAAGGAGATTAT‑AACCTTGGATATAATg                                                                                                                                                                                                                    >  1:688762/1‑239 (MQ=255)                           ccTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGATGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCCGGCAAAGCTTGACCAGACCGGGGATTAT‑TACCCTgactacacttgtcc                                                                                                                                                                                                                >  1:162115/1‑226 (MQ=255)                               ggATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAAATTTTTGCGG‑GTGTTGGGCAAGAAGGTAGTTGGAACCCGGACAAAAATTGTCCA‑ACAAGGGGTTAT‑AAACTTGCATATAATTATCaaggg                                                                                                                                                                                                            >  2:854352/1‑235 (MQ=255)                                  ttGCAGCCGCAGCAGGTTATGATTGAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCggg                                                                                                                                                                                                          <  1:590784/238‑1 (MQ=255)                                       gCCGCAGCAGGTTTGGATTTAGGTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGGCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTg                                                                                                                                                                                                     <  2:225899/238‑1 (MQ=255)                                          gcagcaGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGCGAGGCTCAAAACTTTTGGCGG‑TTGTTGTGCCAGAAGGGAGAGGCAA‑CCGGGCAACGATTTACCAGAGAGGAGATTAA‑AAACATCCCTAAATTGaaaaaacaggaaggggcaa                                                                                                                                                                                                 >  2:587269/1‑221 (MQ=255)                                                  ttATGAATTAGCTAAATCAGAATCTAACTTCGCGGTAAATGAATTGAGTAAGTCTTAATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGGAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGGCAGGGGATTAT‑AACCTTGCATATATTGATCAGGCGGGCATTGCCAACGATGc                                                                                                                                                                                           <  1:820117/237‑1 (MQ=255)                                                          tAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGGAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGGGATTAT‑AACCTTGCATATTATGATCCGGCGGGCAGGGCCAAAGATGCCAGTAtttt                                                                                                                                                                                  >  1:991890/1‑237 (MQ=255)                                                           aGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCACAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTTGCGG‑GTGTTGCGCAAGAAGGTAGTAGGAA‑ACGGGCAAAGATTGTCCAGACAGGAGATTAT‑AAACTTTCATCTATTTATCAGGCGGGGCGTGCCAACGAAGCCAGAATTTca                                                                                                                                                                                 >  1:991351/1‑237 (MQ=255)                                                              tAATTCAGAATATAACTTCGCGTGAAATGAATTTAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAAGAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAg                                                                                                                                                                             <  1:1080352/239‑1 (MQ=255)                                                              tAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCaa                                                                                                                                                                              <  1:774344/238‑1 (MQ=255)                                                                 ttCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCCCAAAACTTTTGGCGG‑TTGTTTGGCAAGAAGGGAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGGGATTAT‑AACCTTGCATATTTTGATCAGGCGGGCAGGGCCAAAGAGGCCAGAATTTCCCAAGGTg                                                                                                                                                                          >  1:581691/1‑239 (MQ=255)                                                                      aaTATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATta                                                                                                                                                                                                                                      <  1:386444/175‑1 (MQ=255)                                                                       atatAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTa                                                                                                                                                                     <  1:718005/237‑1 (MQ=255)                                                                                cgcgGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AAACTTGCATATATTGGTCAGGCGGGCAGTGCCAACCATGCCAGGAATTGCCAAGGTGC‑TTATGGTGATACtt                                                                                                                                                          >  1:740759/1‑238 (MQ=255)                                                                                     ttaaTGATTTGTGTAAGTTTTCATTTAATCAGGCAGCCATAATTGGGCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACTTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGc                                                                                                                                                         <  2:416500/233‑1 (MQ=255)                                                                                               gAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGGAAGAAGGGAGTAGCAA‑ACAGGCAAAGATTGACCAGACAGGGGATTAT‑AAACTTTCATATATGGATAAGGCGGGAGGGGCCAACGATGCCAGTATTTCGCAAGGGGC‑TTATGGTATTaatcggttggttttccctt                                                                                                                                           >  1:416500/1‑221 (MQ=255)                                                                                               gAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGa                                                                                                                                          <  1:854352/239‑1 (MQ=255)                                                                                                      tCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTACCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGc                                                                                                                                                         <  1:772992/218‑1 (MQ=255)                                                                                                        ttCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGGGG‑TTGTTTCGCCAGAAGGTAGGAGCAA‑ACGGGCAAAGATTTGCCAGACAGGGAGTTAT‑AAACATGCATATATTGAACAGGGGGGCAGAGGCAACGATGTCCGTATTTTGCCAAGGGG‑TTATGGTAATAATGA‑GATAAAAATCCAaaaaagatatg                                                                                                                                 >  2:453702/1‑232 (MQ=255)                                                                                                                aTCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTTGCGG‑GTGTTTGGCAAAAAGGGAGTAGCCA‑CCGGGGAAAGATTGACCAGAAAGGGGATTATAAACCTGGCATA‑ATTGATCCAGCGGGCAGGGGCAACGATTCCCGCATTTCGTCAGTGGCTTTATGGGAA‑AAGCC‑GATGATTATCCACAAA‑gggtttcgggaaaaa                                                                                                                        >  2:42579/1‑226 (MQ=255)                                                                                                                aTCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑GTGTTGCGCAAGAAGGGAGGAGGAA‑CCGGGCAAAAATTTACCAGACAGGGGGTTAT‑AACCATGCATATATTTATCAGGCGGGGAGTGGCCACCATGCCAGCATTTCGCCAAGGGG‑TTATGGGGAATCTGC‑GCGGATTATccccgtcttttttaaacaacta                                                                                                                         >  2:772992/1‑219 (MQ=255)                                                                                                                  cAGGGAGCCATAATTGGGCAAGCGGGGACTAATAATAGGGCTCAGTTACGGCAGGGAGGGTCAAAACTTTGGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGc                                                                                                                      <  2:162115/239‑1 (MQ=255)                                                                                                                       acccATAATGGGTCAAGCTGGGCCTAATAATAGTGCGCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGT‑ATAAAGCAAata                                                                                                                 <  2:991351/236‑1 (MQ=255)                                                                                                                          cATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAata                                                                                                                 <  1:453702/236‑1 (MQ=255)                                                                                                                           aTAATTGGTCAAGGTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTa                                                                                                              <  1:316927/238‑1 (MQ=255)                                                                                                                            ttaTTTGTCAAGGTGGGACTAATAATAGGGCCCAGGTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTaca                                                                                                            <  2:675112/237‑1 (MQ=255)                                                                                                                            ttaTTGGTCAATGTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTAAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTaca                                                                                                            <  2:255173/237‑1 (MQ=255)                                                                                                                            tAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTaca                                                                                                            <  1:42579/239‑1 (MQ=255)                                                                                                                            tAATTGGTCAAGCGGGGATTAAAAATGGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTaca                                                                                                            <  2:333720/239‑1 (MQ=255)                                                                                                                            tAAGTGGGCAAGCTGGGACTAATAATAGTGCTCAGTGACGGCAGGGAGGCTCAAAACGTATGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCTAGAAA‑GGTTCTGGTAATAAAGCAAATATTaca                                                                                                            <  1:215688/239‑1 (MQ=255)                                                                                                                                   tCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTTC‑GATGATTTTCCAGAAA‑GGGTCTGGTAATAAAggcaatttttacccgtaat                                                                                                     >  1:326573/1‑221 (MQ=255)                                                                                                                                      aGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGGGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGGATTTTGCAAGGGGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGGAATAAAGCAAATATTACACAGT‑ATGGTa                                                                                                 >  1:625596/1‑239 (MQ=255)                                                                                                                                                 ataataGTGCTCAGTTACGGCAGGGAGGCTCAAAACATTTTGCGG‑GTGGTTCGCCCAGAGGGAATAGCAA‑CCGGGCAAAAATTGACCAGACAGGCGGATTT‑GACACTGCATATATTTAACTAGGGGGCAGGGGCCACAAGTCCCGCCTTTCACACTGTGG‑TTTTGGTCAATCTGCAGCAGATTA‑CCAGAAA‑AGTTTTAGGAATAAAAGAAAAAATACACAGT‑Attagaaaaaaaaaaag                                                                                      >  2:30053/1‑224 (MQ=255)                                                                                                                                                     taGTGCTCAGTTACGGCAGGGAGGATCAAAACATTTGGCAG‑TTGTTGCGCAAGAAGGTAGGAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCaa                                                                                  <  1:203881/239‑1 (MQ=255)                                                                                                                                                        tGCGCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATAtt                                                                                                               <  1:991540/208‑1 (MQ=255)                                                                                                                                                         gCTCAGTTTCGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTTGTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGt                                                                              <  2:154173/239‑1 (MQ=255)                                                                                                                                                          cTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGGGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTTCACAGT‑ATGGTAACCCAAAAACGGCAATTGTa                                                                             >  1:729471/1‑239 (MQ=255)                                                                                                                                                                   ggCAGGAGGGCTTAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAgaga                                                                    <  2:326573/239‑1 (MQ=255)                                                                                                                                                                       gggAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGGAGGAGCAA‑CCGGGCAAAGATTGACCCGACAGGAGATTAT‑AACCTTGCATATATTGATCAAGCGGGCAGGGCCAACGATTCCACTATTTTCCCAGGGGC‑TTATGGGAAAACTCG‑GCTGGTTTTCCACAAA‑AGTTTTGGGAATAAAACAAAAATTACACCAT‑ATGTGACTCCAAAAACCGCAAATTTAGAGcgaaaaaaaa                                                                >  2:698502/1‑230 (MQ=255)                                                                                                                                                                         gAGGCTAAAAACTTTTGGCGGTTTGTTG‑GCAAGAAGGGAGTAGCAA‑CCGGGCAAAAATTGGACCGACCGGGGGTTTT‑AAACTTGCATATATTGATCAGGCGGGCAGTGGCAAAGATGGCAAGATTTTGCCAGGGGC‑TTTTTGGTATAATTG‑GATTGGTATCCCGAAAGGGTTCTGGAAAAAAACCAAAAATTAC‑CAGTAATGTTCCTAAAAAAAAggaatgttaaggaagaaaaaaaac                                                               >  2:835602/1‑217 (MQ=255)                                                                                                                                                                             cttaaaaCTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATGAGGCGGGCAGTGCCAACGATGCCAGTATTGCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCaa                                                           <  2:1072418/235‑1 (MQ=255)                                                                                                                                                                             cTCAAGACTGTTGGCGG‑TTGTTGGGCAATAAGGTAGTAGCAA‑GCGGGAAAAGATTGACCAGACGGGGGATTAT‑AACGTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGTTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAACAAATATTACACAGT‑ATGGTACTCAAAAAACGGc                                                                                    <  1:1055342/213‑1 (MQ=255)                                                                                                                                                                             cTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATTGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTTGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGt                                                                              >  1:155679/1‑219 (MQ=255)                                                                                                                                                                             cTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTAGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGt                                                                              <  2:155679/219‑1 (MQ=255)                                                                                                                                                                                       ttGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGACGCAAATGGCTATTcg                                                >  1:441463/1‑239 (MQ=255)                                                                                                                                                                                       ttGGCGG‑GTGGTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGGCAGGGGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTAATTCGCAAGGTGC‑TTATGTGAATACTGC‑GATGATTATCCAGAAG‑AGTTCTGGGAATTAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGCAGGCAAGATACAAACACAAAAGGCTAATcg                                                >  2:1044942/1‑239 (MQ=255)                                                                                                                                                                                                tgttgCGCAAGAAGGGTATAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCAGATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATTCTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGaca                                         <  2:688762/238‑1 (MQ=255)                                                                                                                                                                                                 gtttCGCAAAAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGAcaca                                       <  2:625596/239‑1 (MQ=255)                                                                                                                                                                                                         agaagGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGc                                                                                                                      <  2:466591/153‑1 (MQ=255)                                                                                                                                                                                                         agaagGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGc                                                                                                                      >  1:466591/1‑153 (MQ=255)                                                                                                                                                                                                                   agCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCg                      <  1:587269/238‑1 (MQ=255)                                                                                                                                                                                                                      aa‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTTGACtt                  >  1:70634/1‑239 (MQ=255)                                                                                                                                                                                                                        aacGGGGAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCGACttt                 <  1:698502/237‑1 (MQ=255)                                                                                                                                                                                                                         ccGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCaaaaa                                                                                         >  2:362751/1‑166 (MQ=255)                                                                                                                                                                                                                         ccGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCaaaaa                                                                                         <  1:362751/166‑1 (MQ=255)                                                                                                                                                                                                                         ccGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTAACCAGGAA‑GGGTCTGGTAATAAAGCCAATATTACACAGG‑ATGGTACTCAGAAAACGGCAATTGAAGTGCAGAGACAATCGCACAAAGAACATCGCGGGACACAACGATAATTTCCATTCGACtttt                >  2:380036/1‑239 (MQ=255)                                                                                                                                                                                                                         ccGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATACCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGC‑ATGGTATTCAAAAAACGGCAATTGTAATGCCGAGACAGTCGCAAATGGCCACTCACGTGAAAAAAGGATAATTTCCATTGGaactta                >  2:441028/1‑234 (MQ=255)                                                                                                                                                                                                                            ggCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCGACTTTTaaa             >  1:241328/1‑239 (MQ=255)                                                                                                                                                                                                                             gCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACCCAGT‑ATGGTACTCAAAAAACGGCAATTGTAGCGCAGAGACAGTCGCCAATGGCTATTCGCGTGCCACAACGTTAATTTCCATTCGACTTTTaaa             >  1:1015033/1‑238 (MQ=255)                                                                                                                                                                                                                              aaaaGGTTGACCCGACAGGAATTTAT‑AACCTTTCATATTTTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCGACTTTTAaat            <  2:441463/237‑1 (MQ=255)                                                                                                                                                                                                                                      gACCGGAAAGGAGATTAT‑AACCTTGCATATTTTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCGACTTTTAAATCAATCCGAt   <  2:991890/239‑1 (MQ=255)                                                                                                                                                                                                                                       aCCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGCTGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCACATGGCTATTCCCGAGACACAACAATAATTTCCATTCGCCTTTTCACTCAATCCGATg  >  1:930723/1‑239 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          TTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAA‑CCGGGCAAAGATTGACCAGACAGGAGATTAT‑AACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGC‑TTATGGTAATACTGC‑GATGATTATCCAGAAA‑GGTTCTGGTAATAAAGCAAATATTACACAGT‑ATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCGACTTTTAAATCAATCCGATG  >  NC_000913/1103970‑1104436

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