BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,117,004	T→A	D73V (GAT→GTT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,004	0	T	A	98.0%	292.5 / NA	99	D73V (GAT→GTT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base T (0/1); major base A (46/50); minor base C (0/2); total (46/53)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGG  >  NC_000913/4116859‑4117145

tAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACgg                                                                                                                                            <  2:1199263/149‑1 (MQ=255)

tAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAAAAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACgg                                                                                                                                            >  2:1210859/1‑149 (MQ=255)

  aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCTGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGaa                                                                                                                                          <  1:617915/149‑1 (MQ=255)

  aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGaa                                                                                                                                          <  2:526702/149‑1 (MQ=255)

      gTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGa                                                                                                                                      >  2:1182268/1‑149 (MQ=255)

           gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTGCTTTTTCCCAGACAATAGTGGTTTTTCAGTTGGGCGTAATCCCGTTAGCTGCAATTTGAACGGAACTGATTTCg                                                                                                                                 >  1:509406/1‑149 (MQ=255)

           gTTGTAAATGCTCGCAGAGTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCg                                                                                                                                 <  1:113926/149‑1 (MQ=255)

            ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGAATCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAATTGAACGGAACTGAGATCgg                                                                                                                                >  1:663166/1‑149 (MQ=255)

              ttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCt                                                                                                                              >  2:970379/1‑149 (MQ=255)

               aaaaaacTCGCAGATTCCTCCGGTACGACGGCACCGTCAGACAATGGCGTTATAGATAGGCATGCCGGTTTCTTTTTCGCAGACAATAGTGGTTTCACGCTGGTACGCAATACCGACAACTCTACTGTGAACGGAACTGATATCGGCtt                                                                                                                             >  2:322229/7‑149 (MQ=255)

                aaaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCttt                                                                                                                            >  2:991137/1‑149 (MQ=255)

                aaaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCttt                                                                                                                            >  2:991070/1‑149 (MQ=255)

                   tgtccggaGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGATAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGc                                                                                                                         <  2:30541/142‑1 (MQ=255)

                   tGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGc                                                                                                                         <  1:711422/149‑1 (MQ=255)

                   cgCTCGCAGGTTTCTGCGGTACGACGGCTCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACGATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGc                                                                                                                         <  2:812450/148‑1 (MQ=255)

                       cGCCGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGc                                                                                                                     <  2:1285764/149‑1 (MQ=255)

                       cGCAGATTTCTGCGGTACGACGGCACTGCCAGACACTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGc                                                                                                                     <  2:1285531/149‑1 (MQ=255)

                        gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCa                                                                                                                    <  1:670224/149‑1 (MQ=255)

                        gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCa                                                                                                                    <  1:931122/149‑1 (MQ=255)

                        gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCa                                                                                                                    >  2:502844/1‑149 (MQ=255)

                        gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATACCGGCTTTCGCCAGCa                                                                                                                    <  2:862225/149‑1 (MQ=255)

                           gATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACtt                                                                                                                 <  2:1081727/149‑1 (MQ=255)

                            aTTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTc                                                                                                                <  2:1320839/149‑1 (MQ=255)

                             tttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCt                                                                                                               <  2:825653/149‑1 (MQ=255)

                               tCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATCCCGCCAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTAc                                                                                                             <  1:1215866/149‑1 (MQ=255)

                                 tGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCa                                                                                                           >  2:282516/1‑149 (MQ=255)

                                  gCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCCCTTCTACAAg                                                                                                          >  2:656310/1‑149 (MQ=255)

                                      tACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTg                                                                                                      <  2:290937/149‑1 (MQ=255)

                                         gacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAg                                                                                                   <  1:723772/149‑1 (MQ=255)

                                                cTGCAAGACAATGGCGGTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGCAAGCCCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTggg                                                                                            <  1:577185/149‑1 (MQ=255)

                                                 gtccggACAACGGCGTTCTCGCTAGCCCTGTCGTGTTCTTTTTCCCAAACAATAGTGGTTTCAAGCTGGTTCGTAATTCCGATAGTTGCAAGTCCAACGGAACTGATATCGGCTTTCGCCATCACTTCTACCAGCGTAGAGCTTGGGGt                                                                                           <  2:324511/144‑1 (MQ=255)

                                                   ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGGGCTTTGGGTgg                                                                                         >  1:221027/1‑149 (MQ=255)

                                                   ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCGGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGGGCTTTgggggg                                                                                         >  1:234301/1‑149 (MQ=255)

                                                     aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAAGTGATATCGGCTTTCGCCAGCACTTCTACCAGCGGGGAGCTTTGGGTgggc                                                                                       >  2:1000508/1‑147 (MQ=255)

                                                     aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGcc                                                                                       >  2:451002/1‑149 (MQ=255)

                                                     aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGcc                                                                                       <  1:451002/149‑1 (MQ=255)

                                                            ggCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGGACGGAACTGATATTGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGGttt                                                                                >  2:589777/1‑149 (MQ=255)

                                                             gCGTTATAGATAGGCTTGCCGGTGTCTTCTTCCCCGACAATAGTGGTTTCACGGTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCCGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTc                                                                               >  1:855606/1‑149 (MQ=255)

                                                              cGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                                                                              <  1:1404871/149‑1 (MQ=255)

                                                                ttATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAt                                                                            >  2:1331904/1‑149 (MQ=255)

                                                                   tagatagGCTTGCCGGTTTCTTTTTCCCCGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg                                                                         <  1:381517/149‑1 (MQ=255)

                                                                      atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                      <  1:1529669/149‑1 (MQ=255)

                                                                      atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                      <  1:1275229/149‑1 (MQ=255)

                                                                      atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                      <  2:110948/149‑1 (MQ=255)

                                                                       ttggCTTGCCGGTTTCTTTTTCCCAGTCAGTAGTGGTTTCACGCTGGTTCGTAATACCGATGGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCg                                                                     <  2:147798/147‑1 (MQ=255)

                                                                        agGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGCTATCGGCTTTCGCCAGCCCTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATGTCCGTTGGGGCgt                                                                    >  2:1401557/1‑149 (MQ=255)

                                                                        agGCTTGCCGGTTTCGTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACAAGTTTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                                                    <  2:718904/149‑1 (MQ=255)

                                                                        agGCTTCCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                                                    <  2:846918/149‑1 (MQ=255)

                                                                          gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCCCTTCTACCAGCGTGGGGCTTTGGGTGGCCCCGATTTCCATTGGGCCgtgt                                                                  >  2:655421/1‑149 (MQ=255)

                                                                          gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                                                  >  2:1094533/1‑149 (MQ=255)

                                                                             tGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGATCTCGCCAGCAATTCTACCAGTTGGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCt                                                               <  1:1140849/149‑1 (MQ=255)

                                                                              gCCGGTTTCTCTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTCGAACGGAACTGATATCCGCTGTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                                              <  1:454499/149‑1 (MQ=255)

                                                                               ccGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACGTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCtca                                                             >  2:865087/1‑147 (MQ=255)

                                                                                       tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACGGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                     >  1:1384602/1‑149 (MQ=255)

                                                                                       tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCACCCTg                                                     >  1:1343757/1‑149 (MQ=255)

                                                                                            cccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                >  2:554931/1‑149 (MQ=255)

                                                                                            cccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGCCGTGTTCTACCCAACCTGGtttt                                                >  2:102875/1‑149 (MQ=255)

                                                                                             ccAGAGAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCCTTGGGTCGTGGGCTACCCAACCTGGGTTTg                                               >  2:1089385/1‑149 (MQ=255)

                                                                                              cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCGGCGCTTCTACCAGCGTGGCGCTTTGGGTGGCCCAGATTTCCAATGGGGCGTGTTCTACCCACCCTGGTGTTgg                                              >  1:293855/1‑149 (MQ=255)

                                                                                              cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTGCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                              >  1:191546/1‑149 (MQ=255)

                                                                                              cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                              >  1:908124/1‑149 (MQ=255)

                                                                                              cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                              <  1:502844/149‑1 (MQ=255)

                                                                                              cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                              >  1:293650/1‑149 (MQ=255)

                                                                                               aGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTggg                                             >  1:719365/1‑149 (MQ=255)

                                                                                                gACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACACAACCTGGTTTTGGGt                                            >  2:629364/1‑149 (MQ=255)

                                                                                                 cacatggTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTa                                           <  2:509406/143‑1 (MQ=255)

                                                                                                     tAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAtt                                       >  2:1352320/1‑149 (MQ=255)

                                                                                                       gTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAACACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTACCCTTGGGTCGTGTTCTACCCAACCTGGTTTTGGGttagtttg                                     >  2:354500/1‑142 (MQ=255)

                                                                                                        tGGTGTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc                                    <  2:927212/149‑1 (MQ=255)

                                                                                                        ggtcttCACTGTGGTCCGTCATCCCGATCGCTCCAATTTGAACGGACCTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTCTTCTACCCAACCTGGTTTTGGGTAGATTTGc                                    <  1:877201/145‑1 (MQ=255)

                                                                                                          gTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTc                                  >  2:1436473/1‑149 (MQ=255)

                                                                                                          gTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTc                                  >  2:706709/1‑149 (MQ=255)

                                                                                                          gTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTGGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTc                                  <  1:676871/149‑1 (MQ=255)

                                                                                                                  cTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgc                          <  1:991070/149‑1 (MQ=255)

                                                                                                                  cTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgc                          <  1:991137/149‑1 (MQ=255)

                                                                                                                    ggTTCGTAGTACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCgcgc                        <  1:282516/149‑1 (MQ=255)

                                                                                                                     ggtCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCt                       <  1:654624/147‑1 (MQ=255)

                                                                                                                      ttCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGCTTTGGGTAGAGTTGCTCAAATTCGCGCTg                      >  1:691492/1‑149 (MQ=255)

                                                                                                                       tCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGc                     <  1:1182268/149‑1 (MQ=255)

                                                                                                                         ttAATACCGATAGCCGCAATTTGAACGGAGCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGa                   <  2:372495/148‑1 (MQ=255)

                                                                                                                          tAATACCGATAGCTGCAATTTGAACGGAATTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                              <  2:590297/121‑1 (MQ=255)

                                                                                                                          tAATACCGATAGCTGCAATTTGAACGGAATTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                              >  1:590297/1‑121 (MQ=255)

                                                                                                                           aaTACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGaca                 >  1:928683/1‑149 (MQ=255)

                                                                                                                             tACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACg               <  1:701789/149‑1 (MQ=255)

                                                                                                                              aCCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCCGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTGTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGc              >  2:8131/1‑149 (MQ=255)

                                                                                                                               ccGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCt             >  2:1196352/1‑149 (MQ=255)

                                                                                                                               ccGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCt             <  1:729449/149‑1 (MQ=255)

                                                                                                                                cGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTa            >  1:1196688/1‑149 (MQ=255)

                                                                                                                                cGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATAGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTa            >  1:424830/1‑149 (MQ=255)

                                                                                                                                 gATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTaa           <  1:1094533/149‑1 (MQ=255)

                                                                                                                                  taaGCGGCAATTTGAACGGAACTGGTATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAAt          <  2:440843/147‑1 (MQ=255)

                                                                                                                                     gCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGat       <  2:1143769/149‑1 (MQ=255)

                                                                                                                                      cagCAATTTGAACGGATCTGATATCGGCTTGCGCCAGCACTTCTACCAGCGTGGAGCTTAGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGata      <  1:364196/147‑1 (MQ=255)

                                                                                                                                      cTGCAATTTGACCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGata      <  2:790382/149‑1 (MQ=255)

                                                                                                                                      cTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGata      <  2:790468/149‑1 (MQ=255)

                                                                                                                                        gCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTGTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt    >  1:1170119/1‑149 (MQ=255)

                                                                                                                                         cAATTTGACCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTTGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTg   <  2:663166/149‑1 (MQ=255)

                                                                                                                                         cAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTg   <  2:570911/149‑1 (MQ=255)

                                                                                                                                          aaTTTTAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGATGCGACACGCTAATGATATTgg  >  2:658325/1‑149 (MQ=255)

TAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGG  >  NC_000913/4116859‑4117145

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