Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,101,635	A→T	22.8%	L17* (TTA→TAA)	csgG ←	curli production assembly/transport outer membrane lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,101,635	0	A	T	22.8%	51.8 / 2.4	38	L17* (TTA→TAA)	csgG	curli production assembly/transport outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (8/21);  new base T (7/2);  total (15/23)
Fisher's exact test for biased strand distribution p-value = 1.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

CGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATATTC  >  NC_000913/1101413‑1101873                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    cGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTTAGCAATCGCCCCGTaa                                                                                                                                                                                                                                     >  2:15072/1‑239 (MQ=255) cGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGGTTGGCGGCTTTTTTAGGGGGGGGGGGTatagaaccccccagcaa                                                                                                                                                                                                                                     >  1:999489/1‑223 (MQ=255)                    tGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGGTTAGCATCCGCTCCGGAAACACACGGC‑AACC‑caaaaa                                                                                                                                                                                                                >  2:401215/1‑237 (MQ=255)                    tGACCAGCATTGCCGCGGCGCTTTGCGGAACAGCAGTGGCGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCTTTAATGTCGGACTGGCGGCTTCTTTAGGCGGGGCGGGTTAGAATCCCGCCCGGAACATTACGGC‑AAAC‑aaaaac                                                                                                                                                                                                                >  2:265736/1‑237 (MQ=255)                      aCCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGGTTAGGATCCCCTCCAGAAAATGAAGGG‑CACC‑CAAAGaaa                                                                                                                                                                                                              >  2:707882/1‑239 (MQ=255)                                ccGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCACACCCCCGCTCACATTGCGGC‑AACa‑aaaaaaaaaaaaaccgcg                                                                                                                                                                                                    >  1:366106/1‑232 (MQ=255)                                     gcgcTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTTAACATACGCCCCATAAACAGTGGGC‑AACCAAAAAGAAAAAAG‑GCTGCATgtg                                                                                                                                                                                               >  2:1004871/1‑237 (MQ=255)                                        gtttGCGTTACAGCATGGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGGTTAGTCCTGAATGTTGTATTCC‑GAAACAAAGATTGTACCCGTCGGCGTTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGTATTAATGTCGGTCTGGCGGATTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGatta                                                                                                                                                                                             <  2:345379/237‑1 (MQ=255)                                                 accgcgtttgataattttttttccGGTTTGGGTTTAAATTGCCCGGTTTCGGCCTGAACTTTGTATCCC‑GAAACAAAGCTTTTACCCGTTGGCTGTGGCAGATGGTTCAAATTTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCGTCTTTATGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTa                                                                                                                                                                                    <  1:326710/216‑1 (MQ=255)                                                     cagTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTTAGCATTCGGTCCGTAAACAGAAGGC‑AAAC‑ACAAGAAAAAAGCGCTGCCTGATATTTTTTCTTATGA‑Ag                                                                                                                                                                              >  2:683937/1‑239 (MQ=255)                                                     cagTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAACATCCGCTCCGTAACATGACGGC‑CACC‑CAAAGAAATAAGGGCTGGCTGATTTTTTTCCTTTTGA‑Ag                                                                                                                                                                              >  2:387518/1‑239 (MQ=255)                                                      tggtggtaattttcTTGCCGGGAAGGTTTTAATTTGCCCGGTTTTGTCC‑GGATGTTTGTTACCGGAATAAAAGTGTTTCCCCGTGGGCGCTGGCGGTTGGGCCAAATCTTTGTAGCTCTGAGCCCTGGCAATTATGGTCGTTCGGGCGTCTTCTTTA‑GCGGGGCGGTTTAGCCTCCGCCTCGTAACAAGACGGC‑AACC‑AAAAGAAATAAGCGCTGCAATATTATTTTCCTTATGA‑AGCt                                                                                                                                                                            <  1:623701/226‑1 (MQ=255)                                                          gagaAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGGTAAGCATCCCCTCAGTAACAAGACGGC‑GACC‑CAAAGAAAAAAAAGACGCCGGATTAtttttccttttataaaggg                                                                                                                                                                          >  1:854825/1‑224 (MQ=255)                                                               gTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCCCCCAGTAACATGACGGC‑ACCC‑AAAAGAAAAAAACGGTGCCTAATTATTTTTTTTATTa‑aaaggggggtta                                                                                                                                                                    >  1:318731/1‑228 (MQ=255)                                                                                          cccGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTTATCATCCGCCCCTGAAAAATACGGCGAACC‑AAAAGAAAAAAGAG‑TCCAATATTTTTTTTCTTTTTA‑ACAACGGGCGTTAAAAAAAGGATGAG‑TTTTTTTTTaaa                                                                                                                                         >  1:673692/1‑238 (MQ=255)                                                                                          cccGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCCCAGTAACATTACGGG‑AAAC‑CAAAAAAAAAAACGCTGCATGGTTATTTTCCTTTTGA‑AGCTGGGGC‑TTTAAAAACGGTTTGG‑TTATTTTTTAAAc                                                                                                                                        >  2:294408/1‑238 (MQ=255)                                                                                            cGGTTTCGTCCTGAATGTTTTATACC‑GGAACAAAGGTTTTCCCCGTCGGGGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCg                                                                                                                                     <  1:401215/239‑1 (MQ=255)                                                                                                ttctcccTGAATTTTGTAATCC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTTAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAAc                                                                                                                                 <  1:294408/234‑1 (MQ=255)                                                                                                     ccTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCATGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGa                                                                                                                            <  2:879225/239‑1 (MQ=255)                                                                                                     ccTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTgg                                                                                                                             <  2:854825/238‑1 (MQ=255)                                                                                                             ttgtATACC‑GAAACAAAGATTTTACCCGGCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGtt                                                                                                                                                         <  2:42642/203‑1 (MQ=255)                                                                                                                   aCC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGGTAAACATCCCCTCAGTAAAATGGCGGC‑ACAC‑AAAAGAAATAAGAGATGCAGGAGTTTTTTTCTTTTGA‑AACTGGGGG‑TTAAAAAACGGGTGAGATTTTTTTTTAAACCCAAACCCTgaggggtcagggtttt                                                                                                               >  1:169644/1‑224 (MQ=255)                                                                                                                        aaccAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGGTCCGGGACATGTCGGC‑ACAC‑CAAAGAAAAAAGCGCTTGCTCACTATTTTCCTTTTGACAGCAGGGAC‑TTAAAAAAGGGTGGGG‑TT‑TTTTTTTAACACGAAACCAgggaggttagggatgtgtggg                                                                                                          >  1:431878/4‑220 (MQ=255)                                                                                                                        aaaCAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCAACCGCTCAGTAACATGACGGG‑CACC‑AAAAGAAAAAAACGCTGCATGATTATTTTCCTTATGA‑AGCCGGGGG‑TTAAAAAACGGTTGAG‑TTATTTTTTTAACCCGAAACCTGGATGGGCGGGGGTTGTCCgg                                                                                                          >  2:291488/1‑239 (MQ=255)                                                                                                                                                                          ctctGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAat                                                         <  2:366106/238‑1 (MQ=255)                                                                                                                                                                            ctGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATc                                                      <  2:431878/239‑1 (MQ=255)                                                                                                                                                                                    ggggCATTTATTTCGGTCTGGCGGCTTTTTTAGGCGGGGCGGTTAAGCATCCGCTAAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTTCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAatat                                              <  1:387518/237‑1 (MQ=255)                                                                                                                                                                                    gAGGCATTAAGGTCGGTCTGGAGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAatat                                              <  2:684369/239‑1 (MQ=255)                                                                                                                                                                                        cATTAATGTCGGTCTGGCGGCTTCTGTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATc                                                      <  2:239235/227‑1 (MQ=255)                                                                                                                                                                                          ttAAGGTCGGTCTGGAGGGTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACAGGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTGTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATGGACCATCGCGGTTGGCAATATCGACAATATAATCg                                         <  2:158712/238‑1 (MQ=255)                                                                                                                                                                                                     cttgcgggtttgtgaGGGGGGGCGGTTAAGCATCCTGTCAGTAACAGTACGGA‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑GTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGATCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAAACGTTGGTCACCATg                             <  1:457800/225‑1 (MQ=255)                                                                                                                                                                                                        ggggcTTCTTTAGGCGGTGGGGTTAGGCTTCCGTTCGGTAACATGGCGGC‑AACC‑AAACTGAATAAGCGCTTCAGTCTTTCTTTCCTTATGA‑AGCTGGGTA‑TTTAAAATCTGTTGAT‑TTATTTTGTAAACCCGAAGCCTGGATGGTCGAGGTTTTGCCGGTTTTATGATCGGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTCCTCACCATGCgg                          <  1:334844/237‑1 (MQ=255)                                                                                                                                                                                                          ggcTTTTTTGGGGGGGGGGGTTAAGCATCCGGTTAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGcc                        <  1:19335/239‑1 (MQ=255)                                                                                                                                                                                                                 ttAGGGGGGGGTGTTAAGCATCCGCACAGTAACATGACGGA‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTa                 <  2:389109/239‑1 (MQ=255)                                                                                                                                                                                                                  tAGGGGGGGCGGTTTAGCATCCGTTCAGTAAACTTACGGC‑AACC‑AAAAGAAATAAACCGTTCATGATAGTTTTCCTTTGGA‑AGGTGGGGC‑GTAAAAAGCGGTTTAG‑TTCTTTTGTAATCCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTAc                <  1:265736/239‑1 (MQ=255)                                                                                                                                                                                                                        gggcggTTAATCATCCGCTTAGTAACAAGACGGC‑AAAC‑AAAAGAAATAAGCGCTGCATGAATTTTTTCCGTATGA‑AACTGGGGC‑TTAAAAATCGGTTTAG‑TTTTTTTGTAAACCTGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTAt          <  1:804900/239‑1 (MQ=255)                                                                                                                                                                                                                            ggTTAAGCATCCGCTCAGTAACATGACGGA‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTGCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTaa       <  2:608859/238‑1 (MQ=255)                                                                                                                                                                                                                               tAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATAtt   <  2:547875/239‑1 (MQ=255)                                                                                                                                                                                                                               tAAGCACCCCCTTAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATAtt   <  1:291488/239‑1 (MQ=255)                                                                                                                                                                                                                                aaGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATATTc  <  2:267169/239‑1 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    CGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACC‑GAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGC‑AACC‑AAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGA‑AGCTGGGGC‑TTAAAAATCGGTTGAG‑TTATTTTGTAAACCCGAAACCTGGATGGTCGAGGTTTGTCCGGTTTTACGATCTGTCACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATATTC  >  NC_000913/1101413‑1101873

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