BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,823,335	G→C	9.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	9.8%	620.1 / 47.7	183	intergenic (+50/‑180)	nadE/cho	NAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand): ref base G (83/82); new base C (8/10); total (91/92)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-02

Only 100 of 185 total aligned reads displayed.

TAAGAAAGCGCCAACGGCCGATCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCAT  >  NC_000913/1823075‑1823593

taAGAAAGCGCCAACGGCCGATCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACg                                                                                                                                                                                                                                                      <  2:203891/275‑1 (MQ=255)

     aaGCGCCAACGGCCGATCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                                                                                                                 >  4:47597/1‑275 (MQ=255)

     aaGCGCCAACGGCCGATCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                                                                                                                 >  4:188063/1‑275 (MQ=255)

                         caagcagatCCCCCTTCTTTCCNNGTTGANGGTGGACCCGGCGTGACCCACGACAATATCCAAGACTATCTGGAAGGGAAAAAAGTCCCTCAACAGGTCGCCAGAACAATAGAGAACTGGTAGCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTAGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGc                                                                                                                                                                                                                             <  3:8176/269‑1 (MQ=255)

                                   gCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCg                                                                                                                                                                                                                                             <  2:222549/249‑1 (MQ=255)

                                   gCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATa                                                                                                                                                                                                                   <  2:15885/275‑1 (MQ=255)

                                                                             aCAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAg                                                                                                                                                                         <  2:188063/275‑1 (MQ=255)

                                                                                    cgacgaTGATATGGAAGGGAAACACGTACCTCAAAAGGTCGCCAGAACAATAGAAAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgtgt                                                                                                                                                                  <  1:41279/275‑1 (MQ=255)

                                                                                             tCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                                                                                                                 >  4:218419/1‑187 (MQ=255)

                                                                                             tCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCttt                                                                                                                                                           >  2:170442/1‑273 (MQ=255)

                                                                                                       aaaaaCGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGNGTGAACGCCTTATTCGACCTATAANACNATGCAAAt                                                                                                                                                                                                                         >  4:379029/1‑201 (MQ=255)

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

                                                                                                                ccTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCgg                                                                                                                                                                                                                                                                    >  8:1223062/1‑149 (MQ=255)

                                                                                                                   cAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATg                                                                                                                                                                                                                                                                 >  7:1320183/1‑149 (MQ=255)

                                                                                                                     aCAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcg                                                                                                                                                                                                                                                               >  6:680037/1‑149 (MQ=255)

                                                                                                                     aCAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcg                                                                                                                                                                                                                                                               <  8:1320183/149‑1 (MQ=255)

                                                                                                                       aGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcggc                                                                                                                                                                                                                                                             <  8:563258/149‑1 (MQ=255)

                                                                                                                          tCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTg                                                                                                                                                                                                                                                          <  5:500636/149‑1 (MQ=255)

                                                                                                                           cGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGa                                                                                                                                                                                                                                                         <  8:775042/149‑1 (MQ=255)

                                                                                                                            gCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGa                                                                                                                                                                                                                                                         <  7:1401841/148‑1 (MQ=255)

                                                                                                                              cAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACg                                                                                                                                                                                                                                                      >  7:927337/1‑149 (MQ=255)

                                                                                                                                gAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGcc                                                                                                                                                                                                                                                    >  7:412835/1‑149 (MQ=255)

                                                                                                                                 aaCAATAGAGAGCTGCTATCTGAAAACCGAACATAAACGCCGTCCGCCATTTACCTTTTTCGAGGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGAAcaaat                                                                                                                                                                                                                                                   <  5:1148062/149‑5 (MQ=255)

                                                                                                                                 aaCAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTCGTTCACGCCGCATCCGGCGTGAAcaaat                                                                                                                                                                                                                                                   <  6:1262882/149‑5 (MQ=255)

                                                                                                                                 aaCAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTGCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTg                                                                                                                     <  2:47597/275‑1 (MQ=255)

                                                                                                                                  aCAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                                                                                                                  <  6:1311361/149‑1 (MQ=255)

                                                                                                                                    aaTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGcc                                                                                                                                                                                                                                                    >  5:379605/1‑145 (MQ=255)

                                                                                                                                    aaTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                                                                                                                                >  7:215095/1‑149 (MQ=255)

                                                                                                                                        gagaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                                                                                                                                                                                            <  6:1053720/149‑1 (MQ=255)

                                                                                                                                        aagaACTTGTATCTACAAACCGAACATAAACGCCTTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGGAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCTTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGCGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTaca                                                                                                                <  1:10584/272‑1 (MQ=255)

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

                                                                                                                                           aaCTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCAAGCCGGATGCGGCGTGAAACCCTTATTCGACCTATAAAACTATGCAAATTc                                                                                                                                                                                                                       >  2:259918/1‑167 (MQ=255)

                                                                                                                                              tGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCtata                                                                                                                                                                                                                                      >  6:117190/1‑149 (MQ=255)

                                                                                                                                                gTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGTGAAAAAGAGTAATTTGTTCACGCCGGATGCTGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCTTAGCACGTCAGAAGGTGTAACCTTTGTAATCGACCCCCCTCTTTTTTAATCGCTTCCGGCCTGTTAAACGGgatagaa                                                                                                      >  4:279450/1‑274 (MQ=255)

                                                                                                                                                gTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGt                                                                                                                    >  4:124291/1‑261 (MQ=255)

                                                                                                                                                gTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCAATAAAACTATGAGAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGgatagat                                                                                                      >  1:61497/1‑275 (MQ=255)

                                                                                                                                                  aTCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAAc                                                                                                                                                                                                                                  >  7:1361214/1‑149 (MQ=255)

                                                                                                                                                    cTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCg                                                                                                                                                                                                                                             <  2:29302/136‑1 (MQ=255)

                                                                                                                                                    cTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTaca                                                                                                                <  2:124291/261‑1 (MQ=255)

                                                                                                                                                     tGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTAt                                                                                                                                                                                                                               >  5:130833/1‑149 (MQ=255)

                                                                                                                                                      gAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATg                                                                                                                                                                                                                              >  5:75289/1‑149 (MQ=255)

                                                                                                                                                       aaaaCCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGACAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGTCCTATAAAACTATGc                                                                                                                                                                                                                             >  6:432650/1‑149 (MQ=255)

                                                                                                                                                         aaCCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCaa                                                                                                                                                                                                                           >  6:980367/1‑149 (MQ=255)

                                                                                                                                                          aCCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCaaa                                                                                                                                                                                                                          <  7:1093310/149‑1 (MQ=255)

                                                                                                                                                              aaCATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGa                                                                                        >  2:143503/1‑275 (MQ=255)

                                                                                                                                                              aaCATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGa                                                                                                                                                                                                                                                         <  8:521812/114‑1 (MQ=255)

                                                                                                                                                               aCATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCaa                                                                                                                                                                                                                     <  6:75289/149‑1 (MQ=255)

                                                                                                                                                                cATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAAt                                                                                                                                                                                                                    <  8:636297/149‑1 (MQ=255)

                                                                                                                                                                 aTAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATa                                                                                                                                                                                                                   <  5:1050181/149‑1 (MQ=255)

                                                                                                                                                                 aTAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGaa                                                                                                                                                                                                                                                        <  8:721030/112‑1 (MQ=255)

                                                                                                                                                                    aaCGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACAt                                                                                                                                                                                                                >  6:355214/1‑149 (MQ=255)

                                                                                                                                                                     aCGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACAtt                                                                                                                                                                                                               <  7:328662/149‑1 (MQ=255)

                                                                                                                                                                           tCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTAAGCAAATTCAATACATTGCAGGa                                                                                                                                                                                                         <  5:897928/149‑1 (MQ=255)

                                                                                                                                                                              gCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTc                                                                                                                                                                                                      <  6:417677/149‑1 (MQ=255)

                                                                                                                                                                                   ttACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATa                                                                                                                                                                                                                   >  6:1278412/1‑131 (MQ=255)

                                                                                                                                                                                     aCCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAgg                                                                                                                                                                                               >  5:1112527/1‑149 (MQ=255)

                                                                                                                                                                                       cGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCaaa                                                                                                                                                                                                                          >  6:57206/1‑120 (MQ=255)

                                                                                                                                                                                       cGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCaaa                                                                                                                                                                                                                          <  5:57206/120‑1 (MQ=255)

                                                                                                                                                                                        gTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAAt                                                                                                                                                                                                  >  8:763817/1‑143 (MQ=255)

                                                                                                                                                                                        gTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCt                                                                                                                                                                                            <  7:1223062/149‑1 (MQ=255)

                                                                                                                                                                                        gTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCt                                                                                                                                                                                            <  5:528511/149‑1 (MQ=255)

                                                                                                                                                                                          tttCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGa                                                                                                                                                                                          <  8:881922/149‑1 (MQ=255)

                                                                                                                                                                                            tCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACa                                                                                                                                                                                        >  7:507011/1‑149 (MQ=255)

                                                                                                                                                                                              gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAg                                                                                                                                                                                                           <  8:666520/128‑1 (MQ=255)

                                                                                                                                                                                              gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAg                                                                                                                                                                                                           >  7:666520/1‑128 (MQ=255)

                                                                                                                                                                                              gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTc                                                                                                                                                                                                      >  8:990670/1‑133 (MQ=255)

                                                                                                                                                                                              gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTc                                                                                                                                                                                                      <  7:990670/133‑1 (MQ=255)

                                                                                                                                                                                              gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAgg                                                                                                                                                                                      >  7:1019135/1‑149 (MQ=255)

                                                                                                                                                                                              gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGAAc                                                                                                                                                                                                                                                       <  6:525058/84‑1 (MQ=255)

                                                                                                                                                                                              gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGAAc                                                                                                                                                                                                                                                       >  5:525058/1‑84 (MQ=255)

                                                                                                                                                                                                         gAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTc                                                                                                                                                                           >  5:257492/1‑149 (MQ=255)

                                                                                                                                                                                                           aaaaGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAg                                                                                                                                                                         >  7:619683/1‑149 (MQ=255)

                                                                                                                                                                                                                 aataatTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgtg                                                                                                                                                                   <  7:770208/149‑1 (MQ=255)

                                                                                                                                                                                                                 aataatTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgtg                                                                                                                                                                   >  7:823746/1‑149 (MQ=255)

                                                                                                                                                                                                                    aatTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACgg                                                                                                                                                                     >  8:337961/1‑144 (MQ=255)

                                                                                                                                                                                                                    aatTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTaa                                                                                                                                                                <  7:542017/149‑1 (MQ=255)

                                                                                                                                                                                                                     atTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAAc                                                                                                                                                               >  8:270618/1‑149 (MQ=255)

                                                                                                                                                                                                                        tGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGt                                                                                                                                                                            >  4:228946/1‑133 (MQ=255)

                                                                                                                                                                                                                        tGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCtt                                                                                                                                                            <  8:1368687/149‑1 (MQ=255)

                                                                                                                                                                                                                         gCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCttt                                                                                                                                                           >  8:1101523/1‑149 (MQ=255)

                                                                                                                                                                                                                          cacaTATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTg                                                                                                                                                          <  6:1081314/149‑1 (MQ=255)

                                                                                                                                                                                                                           acaTATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGCGTAACCTTTGt                                                                                                                                                         <  7:1166409/149‑1 (MQ=255)

                                                                                                                                                                                                                              tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCAt                                                                                                                                                      <  8:215095/149‑1 (MQ=255)

                                                                                                                                                                                                                                ttGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCg                                                                                                                                                    >  6:84104/1‑149 (MQ=255)

                                                                                                                                                                                                                                ttGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCg                                                                                                                                                    >  8:1014313/1‑149 (MQ=255)

                                                                                                                                                                                                                                    aTTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAccc                                                                                                                                                >  7:378775/1‑149 (MQ=255)

                                                                                                                                                                                                                                    aTTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAccc                                                                                                                                                <  5:1407807/149‑1 (MQ=255)

                                                                                                                                                                                                                                    aTTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAccc                                                                                                                                                <  6:257492/149‑1 (MQ=255)

                                                                                                                                                                                                                                            aaaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTg                                                                                                                     <  4:152801/168‑1 (MQ=255)

                                                                                                                                                                                                                                             aaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCg                                                                                                                                                                                                     >  5:796369/1‑87 (MQ=255)

                                                                                                                                                                                                                                                 gagTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTaa                                                                                                                                   <  8:192565/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    tAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCg                                                                                                                                >  6:18348/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      aTTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCtt                                                                                                                                         >  8:834962/1‑138 (MQ=255)

                                                                                                                                                                                                                                                       tttGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTg                                                                                                                                                          >  7:29975/1‑120 (MQ=255)

                                                                                                                                                                                                                                                          gTTCATGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTccc                                                                                                                          >  5:1459978/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           ttCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCttttt                                                                                                                                      >  7:854784/1‑136 (MQ=255)

                                                                                                                                                                                                                                                           ttCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAAt                                                                                                                                  <  7:652702/140‑1 (MQ=255)

                                                                                                                                                                                                                                                              aCGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCt                                                                                                                      >  5:1273350/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                gCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAg                                                                                                                                                                                       >  6:313106/1‑82 (MQ=255)

                                                                                                                                                                                                                                                                gCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAg                                                                                                                                                                                       <  5:313106/82‑1 (MQ=255)

TAAGAAAGCGCCAACGGCCGATCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCAT  >  NC_000913/1823075‑1823593

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