BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,335,655	T→A	5.1%	intergenic (‑65/+39)	eptA ← / ← adiC	lipid A phosphoethanolamine transferase/arginine:agmatine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,335,655	0	T	A	5.1%	1039.1 / 32.3	292	intergenic (‑65/+39)	eptA/adiC	lipid A phosphoethanolamine transferase/arginine:agmatine antiporter
Reads supporting (aligned to +/- strand): ref base T (126/150); new base A (9/6); total (135/157)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.41e-01

Only 100 of 296 total aligned reads displayed.

CTTAGTGTCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCAC  >  NC_000913/4335385‑4335926

cTTAGTGTCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTa                                                                                                                                                                                                                                                                             >  2:509502/1‑275 (MQ=255)

                    aaTCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTg                                                                                                                                                                                                                                                          <  2:239235/274‑1 (MQ=255)

                                                                          cGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGACAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTTTTGGTGCa                                                                                                                                                                                                                     >  2:7194/1‑257 (MQ=255)

                                                                                                                  gCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTg                                                                                                                                                                                                                                                          >  1:64166/1‑180 (MQ=255)

                                                                                                                  gCAATATTANGGCAGATAGAGATATAAAAAGCCGGCAACAANGGCCATGCNNGTAAATTCAAAGCGTGTCTTTTTNGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCTCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTGTTATGTTTg                                                                                                                                                                                                                                                          <  3:64166/180‑1 (MQ=255)

                                                                                                                                 aTAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAgg                                                                                                                                                                                                                                                                          >  6:1418167/1‑149 (MQ=255)

                                                                                                                                  tAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGt                                                                                                                                                                                                                                                                         >  7:1357478/1‑149 (MQ=255)

                                                                                                                                   agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGtt                                                                                                                                                                                                                                                                        <  5:1577526/149‑1 (MQ=255)

                                                                                                                                   agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGtt                                                                                                                                                                                                                                                                        <  7:379868/149‑1 (MQ=255)

                                                                                                                                     agaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTgg                                                                                                                                                                                                                                                                      <  6:1222964/149‑1 (MQ=255)

                                                                                                                                     agaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTC‑TTTTAGTAGGCGCTTCAACATCACGGTGGTTCCGTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGc                                                                                                                                                                                                                                                                     >  7:1593935/1‑149 (MQ=255)

                                                                                                                                         ataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCttt                                                                                                                                                                                                                                                                  >  5:1255100/1‑149 (MQ=255)

                                                                                                                                           aaaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTCGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTATGCCAACCTTAAGTTCTTAAGGTTGGCTTTTa                                                                                                                                                                                                                                                                <  8:1033938/149‑1 (MQ=255)

                                                                                                                                              aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGt                                                                                                                                                                                                                                                             <  7:1653275/149‑1 (MQ=255)

                                                                                                                                               aGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGtt                                                                                                                                                                                                                                                            >  5:511280/1‑149 (MQ=255)

                                                                                                                                                  ggacacaaaTATCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGAAGGCGCTTTAACCTCACGGTGTTTCCATCTAACAAAGTGCGCAGATGCTTGAAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCGTTTATGTTTGCTGGATTAAGGGGAATTAAGCTTGGCTTAGTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAAc                                                                                                                            <  4:116155/267‑1 (MQ=255)

                                                                                                                                                  ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACctct                                                                                                                        >  1:277926/1‑278 (MQ=255)

                                                                                                                                                   gCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATg                                                                                                                                                                                                                                                              >  4:371637/1‑143 (MQ=255)

                                                                                                                                                       acaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTg                                                                                                                                                                                                                                                          <  2:371641/143‑1 (MQ=255)

                                                                                                                                                           tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAg                                                                                                                                                                                                                                                <  5:52560/149‑1 (MQ=255)

                                                                                                                                                             gCCACGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCACTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAATGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAggg                                                                                                                                                                                                                                              <  8:1421186/149‑1 (MQ=255)

                                                                                                                                                                aGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGaa                                                                                                                                                                                                                                           >  8:404605/1‑149 (MQ=255)

                                                                                                                                                                aGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGaa                                                                                                                                                                                                                                           >  7:679384/1‑149 (MQ=255)

                                                                                                                                                                  gCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAg                                                                                                                                                                                                                                                <  1:170558/142‑1 (MQ=255)

                                                                                                                                                                  gCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTg                                                                                                                                                                                                                                 >  4:277214/1‑157 (MQ=255)

                                                                                                                                                                  gCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGNTTATTGGTGCATCTAAGGGATACGGTTTTTTATGNTNCCGGTTNNAAt                                                                                                                                                                                >  4:734170/1‑206 (MQ=255)

                                                                                                                                                                  gCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAtt                                                                                                                                                                               <  3:84582/207‑1 (MQ=255)

                                                                                                                                                                  gCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAtt                                                                                                                                                                               >  1:84582/1‑207 (MQ=255)

                                                                                                                                                                    gAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGc                                                                                                                                                                                                                                                         <  2:644677/131‑1 (MQ=255)

                                                                                                                                                                      gTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTg                                                                                                                                                                                                                                                          <  4:499127/128‑1 (MQ=255)

                                                                                                                                                                         caTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCttt                                                                                                                                                                                                                                  <  5:46184/148‑1 (MQ=255)

                                                                                                                                                                         aaTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCttt                                                                                                                                                                                                                                  >  7:528163/1‑149 (MQ=255)

                                                                                                                                                                           ttCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGc                                                                                                                                                                                                                                <  5:1468598/149‑1 (MQ=255)

                                                                                                                                                                           ttCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGc                                                                                                                                                                                                                                >  7:90519/1‑149 (MQ=255)

                                                                                                                                                                           ttCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGc                                                                                                                                                                                                                                >  6:834716/1‑149 (MQ=255)

                                                                                                                                                                                agagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAtt                                                                                                                                                                                                                           <  8:528163/149‑1 (MQ=255)

                                                                                                                                                                                agagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAg                                                                                                                                                                                           <  4:602912/181‑1 (MQ=255)

                                                                                                                                                                                 gagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                                                                                                                                                                                                          >  8:120934/1‑149 (MQ=255)

                                                                                                                                                                                  agGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAAt                                                                                                                                                                                                                                      <  5:968040/136‑1 (MQ=255)

                                                                                                                                                                                   gggTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                                                                                                                                                        <  6:1152502/149‑1 (MQ=255)

                                                                                                                                                                                      tCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGtt                                                                                                                                                                                                                                                                        <  7:44404/98‑1 (MQ=255)

                                                                                                                                                                                       cTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                                                                                                                                                        <  8:1284954/145‑1 (MQ=255)

                                                                                                                                                                                        tttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATc                                                                                                                                                                                                                   >  6:1557780/1‑149 (MQ=255)

                                                                                                                                                                                          tttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCtt                                                                                                                                                                                                                                   <  5:733052/131‑1 (MQ=255)

                                                                                                                                                                                            tagtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAg                                                                                                                                                                                                               <  5:1685243/149‑1 (MQ=255)

                                                                                                                                                                                                agGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAtt                                                                                                                                                                                                                           >  4:713369/1‑133 (MQ=255)

                                                                                                                                                                                                agGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGt                                                                                                                                                                                   >  2:781107/1‑173 (MQ=255)

                                                                                                                                                                                                  gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAAt                                                                                                                                                                                                                                          >  6:469700/1‑116 (MQ=255)

                                                                                                                                                                                                  gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGc                                                                                                                                                                                                                      >  8:796863/1‑136 (MQ=255)

                                                                                                                                                                                                    gcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACgg                                                                                                                                                                                                       >  7:185540/1‑149 (MQ=255)

                                                                                                                                                                                                    gcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAtt                                                                                                                                                                               <  4:781101/173‑1 (MQ=255)

                                                                                                                                                                                                     cTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATc                                                                                                                                                                                                                                     <  8:1569116/118‑1 (MQ=255)

                                                                                                                                                                                                      ttCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGt                                                                                                                                                                                                     <  7:585664/149‑1 (MQ=255)

                                                                                                                                                                                                        cAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTc                                                                                                                                                                                                                                                                                >  6:992922/1‑72 (MQ=255)

                                                                                                                                                                                                           cATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                                                                                                                                                                                                >  5:1530044/1‑149 (MQ=255)

                                                                                                                                                                                                           cATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                                                                                                                                                                                                <  8:233563/149‑1 (MQ=255)

                                                                                                                                                                                                           cATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                                                                                                                                                                                                <  5:637678/149‑1 (MQ=255)

                                                                                                                                                                                                            aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                                                                                                                                                                                               >  6:698/1‑149 (MQ=255)

                                                                                                                                                                                                            aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                                                                                                                                                                                               >  8:684567/1‑149 (MQ=255)

                                                                                                                                                                                                             tCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                                                                                                                                                        <  5:1375370/123‑1 (MQ=255)

                                                                                                                                                                                                             tCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                                                                                                                                                        >  6:1375370/1‑123 (MQ=255)

                                                                                                                                                                                                              cACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGt                                                                                                                                                                                             >  8:751620/1‑149 (MQ=255)

                                                                                                                                                                                                               aCGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTa                                                                                                                                                                                            <  5:85211/149‑1 (MQ=255)

                                                                                                                                                                                                                cGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAgg                                                                                                                                                                                                              >  2:300105/1‑130 (MQ=255)

                                                                                                                                                                                                                      ttCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGc                                                                                                                                                                                                                      <  6:286277/116‑1 (MQ=255)

                                                                                                                                                                                                                          aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                                                                                                                                                                                                <  5:648193/134‑1 (MQ=255)

                                                                                                                                                                                                                          aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                                                                                                                                                 >  7:565756/1‑149 (MQ=255)

                                                                                                                                                                                                                             gAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTc                                                                                                                                                                              >  8:625644/1‑149 (MQ=255)

                                                                                                                                                                                                                             gAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACCGCCAGATATGCCggg                                                >  2:122273/1‑275 (MQ=255)

                                                                                                                                                                                                                              aaCAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCa                                                                                                                                                                             >  6:739633/1‑149 (MQ=255)

                                                                                                                                                                                                                                cAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAgc                                                                                                                                                                           <  5:200513/149‑1 (MQ=255)

                                                                                                                                                                                                                                             tGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCgg                                                                                                                                                              >  7:794841/1‑149 (MQ=255)

                                                                                                                                                                                                                                                tCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtga                                                                                                                                                           <  7:29379/149‑1 (MQ=255)

                                                                                                                                                                                                                                                 cGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgat                                                                                                                                                          <  8:575319/149‑1 (MQ=255)

                                                                                                                                                                                                                                                 cGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGa                                                                                            <  1:240552/211‑1 (MQ=255)

                                                                                                                                                                                                                                                 cGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGa                                                                                            >  3:240552/1‑211 (MQ=255)

                                                                                                                                                                                                                                                  gcaaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatg                                                                                                                                                         >  5:404061/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       gcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGAACAt                                                                                                                                                    >  6:927841/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        caaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATc                                                                                                                                                   <  8:888062/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         ttataaaaGACAGCCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCg                     <  2:552698/268‑1 (MQ=255)

                                                                                                                                                                                                                                                           aaTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAgg                                                                                                                                                <  6:226452/149‑1 (MQ=255)

                                                                                                                                                                                                                                                            aTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGATTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAggg                                                                                                                                               <  6:1473882/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             ttAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGt                                                                                                                                              >  6:373186/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                aGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGAc                                                                                                                                           <  7:817031/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 gCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCgg                                                                                                                                                              >  8:319060/1‑129 (MQ=255)

                                                                                                                                                                                                                                                                 gCCAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACa                                                                                                                                          >  8:1628738/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  ccAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaa                                                                                                                                         >  7:782237/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   cAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCgg                                                                                                                                                                                       <  7:1270051/102‑1 (MQ=255)

                                                                                                                                                                                                                                                                   cAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaaa                                                                                                                                        >  5:355278/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   cAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaaa                                                                                                                                        >  7:1192268/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    aaCCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGc                                                                                         <  4:533536/195‑1 (MQ=255)

                                                                                                                                                                                                                                                                    aaCCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAAt                                                                                                                                       >  7:1298330/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     aCCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTa                                                                                                                                                                                  <  5:1157990/105‑1 (MQ=255)

                                                                                                                                                                                                                                                                      ccTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGa                                                                                                                                     >  6:159956/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      ccTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGa                                                                                                                                     <  6:521394/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                        ttAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGAcc                                                                                                                                   <  8:1357478/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                        ttAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGAcc                                                                                                                                   <  5:748546/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                          aaGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGca   <  2:591837/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                             ttCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATa                                                                                                                              >  7:339135/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                              tCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATaa                                                                                                                             <  5:373186/149‑1 (MQ=255)

CTTAGTGTCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCAC  >  NC_000913/4335385‑4335926

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