BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,060	C→T	86.2%	A55T (GCC→ACC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,060	0	C	T	86.2%	2182.9 / 276.3	814	A55T (GCC→ACC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (59/54); new base T (354/347); total (413/401)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 819 total aligned reads displayed.

AGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA  >  NC_000913/4116913‑4117207

     aTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAt                                                                                                                                               <  2:3519508/149‑1 (MQ=255)

      tGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGAtt                                                                                                                                              >  4:3324171/1‑149 (MQ=255)

        gCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTc                                                                                                                                            <  3:1418593/149‑1 (MQ=255)

           ttATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCAt                                                                                                                                         >  3:3774078/1‑149 (MQ=255)

             atagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTg                                                                                                                                       <  1:4014561/149‑1 (MQ=255)

                 atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTc                                                                                                                                   >  3:3390922/1‑149 (MQ=255)

                 atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAATTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTc                                                                                                                                   >  3:6418971/1‑149 (MQ=255)

                   agGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCgt                                                                                                                                 >  4:2663466/1‑149 (MQ=255)

                     gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCgtgt                                                                                                                               <  4:4186038/149‑1 (MQ=255)

                      cTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGtt                                                                                                                              <  3:4736499/149‑1 (MQ=255)

                       ttGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTc                                                                                                                             >  1:5090447/1‑149 (MQ=255)

                         tctggTTTGTTTTTCCCAGACAATAGTGGTTTGACGCTGGTTCGTAATACCGATAGCTGCACTTTGATCGGAACTGATATCGGCTATCGCCAGCACTTCTACCATCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTa                                                                                                                           <  3:4449386/146‑1 (MQ=255)

                          ccGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTAc                                                                                                                          <  2:1701554/149‑1 (MQ=255)

                           cGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTAcc                                                                                                                         >  2:6547134/1‑149 (MQ=255)

                            ggTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTAccc                                                                                                                        >  2:4742755/1‑149 (MQ=255)

                             gTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCa                                                                                                                       >  4:3537052/1‑149 (MQ=255)

                             gTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTggg                                                                                                                                     >  2:1476742/1‑135 (MQ=255)

                              tttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTAcc                                                                                                                         >  3:2933621/1‑146 (MQ=255)

                              tttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTAcc                                                                                                                         <  4:2933621/146‑1 (MQ=255)

                              attCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCa                                                                                                                                          >  2:5066783/2‑129 (MQ=255)

                                 cTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCt                                                                                                                   <  1:7018980/149‑1 (MQ=255)

                                  tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                                                  >  1:4587514/1‑149 (MQ=255)

                                  tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                                                  <  2:6774575/149‑1 (MQ=255)

                                  tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                                                  >  3:5321663/1‑149 (MQ=255)

                                  tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                                                  >  1:1007795/1‑149 (MQ=255)

                                   ttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTgg                                                                                                                 >  1:7193582/1‑149 (MQ=255)

                                     ttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtt                                                                                                               <  1:3265819/149‑1 (MQ=255)

                                        ccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                                                                            >  1:1202758/1‑149 (MQ=255)

                                        ccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                                                                            <  1:6488784/149‑1 (MQ=255)

                                         cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCa                                                                                                                                          >  2:6181627/1‑118 (MQ=255)

                                         cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTgg                                                                                                                 <  1:992426/143‑1 (MQ=255)

                                         cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                           <  4:6321195/149‑1 (MQ=255)

                                         cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                           <  1:570133/149‑1 (MQ=255)

                                         cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                           >  2:5510194/1‑149 (MQ=255)

                                         cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                           <  1:1467968/149‑1 (MQ=255)

                                          aGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTggg                                                                                                          >  3:1659214/1‑149 (MQ=255)

                                              aaTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGtt                                                                                                                              <  4:525804/125‑1 (MQ=255)

                                                 aGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                           <  1:5281635/141‑1 (MQ=255)

                                                 aGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAttt                                                                                                   <  4:1149294/149‑1 (MQ=255)

                                                      tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCa                                                                                              >  3:6196829/1‑149 (MQ=255)

                                                                 ttCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTgg                                                                                                                                      >  3:6750105/1‑98 (MQ=255)

                                                                    gTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGa                                                                                >  4:714208/1‑149 (MQ=255)

                                                                     tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGac                                                                               >  3:5606100/1‑149 (MQ=255)

                                                                     tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGac                                                                               <  1:218951/149‑1 (MQ=255)

                                                                     tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGac                                                                               >  4:3907115/1‑149 (MQ=255)

                                                                     tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGa                                                                                <  2:3398725/148‑1 (MQ=255)

                                                                      aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTc                                                                                                                                   >  1:378693/1‑96 (MQ=255)

                                                                        tACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACg                                                                            >  2:5266443/1‑149 (MQ=255)

                                                                         aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGc                                                                           >  4:4558752/1‑149 (MQ=255)

                                                                           cGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgcg                                                                                      >  3:5184116/1‑136 (MQ=255)

                                                                            gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTaa                                                                        <  3:4825485/149‑1 (MQ=255)

                                                                            gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTaa                                                                        <  3:3537052/149‑1 (MQ=255)

                                                                             aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAAt                                                                       <  1:3430272/149‑1 (MQ=255)

                                                                                 cTGCGATTTGATCGGAACTGATATAGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGata                                                                   >  1:6180928/1‑149 (MQ=255)

                                                                                   gCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt                                                                 <  3:6886306/149‑1 (MQ=255)

                                                                                     aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCAGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTgg                                                               <  3:966766/149‑1 (MQ=255)

                                                                                       tttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCa                                                             >  1:2031084/1‑149 (MQ=255)

                                                                                       tttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCa                                                             >  4:1848335/1‑149 (MQ=255)

                                                                                             cGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCgt                                                                                                                                 >  1:2558879/1‑75 (MQ=255)

                                                                                             cGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGa                                                       <  3:5958207/149‑1 (MQ=255)

                                                                                              ggAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGAt                                                      <  3:7123374/149‑1 (MQ=255)

                                                                                              ggAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGAt                                                      <  4:513087/149‑1 (MQ=255)

                                                                                              ggAACTGATATCGGCTGTCGCCAGCACTTCTACCAGCGTGGAGCTTTGTGTGGTCCAGATTTCCATTGGGTCGTGTTCTACTCAACCTCGTTTTGGGGAGTTTTGCTCAAATTCGCGCTGCGACACGCTAATGCTATTGGCATCGTGAt                                                      <  2:4247386/149‑1 (MQ=255)

                                                                                                aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTg                                                                                                                                       >  3:103414/1‑66 (MQ=255)

                                                                                                aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt                                                                 >  4:5699573/1‑136 (MQ=255)

                                                                                                aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATcc                                                    >  2:2329168/1‑149 (MQ=255)

                                                                                                aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATcc                                                    >  4:6603369/1‑149 (MQ=255)

                                                                                                 aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGGATCGTGATCCa                                                   <  3:4056382/149‑1 (MQ=255)

                                                                                                   tGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGGCACGCTAATGATATTGGCATCGTGATCCAtt                                                 <  2:2172125/149‑1 (MQ=255)

                                                                                                    gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTa                                                <  4:2433216/149‑1 (MQ=255)

                                                                                                     atatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTAc                                               <  2:4459548/149‑1 (MQ=255)

                                                                                                      tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                           <  2:796038/88‑1 (MQ=255)

                                                                                                        tCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTg                                                                                                                                       <  4:1564668/58‑1 (MQ=255)

                                                                                                        tCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGc                                                                                  <  2:288000/111‑1 (MQ=255)

                                                                                                             tttCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAAt                                                                       <  3:5705096/117‑1 (MQ=255)

                                                                                                             tttCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcgcgc                                       <  2:1202307/149‑1 (MQ=255)

                                                                                                               tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCgg                                     <  1:6779544/149‑1 (MQ=255)

                                                                                                               tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCgg                                     <  2:1841002/149‑1 (MQ=255)

                                                                                                               tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCgg                                     >  2:3719181/1‑149 (MQ=255)

                                                                                                                 gCCAGCACTTCTACGAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGa                                   >  1:5454045/1‑149 (MQ=255)

                                                                                                                 gCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGa                                   <  3:2684646/149‑1 (MQ=255)

                                                                                                                  ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCACCCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCgg                                     <  4:2081513/146‑1 (MQ=255)

                                                                                                                     gCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtg                               >  1:4270417/1‑149 (MQ=255)

                                                                                                                        cctCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                            <  1:4836579/147‑1 (MQ=255)

                                                                                                                        cTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                            >  2:725497/1‑149 (MQ=255)

                                                                                                                         ttCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTg                                                        >  2:3099635/1‑120 (MQ=255)

                                                                                                                         ttCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtgg                           >  4:677034/1‑149 (MQ=255)

                                                                                                                         ttCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtgg                           >  2:4846136/1‑149 (MQ=255)

                                                                                                                          tCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTggtggt                          <  1:172783/149‑1 (MQ=255)

                                                                                                                           cTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGgtggtg                         >  1:2926035/1‑149 (MQ=255)

                                                                                                                            tACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGAGc                        >  1:2740176/1‑149 (MQ=255)

                                                                                                                              ccAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGccc                      <  3:3255609/149‑1 (MQ=255)

                                                                                                                              ccAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGccc                      <  4:5991654/149‑1 (MQ=255)

                                                                                                                                  cGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGt                  >  1:5891631/1‑149 (MQ=255)

                                                                                                                                   gtagAGCTTTGGGTGGTCTAGATTTCCATTGGGTCGTGTTCTACCCACCCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATCGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTc                 >  4:4781188/4‑149 (MQ=255)

                                                                                                                                       aGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgc             <  4:3079164/149‑1 (MQ=255)

                                                                                                                                        gCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgcg            >  1:3614534/1‑149 (MQ=255)

                                                                                                                                              ggtggtCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGa      >  3:2242216/1‑149 (MQ=255)

                                                                                                                                                tGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgc                                                                                       <  2:4006498/66‑1 (MQ=255)

                                                                                                                                                tGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCagaaatc    <  2:5891290/149‑7 (MQ=255)

AGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA  >  NC_000913/4116913‑4117207

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