BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,774,207	C→T	6.9%	intergenic (+13/‑217)	mtlA → / → mtlD	mannitol‑specific PTS enzyme: IIA, IIB and IIC components/mannitol‑1‑phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,774,207	0	C	T	6.9%	2563.0 / 116.2	786	intergenic (+13/‑217)	mtlA/mtlD	mannitol‑specific PTS enzyme: IIA, IIB and IIC components/mannitol‑1‑phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand): ref base C (420/310); new base T (33/21); total (455/331)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-01

Only 100 of 820 total aligned reads displayed.

TGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAACT  >  NC_000913/3773935‑3774479

                                               gtcgtGGTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTTCCCGTCTGGGGAGTGGTTTTGCTGCCCCTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                                  <  2:132538/274‑1 (MQ=255)

                                                                           tGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttg                                                                                                                                                                                                       >  2:306463/1‑273 (MQ=255)

                                                                                               gaTGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGtcg                                                                                                                                                                                 >  2:21202/1‑273 (MQ=255)

                                                                                                   aCATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGcc                                                                                                                                                                              <  2:112770/274‑1 (MQ=255)

                                                                                                           cGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCAccc                                                                                                                                                                                                                                                                                  <  2:79856/166‑1 (MQ=255)

                                                                                                                               tgcCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCctct                                                                                                                                                                                                                                                                               >  8:1289496/1‑149 (MQ=255)

                                                                                                                               tgcCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCctct                                                                                                                                                                                                                                                                               >  7:195756/1‑149 (MQ=255)

                                                                                                                                gcCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCcctc                                                                                                                                                                                                                                                                              <  7:736644/149‑1 (MQ=255)

                                                                                                                                  ccGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCca                                                                                                                                                                                                                                                                            >  6:358468/1‑149 (MQ=255)

                                                                                                                                   cGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGc                                                                                                                                                                                                                                      >  4:99456/1‑186 (MQ=255)

                                                                                                                                   cGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCAATGGATGATGAGT‑CCGTCATCGAGCGTCTGGCCCACACCACCAGCGTGGATGAAGTGCTGGAACTACTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCAACATGGAGAAGGTGGGGTTACTTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCCCAATTTGTTGAATTTGCACGTTCTTGTAGGCAGGATGAGGCGCTTACGCCGCAGTCGgcgtcgcc                                                                                                                                             >  4:107756/1‑270 (MQ=255)

                                                                                                                                     tAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCtctc                                                                                                                                                                                                                                                                              >  4:379985/1‑144 (MQ=255)

                                                                                                                                           cGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAtt                                                                                                                                                                                                >  4:289595/1‑216 (MQ=255)

                                                                                                                                            gAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGCAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAg                                                                                                                                                                                                                                                                  >  7:1298977/1‑149 (MQ=255)

                                                                                                                                              gCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGt                                                                                                                                                                                                                                                                >  8:358150/1‑149 (MQ=255)

                                                                                                                                                acaTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGGAGAGGGTgg                                                                                                                                                                                                                                                              <  5:151246/149‑1 (MQ=255)

                                                                                                                                                  aTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgggg                                                                                                                                                                                                                                                            <  7:881669/149‑1 (MQ=255)

                                                                                                                                                   ttacggtattacccagcctaaccaatGCACTGGGTAATGAGT‑CCGTCATCGGGCGTCTGGCACACACCCCCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGA                                                                                                                             <  4:13522‑M1/263‑1 (MQ=255)

                                                                                                                                                         ttATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTg                                                                                                                                                                                                                                                     >  5:973677/1‑149 (MQ=255)

                                                                                                                                                            tCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCt                                                                                                                                                                                                                                                  <  6:75198/149‑1 (MQ=255)

                                                                                                                                                            tCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCAAACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATGTGCAGGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGCGGGTTTgg                                                                                                                     >  2:335784/1‑274 (MQ=255)

                                                                                                                                                               ccAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGAt                                                                                                                                                                                                                                               >  5:435563/1‑149 (MQ=255)

                                                                                                                                                               ccAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGAGGGGTTAATTGCCTGAt                                                                                                                                                                                                                                               >  8:674965/1‑149 (MQ=255)

                                                                                                                                                                 aGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATgc                                                                                                                                                                                                                                             <  8:497121/149‑1 (MQ=255)

                                                                                                                                                                   catGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCATGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATgcgc                                                                                                                                                                                                                                           <  8:1253318/147‑1 (MQ=255)

                                                                                                                                                                    cTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCt                                                                                                                                                                                                                                          >  5:1380907/1‑149 (MQ=255)

                                                                                                                                                                         cAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCACCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCt                                                                                                                                                                                                                                     >  5:1207655/1‑149 (MQ=255)

                                                                                                                                                                          aaTGCTCTGAATGTTGGGT‑CCGTCATCGGGCGTCTGGCAAACACCACCAGCGTGGAGGAAGTGCTGGAACTGCTGGCAGGTCTTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCt                                                                                                      <  2:147125/275‑1 (MQ=255)

                                                                                                                                                                          aaTGCACGGGATGTGGTGTCCCGTC‑TGGAGCGTCGGGCCAACACCACCCGCGTGGATGAAGTGCTGGAACTGCTTGCAGGTCGTAAGTAATCCAATCCCCCCCTCTCCACATGGGGAAGGTGGGGTTAATTGCCTGATGCGCTACTCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGCAAAGCCt                                                                                                      <  4:335784/275‑1 (MQ=255)

                                                                                                                                                                           aTGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTa                                                                                                                                                                                                                                   >  5:1074724/1‑149 (MQ=255)

                                                                                                                                                                            tGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                                  <  5:1277625/149‑1 (MQ=255)

                                                                                                                                                                              cACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                                                <  5:453307/149‑1 (MQ=255)

                                                                                                                                                                              cACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                                                <  7:277117/149‑1 (MQ=255)

                                                                                                                                                                               aCTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAg                                                                                                                                                                                                                               >  6:397119/1‑149 (MQ=255)

                                                                                                                                                                                  ggATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGcc                                                                                                                                                                                                                            >  6:947074/1‑149 (MQ=255)

                                                                                                                                                                                   gatgatGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAg                                                                                                                                                                                                                                                                  <  7:1355587/110‑1 (MQ=255)

                                                                                                                                                                                     tgatgaGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                                                                                                                                                         >  8:1086131/1‑149 (MQ=255)

                                                                                                                                                                                       atgaGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAg                                                                                                                                                                                                                       <  7:512480/149‑1 (MQ=255)

                                                                                                                                                                                         gaGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGa                                                                                                                                                                                                                     >  6:90385/1‑149 (MQ=255)

                                                                                                                                                                                         gaGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGa                                                                                                                                                                                                                     <  6:505297/149‑1 (MQ=255)

                                                                                                                                                                                           gT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGc                                                                                                                                                                                                                                      >  2:273685/1‑130 (MQ=255)

                                                                                                                                                                                              ccGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACAt                                                                                                                                                                                                                                                                         >  2:379602/1‑93 (MQ=255)

                                                                                                                                                                                                   aTCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGa                                                                                                                                                                                                                     <  1:11246/140‑1 (MQ=255)

                                                                                                                                                                                                   aTCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATcaca                                                                                                                                                                                                            >  8:569316/1‑149 (MQ=255)

                                                                                                                                                                                                   aTCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATcaca                                                                                                                                                                                                            >  7:670005/1‑149 (MQ=255)

                                                                                                                                                                                                    tCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACaa                                                                                                                                                                                                           >  8:1031972/1‑149 (MQ=255)

                                                                                                                                                                                                    tCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACaa                                                                                                                                                                                                           <  6:238521/149‑1 (MQ=255)

                                                                                                                                                                                                      gAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAAtt                                                                                                                                                                                                         >  8:1343063/1‑149 (MQ=255)

                                                                                                                                                                                                         cGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgt                                                                                                                                                                                                      <  5:515462/149‑1 (MQ=255)

                                                                                                                                                                                                         cGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgt                                                                                                                                                                                                      <  6:1263838/149‑1 (MQ=255)

                                                                                                                                                                                                           tCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgttg                                                                                                                                                                                                    <  5:578321/149‑1 (MQ=255)

                                                                                                                                                                                                             tGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGaa                                                                                                                                                                                                  <  6:1019697/149‑1 (MQ=255)

                                                                                                                                                                                                              ggCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgggg                                                                                                                                                                                                                                                            <  6:1353790/90‑1 (MQ=255)

                                                                                                                                                                                                              ggCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAt                                                                                                                                                                                                 <  5:637316/149‑1 (MQ=255)

                                                                                                                                                                                                                cacacaCCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAttt                                                                                                                                                                                               <  8:1454563/149‑1 (MQ=255)

                                                                                                                                                                                                                aacacacacCCCGGGTTGATGAAGGGCGGNNACTGCTGGGAGGTCTTAATTAATCCAATCCNCCCCTCTCCAAAGGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACANTTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGt                                                                 <  3:44859/274‑1 (MQ=38)

                                                                                                                                                                                                                 acacacCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCtctc                                                                                                                                                                                                                                                                              >  7:463353/1‑69 (MQ=255)

                                                                                                                                                                                                                 acacacCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTg                                                                                                                                                                                              >  8:1152610/1‑149 (MQ=255)

                                                                                                                                                                                                                   acacCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCa                                                                                                                                                                                            <  5:930375/149‑1 (MQ=255)

                                                                                                                                                                                                                     accaccAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgggg                                                                                                                                                                                                                                                            >  5:1170788/1‑83 (MQ=255)

                                                                                                                                                                                                                     accaccAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgggg                                                                                                                                                                                                                                                            <  6:1170788/83‑1 (MQ=255)

                                                                                                                                                                                                                         ccaGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCt                                                                                                                                                                                      >  5:702179/1‑149 (MQ=255)

                                                                                                                                                                                                                                ggATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGc                                                                                                                                                                               >  6:615694/1‑149 (MQ=255)

                                                                                                                                                                                                                                   tGAAGTGCTGGAACTGCTGGCAGTTCTTAAGTAAACCAATCCCACACTCTCAACATGAAGAAGGTGTGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGAATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTa                                              <  1:44054/275‑1 (MQ=255)

                                                                                                                                                                                                                                        tGCTGGAACTGTTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAg                                                                                                                                                                       <  7:946773/149‑1 (MQ=255)

                                                                                                                                                                                                                                        tGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAg                                                                                                                                                                       <  7:1237802/149‑1 (MQ=255)

                                                                                                                                                                                                                                         gCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAgg                                                                                                                                                                      <  7:1153514/149‑1 (MQ=255)

                                                                                                                                                                                                                                           tGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGgcg                                                                                                                                                                    >  7:1420433/1‑149 (MQ=255)

                                                                                                                                                                                                                                           tGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGgcg                                                                                                                                                                    >  8:296333/1‑149 (MQ=255)

                                                                                                                                                                                                                                            ggAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGgcgc                                                                                                                                                                   >  7:331532/1‑149 (MQ=255)

                                                                                                                                                                                                                                            ggAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGgcgc                                                                                                                                                                   >  6:172967/1‑149 (MQ=255)

                                                                                                                                                                                                                                            ggAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGgc                                     >  4:202727/1‑275 (MQ=255)

                                                                                                                                                                                                                                              aaCTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCtt                                                                                                                                                                 <  7:983352/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  gctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcgc                                                                                                                                                             <  5:283138/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  gctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACAGGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcgc                                                                                                                                                             <  5:864661/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   ctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcgcc                                                                                                                                                            <  8:823044/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   ctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcgcc                                                                                                                                                            >  6:83281/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        aGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATc                                                                                                                                                       >  5:1430079/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           tCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCgg                                                                                                                                                    >  5:496821/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           tCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCgg                                                                                                                                                    >  8:61141/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            cGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGgc                                                                                                                                                   <  7:1067569/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             gtaagtaaTCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGgcg                                                                                                                                                  >  6:1417032/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               aagtaaTCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCt                                                                                                                                                <  6:1063821/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                agtaaTCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCAAAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTg                                                                                                                                               >  5:1324630/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 gtaaTCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACTTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGc                                                                                                                                              >  8:91204/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    atccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGcccc                                                                                                                                           >  6:685014/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCct                                                                                                                                          >  8:759610/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAACAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGgc                                     >  4:245771/1‑250 (MQ=255)

                                                                                                                                                                                                                                                                      ccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCctc                                                                                                                                         >  5:376829/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       caatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCctct                                                                                                                                        >  6:446473/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                         accccACCTTCTCCATGTGGAGAGGGTGGGATTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTctcc                                                                                                                                      >  6:301615/3‑149 (MQ=25)

                                                                                                                                                                                                                                                                            ccACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTctccgtg                                                                                                                                   >  6:1455296/1‑146 (MQ=255)

                                                                                                                                                                                                                                                                             cACCTTCTCCATGTGGAGAGGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCAc                                                                                                                                  >  5:23038/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                             cACCTTCTCCATGTGGAGAGGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCAc                                                                                                                                  >  7:595993/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                             cACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTctccgtga                                                                                                                                  <  5:817914/149‑5 (MQ=255)

                                                                                                                                                                                                                                                                             cACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCAc                                                                                                                                  >  7:791963/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                             cACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCAc                                                                                                                                  >  8:139738/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                               cccTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGG                                                                                                                                >  8:458081‑M1/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                                ccTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGA                                                                                                                               >  7:1014352‑M1/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                                 ctctCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAG                                                                                                                              >  6:548601‑M1/1‑149 (MQ=255)

TGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGT‑CCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAACT  >  NC_000913/3773935‑3774479

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`