BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,169,380	T→C	5.3%	intergenic (+50/+32)	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,380	0	T	C	5.3%	576.3 / 14.2	169	intergenic (+50/+32)	bhsA/ldtC	biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand): ref base T (92/67); new base C (7/2); total (99/70)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 172 total aligned reads displayed.

CATGTCATTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGT  >  NC_000913/1169117‑1169650

                                                    caacaaAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTAACCGGGCTTTAAAGCGTTTGGGGGCTCCATTGATGAcgg                                                                                                                                                                                                                  >  4:203409/1‑274 (MQ=255)

                                                                                             tCTGGGATGGCTGGAAAAGCAGCTGGCGCAAAAAGCGGATGGGATGGGCGCAAAAGCTTTCCGTATTACTTCTGTAACCGGTCCTAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGt                                                                                                                                                                        <  2:203409/275‑1 (MQ=255)

                                                                                                       cTGGAAGAGCAGATTGCGCAAAAAGCGGATGAGATGGGCGCAAAAACTTTCCGTAGTACTTCTGTAACCGGTACGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATGAATGCCTGCTACTTCTGATTTTTTCCCCGCGACATGCCGTGTCGAGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACtt                                                                                                                                                              <  4:259704/275‑1 (MQ=255)

                                                                                                                cagcTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCtt                                                                                                                                                                                                                                                >  4:327230/1‑184 (MQ=255)

                                                                                                                   cTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACAt                                                                                                                                                                                                                                                                                <  5:1159271/149‑1 (MQ=255)

                                                                                                                    tGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACa                                                                                                                                                                                                                                            <  2:327230/184‑1 (MQ=255)

                                                                                                                     ggCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGc                                                                                                                                                                                                                                                                              >  5:732607/1‑149 (MQ=255)

                                                                                                                     ggCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGc                                                                                                                                                                                                                                                                              <  7:743181/149‑1 (MQ=255)

                                                                                                                      gcgcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGcc                                                                                                                                                                                                                                                                             <  6:348489/149‑1 (MQ=255)

                                                                                                                        gcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgt                                                                                                                                                                                                                                                                           >  6:628866/1‑149 (MQ=255)

                                                                                                                           aaaaGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTc                                                                                                                                                                                                                                                                        >  8:162970/1‑149 (MQ=255)

                                                                                                                                  gATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCTTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGa                                                                                                                                                                                                                                                                 <  5:1291367/149‑1 (MQ=255)

                                                                                                                                  gATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGa                                                                                                                                                                                                                                                                 <  5:628866/149‑1 (MQ=255)

                                                                                                                                  gATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa                                                                                                                                    <  2:244655/274‑1 (MQ=255)

                                                                                                                                      agaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttt                                                                                                                                                                                                                                                             <  8:488800/149‑1 (MQ=255)

                                                                                                                                      agaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGa                                                                                                                               <  2:156520/275‑1 (MQ=255)

                                                                                                                                       gaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttt                                                                                                                                                                                                                                                            >  5:353700/1‑149 (MQ=255)

                                                                                                                                       gaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCATCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCGGGGaaaaaa                                                                                                                                                                                                                                                            <  6:701661/149‑6 (MQ=255)

                                                                                                                                        aTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                                                                                                                                                                                           >  7:1431251/1‑149 (MQ=255)

                                                                                                                                         tGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttttt                                                                                                                                                                                                                                                          <  5:1139428/149‑1 (MQ=255)

                                                                                                                                                   aTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttt                                                                                                                                                                                                                                                            <  7:223741/137‑1 (MQ=255)

                                                                                                                                                   aTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttt                                                                                                                                                                                                                                                            >  8:223741/1‑137 (MQ=255)

                                                                                                                                                           gTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATACCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGt                                                                                                                                                                                                                                        >  7:387690/1‑149 (MQ=255)

                                                                                                                                                             aTTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGttt                                                                                                                                                                                                                                      >  7:371861/1‑149 (MQ=255)

                                                                                                                                                                  ttCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                                                                                                                                                                                                 >  7:1057100/1‑149 (MQ=255)

                                                                                                                                                                  ttCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                                                                                                                                                                                                 <  5:421647/149‑1 (MQ=255)

                                                                                                                                                                  ttCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                                                                                                                                                                                                 <  8:1350759/149‑1 (MQ=255)

                                                                                                                                                                     tGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtgtg                                                                                                                                                                                                                                   >  5:592607/1‑144 (MQ=255)

                                                                                                                                                                     tGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtgtg                                                                                                                                                                                                                                   <  6:592607/144‑1 (MQ=255)

                                                                                                                                                                     tGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTc                                                                                                                                                                                                                              >  8:767134/1‑149 (MQ=255)

                                                                                                                                                                     tGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTc                                                                                                                                                                                                                              >  7:1287759/1‑149 (MQ=255)

                                                                                                                                                                      gTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCa                                                                                                                                                                                                                             >  6:840628/1‑149 (MQ=255)

                                                                                                                                                                        aaCCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACt                                                                                                                                                                                                                           <  7:162970/149‑1 (MQ=255)

                                                                                                                                                                         aCCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGCTTGTGGGCTCACTTGATGACGGCGGa                                                                                                                                                                                                              >  4:135330/1‑161 (MQ=255)

                                                                                                                                                                            ggTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                                                                                                                                                       <  8:1187979/149‑1 (MQ=255)

                                                                                                                                                                            ggTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                                                                                                                                                       >  5:463122/1‑149 (MQ=255)

                                                                                                                                                                                 gAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttt                                                                                                                                                                                                                                                            >  1:22435/1‑107 (MQ=255)

                                                                                                                                                                                  aaTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGgcag                                                                                                 <  4:84348/261‑1 (MQ=255)

                                                                                                                                                                                    tACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcggcgg                                                                                                                                                                                                               >  7:601188/1‑149 (MQ=255)

                                                                                                                                                                                     aCCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTAt                                                                                                                                                                                                                                                        <  3:22435/107‑1 (MQ=255)

                                                                                                                                                                                     aCCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                                                                                                                                                       <  8:645797/140‑1 (MQ=255)

                                                                                                                                                                                     aCCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                                                                                                                                                       >  7:645797/1‑140 (MQ=255)

                                                                                                                                                                                        cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACAt                                                                                                                                                                                                           >  8:29484/1‑149 (MQ=255)

                                                                                                                                                                                         tCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATc                                                                                                                                                                                                          <  8:252473/149‑1 (MQ=255)

                                                                                                                                                                                           cATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCa                                                                                                                                                                                                        >  7:326159/1‑149 (MQ=255)

                                                                                                                                                                                               gAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCgg                                                                                                                                                                                                    >  7:892473/1‑149 (MQ=255)

                                                                                                                                                                                                aaCAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGc                                                                                                                                                                                                   >  8:470352/1‑149 (MQ=255)

                                                                                                                                                                                                aaCAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATc                                                                                                                                                                  >  4:372062/1‑182 (MQ=255)

                                                                                                                                                                                                    gcagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACa                                                                 <  4:50712/275‑1 (MQ=255)

                                                                                                                                                                                                     cagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATccc                                                                                                                                                                                              <  7:29484/149‑1 (MQ=255)

                                                                                                                                                                                                     cagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATccc                                                                                                                                                                                              <  7:470352/149‑1 (MQ=255)

                                                                                                                                                                                                     cagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATAATTGCATTGCGCTATTCAGCGTACTTGGCAGCAGCTGGGGATCGTCCTCTGtcttttttccaca                                                                >  1:10751/1‑266 (MQ=255)

                                                                                                                                                                                                         aaTTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCAGGaa                                                                                                                                                                                          <  7:784555/149‑1 (MQ=255)

                                                                                                                                                                                                             tataAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa                                                                                                                                                                                      <  5:1043604/149‑1 (MQ=255)

                                                                                                                                                                                                             tataAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa                                                                                                                                                                                      >  5:562707/1‑149 (MQ=255)

                                                                                                                                                                                                                    aaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccat                                                                                                                                                                               <  8:1250060/149‑1 (MQ=255)

                                                                                                                                                                                                                      gCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATca                                                                                                                                                                             <  8:326159/149‑1 (MQ=255)

                                                                                                                                                                                                                          tAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACACCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattgca                                                                                                                          <  4:301006/196‑1 (MQ=255)

                                                                                                                                                                                                                               cTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGGACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGAcgcg                                      <  2:62304/275‑1 (MQ=255)

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

                                                                                                                                                                                                                                 cATTAATGCCTGTTACTGCTGATTTTTTCCCCGCGACAGGCCGGGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCa                                    <  2:57978/275‑1 (MQ=255)

                                                                                                                                                                                                                                  aTTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTTCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCAc                                   >  2:106578/1‑275 (MQ=255)

                                                                                                                                                                                                                                   ttAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCAc                                                                                                                                                                <  5:230746/149‑1 (MQ=255)

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

                                                                                                                                                                                                                                       tGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGAAAGCGGCTGATGTACTTCAACCAGACGACCACCGt                                >  2:341804/1‑273 (MQ=255)

                                                                                                                                                                                                                                           tGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGtt                          <  2:40448/275‑1 (MQ=255)

                                                                                                                                                                                                                                              tACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAg                                                                                                                                                     >  5:1339931/1‑149 (MQ=255)

                                                                                                                                                                                                                                                ctgctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGtt                                                                                                                                                   <  7:1175889/149‑1 (MQ=255)

                                                                                                                                                                                                                                                 tgctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGc                                                                                                                                                          >  7:400493/1‑141 (MQ=255)

                                                                                                                                                                                                                                                 tgctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGc                                                                                                                                                          <  8:400493/141‑1 (MQ=255)

                                                                                                                                                                                                                                                 tgctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGttt                                                                                                                                                  <  8:1057100/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  gctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtg                                                                                                                                                 >  8:1051476/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   ctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                                                                                                                                                >  8:78750/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   ctgATTTGGTCCCCGCAAAATGCCGAGACGCGGGGGATATATTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                                                                                                                                                >  6:698839/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    tgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgtg                                                                                                                                               >  7:84357/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    tgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgtg                                                                                                                                               >  7:986391/1‑149 (MQ=255)

                                                                                                                                                                                                                                                         ttttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCa                                                                                                                                          >  8:93503/1‑149 (MQ=255)

                                                                                                                                                                                                                                                         ttttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCa                                                                                                                                          >  6:1466824/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          tttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCAt                                                                                                                                         >  6:1087374/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCttt                                                                                                                                     >  5:1462554/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCttt                                                                                                                                     >  7:588783/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              cccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCttt                                                                                                                                     >  8:372779/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               ccGCGTCATGCCGTGTCGCGGGGA‑TTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaa                                                                                                                                   >  8:799051/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               ccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCa                                                                                                                                                            >  6:1237399/1‑125 (MQ=255)

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

                                                                                                                                                                                                                                                               ccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa                                                                                                                                    >  5:1457363/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                                               ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa                                                                                                                                    >  5:866875/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa                                                                                                                                    >  5:913373/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa                                                                                                                                    >  6:774520/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                cgcgACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGt                                                                                                                                                    <  6:169416/132‑1 (MQ=255)

                                                                                                                                                                                                                                                                cgcgACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGt                                                                                                                                                    >  5:169416/1‑132 (MQ=255)

                                                                                                                                                                                                                                                                 gcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaaa                                                                                                                                  <  6:1462554/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 gcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaaa                                                                                                                                  <  6:371806/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 gcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaaa                                                                                                                                  >  8:834267/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   gACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGATGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAgtgt  <  2:3152/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                    aCATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGa                                                                                                                               <  6:58735/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      aTGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGatt                                                                                                                             <  8:371861/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                       tGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatga                                                                                                                                                                                                                     >  8:1393000/1‑60 (MQ=255)

                                                                                                                                                                                                                                                                       tGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatga                                                                                                                                                                                                                     <  7:1393000/60‑1 (MQ=255)

CATGTCATTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGT  >  NC_000913/1169117‑1169650

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