BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,318,095	C→G	100%	L594F (TTG→TTC)	rcsC ←	hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,318,095	0	C	G	99.7%	976.6 / ‑6.2	266	L594F (TTG→TTC)	rcsC	hybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand): ref base C (1/1); major base G (122/139); minor base A (0/2); total (124/142)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 266 total aligned reads displayed.

AAACTGGCTGCCCATTCCCGGTTCTGAATCTACCGAGATATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGTCACCAGCCGATCGACG  >  NC_000913/2317825‑2318368

aaaCTGGCTGCCCTATCCCGGTCTTGAATCTAACGGGATAACGCCGCCCACAATGCAGGTCAGTGTTGCACAAATCGCCCGACCCAGACCGTTCCCCTGGAAATGACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGAGCAAACAAGCGAACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGAGAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGaa                                                                                                                                                                                                                                                                                 <  4:544706/273‑1 (MQ=255)

         gCCCATTCCCGGTTCTGAATCTACCGAGATATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCACATGGAAATTACGCTGTACGCCCGTTCTGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGTTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAAACGCCATCCGCGCGAACATGCAAAACTTTACAGCCGGTATCGGTGCATTTTTTGGCGTTTCTGAACTGGTTggca                                                                                                                                                                                                                                                                       >  4:186160/1‑272 (MQ=255)

                                                          aTCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTCCTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGt                                                                                                                                                                                                                 >  1:162317/1‑279 (MQ=255)

                                                                                            tCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCTCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGGTACTGAAAAGGTTGGAGATTACCTGCGGTAAACGCATCGGGGCGCCATTTAAGGCCACTGGCCCATCCGGTTCAAt                                                                                                                                                                                                            >  4:66280/1‑250 (MQ=255)

                                                                                            tCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGCATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACAt                                                                                                                                                                                                                      >  2:671196/1‑240 (MQ=255)

                                                                                                cTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTa                                                                                                                                                                               >  2:210768/1‑275 (MQ=255)

                                                                                                  ggAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCCCTTCTTTCGCCCGTATCCCCACGCCGGTATCGCGAACACGGATAGAAAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTAAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTaccac                                                                                                                                                                           <  1:20353/277‑1 (MQ=255)

                                                                                                           ggcgTACGCCCGTTCCGACATGGAAGAAGGGGTCAAAAAAGCGCAACACTTCGTTGGCCGGTATCCCCACGCGGGTATCGCGAACAAGGATATGGAGGTAATGGCCAACCGCGCAAACATGCAAAACTATACAGCCGGTATCGATGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCGGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGt                                                                                                                                                                    <  4:62714/272‑1 (MQ=255)

                                                                                                              gTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATTCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCAACGGGTCACCATTTAAGGCCACGTGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTTCTTGCGTAacaccacgggttat                                                                                                                                                                 >  4:33178/1‑267 (MQ=255)

                                                                                                                 cGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTGTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGACGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTgg                                                                                                                                                          <  1:108850/279‑1 (MQ=255)

                                                                                                                            aCCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGaa                                                                                                                                                                                                                                                                                 <  7:606657/149‑1 (MQ=255)

                                                                                                                            aCCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGaa                                                                                                                                                                                                                                                                                 >  6:79292/1‑149 (MQ=255)

                                                                                                                            aCCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGAc                                                                                                                                                                                                                                                                  >  3:21780/1‑164 (MQ=255)

                                                                                                                                 gaagaaGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGt                                                                                                                                                                                                                                                                            <  6:1352074/149‑1 (MQ=255)

                                                                                                                                    gaagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTg                                                                                                                                                                                                                                                                          >  6:33079/1‑148 (MQ=255)

                                                                                                                                    gaagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGATACTGAACAGGTTgg                                                                                                                                                                                                                                                                         <  8:1460058/149‑1 (MQ=255)

                                                                                                                                     aagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGga                                                                                                                                                                                                                                                                        <  7:1181909/149‑1 (MQ=255)

                                                                                                                                       ggaaTCAAACAAGCGAACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGc                                                                                                                                                                                                                          <  1:99564/190‑1 (MQ=255)

                                                                                                                                          aTCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGgag                                                                                                                                                                                                                                                                       >  5:168216/1‑145 (MQ=255)

                                                                                                                                              aaCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACctg                                                                                                                                                                                                                                                               >  5:1481587/1‑149 (MQ=255)

                                                                                                                                               aCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGa                                                                                                                                                                                                                                                                   <  8:1242524/144‑1 (MQ=255)

                                                                                                                                               aCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACctgc                                                                                                                                                                                                                                                              >  5:351693/1‑149 (MQ=255)

                                                                                                                                               aCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACctgc                                                                                                                                                                                                                                                              <  6:1309514/149‑1 (MQ=255)

                                                                                                                                                cacgcgcACCCCTGCGTTCCACGGTATCCCCACGCCGGCAGCGCGAACACGGATATAGGGATAATCGCCATCCGCGCGAAAATGCAAAACTGTACAGCCGGTATCGGTGAATTTTCTGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCAGCGCGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACTAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCAcgc                                                                                                                                <  4:242075/271‑1 (MQ=255)

                                                                                                                                                   gcgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGga                                                                                                                                                                                                                                                                        >  5:749591/1‑135 (MQ=255)

                                                                                                                                                   gcgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTa                                                                                                                                                                                                                                                          <  6:274837/149‑1 (MQ=255)

                                                                                                                                                    cgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTaa                                                                                                                                                                                                                                                         >  7:1509851/1‑149 (MQ=255)

                                                                                                                                                          acTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCAt                                                                                                                                                                                                                                                   >  8:1230953/1‑149 (MQ=255)

                                                                                                                                                           cttcttTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATc                                                                                                                                                                                                                                                  >  7:44095/1‑149 (MQ=255)

                                                                                                                                                                ttCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTc                                                                                                                                                                                                                                             <  5:1255582/149‑1 (MQ=255)

                                                                                                                                                                ttCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTc                                                                                                                                                                                                                                             <  8:802541/149‑1 (MQ=255)

                                                                                                                                                                 tCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCAcg                                                                                                                  >  2:161052/1‑271 (MQ=255)

                                                                                                                                                                  cGCCGGTACCCCCACGCCGGTATCGGGAACACGGATAGAGAGGTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGtt                                                                                                             <  2:588998/275‑1 (MQ=255)

                                                                                                                                                                   gCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGga                                                                                                                                                                                                                                                                        >  8:971458/1‑119 (MQ=255)

                                                                                                                                                                   gCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAAc                                                                                                                                                                                                                                                       <  7:402954/136‑1 (MQ=255)

                                                                                                                                                                   gCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAAc                                                                                                                                                                                                                                                       >  8:402954/1‑136 (MQ=255)

                                                                                                                                                                      ggTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCAtt                                                                                                                                                                                                                                       <  7:688158/149‑1 (MQ=255)

                                                                                                                                                                       gTATCCCCCCGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATc                                                                                                        <  2:369022/275‑1 (MQ=255)

                                                                                                                                                                         aTCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTaa                                                                                                                                                                                                                                    <  6:99684/149‑1 (MQ=255)

                                                                                                                                                                          tCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAg                                                                                                                                                                                                                                   <  6:740524/149‑1 (MQ=255)

                                                                                                                                                                            cccACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGc                                                                                                                                                                                                                                 >  6:1497819/1‑149 (MQ=255)

                                                                                                                                                                               aCGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCAc                                                                                                                                                                                                                              <  6:928709/149‑1 (MQ=255)

                                                                                                                                                                                 gCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTg                                                                                                                                                                                                                            <  5:226541/149‑1 (MQ=255)

                                                                                                                                                                                 gCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTg                                                                                                                                                                                                                            >  7:165210/1‑149 (MQ=255)

                                                                                                                                                                                  ccGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTgg                                                                                                                                                                                                                           <  5:79292/149‑1 (MQ=255)

                                                                                                                                                                                   cGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGc                                                                                                                                                                                                                          <  5:1166113/149‑1 (MQ=255)

                                                                                                                                                                                        tCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATc                                                                                                                                                                                                                     <  6:587507/149‑1 (MQ=255)

                                                                                                                                                                                        tCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATc                                                                                                                                                                                                                     >  7:1603951/1‑149 (MQ=255)

                                                                                                                                                                                           cgAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCgg                                                                                                                                                                                                                  <  5:562227/149‑1 (MQ=255)

                                                                                                                                                                                           cgAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCgg                                                                                                                                                                                                                  >  7:1234332/1‑149 (MQ=255)

                                                                                                                                                                                           cgAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCgg                                                                                                                                                                                                                  <  7:616429/149‑1 (MQ=255)

                                                                                                                                                                                            gAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGt                                                                                                                                                                                                                 >  7:397869/1‑149 (MQ=255)

                                                                                                                                                                                                acGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCaa                                                                                                                                                                                                             <  6:1702114/149‑1 (MQ=255)

                                                                                                                                                                                                  ggATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACGGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAg                                                                                                                                                                                                                                   >  7:965432/1‑125 (MQ=255)

                                                                                                                                                                                                     tAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAg                                                                                                                                                                                                        <  7:898629/149‑1 (MQ=255)

                                                                                                                                                                                                       gagagaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCa                                                                                                                                                                                                      <  8:1603951/149‑1 (MQ=255)

                                                                                                                                                                                                       gagagaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCAACGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCa                                                                                                                                                                                                      <  6:1591379/149‑1 (MQ=255)

                                                                                                                                                                                                         gagaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGt                                                                                                                                                                                                    <  7:1300727/149‑1 (MQ=255)

                                                                                                                                                                                                             tAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACaa                                                                                                                                                                                                <  6:1566503/149‑1 (MQ=255)

                                                                                                                                                                                                              aaTCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCAt                                                                                                                                                                                                                                        >  4:148174/1‑108 (MQ=255)

                                                                                                                                                                                                               aTCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGTATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGc                                                                                                                                                                                              >  8:148629/1‑149 (MQ=255)

                                                                                                                                                                                                                tCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGc                                                                                                                                                                                                                                           <  6:572552/103‑1 (MQ=255)

                                                                                                                                                                                                                  gCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTaa                                                                                                                                                                                                                                    <  2:148176/108‑1 (MQ=255)

                                                                                                                                                                                                                  gCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTa                                                                                                                                                                                           <  5:279905/149‑1 (MQ=255)

                                                                                                                                                                                                                   ccATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTaa                                                                                                                                                                                          <  5:1419025/149‑1 (MQ=255)

                                                                                                                                                                                                                   ccATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACAAGCGTTAAATAGTTGGCGGTGATCTGGTTCATCACTTCCCGCGGTGAAAACTCCcgc                                                                                                                 >  4:309419/1‑222 (MQ=255)

                                                                                                                                                                                                                    cATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAAc                                                                                                                                                                                         <  8:268512/149‑1 (MQ=255)

                                                                                                                                                                                                                         gcgcgAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCtt                                                                                                                                                                                    <  7:203828/149‑1 (MQ=255)

                                                                                                                                                                                                                         gcgcgAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCtt                                                                                                                                                                                    <  6:753810/149‑1 (MQ=255)

                                                                                                                                                                                                                          cgcgAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACt                                                                                                                      >  4:377710/1‑210 (MQ=255)

                                                                                                                                                                                                                             gAACATGCAAAACTATACTGCCGGTATCTGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTGAGGCCACTGGCACATCCGGTTCaa                                                                                                                                                                                                             >  8:1035350/1‑120 (MQ=255)

                                                                                                                                                                                                                              aaCATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCAcg                                                                                                                  <  2:377714/210‑1 (MQ=255)

                                                                                                                                                                                                                                   gCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGtt                                                                                                                                                                                                                                                                           <  6:1160390/52‑1 (MQ=255)

                                                                                                                                                                                                                                   gCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCAt                                                                                                                                                                                                                                        <  8:1539179/87‑1 (MQ=255)

                                                                                                                                                                                                                                   gCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTaccacc                                                                                                                                                                          >  6:1021651/1‑149 (MQ=255)

                                                                                                                                                                                                                                     aaaaCTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAg                                                                                                                                                                        >  8:945561/1‑149 (MQ=255)

                                                                                                                                                                                                                                       aaCTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCtt                                                                                                                                                                                    >  5:702380/1‑135 (MQ=255)

                                                                                                                                                                                                                                               aGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCa                                                                                                                                                                                                                               <  4:231837/84‑1 (MQ=255)

                                                                                                                                                                                                                                               aGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAg                                                                                                                                                              >  8:1673150/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 ccGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGAc                                                                                                                                                                                                                                                                  <  5:1555625/47‑1 (MQ=255)

                                                                                                                                                                                                                                                 ccGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAg                                                                                                                                                              >  2:750931/1‑147 (MQ=255)

                                                                                                                                                                                                                                                  cGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGTTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCTCGGATAATTTTCAACAACAGGCTGGAAGAgtt                              >  2:304416/1‑274 (MQ=255)

                                                                                                                                                                                                                                                   ggTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTgg                                                                                                                                                          <  8:1343915/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   ggTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTgg                                                                                                                                                          <  6:640516/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   ggTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCg                                                                                                           <  2:638579/196‑1 (MQ=255)

                                                                                                                                                                                                                                                    gTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGc                                                                                                                                                         <  8:397869/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     tATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTgg                                                                                                                                                          <  4:750925/147‑1 (MQ=255)

                                                                                                                                                                                                                                                     tATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGtt                                                                                                             <  2:569998/192‑1 (MQ=255)

                                                                                                                                                                                                                                                          gTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAg                                                                                                                                                                                                        <  8:607691/96‑1 (MQ=255)

                                                                                                                                                                                                                                                          gTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGAt                                                                                                                                                   >  6:283209/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          gTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGAt                                                                                                                                                   <  6:1481587/149‑1 (MQ=255)

                                                                                                                                                                                                                                                           tGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCCCTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTCCCACCAGCGGTAAATAGTTGGCGGTGAtg                                                                                                                                                  <  5:1364336/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              aTTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGCTTCGCTCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGTTGGAAGAGTTTTTCATTcccat                  >  4:121822/1‑272 (MQ=255)

                                                                                                                                                                                                                                                               ttttATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGtt                                                                                                             <  2:298212/182‑1 (MQ=255)

                                                                                                                                                                                                                                                                     ggCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTtcatca                                                                                                                                        >  5:1154039/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     ggCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTtcatca                                                                                                                                        >  6:1376978/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                         ttACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTc                                                                                                                                    >  6:1535145/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                         ttACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCAcg                                                                                                                                 >  1:258298/1‑152 (MQ=255)

                                                                                                                                                                                                                                                                         ttACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCtt                                                                                 >  2:429606/1‑200 (MQ=255)

                                                                                                                                                                                                                                                                         ctACTGAACAGGTCGGTGATGACCTGCTGTAAACGCCTCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTTCTTGCGTACCACTAGCGGTAAAAAGTTGGCGTTGATGTGGNTCATCCCTTCCCTCGGTGTAAACTCACGCGGTTCGATCTCTAACTGTTCTGATTCAATTTTCGATAATACGATAATATTCNCGTTATATTACAACAACAAGCTGGGNAAGGTTTTTATTTCTGTCACCAcccacacacac  >  3:169267/2‑271 (MQ=255)

AAACTGGCTGCCCATTCCCGGTTCTGAATCTACCGAGATATCGCCGTCCATCATGCTGATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGTCACCAGCCGATCGACG  >  NC_000913/2317825‑2318368

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