BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,113,548	T→C	5.1%	intergenic (+142/‑31)	opgH → / → yceK	OPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist/outer membrane integrity lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,113,548	0	T	C	5.1%	708.8 / 18.8	216	intergenic (+142/‑31)	opgH/yceK	OPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist/outer membrane integrity lipoprotein
Reads supporting (aligned to +/- strand): ref base T (111/93); new base C (5/6); total (116/100)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.66e-02

Only 100 of 217 total aligned reads displayed.

GGCCCGTCTGCATTTCCGTGTCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATT  >  NC_000913/1113286‑1113808

                    tCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCACTTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGAATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAAt                                                                                                                                                                                                                                       >  2:291995/1‑275 (MQ=255)

                                   gtggtttttcttaATGGGTGAGTTATTACGAAGGGATAAAGCTTAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAACCGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtg‑gt                                                                                                                                                                                                                       <  4:291995/268‑1 (MQ=255)

                                                       aGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATAGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTAACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCACATGAGTTAAATTGTG‑TTGACCATCATGGTGACCTTGc                                                                                                                                                                                                   >  2:55405/1‑275 (MQ=255)

                                                                   gggATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGAGAAGGCGTTCACGACGCATCCGGCAAAAACAACGAGCCAATAAAAATACAGGCGTAATGCGGGTATTTTTTTACGCAAGAGGATTCCAATGTGTTTAATTGTG‑TTTAGCATGATGGTGACCTTCATGAGTGGAtgtg                                                                                                                                                                                       >  2:228100/1‑275 (MQ=38)

                                                                      aTAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTAACaaa                                                                                                                                                                                                                                        >  4:330848/1‑224 (MQ=255)

                                                                               aaTCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGc                                                                                                                                                                                                                                                                     >  2:345518/1‑186 (MQ=255)

                                                                                   cccTGGCATTGCGTAAACCGTATTGGGCTTCTCAATAAAAACTTAGTTGCCGGATGCGTTACGCTTATAAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGt                                                                                                                                                                         <  4:43170/271‑1 (MQ=255)

                                                                                   cACTGGCATTGCGTAAACCGGATGCGGCTTCGANATAANAACGTAGTTGCCTGATGCGCTACNCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGt                                                                                                                                                                                                                                                                 <  4:345518/186‑1 (MQ=255)

                                                                                                              cTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCa                                                                                                                                             <  2:169410/274‑1 (MQ=255)

                                                                                                                    aaTAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGc                                                                                                                                                                                                                                                                     >  7:865455/1‑149 (MQ=255)

                                                                                                                     aTAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGcc                                                                                                                                                                                                                                                                    >  6:575096/1‑149 (MQ=255)

                                                                                                                      tAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCATTTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGCTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCg                                                                                                                                                                                                                                                                   >  5:519465/1‑149 (MQ=255)

                                                                                                                      tAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCg                                                                                                                                                                                                                                                                   >  5:276793/1‑149 (MQ=255)

                                                                                                                         aaaCGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTa                                                                                                                                                                                                                                                                <  7:840141/149‑1 (MQ=255)

                                                                                                                          aaCGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAt                                                                                                                                                                                                                                                               <  6:171236/149‑1 (MQ=255)

                                                                                                                          aaCGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAt                                                                                                                                                                                                                                                               <  8:453289/149‑1 (MQ=255)

                                                                                                                          aaCGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAt                                                                                                                                                                                                                                                               >  6:695315/1‑149 (MQ=255)

                                                                                                                           aCGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAtt                                                                                                                                                                                                                                                              >  6:1328532/1‑149 (MQ=255)

                                                                                                                            cGTAGTTGCCTTATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAACAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAttt                                                                                                                                                                                                                                                             <  7:952089/149‑1 (MQ=25)

                                                                                                                             gTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAtttt                                                                                                                                                                                                                                                            >  7:343338/1‑149 (MQ=255)

                                                                                                                              tAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAttttt                                                                                                                                                                                                                                                           <  5:575096/149‑1 (MQ=255)

                                                                                                                                 ttGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACACCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTa                                                                                                                                                                                                                                                        >  5:1054499/1‑149 (MQ=255)

                                                                                                                                   gCCTGATTCGCTACGCTTATCAGCCCTAAATCGTTCCTGCAATTTATTGATTTTGCAATACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTGTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCa                                                                                                                       <  1:65148/275‑1 (MQ=255)

                                                                                                                                   gCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAttt                                                                                                                                                                                                                                                             >  7:70193/1‑142 (MQ=255)

                                                                                                                                   gCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAttt                                                                                                                                                                                                                                                             <  8:70193/142‑1 (MQ=255)

                                                                                                                                    ccTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGa                                                                                                                                                                                                                                                     >  8:248170/1‑149 (MQ=255)

                                                                                                                                         tGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAagag                                                                                                                                                                                                                                                <  8:804111/149‑1 (MQ=255)

                                                                                                                                             cTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAt                                                                                                                                                                                                                                                               >  5:1364985/1‑130 (MQ=255)

                                                                                                                                             cTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAt                                                                                                                                                                                                                                                               <  6:1364985/130‑1 (MQ=255)

                                                                                                                                              tACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGgtgt                                                                                                                          >  2:281542/1‑261 (MQ=255)

                                                                                                                                              tACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGgtgt                                                                                                                          <  4:281542/261‑1 (MQ=255)

                                                                                                                                                 gCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCaaa                                                                                                                                                                                                                                        >  7:151401/1‑149 (MQ=255)

                                                                                                                                                 gCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTAGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCaaa                                                                                                                                                                                                                                        <  7:649139/149‑1 (MQ=255)

                                                                                                                                                  cTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAAt                                                                                                                                                                                                                                       >  6:264625/1‑149 (MQ=255)

                                                                                                                                                   ttATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATg                                                                                                                                                                                                                                      >  7:1080199/1‑149 (MQ=255)

                                                                                                                                                    tATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGcc                                                                                                                                                                                                                                                                    >  7:214379/1‑118 (MQ=255)

                                                                                                                                                    tATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGc                                                                                                                                                                                                                                     <  6:1477482/149‑1 (MQ=255)

                                                                                                                                                     aTCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCg                                                                                                                                                                                                                                    <  8:769673/149‑1 (MQ=255)

                                                                                                                                                     aTCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTCTAATTgtg‑gt                                                                                                                                                                                                                       >  2:298272/1‑161 (MQ=255)

                                                                                                                                                      tCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGt                                                                                                                                                                                                                                   >  5:150669/1‑149 (MQ=255)

                                                                                                                                                      tCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGt                                                                                                                                                                                                                                   >  6:995627/1‑149 (MQ=255)

                                                                                                                                                      tCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGt                                                                                                                                                                                                                                   >  8:1355434/1‑149 (MQ=255)

                                                                                                                                                       cAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGtt                                                                                                                                                                                                                                  <  7:614815/149‑1 (MQ=255)

                                                                                                                                                         ggCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTa                                                                                                                                                                                                                                >  8:1398023/1‑149 (MQ=255)

                                                                                                                                                           ccTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAAt                                                                                                                                                                                                                              <  5:844893/149‑1 (MQ=255)

                                                                                                                                                            cTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAGTGCGTTTAAtt                                                                                                                                                                                                                             <  5:1096900/149‑1 (MQ=255)

                                                                                                                                                            cTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAAtt                                                                                                                                                                                                                             >  6:190235/1‑149 (MQ=255)

                                                                                                                                                               cATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgtg                                                                                                                                                                                                                          <  5:92606/149‑1 (MQ=255)

                                                                                                                                                               cATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgtg                                                                                                                                                                                                                          <  7:493535/149‑1 (MQ=255)

                                                                                                                                                               cATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAttt                                                                                                                                                                                                                                                             <  6:1083652/114‑1 (MQ=37)

                                                                                                                                                                    ttCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGt                                                                                                                                                                         >  1:52778/1‑192 (MQ=255)

                                                                                                                                                                     tCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGc                                                                                                                                                                                                                   <  5:823667/149‑1 (MQ=255)

                                                                                                                                                                        tGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGcatc                                                                                                                                                                                                                >  7:1211472/1‑149 (MQ=255)

                                                                                                                                                                         gCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTc                                                                                                                                                                        >  1:31547/1‑188 (MQ=255)

                                                                                                                                                                           aaTTTATTGATTTTACAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATg                                                                                                                                                                                                             <  6:45854/149‑1 (MQ=255)

                                                                                                                                                                            aTTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAaga                                                                                                                                                                                                                                                 <  8:835418/113‑1 (MQ=255)

                                                                                                                                                                             tttATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCAc                                                                                                                                                                    <  3:31547/188‑1 (MQ=255)

                                                                                                                                                                                aTTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGAc                                                                                                                                                                                                        >  5:1193679/1‑149 (MQ=255)

                                                                                                                                                                                aTTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGAAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGAc                                                                                                                                                                                                        <  5:1018899/149‑1 (MQ=255)

                                                                                                                                                                                 ttGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGAcc                                                                                                                                                                                                       >  5:426862/1‑149 (MQ=255)

                                                                                                                                                                                 ttGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATGGTG‑GTGAGCATCATGGTGAcc                                                                                                                                                                                                       >  6:517010/1‑149 (MQ=255)

                                                                                                                                                                                  tGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCt                                                                                                                                                                                                      <  5:410673/149‑1 (MQ=255)

                                                                                                                                                                                  tGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCt                                                                                                                                                                                                      <  8:361082/149‑1 (MQ=255)

                                                                                                                                                                                    aTTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTg                                                                                                                                                                                                    >  5:335339/1‑149 (MQ=255)

                                                                                                                                                                                     tgttGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGc                                                                                                                                                                                                   <  8:1253042/147‑1 (MQ=255)

                                                                                                                                                                                            aGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGa                                                                                                                                                                                                <  8:877971/145‑1 (MQ=255)

                                                                                                                                                                                            aGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGa                                                                                                                                                                                                >  7:877971/1‑145 (MQ=255)

                                                                                                                                                                                                tttGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAagag                                                                                                                                                                                                                                                >  6:86378/1‑94 (MQ=255)

                                                                                                                                                                                                  tgtgGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATg                                                                                       >  2:45289/1‑244 (MQ=255)

                                                                                                                                                                                                       gTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgcagca                                                                                                                                                                                 >  6:1144382/1‑149 (MQ=255)

                                                                                                                                                                                                           gATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCattat                                                                                                                                                                             >  7:428086/1‑149 (MQ=255)

                                                                                                                                                                                                             tAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATTATGCCGCTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTTACGTTGCTGAGTGGCTGTGGCAGCAttatta                                                                                                                                                                           <  8:1211472/149‑1 (MQ=255)

                                                                                                                                                                                                              aaGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAg                                                                                                                                                                          <  8:885723/149‑1 (MQ=255)

                                                                                                                                                                                                                        cgccgcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACtata                                                                                                                                                                <  5:305053/149‑1 (MQ=255)

                                                                                                                                                                                                                           cgcATCCGGTAAAAACAACGGGCCAATAAAAATACCGGCGTTAGGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GGGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGAt                                <  4:261886/274‑1 (MQ=255)

                                                                                                                                                                                                                             cATCCGGCAGAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggg                                                                                                                                                           <  7:1230161/149‑1 (MQ=255)

                                                                                                                                                                                                                               tCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgt                                                                                                                                                                                                                           <  7:196432/84‑1 (MQ=255)

                                                                                                                                                                                                                               tCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgt                                                                                                                                                                                                                           >  8:196432/1‑84 (MQ=255)

                                                                                                                                                                                                                                    cAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggggcagggg                                                                                                                                                    <  6:303585/149‑1 (MQ=255)

                                                                                                                                                                                                                                      aaaaaaACTGGCCAATAAAAAATACGTGGTTATGCGGGTAGTTGTTGACGAAAGGGGTATCAAATGCGTTTAATGGTG‑GTGGGCATCAGTGGGACCTTGCTGAGTGGCTGTGGCAGGATGATTAGGCGCACTATACCGGGGCAGGGTCATGGCAACCAATATTATCCCGGTGTGCAATTGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATtcatca                    <  4:140487/275‑1 (MQ=255)

                                                                                                                                                                                                                                      aaaaCAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggggcaggggca                                                                                                                                                  >  5:633239/1‑149 (MQ=255)

                                                                                                                                                                                                                                       aaaCAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCAt                                                                                                                                                 <  7:1336516/149‑1 (MQ=255)

                                                                                                                                                                                                                                        aacaacGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATg                                                                                                                                                >  7:1367856/1‑149 (MQ=255)

                                                                                                                                                                                                                                           aacGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCa                                                                                                                                             >  8:1132869/1‑149 (MQ=255)

                                                                                                                                                                                                                                              gAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAAcc                                                                                                                                          >  7:942388/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   aaTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggggcag                                                                                                                                                       >  6:1265379/1‑131 (MQ=255)

                                                                                                                                                                                                                                                   aaTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggggcag                                                                                                                                                       <  5:1265379/131‑1 (MQ=255)

                                                                                                                                                                                                                                                   aaTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATg                                                                                                                                                >  8:60324/1‑138 (MQ=255)

                                                                                                                                                                                                                                                    aTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAAtatt                                                                                                                                    >  7:1237059/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAAtatt                                                                                                                                    >  6:546827/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     tAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAAtatta                                                                                                                                   >  5:679634/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     tAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAAtatta                                                                                                                                   >  6:1155966/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        aaaTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGt                                                                                                                                                                                              <  8:531363/87‑1 (MQ=38)

                                                                                                                                                                                                                                                        aaaTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGt                                                                                                                                                                                              >  7:531363/1‑87 (MQ=38)

                                                                                                                                                                                                                                                        aaaTACCCGCGTTATGCCGGTTTGTTTTTTCGAAAGAGGTATCAAATTCGTTTAATGGTGTGTG‑GCATCATGTGTACCTTGGTTGTGGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGAATGGCAACCAATATTATCCCGGTGTGAAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATGCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGAtt  <  4:55405/275‑1 (MQ=255)

                                                                                                                                                                                                                                                          aTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCt                                                                                                                                                                                                  >  5:1411836/1‑81 (MQ=255)

                                                                                                                                                                                                                                                            aCCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGgt                                                                                                                            <  7:1398023/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                    tATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATg                                                                                                                    >  5:188005/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     aTGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATgg                                                                                                                   <  7:229666/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      tGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATggg                                                                                                                  <  5:1144382/149‑1 (MQ=255)

GGCCCGTCTGCATTTCCGTGTCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTG‑GTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATT  >  NC_000913/1113286‑1113808

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