BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,169,382	C→G	6.3%	intergenic (+52/+30)	bhsA → / ← ldtC	biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,169,382	0	C	G	6.3%	791.0 / 33.2	227	intergenic (+52/+30)	bhsA/ldtC	biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand): ref base C (113/99); new base G (8/7); total (121/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

Only 100 of 231 total aligned reads displayed.

TTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGAT  >  NC_000913/1169124‑1169654

tttGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGCttt                                                                                                                                                                                                                                                                  >  2:546160/1‑275 (MQ=255)

        cTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATcc                                                                                                                                                                                                                                                          <  2:558565/275‑1 (MQ=255)

                  gTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCCTGTCCAGGGGATTTTTTTATTCGTGCTTTACATCg                                                                                                                                                                                                                                             >  1:32860/1‑278 (MQ=255)

                   tCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAg                                                                                                                                                                                                                                               <  2:664505/275‑1 (MQ=255)

                                             caacaaAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATcc                                                                                                                                                                                                                                                          <  1:175102/238‑1 (MQ=255)

                                                                                     aTCTGGGATCGCTGGAAGAGCCGCTGGCGCAAACAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGTTTTTTTTATCCGGGCTTTCCAGCGTTTGTGGGCTCCCTTTCGGACGGCGGACACCCCACGGCATCCCGTACCGGCCATCCCTCACAt                                                                                                                                                                              >  4:261390/1‑274 (MQ=255)

                                                                                                         cagcTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatga                                                                                                                                                                                                                         >  3:258347/1‑211 (MQ=255)

                                                                                                                gggaaaaaGCGGAGTAGTGGGGGCCAAAACCTTACCGTATTACTTCTTTAACCGGTCCGAATACCCTCCAGGGAACACAAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCAGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgtg                                                                                                                                                   <  4:546156/271‑1 (MQ=255)

                                                                                                                 gcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgt                                                                                                                                                                                                                                                                               >  7:735123/1‑149 (MQ=255)

                                                                                                                 gcAAAAAGCGGATGAGATCGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgt                                                                                                                                                                                                                                                                               >  6:1588676/1‑149 (MQ=255)

                                                                                                                   aaaaaGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgtgt                                                                                                                                                                                                                                                                             >  7:491379/1‑149 (MQ=255)

                                                                                                                   aaaaaGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgtgt                                                                                                                                                                                                                                                                             <  8:1577584/149‑1 (MQ=255)

                                                                                                                     aaaGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCAGTGTcg                                                                                                                                                                                                                                                                           >  7:1122600/1‑149 (MQ=255)

                                                                                                                       aGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa                                                                                                                                        <  1:103164/278‑1 (MQ=255)

                                                                                                                        gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                                                                                                                                                        <  6:479575/149‑1 (MQ=255)

                                                                                                                           gATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGa                                                                                                                                                                                                                                                                     >  6:610721/1‑149 (MQ=255)

                                                                                                                               agaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttt                                                                                                                                                                                                                                                                 <  5:290903/149‑1 (MQ=255)

                                                                                                                                gaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttt                                                                                                                                                                                                                                                                <  8:735123/149‑1 (MQ=255)

                                                                                                                                 aTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                                                                                                                                                                                               <  5:992726/149‑1 (MQ=255)

                                                                                                                                 aTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                                                                                                                                                                                               >  5:891832/1‑149 (MQ=255)

                                                                                                                                  tGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATCAATAAGCATTAACACTCATTAATGACTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                                                                                                                                                        <  6:1129123/139‑1 (MQ=255)

                                                                                                                                  tGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttttt                                                                                                                                                                                                                                                              <  6:891832/149‑1 (MQ=255)

                                                                                                                                  tGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttttt                                                                                                                                                                                                                                                              >  7:328988/1‑149 (MQ=255)

                                                                                                                                  tGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                                                                                                                                                        >  5:1129123/1‑139 (MQ=255)

                                                                                                                                      cgcAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATcc                                                                                                                                                                                                                                                          <  6:396795/149‑1 (MQ=255)

                                                                                                                                      cgcAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGACAGGACATCCATCACATTTTGCATCACTTCAGCGTCAATTTGTGCTGCATCTTTAAATGatttcatt                                                                                                                            >  4:506771/1‑275 (MQ=255)

                                                                                                                                         aaaaTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCggg                                                                                                                                                                                                                                                       <  6:284665/149‑1 (MQ=255)

                                                                                                                                                ttCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACa                                                                                                                                                                                                                                                >  7:942435/1‑149 (MQ=255)

                                                                                                                                                     tATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGtt                                                                                                                                                                                                                                           <  5:1491877/149‑1 (MQ=255)

                                                                                                                                                     tATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGtt                                                                                                                                                                                                                                           >  6:583639/1‑149 (MQ=255)

                                                                                                                                                         aCTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTAt                                                                                                                                                                                                                                                            <  4:693544/128‑1 (MQ=255)

                                                                                                                                                           ttCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttt                                                                                                                                                                                                                                                                >  8:1147010/1‑122 (MQ=255)

                                                                                                                                                           ttCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttt                                                                                                                                                                                                                                                                 <  8:302560/121‑1 (MQ=255)

                                                                                                                                                           ttCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                                                                                                                                                                                                     >  8:1370929/1‑149 (MQ=255)

                                                                                                                                                                   aaacagcaaaTACCCCCCATGGAACAGCAGTATNNNTTTAATAAGCCTTATCCCTCATTATTGCCTGCTACGTCTGTTTTTTTCCCCGCGACATGCCTTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGa                                                                                           <  3:296346/271‑1 (MQ=255)

                                                                                                                                                                    cGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                                                                                                                                                            >  8:621925/1‑149 (MQ=255)

                                                                                                                                                                     ggTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                                                                                                                                                           <  5:45802/149‑1 (MQ=255)

                                                                                                                                                                     ggTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                                                                                                                                                           <  7:1547803/149‑1 (MQ=255)

                                                                                                                                                                               cccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatga                                                                                                                                                                                                                         >  6:789460/1‑141 (MQ=255)

                                                                                                                                                                               cccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatga                                                                                                                                                                                                                         <  5:789460/141‑1 (MQ=255)

                                                                                                                                                                                ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACa                                                                                                                                                                                                                >  7:1666020/1‑149 (MQ=255)

                                                                                                                                                                                ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATcaca                                                                                                                                                                               <  1:238482/182‑1 (MQ=255)

                                                                                                                                                                                 cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTc                                                                                                                                                                                                                                                                            >  2:426863/1‑88 (MQ=255)

                                                                                                                                                                                 cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACAt                                                                                                                                                                                                               >  7:103220/1‑149 (MQ=255)

                                                                                                                                                                                 cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGat                                                                                                                                  >  2:3356/1‑226 (MQ=255)

                                                                                                                                                                                    cATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCa                                                                                                                                                                                                            >  6:1612907/1‑149 (MQ=255)

                                                                                                                                                                                     aTGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                                                                                                                                                        <  4:426859/88‑1 (MQ=255)

                                                                                                                                                                                     aTGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCAc                                                                                                                                                                                                           <  8:491379/149‑1 (MQ=255)

                                                                                                                                                                                     aTGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACgg                                                                                                                                                                                           <  2:263859/165‑1 (MQ=255)

                                                                                                                                                                                      tGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtgt                                                                                                                                                                                                                                        <  8:1470655/119‑1 (MQ=255)

                                                                                                                                                                                         aaCAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGtt                                                                                                                                                                                                                                           >  7:1165015/1‑113 (MQ=255)

                                                                                                                                                                                          aCAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGc                                                                                                                                                                                                                                                      <  8:694243/101‑1 (MQ=255)

                                                                                                                                                                                            agcagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATc                                                                                                                                                                                                    <  8:103220/149‑1 (MQ=255)

                                                                                                                                                                                            agcagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCACGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTCTTCAGCGTAATTGACCGACGCTGCGGATCGTCATCAATCTTCTCTg                                                                        >  2:492817/1‑273 (MQ=255)

                                                                                                                                                                                               agTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCg                                                                                                                                                                                                 >  5:1260105/1‑149 (MQ=255)

                                                                                                                                                                                                           ataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcc                                                                                                                                                                                     >  5:581796/1‑149 (MQ=255)

                                                                                                                                                                                                            taaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTAc                                                                                                                                                                                                                                                 >  8:117034/1‑88 (MQ=255)

                                                                                                                                                                                                             aaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccat                                                                                                                                                                                   <  7:1370929/149‑1 (MQ=255)

                                                                                                                                                                                                               gCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATca                                                                                                                                                                                 <  7:621925/149‑1 (MQ=255)

                                                                                                                                                                                                                       cccacttAATGCCGGCGACTGCTGTTTTTTTCCCCGCGACAGGGCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCCTCCCGGAACTGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGAcgc                                           <  4:20016/270‑1 (MQ=255)

                                                                                                                                                                                                                       ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACa                                                                                                                                                                                                                                                >  6:1490156/1‑78 (MQ=255)

                                                                                                                                                                                                                        cTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACg                                                                                                                                                                                            <  6:301871/129‑1 (MQ=255)

                                                                                                                                                                                                                        cTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACg                                                                                                                                                                                            >  5:301871/1‑129 (MQ=255)

                                                                                                                                                                                                                        cTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCa                                                                                                                                                                        >  8:642582/1‑149 (MQ=255)

                                                                                                                                                                                                                         tCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCAt                                                                                                                                                                       >  7:816818/1‑149 (MQ=255)

                                                                                                                                                                                                                          cATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATc                                                                                                                                                                      <  6:730095/149‑1 (MQ=255)

                                                                                                                                                                                                                          cATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCCGCTGCGGTTCGTCATCAATTTTCTCTTGCAGCGGCTGATGTACTTCAACCAGCCacgc                                         >  4:655110/1‑274 (MQ=255)

                                                                                                                                                                                                                          cATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTCGGCTCACTTGATGACGGCGGACATCCACCGGCACCCCGGAACGGACATCCATTACATGTTGCATCCCTTCAGCGTCCGTTTGTGCTGCATCTTTAAAGTATTGCATTGCGCTTTTCAGCGTAATTGGTCGAAGCTGCGGATCTGCATCAATTTTCCCTGACAGCCGCCGAGTTACTTCAACCAGAccgtca                                        >  4:415531/1‑270 (MQ=255)

                                                                                                                                                                                                                           aTTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa                                                                                                                                                                                          >  8:666432/1‑128 (MQ=255)

                                                                                                                                                                                                                           aTTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa                                                                                                                                                                                          <  7:666432/128‑1 (MQ=255)

                                                                                                                                                                                                                           aTTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGAGGTACGTCAACCAGCAgcgcca                                       >  2:145628/1‑273 (MQ=255)

                                                                                                                                                                                                                               aTGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTc                                                                                                                                                                 >  5:375279/1‑149 (MQ=255)

                                                                                                                                                                                                                                tGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa                                                                                                                                                                                          >  3:180503/1‑123 (MQ=255)

                                                                                                                                                                                                                                   cTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTc                                                                                                                                                                                                                                  <  5:97332/80‑1 (MQ=255)

                                                                                                                                                                                                                                   cTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTc                                                                                                                                                                                                                                  >  6:97332/1‑80 (MQ=255)

                                                                                                                                                                                                                                   cTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCg                                                                                                                                                             >  5:247550/1‑149 (MQ=255)

                                                                                                                                                                                                                                   cTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCg                                                                                                                                                             >  7:714710/1‑149 (MQ=255)

                                                                                                                                                                                                                                    tGCTACTGCTGATTTTTTCCCCGCGACATGCTGTGTCGCGGGGGCTTTTTTTTCCGTGGTTTACATCGTTTGTGGGCTAACTTGATGACGGCGGACATCCACGCGCATCCCGGAACGTACTCCCATCACATGTGGCATAACTACTGCGt                                                                                                                                                            >  8:772281/1‑149 (MQ=255)

                                                                                                                                                                                                                                    tGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGtt                                                                                                                                                                           <  8:1120641/134‑1 (MQ=255)

                                                                                                                                                                                                                                     gCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtg                                                                                                                                                                                                                                         >  4:273317/1‑71 (MQ=255)

                                                                                                                                                                                                                                     gCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAAc                                                  <  1:156319/254‑1 (MQ=255)

                                                                                                                                                                                                                                     gCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATACATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCCGCGTACTTGGCAGAAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCGGATGTACTTCAACCAGACGCCCCCCTTTCGGTTc                             >  4:575828/1‑275 (MQ=255)

                                                                                                                                                                                                                                      cTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCa                                                                                                                                                          <  6:1265362/149‑1 (MQ=255)

                                                                                                                                                                                                                                         ctgctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTa                                                                                                                                                                                                                                                  <  2:596758/58‑1 (MQ=255)

                                                                                                                                                                                                                                            ctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                                                                                                                                                    >  6:1650524/1‑149 (MQ=255)

                                                                                                                                                                                                                                               tttttttCCCCGCGACATGCCGTGTCGCGGGTATTTTTTTATCCAGGGTTTACAGCGTTTGTGGGCTCACTTGGTGAAGGCGGACAACCAACGGCATCCGGGAACGGACTTCCAACAGATGTTGAATCACTTCGGCGTCAGTTTGtgct                                                                                                                                                 <  6:1387664/148‑1 (MQ=255)

                                                                                                                                                                                                                                               aTTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgct                                                                                                                                                 <  8:169857/149‑1 (MQ=255)

                                                                                                                                                                                                                                                ttttttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                                                                                                                                                            <  5:693002/73‑1 (MQ=255)

                                                                                                                                                                                                                                                  aaatCCCCGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                                                                                                                                                            <  5:1668820/68‑1 (MQ=25)

                                                                                                                                                                                                                                                      cgccGCGACAGGCCGGGGCGCGGGTATTTTTTTAGTCGGGCTTTACAGCGTGTGTGGGCTCACTTGATGACTGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACAGGTTGCATCACTTCAGCGGCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTAGTCAGCTTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTGCAACCAGACGCGCACCGTTTGGTTCGTCAGAGACTTTTATCg            <  2:184523/273‑1 (MQ=255)

                                                                                                                                                                                                                                                      ccccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCtt                                                                                                                                          >  5:1386803/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        ccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa                                                                                                                                        >  6:968255/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa                                                                                                                                        >  7:21102/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa                                                                                                                                        >  8:914397/1‑149 (MQ=255)

                                                                                                                                                                                                                                                         cgcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaa                                                                                                                                       >  5:622049/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          gcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaaa                                                                                                                                      <  5:1588676/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              cATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTtg    <  2:286206/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                tGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattg                                                                                                                                >  8:1021269/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                tGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattg                                                                                                                                <  6:197473/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                tGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTtgat  <  2:521244/275‑1 (MQ=255)

TTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGAT  >  NC_000913/1169124‑1169654

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