BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,169,384	G→A	6.4%	intergenic (+54/+28)	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,384	0	G	A	6.4%	723.0 / 25.1	227	intergenic (+54/+28)	bhsA/ldtC	biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand): ref base G (110/99); new base A (8/6); total (120/107)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

Only 100 of 233 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)

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

                                                                          cgcgGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGCTGACGGCGGACATCCCCCGGAATCCCGGAACGGAc                                                                                                                                                                                         >  4:521240/1‑274 (MQ=255)

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

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

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

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

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

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

                                                                                                                        gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                                                                                                                                                        <  7:415333/149‑1 (MQ=255)

                                                                                                                         cGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCggg                                                                                                                                                                                                                                                                       <  6:986858/149‑1 (MQ=255)

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

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

                                                                                                                                 aTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                                                                                                                                                                                               <  6:1523393/149‑1 (MQ=255)

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

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

                                                                                                                                        cAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCgg                                                                                                                                                                                                                                                        <  8:796770/149‑1 (MQ=255)

                                                                                                                                               tttCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTAc                                                                                                                                                                                                                                                 >  6:783754/1‑149 (MQ=255)

                                                                                                                                                     tATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGtt                                                                                                                                                                                                                                           >  8:1442677/1‑149 (MQ=255)

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

                                                                                                                                                      aTTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGttt                                                                                                                                                                                                                                          <  5:72844/149‑1 (MQ=255)

                                                                                                                                                          cTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTgg                                                                                                                                                                                                                                      <  8:634903/149‑1 (MQ=255)

                                                                                                                                                           ttCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttt                                                                                                                                                                                                                                                                <  7:1147010/122‑1 (MQ=255)

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

                                                                                                                                                                       tCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACa                                                                                                                                                                                                                                                <  7:735025/126‑1 (MQ=255)

                                                                                                                                                                            aTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGGGCTGCATCTTTAAATGATTTCATTGCGCTCTTCAGCGTAATTTGCAGCAGCTGCGGATCgc                                                                                       >  4:286203/1‑273 (MQ=255)

                                                                                                                                                                             tACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcggcgg                                                                                                                                                                                                                   <  6:1610839/149‑1 (MQ=255)

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

                                                                                                                                                                                ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACa                                                                                                                                                                                                                >  5:1387664/1‑149 (MQ=255)

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

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

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

                                                                                                                                                                                 cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                                                                                                                                                        >  6:300366/1‑92 (MQ=39)

                                                                                                                                                                                 cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                                                                                                                                                        <  5:300366/92‑1 (MQ=255)

                                                                                                                                                                                  tCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATc                                                                                                                                                                                                              >  5:868869/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                                                                                                                                                                                                                                        >  7:1470655/1‑119 (MQ=255)

                                                                                                                                                                                        gAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTgg                                                                                                                                                                                                                                      >  5:1329519/1‑119 (MQ=255)

                                                                                                                                                                                        gAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTgg                                                                                                                                                                                                                                      <  6:1329519/119‑1 (MQ=255)

                                                                                                                                                                                        gAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                                                                                                                                                           >  7:1609201/1‑130 (MQ=255)

                                                                                                                                                                                         aaCAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTGTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTAATCCGGGCTTTACAGCGtt                                                                                                                                                                                                                                           <  8:1165015/113‑1 (MQ=255)

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

                                                                                                                                                                                          aCAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGc                                                                                                                                                                                                                                                      >  7:694243/1‑101 (MQ=255)

                                                                                                                                                                                            agcagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTgg                                                                                                                                                                                                                                      <  6:202110/115‑1 (MQ=255)

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

                                                                                                                                                                                              cagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGc                                                                                                                                                                                                       >  4:689649/1‑144 (MQ=255)

                                                                                                                                                                                                gTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttttt                                                                                                                                                                                                                                                              >  6:969308/1‑87 (MQ=255)

                                                                                                                                                                                                 tAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                                                                                                                                                        >  6:1640956/1‑76 (MQ=38)

                                                                                                                                                                                                 tAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                                                                                                                                                        <  5:1640956/76‑1 (MQ=38)

                                                                                                                                                                                                  aaTTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATcc                                                                                                                                                                                                   <  2:689654/144‑1 (MQ=255)

                                                                                                                                                                                                     ttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATccc                                                                                                                                                                                                  <  6:360298/142‑1 (MQ=255)

                                                                                                                                                                                                          aataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa                                                                                                                                                                                          >  8:742805/1‑145 (MQ=255)

                                                                                                                                                                                                            taaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTAc                                                                                                                                                                                                                                                 <  7:117034/88‑1 (MQ=255)

                                                                                                                                                                                                              aGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccatc                                                                                                                                                                                  <  6:868869/149‑1 (MQ=255)

                                                                                                                                                                                                              aGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTAGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccatc                                                                                                                                                                                  >  6:1197364/1‑149 (MQ=255)

                                                                                                                                                                                                               gCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCAGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACCTGTTGCATCACTTCACCTTCAGTTTGTGCTGCATCTTTAAATCATTTCATTGCGCTATTAAGCGGAGTTGGCAGCCGCTTCGGATCGTCCTCCATCTTCTCTGACAGCGGCGCATGTGCtttc                                                    >  4:77101/1‑272 (MQ=255)

                                                                                                                                                                                                                    aaCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGt                                                                                                                                                                            <  7:497612/149‑1 (MQ=255)

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

                                                                                                                                                                                                                       ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGc                                                                                                                                                                         <  6:581796/149‑1 (MQ=255)

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

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

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

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

                                                                                                                                                                                                                           aTTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCa                                                                                                                                                                     >  5:990725/1‑149 (MQ=255)

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

                                                                                                                                                                                                                            ttAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccatc                                                                                                                                                                                  <  8:1067843/135‑1 (MQ=255)

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

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

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

                                                                                                                                                                                                                                    tgatACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGtt                                                                                                                                                                           >  7:1120641/4‑134 (MQ=255)

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

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

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

                                                                                                                                                                                                                                        aCTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGt                                                                                                                                                        >  5:197473/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                             tgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgtg                                                                                                                                                   >  8:273584/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                             tgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGACCATCCACCGGCATCCCGGAACGGCCATCCATCACATGTTGCATCCCTTCAGCGCCGGTTtgtg                                                                                                                                                   >  5:510594/1‑149 (MQ=255)

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

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

                                                                                                                                                                                                                                               aTTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgct                                                                                                                                                 >  6:1465477/1‑149 (MQ=255)

                                                                                                                                                                                                                                               aTTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgct                                                                                                                                                 <  5:583639/149‑1 (MQ=255)

                                                                                                                                                                                                                                                ttttttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                                                                                                                                                            >  6:693002/1‑73 (MQ=255)

                                                                                                                                                                                                                                                  aaatCCCCGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                                                                                                                                                            >  6:1668820/4‑71 (MQ=25)

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

                                                                                                                                                                                                                                                      ccccGGGACATGCCGTGTCTCGGGGATTTTTTTATCCGTGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCg            <  4:130319/275‑1 (MQ=255)

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

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

                                                                                                                                                                                                                                                        ccGCGACACGGCATGTCGCGGGGATTATTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa                                                                                                                                        >  6:596408/1‑149 (MQ=255)

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

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

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

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

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

                                                                                                                                                                                                                                                                  ccGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattgca                                                                                                                              >  7:1684306/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    gtgtCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattgcatt                                                                                                                            >  8:720756/1‑149 (MQ=255)

TTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGAT  >  NC_000913/1169124‑1169654

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