BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,074,691	C→G	5.1%	intergenic (+33/‑88)	flu → / → yeeR	CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter/CP4‑44 prophage; putative membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,074,691	0	C	G	5.1%	844.2 / 26.4	237	intergenic (+33/‑88)	flu/yeeR	CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter/CP4‑44 prophage; putative membrane protein
Reads supporting (aligned to +/- strand): ref base C (105/120); new base G (2/10); total (107/130)
Fisher's exact test for biased strand distribution p-value = 7.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

Only 100 of 250 total aligned reads displayed.

TGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACA‑GAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCA  >  NC_000913/2074424‑2074938

tGTTATCCGCACTCTCAGCTCCCGGGGAGATAGGCGTGTGGGGACTTCCACTGCAGGCAGCTGGATGAAGTTCTCTCCCTCACA‑GAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCg                                                                                                                                                                                                                                                  <  4:642041/275‑1 (MQ=255)

         cACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACA‑GAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCCGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCCCCCCTCTGTTGTCCCGGGCATCATGCCCgtcctctta                                                                                                                                                                                                                                          >  4:509667/1‑266 (MQ=255)

                                                 aCTGCAGGCAGCGGGATGACGTTCTCTCCCTCACA‑GAAGGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCAtt                                                                                                                                                                                                  <  2:655280/274‑1 (MQ=255)

                                                   tGCAGGCAGCGGGATGACGTTCTCTCCCTCACA‑GAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCg                                                                                                                                                                                                                                                  <  2:31481/224‑1 (MQ=255)

                                                      aggcagCGGGATGACGTTCTCTCCATCAAA‑GAAGGGCACATCACTGGAACAGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGGCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACaa                                                                                                                                                                                             <  2:540722/274‑1 (MQ=255)

                                                                     gatctctcCC‑CACAGGAATGGCAAATCACTGGACATGCAGGCCGGACGGGAAGCCCGTGTCCGGGAAAATATGACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGGCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtcat                                                                                                                                                                              <  4:194113/272‑1 (MQ=255)

                                                                                   a‑GAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCCCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCCACGCCACATTCGCTGGCATTAACAATAACATGATATTCCTCACGGAGTGatcatg                                                                                                                                                               >  2:404498/1‑270 (MQ=255)

                                                                                                  cTGGCCGTGCAGGCCTCACCGGAAGCCCGTTTCCGTGAATATCTCCCCCTGGGCGTTCAGGCCGCGTTTGCCCACTGCGTAAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGCCAGAACCATCGCCCCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGACgg                                                                                                                                                 <  4:293175/275‑1 (MQ=255)

                                                                                                          gCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCa                                                                                                                                                                                                    >  4:388302/1‑216 (MQ=255)

                                                                                                                                gtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCgg                                                                                                                                                                                                                                                 <  5:1144078/149‑1 (MQ=255)

                                                                                                                                 tCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCGTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCggg                                                                                                                                                                                                                                                <  6:60789/149‑1 (MQ=255)

                                                                                                                                  ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCGGGa                                                                                                                                                                                                                                               <  8:1415001/149‑1 (MQ=255)

                                                                                                                                  ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCGGGa                                                                                                                                                                                                                                               <  5:1325587/149‑1 (MQ=255)

                                                                                                                                  ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTGATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTGATGACAGAACCATCGCCTCTCTGTTGTCCCGGTCATCATGACCGGGa                                                                                                                                                                                                                                               <  8:1672977/149‑1 (MQ=255)

                                                                                                                                   cGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAGCGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGAc                                                                                                                                                                                                                                              <  8:30493/149‑1 (MQ=255)

                                                                                                                                   cGGGAAAATATCACCCGGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGTGTATAACGGTCAGGCCACACTGGATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCGGGAc                                                                                                                                                                                                                                              <  6:485322/149‑1 (MQ=255)

                                                                                                                                      gAAAATATCACCCTGGGCGTTCGGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCg                                                                                                                                                                                                                                           >  5:1420559/1‑149 (MQ=255)

                                                                                                                                         aaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGAc                                                                                                                                                                                                                                        <  6:761870/149‑1 (MQ=255)

                                                                                                                                             tCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGgc                                                                                                                                                                                                                                    <  6:678090/149‑1 (MQ=255)

                                                                                                                                                  cTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAAc                                                                                                                                                                                                                               >  8:1551999/1‑149 (MQ=255)

                                                                                                                                                   tGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACg                                                                                                                                                                                                                              >  5:569332/1‑149 (MQ=255)

                                                                                                                                                      gCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGAt                                                                                                                                                                                                                           >  7:865906/1‑149 (MQ=255)

                                                                                                                                                      gCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGCGACCCTCTGACAGCACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGAt                                                                                                                                                                                                                           <  8:1038209/149‑1 (MQ=255)

                                                                                                                                                       cGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATc                                                                                                                                                                                                                          <  7:1476021/149‑1 (MQ=255)

                                                                                                                                                       cGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATc                                                                                                                                                                                                                          <  5:1200074/149‑1 (MQ=255)

                                                                                                                                                         ttCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCg                                                                                                                                                  <  2:363981/219‑1 (MQ=255)

                                                                                                                                                           cAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGTTAGTCGCCGCGCTGATCCTGCTGATCGCAGGATTTGCCAGTCTTTGCCTTTTGTTCAGAGCCtt                                                                                        >  2:201670/1‑275 (MQ=255)

                                                                                                                                                            aGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCaa                                                                                                                                                                                                                     >  7:783443/1‑149 (MQ=255)

                                                                                                                                                            aGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCaa                                                                                                                                                                                                                     <  6:1627280/149‑1 (MQ=255)

                                                                                                                                                              gCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACg                                                                                                                                                                                                                   >  6:830836/1‑149 (MQ=255)

                                                                                                                                                                cGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCa                                                                                <  1:118646/278‑1 (MQ=255)

                                                                                                                                                                    tATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACAt                                                                                                                                                                                                             <  8:1315013/149‑1 (MQ=255)

                                                                                                                                                                        cccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGc                                                                                                                                                                                                         >  7:513690/1‑149 (MQ=255)

                                                                                                                                                                        cccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGc                                                                                                                                                                                                         >  6:278646/1‑149 (MQ=255)

                                                                                                                                                                        cccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGc                                                                                                                                                                                                         <  6:1507236/149‑1 (MQ=255)

                                                                                                                                                                         ccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCt                                                                                                                                                                                                        >  5:153264/1‑149 (MQ=255)

                                                                                                                                                                         ccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCt                                                                                                                                                                                                        >  6:1372999/1‑149 (MQ=255)

                                                                                                                                                                           acaGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTgg                                                                                                                                                                                                      >  7:845662/1‑149 (MQ=255)

                                                                                                                                                                              gCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCt                                                                                                                                                                                                                         >  8:368141/1‑127 (MQ=255)

                                                                                                                                                                                     cGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAAt                                                                                                                                                                                            >  8:356459/1‑149 (MQ=255)

                                                                                                                                                                                          gcagcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCACCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACAt                                                                                                                                                                                       <  6:1258152/149‑1 (MQ=255)

                                                                                                                                                                                           cagcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcatca                                                                                                                                                                                                                                                        >  1:248288/1‑83 (MQ=255)

                                                                                                                                                                                            agcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCt                                                                                                                                                                                                        <  6:1442833/130‑1 (MQ=255)

                                                                                                                                                                                            agcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGa                                                                                                                                                                                     <  7:1573493/149‑1 (MQ=255)

                                                                                                                                                                                             gcgcTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGAc                                                                                                                                                                                                                                                    >  8:1624693/1‑85 (MQ=255)

                                                                                                                                                                                             gcgcTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGAACGGGACCCGGACCGCCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAAAAATAACATGat                                                                                                                                                                                    >  6:86940/1‑149 (MQ=255)

                                                                                                                                                                                                   aaGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGc                                                                                                                                                                                                         >  7:296092/1‑122 (MQ=255)

                                                                                                                                                                                                      ggTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCAc                                                                                                                                                                           <  7:192902/149‑1 (MQ=255)

                                                                                                                                                                                                          taACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAg                                                                                                                                                                       >  8:507489/1‑149 (MQ=255)

                                                                                                                                                                                                              ggTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGAc                                                                                                                                                                   <  6:386930/149‑1 (MQ=255)

                                                                                                                                                                                                                tCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTa                                                                                                                                                                 <  7:680300/149‑1 (MQ=255)

                                                                                                                                                                                                                 cAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTAt                                                                                                                                                                <  7:925325/149‑1 (MQ=255)

                                                                                                                                                                                                                  aGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCg                                                                                                                                                                                                          <  1:149081/106‑1 (MQ=255)

                                                                                                                                                                                                                  aGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCg                                                                                                                                                                                                          >  3:149081/1‑106 (MQ=255)

                                                                                                                                                                                                                  aGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTTACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCttt                                                                                                     >  2:238553/1‑207 (MQ=255)

                                                                                                                                                                                                                  aGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATcc                                                                                                                                      >  2:334048/1‑174 (MQ=255)

                                                                                                                                                                                                                    gCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTCACTATGtt                                                                                                                                                             >  6:1056170/1‑149 (MQ=255)

                                                                                                                                                                                                                      cacaCTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTTACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGtt                                                                                                 <  4:238550/207‑1 (MQ=255)

                                                                                                                                                                                                                         acTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGa                                                                                                                                                        <  5:1309756/149‑1 (MQ=255)

                                                                                                                                                                                                                           tGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATa                                                                                                                                                      <  8:665427/149‑1 (MQ=255)

                                                                                                                                                                                                                            gAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAg                                                                                                                                                     >  8:499056/1‑149 (MQ=255)

                                                                                                                                                                                                                                     cTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATTTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgcgcgc                                                                                                                                            >  7:112510/1‑149 (MQ=255)

                                                                                                                                                                                                                                      ttCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCt                                                                                                                                           <  6:1585662/149‑1 (MQ=255)

                                                                                                                                                                                                                                      ttCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCt                                                                                                                                           <  7:1551999/149‑1 (MQ=255)

                                                                                                                                                                                                                                          gACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATc                                                                                                                                       >  7:1259464/1‑149 (MQ=255)

                                                                                                                                                                                                                                          gACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATc                                                                                                                                       <  8:911632/149‑1 (MQ=255)

                                                                                                                                                                                                                                          gACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATc                                                                                                                                       >  5:203235/1‑149 (MQ=255)

                                                                                                                                                                                                                                           aCAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTg                                                                                                                                                                     <  5:559385/118‑1 (MQ=255)

                                                                                                                                                                                                                                           aCAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTg                                                                                                                                                                     >  6:559385/1‑118 (MQ=255)

                                                                                                                                                                                                                                           aCAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTAt                                                                                                                                                                <  7:437577/123‑1 (MQ=255)

                                                                                                                                                                                                                                           aCAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATcc                                                                                                                                      >  5:848492/1‑149 (MQ=255)

                                                                                                                                                                                                                                             aGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctg                                                                                                                                    >  7:220598/1‑149 (MQ=255)

                                                                                                                                                                                                                                             aGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctg                                                                                                                                    >  6:121273/1‑149 (MQ=255)

                                                                                                                                                                                                                                               aaCCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCcaca                                                                                                                                                                                                              >  2:574585/1‑73 (MQ=255)

                                                                                                                                                                                                                                                  cATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGAt                                                                                                                               >  5:992865/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   aTCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCg                                                                                                                                                                                                          <  4:574581/73‑1 (MQ=255)

                                                                                                                                                                                                                                                      gCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCa                                                                                                                           <  5:764602/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      gCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCa                                                                                                                           >  5:876897/1‑149 (MQ=255)

                                                                                                                                                                                                                                                         tctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGa                                                                                                                                                                    <  1:43440/105‑1 (MQ=255)

                                                                                                                                                                                                                                                         tctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGa                                                                                                                                                                    >  3:43440/1‑105 (MQ=255)

                                                                                                                                                                                                                                                          ctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGAt                                                                                                                       <  8:1560369/149‑1 (MQ=255)

                                                                                                                                                                                                                                                           tctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTgcgc                                                                    >  2:768884/1‑199 (MQ=255)

                                                                                                                                                                                                                                                           tctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTgtgt                                             >  4:436288/1‑222 (MQ=255)

                                                                                                                                                                                                                                                           gtccgggTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGt                                                                                                                                                                      >  2:707724/6‑101 (MQ=255)

                                                                                                                                                                                                                                                              gtgGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGc                                                                                                                   <  8:845662/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               gggTCCCNNNCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGc                                                                <  4:768878/198‑1 (MQ=255)

                                                                                                                                                                                                                                                               gggTCCCGGTCATCATGACCGGGACCCGGCCCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACt                                                                                                                                                                  <  4:707719/100‑1 (MQ=255)

                                                                                                                                                                                                                                                                 gTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGa                                                                                                                                                                    >  5:1633366/1‑97 (MQ=255)

                                                                                                                                                                                                                                                                   cccGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATTTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttc                                                                                                              <  8:112510/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                   cccGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttc                                                                                                              >  6:1675558/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    ccGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGt                                                                                                                                                              <  5:712226/100‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ccGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct                                                                                                             <  6:969163/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ccGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct                                                                                                             >  5:1042677/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttctt                                                                                                            >  5:1337988/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      ggTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtc                                                                                                                                                                                >  8:340820/1‑80 (MQ=255)

                                                                                                                                                                                                                                                                      ggTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtc                                                                                                                                                                                <  7:340820/80‑1 (MQ=255)

                                                                                                                                                                                                                                                                      ggTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttcttc                                                                                                           <  7:1384292/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      ggTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttcttc                                                                                                           <  5:449930/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                       gTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCg                                                                                                          <  8:1259464/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                           tcatGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGt                                                                                                                                                              <  8:604117/93‑1 (MQ=255)

TGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACA‑GAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCA  >  NC_000913/2074424‑2074938

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