BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,823,335	G→C	6.8%	intergenic (+50/‑180)	nadE → / → cho	NAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,823,335	0	G	C	6.8%	702.4 / 31.4	203	intergenic (+50/‑180)	nadE/cho	NAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand): ref base G (92/97); new base C (9/5); total (101/102)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-01

Only 100 of 206 total aligned reads displayed.

AAGCGCCAACGGCCGATCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAAT  >  NC_000913/1823080‑1823601

aaGCGCCAACGGCCGAGCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCg                                                                                                                                                                                                                                                     <  1:189594/279‑1 (MQ=255)

   cgcCAACGGCCGATCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcggc                                                                                                                                                                                                                                                                     >  2:10197/1‑260 (MQ=255)

    gcCGACGGACGATCTGGAAGATGATCGCCCGTCTCTGCCAGATGAAGTGGCACTCGGCCGGACCGATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGGACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAATAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAAGTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCg                                                                                                                                                                                                                                                     <  2:609976/275‑1 (MQ=255)

                    gaagattacggaccttccttGCCAGATGAAGTGGAAATCGGCTGGACCTATGAAAATATCGACGACTAGCTGGAAGGGAAAATCGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGc                                                                                                                                                                                                                                     <  4:392416/263‑1 (MQ=255)

                    gAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAAGATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGc                                                                                                                                                                                                                                     <  2:529919/275‑1 (MQ=255)

                          gatCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATACTTTGCACATATTGTATTGTGCGAAAAAGAGTAATTTTTTACCGCCGGATGCGGCGTGAACNNNTTATTCCACCTATAAAACTATACAAttttcttt                                                                                                                                                                                                                           >  3:266700/1‑274 (MQ=255)

                          gatCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGAAAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAAc                                                                                                                                                                                                                                          <  1:266700/264‑1 (MQ=255)

                                       tGCCAAATGAAATTGGACTCGGCGTGAACTATGACAATATCGACAACTATCTGGAATGGACAAACGTACCCCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACAGTGCAg                                                                                                                                                                                                                   <  4:280126/274‑1 (MQ=255)

                                          caaaTGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACGATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAAGAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGaa                                                                                                                                                                                                           <  1:157873/276‑1 (MQ=255)

                                                  gTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCg                                                                                                                                                                                                                                                     <  2:44994/229‑1 (MQ=255)

                                                                      tGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGc                                                                                                                                                                                                                                     <  1:231874/225‑1 (MQ=255)

                                                                      tGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGc                                                                                                                                                                                                                                     >  3:231874/1‑225 (MQ=255)

                                                                            tatCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAAt                                                                                                                                                                                                                                 >  2:46305/1‑223 (MQ=255)

                                                                                                           ccTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCgg                                                                                                                                                                                                                                                                            >  7:1136439/1‑149 (MQ=255)

                                                                                                             tCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGAt                                                                                                                                                                                                                                                                          >  7:1279346/1‑149 (MQ=255)

                                                                                                             tCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGAt                                                                                                                                                                                                                                                                          >  8:1467737/1‑149 (MQ=255)

                                                                                                             tCAACAGATCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGAt                                                                                                                                                                                                                                                                          >  5:982031/1‑149 (MQ=255)

                                                                                                               aaCAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgc                                                                                                                                                                                                                                                                        >  6:1124240/1‑149 (MQ=255)

                                                                                                                     tCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTg                                                                                                                                                                                                                                                                  <  8:13630/149‑1 (MQ=255)

                                                                                                                         cAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGc                                                                                                                                >  4:54386/1‑275 (MQ=255)

                                                                                                                          aGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGc                                                                                                                                                                                                                                                             >  8:434573/1‑149 (MQ=255)

                                                                                                                            aaCAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCt                                                                                                                                                                                                                                                           <  6:326208/149‑1 (MQ=255)

                                                                                                                               aaTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                                                                                                                                        <  5:1197367/149‑1 (MQ=255)

                                                                                                                               aaTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                                                                                                                                        >  6:1396419/1‑149 (MQ=255)

                                                                                                                                aTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcg                                                                                                                                                                                                                                                                       <  8:1266873/133‑1 (MQ=255)

                                                                                                                                 tAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTc                                                                                                                                                                                                                                                      <  7:434573/149‑1 (MQ=255)

                                                                                                                                 tAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTc                                                                                                                                                                                                                                                      <  7:931727/149‑1 (MQ=255)

                                                                                                                                   gggaATTGGTAACTGAAAACCAAACATAAACGCCTTCCGCAAATTACCGTTTGCGATGATTTCTGGAAAAATTAATAATTTGCACATATTGGAGTGGGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAAGAAATTGCAGAGGTCGATTAGGCCTGACAGTGGTAGAAAGTAAGAAGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTaca                                                                                                                        <  4:507032/271‑1 (MQ=255)

                                                                                                                                   gagaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAGGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                                                                                                                                                                                                    >  8:1339910/1‑149 (MQ=255)

                                                                                                                                    gaatgtgttATCTGAAAACCGAACATAAACTCCGTCCGCAAATTACGTGTTTCGAGCATTTCTGGAAAAAGAAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTGGTTCACGCCGGATGCGGAGTGAACGCCTTATTCGACCTATAAAACTACGCAAATTCAATACATTGCAGGGGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTg                                                                                                                     <  4:654146/267‑1 (MQ=255)

                                                                                                                                    agaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAc                                                                                                                                                                                                                                                   >  6:740820/1‑149 (MQ=255)

                                                                                                                                    agaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAc                                                                                                                                                                                                                                                   <  6:856432/149‑1 (MQ=255)

                                                                                                                                     gaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAcc                                                                                                                                                                                                                                                  <  7:786286/149‑1 (MQ=255)

                                                                                                                                      aaCTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCt                                                                                                                                                                                                                                                           >  2:732050/1‑139 (MQ=255)

                                                                                                                                           gTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTACTCGCTTACCGCCTTGTCCACTGgatagat                                                                                                              >  2:561768/1‑275 (MQ=255)

                                                                                                                                            tATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATaaaa                                                                                                                                                                                                                                           <  7:1287867/149‑1 (MQ=255)

                                                                                                                                              tCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACt                                                                                                                                                                                                                                         >  6:1117909/1‑149 (MQ=255)

                                                                                                                                               cTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTa                                                                                                                                                                                                                                        <  5:1521278/149‑1 (MQ=255)

                                                                                                                                                tGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTAt                                                                                                                                                                                                                                       >  6:1583246/1‑149 (MQ=255)

                                                                                                                                                 gAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATg                                                                                                                                                                                                                                      >  8:24658/1‑149 (MQ=255)

                                                                                                                                                 gAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGTCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCATCTTTTTTAATCGCTTACCTCATGTTACACTTGCTAGAGAACCTGCAtt                                                                                                    >  1:295908/1‑279 (MQ=255)

                                                                                                                                                      acgaacagagagtccGTCCGGCAATTACCGGATTGGATGATGTCTGGAAAAGGCAATAATTGAGAAAGTATGGATTGGTCGAAAAAGAGTAAATGGTTCACGCCGGATGCGGTGTGAACGCTTTATTAGACCGAGAAAACTACAAAAATTAAATAAGTTGCAGTAGTCGAAAAGGCCTGAAGGGTGTAAGAGGTCAGAGGGTGTAACCTTTGACAGCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTc                                                                                                   <  4:189484/262‑1 (MQ=255)

                                                                                                                                                         aaCATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCa                                                                                                                                                                                                                              <  5:1269951/149‑1 (MQ=255)

                                                                                                                                                          aCATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCaa                                                                                                                                                                                                                             <  5:566895/149‑1 (MQ=255)

                                                                                                                                                            aTAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATa                                                                                                                                                                                                                           <  6:70849/149‑1 (MQ=255)

                                                                                                                                                                   ccGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCa                                                                                                                                                                                                                    <  7:1489061/149‑1 (MQ=255)

                                                                                                                                                                    cGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcggc                                                                                                                                                                                                                                                                     >  4:706999/1‑99 (MQ=255)

                                                                                                                                                                        cGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGt                                                                                                                                                                                                               <  5:1124240/149‑1 (MQ=255)

                                                                                                                                                                            aaTTACCGTTTTCGATGATTTCTGGAAAGAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGaa                                                                                                                                                                                                           >  6:195202/1‑149 (MQ=255)

                                                                                                                                                                            aaTTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGaa                                                                                                                                                                                                           <  8:1136439/149‑1 (MQ=255)

                                                                                                                                                                             aTTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAAt                                                                                                                                                                                                          <  7:1330092/149‑1 (MQ=255)

                                                                                                                                                                             aTTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAAt                                                                                                                                                                                                          >  5:1366623/1‑149 (MQ=255)

                                                                                                                                                                                aCCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAgg                                                                                                                                                                                                       <  5:14502/149‑1 (MQ=255)

                                                                                                                                                                                   gTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCt                                                                                                                                                                                                    <  8:79969/149‑1 (MQ=255)

                                                                                                                                                                                         gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAgg                                                                                                                                                                                              >  5:745270/1‑149 (MQ=255)

                                                                                                                                                                                            gatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcg                                                                                                                                                                                                                                                                       <  4:457010/73‑1 (MQ=255)

                                                                                                                                                                                             atTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTa                                                                                                                                                                                          >  7:1403196/1‑149 (MQ=255)

                                                                                                                                                                                                tCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATa                                                                                                                                                                                                         >  7:889227/1‑131 (MQ=255)

                                                                                                                                                                                                tCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATa                                                                                                                                                                                                         <  8:889227/131‑1 (MQ=255)

                                                                                                                                                                                                 cTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCAc                                                                                                                                                                                      <  8:1255972/149‑1 (MQ=255)

                                                                                                                                                                                                    gAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCATGAGTCGAATAGACCTGACAGGCGTAGCACGTc                                                                                                                                                                                   >  8:1221772/1‑149 (MQ=255)

                                                                                                                                                                                                    gAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTc                                                                                                                                                                                   <  8:215344/149‑1 (MQ=255)

                                                                                                                                                                                                       aaaGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                                                                                                                         >  2:348120/1‑76 (MQ=255)

                                                                                                                                                                                                           taataaTTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCg                                                                                                                                                                                                                                                     <  4:348116/76‑1 (MQ=255)

                                                                                                                                                                                                           taataaTTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                                                                                                                                            >  5:705575/1‑149 (MQ=255)

                                                                                                                                                                                                           taataaTTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                                                                                                                                            >  6:1254538/1‑149 (MQ=255)

                                                                                                                                                                                                            aataatTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgtg                                                                                                                                                                           >  6:538059/1‑149 (MQ=255)

                                                                                                                                                                                                                 tttGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAAtt                                                                                                                                                                                                                                >  8:1119986/1‑91 (MQ=255)

                                                                                                                                                                                                                 tttGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAAcc                                                                                                                                                                      <  8:569457/149‑1 (MQ=255)

                                                                                                                                                                                                                    gCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCttt                                                                                                                                                                   <  6:756191/149‑1 (MQ=255)

                                                                                                                                                                                                                         tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAAt                                                                                                                                                                                                                            >  7:27875/1‑87 (MQ=255)

                                                                                                                                                                                                                         tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAAt                                                                                                                                                                                                                            <  8:27875/87‑1 (MQ=255)

                                                                                                                                                                                                                         tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCAt                                                                                                                                                              >  7:560930/1‑149 (MQ=255)

                                                                                                                                                                                                                           ttGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCg                                                                                                                                                            >  5:1155457/1‑149 (MQ=255)

                                                                                                                                                                                                                           ttGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCg                                                                                                                                                            >  5:25642/1‑149 (MQ=255)

                                                                                                                                                                                                                           ttGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCg                                                                                                                                                            >  7:998938/1‑149 (MQ=255)

                                                                                                                                                                                                                             ggATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAc                                                                                                                                                          <  7:504888/149‑1 (MQ=255)

                                                                                                                                                                                                                             ggATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAc                                                                                                                                                          >  6:960803/1‑149 (MQ=255)

                                                                                                                                                                                                                              gATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAcc                                                                                                                                                         <  5:924065/149‑1 (MQ=255)

                                                                                                                                                                                                                                ttGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCg                                                                                                                                                       >  7:1121011/1‑149 (MQ=255)

                                                                                                                                                                                                                                   tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGc                                                                                                                                                                                             >  2:293467/1‑108 (MQ=255)

                                                                                                                                                                                                                                   tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACg                                                                                                                                                                                     >  2:693488/1‑116 (MQ=255)

                                                                                                                                                                                                                                       aaaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGa                                                                                                                                                                                <  2:237372/117‑1 (MQ=255)

                                                                                                                                                                                                                                       aaaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCttt                                                                                                                                                >  8:98548/1‑149 (MQ=255)

                                                                                                                                                                                                                                             agTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACg                                                                                                                                                                              >  4:117753/1‑113 (MQ=255)

                                                                                                                                                                                                                                             agTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCTTTTAACTTCTCCGGGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACTTTc            >  2:191224/1‑275 (MQ=255)

                                                                                                                                                                                                                                              gTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCa                                                                                                                                                               <  8:1594422/127‑1 (MQ=255)

                                                                                                                                                                                                                                              gTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCa                                                                                                                                                               >  7:1594422/1‑127 (MQ=255)

                                                                                                                                                                                                                                                 aTTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCt                                                                                                                                      >  8:553941/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  tttGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCAc                                                                                                                                                                                      >  6:105657/1‑100 (MQ=255)

                                                                                                                                                                                                                                                  tttGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCAc                                                                                                                                                                                      <  5:105657/100‑1 (MQ=255)

                                                                                                                                                                                                                                                  tttGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCtt                                                                                                                                     >  5:157539/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      ttCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCg                                                                                                                                 >  6:469672/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       tCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGTAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGAAGGTGTAACCTTTGTCATCGACTCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAAACAACATTCGGAGTAAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAAATTGAAGCGGGGGCAATTTATGAATATCCGGAACATTTACGTTCATTACTTAAt  >  4:382205/1‑275 (MQ=255)

                                                                                                                                                                                                                                                         aCGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCt                                                                                                                              >  8:844722/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          cGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                                                                                                                                            >  7:671144/1‑102 (MQ=255)

                                                                                                                                                                                                                                                          cGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                                                                                                                                            <  8:671144/102‑1 (MQ=255)

                                                                                                                                                                                                                                                           gCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                                                                                                                                            >  8:1688154/1‑101 (MQ=255)

                                                                                                                                                                                                                                                           gCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACAGTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                                                                                                                                            <  7:1688154/101‑1 (MQ=255)

                                                                                                                                                                                                                                                             cGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTa                                                                                                                          <  5:195202/149‑1 (MQ=255)

AAGCGCCAACGGCCGATCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAAT  >  NC_000913/1823080‑1823601

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