BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,005	T→A	22.7%	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,005	0	T	A	22.7%	236.0 / 71.0	110	D73V (GAT→GTT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base T (44/41); new base A (15/10); total (59/51)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 116 total aligned reads displayed.

TTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGT  >  NC_000913/4116857‑4117152

ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                                                                                     <  3:386095/149‑1 (MQ=255)

ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                                                                                     <  4:512755/149‑1 (MQ=255)

 tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                                                                                                    <  4:52537/149‑1 (MQ=255)

 tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                                                                                                    <  4:57947/149‑1 (MQ=255)

   tAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACg                                                                                                                                                   >  1:401824/1‑148 (MQ=255)

    aaaCCGTCACGTTTTAAATGCTCGCAGCTTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGa                                                                                                                                                 >  4:845706/1‑149 (MQ=255)

      aCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTGCCGGACACCAGTGGTTTCACGCTGGTTCGTAATACCGATAACTGCAATTTGATCGGaa                                                                                                                                                <  2:235572/148‑1 (MQ=255)

       ccGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACt                                                                                                                                              <  4:673107/149‑1 (MQ=255)

           cACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGAAAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGat                                                                                                                                           >  1:69406/1‑148 (MQ=255)

             cGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGatat                                                                                                                                         >  2:305710/1‑148 (MQ=255)

              gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCg                                                                                                                                       >  3:204740/1‑149 (MQ=255)

               ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCgg                                                                                                                                      >  3:614390/1‑149 (MQ=255)

                tttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCgg                                                                                                                                      >  2:766134/1‑148 (MQ=255)

                tttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCgg                                                                                                                                      >  2:931224/1‑148 (MQ=255)

                tttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCgg                                                                                                                                      <  1:372149/148‑1 (MQ=255)

                tttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCgg                                                                                                                                      >  2:77769/1‑148 (MQ=255)

                 ttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGc                                                                                                                                     <  2:623770/148‑1 (MQ=255)

                    aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCttt                                                                                                                                  <  2:47840/148‑1 (MQ=255)

                    aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTc                                                                                                                                 <  4:384803/149‑1 (MQ=255)

                    aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCttt                                                                                                                                  <  2:367571/148‑1 (MQ=255)

                     aTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCg                                                                                                                                >  3:334551/1‑149 (MQ=255)

                     aTGCTCGCAGAGTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                                                                                     >  1:455886/1‑128 (MQ=255)

                      tGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCg                                                                                                                                >  1:727850/1‑148 (MQ=255)

                       gCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCGGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGcc                                                                                                                              <  3:862910/149‑1 (MQ=255)

                         tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAg                                                                                                                            >  4:436657/1‑149 (MQ=255)

                                 ttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTa                                                                                                                    >  3:885837/1‑149 (MQ=255)

                                  tCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAc                                                                                                                   <  4:65735/149‑1 (MQ=255)

                                    tGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCa                                                                                                                 >  4:533115/1‑149 (MQ=255)

                                      cGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATTGGCTTTCGCCAGCACTTCTACCAg                                                                                                                >  1:365646/1‑148 (MQ=255)

                                      cGGTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAg                                                                                                                >  1:812308/1‑148 (MQ=255)

                                        gTACGACGGCACTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGt                                                                                                             <  3:70533/149‑1 (MQ=255)

                                        gTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTGGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCg                                                                                                              >  2:102808/1‑148 (MQ=255)

                                                gCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCtt                                                                                                      <  1:151677/148‑1 (MQ=255)

                                                 cACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTg                                                                                                    <  3:436657/149‑1 (MQ=255)

                                                  aCTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                                                                                   >  3:737537/1‑149 (MQ=255)

                                                  aCTGCCAGACAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTg                                                                                                    <  2:69406/148‑1 (MQ=255)

                                                  aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTgg                                                                                                   <  4:102034/149‑1 (MQ=255)

                                                  aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTgg                                                                                                   >  4:388006/1‑149 (MQ=255)

                                                  aCTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTgg                                                                                                   <  4:614390/149‑1 (MQ=255)

                                                        aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGc                                                                                              >  2:123927/1‑148 (MQ=255)

                                                          aCAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTggg                                                                                                  >  2:256910/1‑142 (MQ=255)

                                                          aCAATGGCGTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTggg                                                                                                  <  1:256910/142‑1 (MQ=255)

                                                                  gTTATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                                                                                    <  2:759543/148‑1 (MQ=255)

                                                                  gTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCa                                                                                                                 <  4:881476/119‑1 (MQ=255)

                                                                  gTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCa                                                                                                                 >  3:881476/1‑119 (MQ=255)

                                                                   ttATAGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAt                                                                                  >  3:426284/1‑149 (MQ=255)

                                                                    tataGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                                                         >  1:704068/1‑141 (MQ=255)

                                                                    tataGATATGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAtt                                                                                 <  3:569528/149‑1 (MQ=255)

                                                                     atagatatGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAtt                                                                                 <  1:226722/148‑1 (MQ=255)

                                                                      tagatatGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg                                                                               <  4:885837/149‑1 (MQ=255)

                                                                          tagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACtt                                                                                                                       >  4:451977/1‑105 (MQ=255)

                                                                          tagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACtt                                                                                                                       <  3:451977/105‑1 (MQ=255)

                                                                            ggCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCGGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                                                          >  1:238244/1‑148 (MQ=255)

                                                                             gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGCGCTTTGGGTGGCCCAGATTTCCATTGGGTAgtg                                                                         >  1:172033/1‑148 (MQ=255)

                                                                             gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                                                        >  3:616545/1‑149 (MQ=255)

                                                                              cTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGtt                                                                       <  3:600279/149‑1 (MQ=255)

                                                                                 gCCGGTTTCTTTTTCCCAGACCATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGaa                                                                                                                                                <  4:475609/73‑1 (MQ=255)

                                                                                 gCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGaa                                                                                                                                                >  3:475609/1‑73 (MQ=255)

                                                                                 gCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                                                    <  3:423386/149‑1 (MQ=255)

                                                                                      tttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCa                                                                                   <  3:835837/129‑1 (MQ=255)

                                                                                      tttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCa                                                                                   >  4:835837/1‑129 (MQ=255)

                                                                                        tCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCGGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAcc                                                             >  4:10212/1‑149 (MQ=255)

                                                                                          tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATTGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTGCTACCAGCGTGGAGCTTTGGGTGGCCCAGTTTTCCATTGGGTCGTGTTCTACCCAACCTg                                                           >  4:86025/1‑149 (MQ=255)

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

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

                                                                                           ttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCg                                                                                                                                <  1:647516/79‑1 (MQ=255)

                                                                                           ttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCg                                                                                                                                >  2:647516/1‑79 (MQ=255)

                                                                                           ttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                           >  1:278387/1‑148 (MQ=255)

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

                                                                                             ttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGt                                                         <  1:924570/148‑1 (MQ=255)

                                                                                                ccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                     >  3:766181/1‑149 (MQ=255)

                                                                                                 cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGTTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                    >  4:779721/1‑149 (MQ=255)

                                                                                                 cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCGGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                     <  2:727850/148‑1 (MQ=255)

                                                                                                 cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                     >  2:665135/1‑148 (MQ=255)

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

                                                                                                 cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                     >  2:321573/1‑148 (MQ=255)

                                                                                                          gTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTTGGTGGCCCAGATTTCCATTGGGGCGTGTTCTACCCAACCTGGTTTTGGGGAGATTTg                                           >  4:256653/1‑149 (MQ=255)

                                                                                                           tGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTg                                           <  1:810621/148‑1 (MQ=255)

                                                                                                               ttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaa                                      >  3:715645/1‑149 (MQ=255)

                                                                                                                  aCGCTGGTTTGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAAta                                   >  4:822903/1‑148 (MQ=255)

                                                                                                                    gCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcg                                 <  3:533115/149‑1 (MQ=255)

                                                                                                                       ggTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCgcgc                              >  3:381125/1‑149 (MQ=255)

                                                                                                                       ggTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCCTGGAGCATTGGGTGGCCCAGATTTCCATTGGGTCGTGATCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgca                               >  1:692985/1‑147 (MQ=255)

                                                                                                                         tccGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTg                            >  4:655812/3‑149 (MQ=255)

                                                                                                                             tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCGGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTc                                  <  4:409413/139‑1 (MQ=255)

                                                                                                                             tAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGa                         >  1:612607/1‑148 (MQ=255)

                                                                                                                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                      >  1:194627/1‑118 (MQ=255)

                                                                                                                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                      <  2:194627/118‑1 (MQ=255)

                                                                                                                                   cGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCt                   >  2:226334/1‑148 (MQ=255)

                                                                                                                                    gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCGGCACTTGTAc                                                                                                                   >  3:653572/1‑51 (MQ=255)

                                                                                                                                    gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCGGCACTTCTAc                                                                                                                   <  4:653572/51‑1 (MQ=255)

                                                                                                                                       aGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGa              >  3:485744/1‑149 (MQ=255)

                                                                                                                                          tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat           >  4:328514/1‑149 (MQ=255)

                                                                                                                                            cAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTg         <  3:734894/149‑1 (MQ=255)

                                                                                                                                            cAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt          >  2:647132/1‑148 (MQ=255)

                                                                                                                                               tttGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAAtt                                   >  1:760586/1‑120 (MQ=255)

                                                                                                                                               tttGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAAtt                                   <  2:760586/120‑1 (MQ=255)

                                                                                                                                                   aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGt  >  3:108776/1‑149 (MQ=255)

                                                                                                                                                   aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATct   <  2:135926/148‑2 (MQ=255)

                                                                                                                                                    tCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGt  <  1:379707/148‑1 (MQ=255)

TTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGT  >  NC_000913/4116857‑4117152

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