BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	858,924	T→G	5.4%	R46R (CGA→CGC)	ybiU ←	DUF1479 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	858,924	0	T	G	5.4%	454.4 / 13.9	129	R46R (CGA→CGC)	ybiU	DUF1479 family protein
Reads supporting (aligned to +/- strand): ref base T (60/62); new base G (5/2); total (65/64)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 130 total aligned reads displayed.

GGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTC  >  NC_000913/858777‑859070

ggggAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACt                                                                                                                                                    >  4:960812/1‑148 (MQ=255)

          ccTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCaa                                                                                                                                          <  4:759112/148‑1 (MQ=255)

          ccTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCaa                                                                                                                                          <  3:351878/148‑1 (MQ=255)

            ttttATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATg                                                                                                                                        >  3:565130/1‑148 (MQ=255)

            ttttATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGt                                                                                                                                       >  1:647150/1‑149 (MQ=255)

                aTCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCAt                                                                                                                                    <  3:808186/148‑1 (MQ=255)

                 tCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATc                                                                                                                                   >  4:834990/1‑148 (MQ=255)

                 tCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCg                                                                                                                                  >  1:1012330/1‑149 (MQ=255)

                  cACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGc                                                                                                                                 >  1:18224/1‑149 (MQ=255)

                    ccGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATg                                                                                                                                        >  3:207663/1‑140 (MQ=255)

                    ccGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATg                                                                                                                                        <  4:207663/140‑1 (MQ=255)

                     cGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTa                                                                                                                              <  1:338142/149‑1 (MQ=255)

                      gCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTa                                                                                                                              >  3:388616/1‑148 (MQ=255)

                               cgACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTaaa                                                                                                                    <  2:595871/149‑1 (MQ=255)

                               cgACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTaa                                                                                                                     <  3:361337/148‑1 (MQ=255)

                               cgACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTaa                                                                                                                     <  4:334261/148‑1 (MQ=255)

                               cgACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCCTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTaa                                                                                                                     <  4:696207/148‑1 (MQ=255)

                                   gTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGAt                                                                                                                 <  3:717448/148‑1 (MQ=255)

                                   gTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATc                                                                                                                <  2:846401/149‑1 (MQ=255)

                                    tttAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATct                                                                                                               <  1:617227/149‑1 (MQ=255)

                                      tAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctgc                                                                                                             <  2:824992/149‑1 (MQ=255)

                                      tAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctg                                                                                                              <  3:130732/148‑1 (MQ=255)

                                        aTCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctgc                                                                                                             <  1:799984/147‑1 (MQ=255)

                                          cTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTgg                                                                                                          >  4:837058/1‑148 (MQ=255)

                                           tGTTCGCGCTGGTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGAc                                                                                                        >  1:798277/1‑149 (MQ=255)

                                            gTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGAc                                                                                                        >  3:932557/1‑148 (MQ=255)

                                             ttCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGt                                                                                                      >  2:135613/1‑149 (MQ=255)

                                                       tctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTg                                                                                            >  2:1010277/1‑149 (MQ=255)

                                                        ctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTgc                                                                                           <  2:1012330/149‑1 (MQ=255)

                                                        ctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTgc                                                                                           <  1:847937/149‑1 (MQ=255)

                                                              tAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGc                                                                                     >  1:135499/1‑149 (MQ=255)

                                                                     atcTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGCCGTCGCCAAGCTGCGCCCGCAGCGCGt                                                                              <  2:458760/147‑1 (MQ=255)

                                                                     aCCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGcgcg                                                                               >  3:1037701/1‑148 (MQ=255)

                                                                     aCCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGt                                                                              <  1:784323/149‑1 (MQ=255)

                                                                       cTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGc                                                                            >  1:474042/1‑149 (MQ=255)

                                                                           tttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTc                                                                         >  4:366801/1‑148 (MQ=255)

                                                                           tttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCa                                                                        >  2:992415/1‑149 (MQ=255)

                                                                            ttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCa                                                                        >  4:293432/1‑148 (MQ=255)

                                                                            ttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCa                                                                        >  4:575351/1‑148 (MQ=255)

                                                                            ttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCa                                                                        >  4:44386/1‑148 (MQ=255)

                                                                             tGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCAt                                                                       <  4:617928/148‑1 (MQ=255)

                                                                             tGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATc                                                                      <  2:881584/149‑1 (MQ=255)

                                                                             tGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATc                                                                      <  2:1046117/149‑1 (MQ=255)

                                                                             tGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATc                                                                      <  2:798277/149‑1 (MQ=255)

                                                                             tGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATc                                                                      <  2:133772/149‑1 (MQ=255)

                                                                                  tCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACg                                                                 >  2:546041/1‑149 (MQ=255)

                                                                                    ggCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGa                                                               <  2:18224/149‑1 (MQ=255)

                                                                                     gCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGAt                                                              >  1:134063/1‑149 (MQ=255)

                                                                                      cATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGAt                                                              <  3:716589/148‑1 (MQ=255)

                                                                                      cATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGAt                                                              <  3:575351/148‑1 (MQ=255)

                                                                                      aaTAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATa                                                             >  2:912091/2‑149 (MQ=255)

                                                                                         aaGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCCGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgct                                                          <  2:836053/149‑1 (MQ=255)

                                                                                           gACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCAGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgctgc                                                        >  2:497758/1‑149 (MQ=255)

                                                                                           gACAGCACCGGCCAGACGGCATCGCCCGGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgctgc                                                        <  2:462061/149‑1 (MQ=255)

                                                                                              aGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGa                                                               >  1:418877/1‑139 (MQ=255)

                                                                                              aGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGa                                                               <  2:418877/139‑1 (MQ=255)

                                                                                              aGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCtt                                                      <  4:360001/148‑1 (MQ=255)

                                                                                              aGCACCGGACAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTTCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCttt                                                     <  2:878993/149‑1 (MQ=255)

                                                                                                cACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTtgt                                                   <  1:1046015/149‑1 (MQ=255)

                                                                                                 aCCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTtgtg                                                  >  1:869877/1‑149 (MQ=255)

                                                                                                  ccGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGa                                                 <  1:927099/149‑1 (MQ=255)

                                                                                                     gCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCg                                              <  2:741153/149‑1 (MQ=255)

                                                                                                        aGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCggca                                           <  1:270991/149‑1 (MQ=255)

                                                                                                        aGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCggca                                           <  1:122771/149‑1 (MQ=255)

                                                                                                         gACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCggca                                           >  4:1025429/1‑148 (MQ=255)

                                                                                                           cGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCggcagg                                         <  4:388616/148‑1 (MQ=255)

                                                                                                             gCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGTCaa                                      <  1:926227/149‑1 (MQ=255)

                                                                                                              cATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCaa                                      <  3:834990/148‑1 (MQ=255)

                                                                                                                tCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAAtg                                    >  4:152236/1‑148 (MQ=255)

                                                                                                                tCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCAGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAAtgt                                   <  1:497758/149‑1 (MQ=255)

                                                                                                                    ccTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCg                               >  1:892467/1‑149 (MQ=255)

                                                                                                                      tgtgCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCt                             >  2:384598/1‑149 (MQ=255)

                                                                                                                      tgtgCTTTGAGTGCGTTGATTTCAGCCACTCGCGTAGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGc                              >  4:1000250/1‑148 (MQ=255)

                                                                                                                       gtgCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCt                             >  4:487291/1‑148 (MQ=255)

                                                                                                                           tttGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTaa                        >  1:61577/1‑149 (MQ=255)

                                                                                                                                gTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTaaaa                      >  1:1028145/1‑146 (MQ=255)

                                                                                                                                gTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTaaaa                      <  2:1028145/146‑1 (MQ=255)

                                                                                                                                   cGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGaa                <  2:134063/149‑1 (MQ=255)

                                                                                                                                    gTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAg               >  1:848321/1‑149 (MQ=255)

                                                                                                                                    gTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGGGTCGCTGGTAAAAGTAGaa                >  3:839960/1‑148 (MQ=255)

                                                                                                                                     ttGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGa                                                                   <  4:181834/96‑1 (MQ=255)

                                                                                                                                     ttGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGa                                                                   >  3:181834/1‑96 (MQ=255)

                                                                                                                                     ttGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAg               >  3:347509/1‑148 (MQ=255)

                                                                                                                                     ttGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGc              >  1:872983/1‑149 (MQ=255)

                                                                                                                                      tGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGAt                                                <  3:58677/114‑1 (MQ=255)

                                                                                                                                      tGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGAt                                                >  4:58677/1‑114 (MQ=255)

                                                                                                                                      tGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGc              >  4:909526/1‑148 (MQ=255)

                                                                                                                                       gATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCa            <  2:61577/149‑1 (MQ=255)

                                                                                                                                           tCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac        >  2:63303/1‑149 (MQ=255)

                                                                                                                                           attgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac        >  1:284506/4‑149 (MQ=255)

                                                                                                                                            ttgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTc                                                                         >  4:253415/3‑83 (MQ=255)

                                                                                                                                            ttgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac        >  4:255919/3‑148 (MQ=255)

                                                                                                                                            cAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac        >  3:227546/1‑148 (MQ=255)

                                                                                                                                             tgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGAt                                                <  3:963927/106‑1 (MQ=255)

                                                                                                                                             tgCCACGCGAGTGGCAAAGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGAt                                                >  4:963927/2‑107 (MQ=255)

                                                                                                                                             aGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAacca      >  2:730106/1‑149 (MQ=255)

                                                                                                                                             aGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAacca      >  1:1061353/1‑149 (MQ=255)

                                                                                                                                             aGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAacc       >  3:148724/1‑148 (MQ=255)

                                                                                                                                                  cTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACctc  >  4:231123/1‑148 (MQ=255)

                                                                                                                                                  cTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACctc  >  3:661225/1‑148 (MQ=255)

GGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTC  >  NC_000913/858777‑859070

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